Estimating conditioning of BVPs for ODEs

Authors:
L. F. Shampine;P. H. Muir
Affiliations:
-;-
Venue:
Mathematical and Computer Modelling: An International Journal
Year:
2004

Citing 7
Cited 4

FORTRAN codes for estimating the one-norm of a real or complex matrix, with applications to condition estimation

ACM Transactions on Mathematical Software (TOMS)
Runge-Kutta Software with Defect Control four Boundary Value ODEs

SIAM Journal on Scientific Computing
LAPACK Users' guide (third ed.)

LAPACK Users' guide (third ed.)
A Block Algorithm for Matrix 1-Norm Estimation, with an Application to 1-Norm Pseudospectra

SIAM Journal on Matrix Analysis and Applications
FORTRAN Packages for Solving Certain Almost Block Diagonal Linear Systems by Modified Alternate Row and Column Elimination

ACM Transactions on Mathematical Software (TOMS)
A BVP solver based on residual control and the Maltab PSE

ACM Transactions on Mathematical Software (TOMS)
Nonlinear Equation Solvers in Boundary Value Problem Codes

Proceedings of a Working Conference on Codes for Boundary-Value Problems in Ordinary Differential Equations

The role of conditioning in mesh selection algorithms for first order systems of linear two point boundary value problems

Journal of Computational and Applied Mathematics - Special issue: International workshop on the technological aspects of mathematics
On the development of effective algorithms for the numerical solution of singularly perturbed two-point boundary value problems

International Journal of Computing Science and Mathematics
The role of conditioning in mesh selection algorithms for first order systems of linear two point boundary value problems

Journal of Computational and Applied Mathematics - Special issue: International workshop on the technological aspects of mathematics
A Runge-Kutta BVODE Solver with Global Error and Defect Control

ACM Transactions on Mathematical Software (TOMS)

Quantified Score

Hi-index	0.98

Visualization

Abstract

An alternative to control of the global error of a numerical solution to a boundary value problem (BVP) for ordinary differential equations (ODEs) is control of its residual, the amount by which it fails to satisfy the ODEs and boundary conditions. Among the methods used by codes that control residuals are collocation, Runge-Kutta methods with continuous extensions, and shooting. Specific codes that concern us are bvp4c of the Matlab problem solving environment and the FORTRAN code MIRKDC for general scientific computation. The residual of a numerical solution is related to its global error by a conditioning constant. In this paper, we investigate a conditioning constant appropriate for BVP solvers that control residuals and show how to estimate it numerically at a modest cost. Codes that control residuals can compute pseudosolutions, numerical solutions to BVPs that do not have solutions. That is, a ''well-behaved'' approximate solution is computed for an ill-posed mathematical problem. The estimate of conditioning is used to improve the robustness of bvp4c and MIRKDC and in particular, help users identify when a pseudosolution may have been computed.