Abstract
The normal order selection strategy used by a variable order, variable stepsize ordinary differential equation solver is based on selecting the formula that can use the largest stepsize while meeting the accuracy constraint. This strategy does not take into account the stiffness or nonstiffness of the differential system. By considering how the stiffness of the differential system influences order selection, a new selection strategy is proposed that is biased in the presence of stiffness toward the lower-order formulas. This new order selection strategy has been incorporated into a version of a widely used ordinary differential equation solver LSODE. The modified solver is shown to be greatly more efficient for problems that the original solver cannot solve efficiently, namely, problems where the Jacobian of the differential system has some eigenvalues of large modulus close to the imaginary axis. On other problems that the unmodified solver solves efficiently, the modified solver still performs with comparable efficiency.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
More From: SIAM Journal on Scientific and Statistical Computing
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.