Efficient numerical approximation of a non-regular Fokker–Planck equation associated with first-passage time distributions

03/08/2021
by   Udo Boehm, et al.
0

In neuroscience, the distribution of a decision time is modelled by means of a one-dimensional Fokker–Planck equation with time-dependent boundaries and space-time-dependent drift. Efficient approximation of the solution to this equation is required, e.g., for model evaluation and parameter fitting. However, the prescribed boundary conditions lead to a strong singularity and thus to slow convergence of numerical approximations. In this article we demonstrate that the solution can be related to the solution of a parabolic PDE on a rectangular space-time domain with homogeneous initial and boundary conditions by transformation and subtraction of a known function. We verify that the solution of the new PDE is indeed more regular than the solution of the original PDE and proceed to discretize the new PDE using a space-time minimal residual method. We also demonstrate that the solution depends analytically on the parameters determining the boundaries as well as the drift. This justifies the use of a sparse tensor product interpolation method to approximate the PDE solution for various parameter ranges. The predicted convergence rates of the minimal residual method and that of the interpolation method are supported by numerical simulations.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

07/19/2020

Dynamic tensor approximation of high-dimensional nonlinear PDEs

We present a new method based on functional tensor decomposition and dyn...
03/21/2017

Controllability to Equilibria of the 1-D Fokker-Planck Equation with Zero-Flux Boundary Condition

We consider the problem of controlling the spatiotemporal probability di...
02/22/2022

Optimal Interpolation Data for PDE-based Compression of Images with Noise

We introduce and discuss shape-based models for finding the best interpo...
02/14/2020

Space-Time Collocation Method: Loop Quantum Hamiltonian Constraints

A space-time collocation method (STCM) using asymptotically-constant bas...
09/17/2011

A KdV-like advection-dispersion equation with some remarkable properties

We discuss a new non-linear PDE, u_t + (2 u_xx/u) u_x = epsilon u_xxx, i...
11/04/2020

Choosing The Best Interpolation Data in Images with Noise

We introduce and discuss shape based models for finding the best interpo...
02/24/2021

Convergence in the maximum norm of ADI-type methods for parabolic problems

Results on unconditional convergence in the Maximum norm for ADI-type me...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.