MADNESS
0.10.1
|
The source is here.
Given , we seek to solve the PDE
for the function at some future time (i.e., for ). is a linear operator that we are able to exponentiate, and is everything else including linear and non-linear parts.
In the semigroup approach the formal solution to the PDE is written
Numerical quadrature of the integral using Gauss-Legendre quadrature points is used resulting in a set of equations that are iteratively solved (presently using simple fixed point iteration from a first-order explicit rule).
The user provides
Have a look in testspectralprop.cc for example use.
With quadrature points, the error is and the number of applications of the exponential operator per time step is where is the number of iterations necessary to solve the equations (typically about 5 but this is problem dependent).