Abstract

Implicit two-step peer methods are introduced for the solution of large stiff systems. Although these methods compute s-stage approximations in each time step one-by-one like diagonally-implicit Runge-Kutta methods the order of all stages is the same due to the two-step structure. The nonlinear stage equations are solved by an inexact Newton method using the Krylov solver FOM (Arnoldi’s method). The methods are zero-stable for arbitrary step size sequences. We construct different methods having order p=s in the multi-implicit case and order p=s−1 in the singly-implicit case with arbitrary step sizes and s≤5. Numerical tests in Matlab for several semi-discretized partial differential equations show the efficiency of the methods compared to other Krylov codes.

Full Text
Paper version not known

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 call

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.