MATRIX_EXPONENTIAL is a FORTRAN77 library which exhibits and compares some algorithms for approximating the matrix exponential function.
Formally, for a square matrix A and scalar t, the matrix exponential exp(A*t) can be defined as the sum:
exp(A*t) = sum ( 0 <= i < oo ) A^i t^i / i!
The simplest form of the matrix exponential problem asks for the value when t = 1. Even for this simple case, and for a matrix of small order, it can be quite difficult to compute the matrix exponential accurately.
MATRIX_EXPONENTIAL needs the R8LIB library. The test code for MATRIX_EXPONENTIAL requires the TEST_MATRIX_EXPONENTIAL library.
The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.
MATRIX_EXPONENTIAL is available in a C version and a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.
EXPOKIT, a FORTRAN77 library which solves various forms of the matrix exponential problem, by Roger Sidje.
R8LIB, a FORTRAN77 library which contains many utility routines using double precision real (R8) arithmetic.
TEST_MAT, a FORTRAN77 library which defines test matrices.
TEST_MATRIX_EXPONENTIAL, a FORTRAN77 library which defines a set of test cases for computing the matrix exponential.
You can go up one level to the FORTRAN77 source codes.