Related papers: An approximate stationary solution for multi-allel…
As a follow up to \cite{Causley2013}, we provide a detailed description of the numerical implementation of an O(N), A-stable, second order accurate solution of the wave equation, constructed from semi-discrete boundary value problems. We…
A generalisation of reaction diffusion systems and their travelling solutions to cases when the productive part of the reaction happens only on a surface in space or on a line on plane but the degradation and the diffusion happen in bulk…
We devise a stabilized method to weakly enforce bound constraints in the discrete solution of advection-dominated diffusion problems. This method combines a nonlinear penalty formulation with a discontinuous Galerkin-based residual…
We consider self-similar approximations of nonlinear hyperbolic systems in one space dimension with Riemann initial data and general diffusion matrix. We assume that the matrix of the system is strictly hyperbolic and the diffusion matrix…
Analytic solutions to the nonlinear radiation diffusion equation with an instantaneous point source for a non-homogeneous medium with a power law spatial density profile, are presented. The solutions are a generalization of the well known…
In this paper we introduce a multilocus diffusion model of a population of $N$ haploid, asexually reproducing individuals. The model includes parent-dependent mutation and interlocus selection, the latter limited to pairwise relationships…
In this paper, we study quasi-stationarity for a large class of Kolmogorov diffusions. The main novelty here is that we allow the drift to go to $- \infty$ at the origin, and the diffusion to have an entrance boundary at $+\infty$. These…
The evolution of the allelic proportion $x$ of a biallelic locus subject to the forces of mutation and drift is investigated in a diffusion model, assuming small scaled mutation rates. The overall scaled mutation rate is parametrized with…
A large class of physically important nonlinear and nonhomogeneous evolution problems, characterized by advection-like and diffusion-like processes, can be usefully studied by a time-differential form of Kolmogorov's solution of the…
We prove the existence and uniqueness of a quasi-stationary distribution for three stochastic processes derived from the model of Muller's ratchet. This model was invented with the aim of evaluating the limitations of an asexual…
The problem of quickest detection of a change in the distribution of a $n\times p$ random matrix based on a sequence of observations having a single unknown change point is considered. The forms of the pre- and post-change distributions of…
This paper establishes the spectral stability of monotone, stationary front solutions for reaction-diffusion equations where the reaction function is of Nagumo (or bistable) type and with diffusion coefficients which are density dependent…
A general system of n ordinary differential equations coupled with one reaction-diffusion equation, considered in a bounded N-dimensional domain, with no-flux boundary condition is studied in a context of pattern formation. Such initial…
This paper investigates the Keller-Segel model with quadratic cellular diffusion over a disk in $\mathbb R^2$ with a focus on the formation of its nontrivial patterns. We obtain explicit formulas of radially symmetric stationary solutions…
We consider a mutation-selection model of a population structured by the spatial variables and a trait variable which is the diffusion rate. Competition for resource is local in spatial variables, but nonlocal in the trait variable. We…
We consider the inverse problem of reconstructing the optical parameters of the radiative transfer equation (RTE) from boundary measurements in the diffusion limit. In the diffusive regime (the Knudsen number $\mathsf{Kn}\ll 1$), the…
The purpose of this paper is to provide a complete description of the eigenvalues of the generator of a neutral multi-type Moran model, and the applications to the study of the speed of convergence to stationarity. The Moran model we…
Building upon recent results obtained in [7,8,9], we describe an efficient second order, A-stable scheme for solving the wave equation, based on the method of lines transpose (MOL$^T$), and the resulting semi-discrete (i.e. continuous in…
The random reshuffling Kaczmarz (RRK) method enjoys the simplicity and efficiency in solving linear systems as a Kaczmarz-type method, whereas it also inherits the practical improvements of the stochastic gradient descent (SGD) with random…
We propose a fixed-point-based numerical framework for computing stationary states of nonlocal Fokker-Planck-type equations. Instead of discretising the differential operators directly, we reformulate the stationary problem as a nonlinear…