Abstract
In the stiff situation, we consider the long-time behavior of the relative error gamma _n in the numerical integration of a linear ordinary differential equation y^{prime }(t)=Ay(t),quad tge 0, where A is a normal matrix. The numerical solution is obtained by using at any step an approximation of the matrix exponential, e.g. a polynomial or a rational approximation. We study the long-time behavior of gamma _n by comparing it to the relative error gamma _n^{mathrm{long}} in the numerical integration of the long-time solution, i.e. the projection of the solution on the eigenspace of the rightmost eigenvalues. The error gamma _n^{mathrm{long}} grows linearly in time, it is small and it remains small in the long-time. We give a condition under which gamma _napprox gamma _n^{mathrm{long}}, i.e. frac{gamma _n}{gamma _n^{mathrm{long}}}approx 1, in the long-time. When this condition does not hold, the ratio frac{gamma _n}{gamma _n^{mathrm{long}}} is large for all time. These results describe the long-time behavior of the relative error gamma _n in the stiff situation.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: 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.