Operator splitting algorithm for structured population models on metric spaces
Fuente:
arXiv
Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866914403203416064 |
|---|---|
| author | Lindow, Carolin Düll, Christian Gwiazda, Piotr Miasojedow, Błażej Marciniak-Czochra, Anna |
| author_facet | Lindow, Carolin Düll, Christian Gwiazda, Piotr Miasojedow, Błażej Marciniak-Czochra, Anna |
| contents | In this paper, we propose a numerical scheme for structured population models defined on a separable and complete metric space. In particular, we consider a generalized version of a transport equation with additional growth and non-local interaction terms in the space of nonnegative Radon measures equipped with the flat metric. The solutions, given by families of Radon measures, are approximated by linear combinations of Dirac measures. For this purpose, we introduce a finite-range approximation of the measure-valued model functions, provided that they are linear. By applying an operator splitting technique, we are able to separate the effects of the transport from those of growth and the non-local interaction. We derive the order of convergence of the numerical scheme, which is linear in the spatial discretization parameters and polynomial of order $α$ in the time step size, assuming that the model functions are $α$ Hölder regular in time. In a second step, we show that our proposed algorithm can approximate the posterior measure of Bayesian inverse models, which will allow us to link model parameters to measured data in the future. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2502_07341 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Operator splitting algorithm for structured population models on metric spaces Lindow, Carolin Düll, Christian Gwiazda, Piotr Miasojedow, Błażej Marciniak-Czochra, Anna Numerical Analysis 28-08, 65J08, 65J22, 65M57 In this paper, we propose a numerical scheme for structured population models defined on a separable and complete metric space. In particular, we consider a generalized version of a transport equation with additional growth and non-local interaction terms in the space of nonnegative Radon measures equipped with the flat metric. The solutions, given by families of Radon measures, are approximated by linear combinations of Dirac measures. For this purpose, we introduce a finite-range approximation of the measure-valued model functions, provided that they are linear. By applying an operator splitting technique, we are able to separate the effects of the transport from those of growth and the non-local interaction. We derive the order of convergence of the numerical scheme, which is linear in the spatial discretization parameters and polynomial of order $α$ in the time step size, assuming that the model functions are $α$ Hölder regular in time. In a second step, we show that our proposed algorithm can approximate the posterior measure of Bayesian inverse models, which will allow us to link model parameters to measured data in the future. |
| title | Operator splitting algorithm for structured population models on metric spaces |
| topic | Numerical Analysis 28-08, 65J08, 65J22, 65M57 |
| url | https://arxiv.org/abs/2502.07341 |