DeepAI
Log In Sign Up

Efficient numerical computation of the basic reproduction number for structured populations

04/27/2020
by   Dimitri Breda, et al.
0

As widely known, the basic reproduction number plays a key role in weighing birth/infection and death/recovery processes in several models of population dynamics. In this general setting, its characterization as the spectral radius of next generation operators is rather elegant, but simultaneously poses serious obstacles to its practical determination. In this work we address the problem numerically by reducing the relevant operators to matrices through a pseudospectral collocation, eventually computing the sought quantity by solving finite-dimensional eigenvalue problems. The approach is illustrated for two classes of models, respectively from ecology and epidemiology. Several numerical tests demonstrate experimentally important features of the method, like fast convergence and influence of the smoothness of the models' coefficients. Examples of robust analysis of instances of specific models are also presented to show potentialities and ease of application.

READ FULL TEXT

page 1

page 2

page 3

page 4

11/27/2021

Asymptotic spectra of large (grid) graphs with a uniform local structure (part II): numerical applications

In the current work we are concerned with sequences of graphs having a g...
11/24/2021

A global quadratic speed-up for computing the principal eigenvalue of Perron-like operators

We consider a new algorithm in light of the min-max Collatz-Wielandt for...
03/22/2019

The basic reproduction number, R_0, in structured populations

In this paper, we provide a straightforward approach to defining and der...
05/05/2020

Dynamics of an Intra-host Diffusive Pathogen Infection Model

In this paper, we first propose a diffusive pathogen infection model wit...
09/12/2019

Numeric Solutions of Eigenvalue Problems for Generic Nonlinear Operators

Numerical methods for solving linear eigenvalue problem are widely studi...
09/07/2021

Bivariate collocation for computing R_0 in epidemic models with two structures

Structured epidemic models can be formulated as first-order hyperbolic P...