timflow.transient.invlapnumba.invlap#

timflow.transient.invlapnumba.invlap(t, tmax, fp, M, alpha=1e-10, tol=1e-09)#

Inverse Laplace tansform with algorithm of De Hoog, Knight and Stokes.

Parameters:

  • t (array) – times for which inverse is computed

  • tmax (float) – maximum time

  • fp (complex array) – Laplace transformed solution

  • M (integer) – number of terms (number of values of fp)

  • alpha (is the real part of the rightmost pole or singularity, which) – is chosen based on the desired accuracy (assuming the rightmost singularity is 0), and tol=10α is the desired tolerance

Returns:

  • result (array) – time domain solution for specified times

  • Reference

  • ———

  • de Hoog, F., J. Knight, A. Stokes (1982). An improved method for

  • numerical inversion of Laplace transforms. SIAM Journal of Scientific

  • and Statistical Computing 3 (357-366, http://dx.doi.org/10.1137/0903022)

  • https (//bitbucket.org/klkuhlm/invlap/src/default/invlap.py)