Collocation for two-point boundary value problems revisited
SIAM Journal on Numerical Analysis
The computational complexity of differential and integral equations: an information-based approach
The computational complexity of differential and integral equations: an information-based approach
Global bounds on numerical error for ordinary differential equations
Journal of Complexity - Festschrift for Joseph F. Traub, Part 1
Algorithm 569: COLSYS: Collocation Software for Boundary-Value ODEs [D2]
ACM Transactions on Mathematical Software (TOMS)
Convergence of Adaptive Finite Element Methods
SIAM Review
Digital filters in adaptive time-stepping
ACM Transactions on Mathematical Software (TOMS)
Explicit, Time Reversible, Adaptive Step Size Control
SIAM Journal on Scientific Computing
Efficient mesh selection for collocation methods applied to singular BVPs
Journal of Computational and Applied Mathematics
Time-step selection algorithms: adaptivity, control, and signal processing
Applied Numerical Mathematics - The third international conference on the numerical solutions of volterra and delay equations, May 2004, Tempe, AZ
Constant coefficient linear multistep methods with step density control
Journal of Computational and Applied Mathematics
Computational complexity of numerical solutions of initial value problems for differential algebraic equations
Automatic grid control in adaptive BVP solvers
Numerical Algorithms
Applied Numerical Mathematics
Hi-index | 0.00 |
In this paper we analyze the problem of adaptivity for one-step numerical methods for solving ODEs, both IVPs and BVPs, with a view to generating grids of minimal computational cost for which the local error is below a prescribed tolerance (optimal grids). The grids are generated by introducing an auxiliary independent variable @t and finding a grid deformation map, t=@Q(@t), that maps an equidistant grid {@t"j} to a non-equidistant grid in the original independent variable, {t"j}. An optimal deformation map @Q is determined by a variational approach. Finally, we investigate the cost of the solution procedure and compare it to the cost of using equidistant grids. We show that if the principal error function is non-constant, an adaptive method is always more efficient than a non-adaptive method.