Abstract
In order to solve a differential problem, the Laplace Transform method, when applicable, replaces the problem with a simpler one; the solution is obtained by solving the new problem and then by computing the inverse Laplace Transform of this function. In a numerical context, since the solution of the transformed problem consists of a sequence of Laplace Transform samples, most of the software for the numerical inversion cannot be used since the transform, among parameters, must be passed as a function.
To fill this gap, we present Talbot Suite DE, a C software collection for Laplace Transform inversions, specifically designed for these problems and based on Talbot’s method. It contains both sequential and parallel implementations; the latter is accomplished by means of OpenMP. We also report some performance results.
Aimed at non-expert users, the software is equipped with several examples and a User Guide that includes the external documentation, explains how to use all the sample code, and reports its results about accuracy and efficiency. Some examples are entirely in C and others combine different programming languages (C/MATLAB, C/FORTRAN). The User Guide also contains useful hints to avoid possible errors issued during the compilation or execution of mixed-language code.
Supplemental Material
Available for Download
Software for Talbot Suite DE: Application of Modified Talbots Method to Solve Differential Problems
- A. Ghizzetti and A. Ossicini. 1943,1972. Trasformate di Laplace e Calcolo Simbolico. UTET.Google Scholar
- A. Murli and M. Rizzardi. 1990. Talbot’s method for the Laplace inversion problem. ACM Trans. Math. Soft. 16, 2 (1990), 158--168. Google ScholarDigital Library
- F. R. de Hoog, J. H. Knight, and A. N. Stokes. 1982. An improved method for numerical inversion of Laplace transforms. SIAM J. Sci. Stat. Comput. 3 (1982), 357--366. Google ScholarDigital Library
- D. G. Duffy. 2004. Transform Methods for Solving Partial Differential Equations. Chapman 8 Hall/CRC.Google Scholar
- J. Burkardt. 2012. C-version of ODE initial-value problem solver. Retrieved from http://people.sc.fsu.edu/∼jburkardt/c_src/c_src.html.Google Scholar
- J. L. Schiff. 1999. The Laplace Transform: Theory and Applications. Springer.Google Scholar
- J. R. Cash, D. Hollevoet, F. Mazzia, and M. Nagy. 2013. Algorithm 927: The MATLAB code bvptwp.m for the numerical solution of two point boundary value problems. ACM Trans. Math. Softw. 39, 2 (2013), 15:1--15:12. Google ScholarDigital Library
- J. R. Cash and M. H. Wright. 1990. Implementation issues in solving nonlinear equations for two-point boundary value problems. Computing 45 (1990), 17--37. Google ScholarDigital Library
- J. R. Cash and M. H. Wright. 1991. A deferred correction method for nonlinear two-point boundary value problems: Implementation and numerical evaluation. SIAM J. Sci. Stat. Comput. 12 (1991), 971--989.Google ScholarDigital Library
- J. R. Cash and M. H. Wright. 1996. twpbvp.f. Retrieved from http://www.netlib.org/ode/.Google Scholar
- L. Antonelli, S. Corsaro, Z. Marino, and M. Rizzardi. 2014. Algorithm 944: Talbot suite: Parallel implementations of Talbot’s method for the numerical inversion of Laplace transforms. ACM Trans. Math. Softw. 40, 4 (2014), 29:1--29:18. Google ScholarDigital Library
- L. F. Shampine and M. K. Gordon. 1975. Computer Solution of Ordinary Differential Equations: The Initial Value Problem. Freeman.Google Scholar
- L. F. Shampine and M. K. Gordon. 1975. Ordinary differential equation initial-value problem solver (alg.: Adam’s methods) double precision FORTRAN. Retrieved from http://www.netlib.org/ode/.Google Scholar
- M. Rizzardi. 1995. A modification of Talbot’s method for the simultaneous approximation of several values of the inverse Laplace transform. ACM Trans. Math. Softw. 21, 4 (1995), 347--371. Google ScholarDigital Library
- M. Rizzardi. 2017. Detection of the singularities of a complex function by numerical approximations of its Laurent coefficients. Numer. Algor. (online: May 2017).Google Scholar
- M. Rizzardi. 2017. LaurentCoeffs: MATLAB Suite to detect the Singularities of a Complex Function. User Guide. Technical Report. http://dist.uniparthenope.it/mariarosaria.rizzardi.Google Scholar
- M. A. de Rosa, G. Giunta, and M. Rizzardi. 1995. Parallel Talbot’s algorithm for distributed memory machines. Parallel Comput. 21, 5 (1995), 783--801. Google ScholarDigital Library
- Maplesoft. 2016. Maple. Retrieved from http://www.maplesoft.com/products/Maple/.Google Scholar
- MathWorks. 2016. MATLAB. Retrieved from http://www.mathworks.com/.Google Scholar
- MathWorks. 2016. MATLAB Symbolic Math Toolbox. Retrieved from http://mathworks.com/help/symbolic/.Google Scholar
- MCS Div., Argonne National Laboratory. 2016. The Message Passing Interface (MPI) standard. Retrieved from http://www.mcs.anl.gov/research/projects/mpi/.Google Scholar
- MPICH Org. 2016. MPICH: A high performance and widely portable implementation of the Message Passing Interface (MPI) standard. Retrieved from http://www.mpich.org/.Google Scholar
- OpenMP Architecture Review Board. 2017. The OpenMP API specification for parallel programming. http://www.openmp.org/. (2017).Google Scholar
- S. Balay, S. Abhyankar, M. F. Adams, J. Brown, P. Brune, K. Buschelman, L. Dalcin, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. Curfman McInnes, K. Rupp, B. F. Smith, S. Zampini, and H. Zhang. 2016. PETSc Web page. Retrieved from https://www.mcs.anl.gov/petsc.Google Scholar
- A. Talbot. 1979. The accurate numerical inversion of Laplace transforms. J. Inst. Math. Appl. 23 (1979), 97--120.Google ScholarCross Ref
- T. Weeks. 1966. Numerical inversion of Laplace transform using Laguerre functions. J. ACM 13 (1966). Google ScholarDigital Library
- WOLFRAM. 2016. MATHEMATICA. Retrieved from https://www.wolfram.com/mathematica/.Google Scholar
Index Terms
- Algorithm 981: Talbot Suite DE: Application of Modified Talbot’s Method to Solve Differential Problems
Recommendations
A modification of Talbot's method for the simultaneous approximation of several values of the inverse Laplace transform
In recent years many results have been obtained in the field of the numerical inversion of Laplace transforms. Among them, a very accurate and general method is due to Talbot: this method approximates the value of the inverse Laplace transform f(t), for ...
Algorithm 944: Talbot Suite: Parallel Implementations of Talbot's Method for the Numerical Inversion of Laplace Transforms
We present Talbot Suite, a C parallel software collection for the numerical inversion of Laplace Transforms, based on Talbot's method. It is designed to fit both single and multiple Laplace inversion problems, which arise in several application and ...
An improved Talbot method for numerical Laplace transform inversion
The classical Talbot method for the computation of the inverse Laplace transform is improved for the case where the transform is analytic in the complex plane except for the negative real axis. First, by using a truncated Talbot contour rather than the ...
Comments