Abstract
FORTRAN 95 software is provided for computing the matrix values of φ-functions required in exponential integrators. The subroutines in the library accept as their argument a full, diagonal, or upper quasitriangular matrix with real or complex entries in one of four precisions. Two different algorithms are implemented, one is the scaling and squaring method, and the other is a modified block Schur--Parlett algorithm. In the latter algorithm, a recursive three-by-three blocking is applied to the argument based on an estimate of the sep-inverse function. The estimation of the sep-inverse function is carried out by Hager--Higham estimator implemented as the subroutine xLACON in LAPACK. Our modifications to the block Schur--Parlett algorithm are described together with the results of numerical experiments.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
