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

by   Udo Boehm, et al.

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.



There are no comments yet.


page 1

page 2

page 3

page 4


Dynamic tensor approximation of high-dimensional nonlinear PDEs

We present a new method based on functional tensor decomposition and dyn...

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...

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

We introduce and discuss shape-based models for finding the best interpo...

Space-Time Collocation Method: Loop Quantum Hamiltonian Constraints

A space-time collocation method (STCM) using asymptotically-constant bas...

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...

Choosing The Best Interpolation Data in Images with Noise

We introduce and discuss shape based models for finding the best interpo...

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.