Mathematics – Numerical Analysis
Scientific paper
2009-07-27
Mathematics
Numerical Analysis
20 pages, 3 colour figures, code available from http://www.maths.leeds.ac.uk/~jitse/software.html . v2: Various changes to i
Scientific paper
We develop an algorithm for computing the solution of a large system of linear ordinary differential equations (ODEs) with polynomial inhomogeneity. This is equivalent to computing the action of a certain matrix function on the vector representing the initial condition. The matrix function is a linear combination of the matrix exponential and other functions related to the exponential (the so-called phi-functions). Such computations are the major computational burden in the implementation of exponential integrators, which can solve general ODEs. Our approach is to compute the action of the matrix function by constructing a Krylov subspace using Arnoldi or Lanczos iteration and projecting the function on this subspace. This is combined with time-stepping to prevent the Krylov subspace from growing too large. The algorithm is fully adaptive: it varies both the size of the time steps and the dimension of the Krylov subspace to reach the required accuracy. We implement this algorithm in the Matlab function phipm and we give instructions on how to obtain and use this function. Various numerical experiments show that the phipm function is often significantly more efficient than the state-of-the-art.
Niesen Jitse
Wright Will M.
No associations
LandOfFree
A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-246077