Algorithm 962

BACOL and BACOLR are (Fortran 77) B-spline adaptive collocation packages for the numerical solution of 1D parabolic Partial Differential Equations (PDEs). The packages have been shown to be superior to other similar packages, especially for problems exhibiting sharp, moving spatial layer regions, where a stringent tolerance is imposed. In addition to providing temporal error control through the timestepping software, BACOL and BACOLR feature control of a high-order estimate of the spatial error of the approximate solution, obtained by computing a second approximate solution of one higher order of accuracy; the cost is substantial—execution time and memory usage are almost doubled. In this article, we discuss BACOLI, a new version of BACOL that computes only one approximate solution and uses efficient interpolation-based schemes to obtain a spatial error estimate. In previous studies these schemes have been shown to provide spatial error estimates of comparable quality to those of BACOL. We describe the substantial modification of BACOL needed to obtain BACOLI, and provide numerical results showing that BACOLI is significantly more efficient than BACOL, in some cases by as much as a factor of 2. We also introduce a Fortran 95 wrapper for BACOLI (called BACOLI95) and discuss its simplified user interface.

[1]  Patrick Keast,et al.  Algorithm 603: COLROW and ARCECO: FORTRAN Packages for Solving Certain Almost Block Diagonal Linear Systems by Modified Alternate Row and Column Elimination , 1983, TOMS.

[2]  Seymour V. Parter,et al.  Collocation methods for parabolic partial differential equations in one space dimension , 1976 .

[3]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[4]  T. Smith,et al.  ASYMPTOTICALLY CORRECT INTERPOLATION-BASED SPATIAL ERROR ESTIMATION FOR 1D PDE SOLVERS , 2013 .

[5]  Paul H. Muir,et al.  Algorithm 874: BACOLR—spatial and temporal error control software for PDEs based on high-order adaptive collocation , 2008, TOMS.

[6]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis: Oden/A Posteriori , 2000 .

[7]  Gary R. Consolazio,et al.  Finite Elements , 2007, Handbook of Dynamic System Modeling.

[8]  Peter K. Moore,et al.  Interpolation error-based a posteriori error estimation for hp -refinement using first and second derivative jumps , 2004 .

[9]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[10]  W. Finden,et al.  An error term and uniqueness for Hermite-Birkhoff interpolation involving only function values and/or first derivative values , 2008 .

[11]  Peter K. Moore,et al.  Applications of Lobatto polynomials to an adaptive finite element method: A posteriori error estimates for HP-Adaptivity and Grid-to-Grid Interpolation , 2003, Numerische Mathematik.

[12]  Peter K. Moore,et al.  An implicit interpolation error-based error estimation strategy for hp-adaptivity in one space dimension , 2004, J. Num. Math..

[13]  Zhi Li,et al.  B -Spline Gaussian Collocation Software for Two-Dimensional Parabolic PDEs , 2013 .

[14]  Roy D. Williams,et al.  Estimating the Error of Numerical Solutions of Systems of Reaction-Diffusion Equations , 2000 .

[15]  C. R. Deboor,et al.  A practical guide to splines , 1978 .

[16]  Rong Wang,et al.  A comparison of adaptive software for 1D parabolic PDEs , 2004 .

[17]  Jim Douglas,et al.  Collocation Methods for Parabolic Equations in a Single Space Variable , 1974 .

[18]  T. Smith,et al.  SUPERCONVERGENT INTERPOLANTS FOR EFFICIENT SPATIAL ERROR ESTIMATION IN 1D PDE COLLOCATION SOLVERS , 2022 .

[19]  Peter K. Moore,et al.  Interpolation error-based a posteriori error estimation for two-point boundary value problems and parabolic equations in one space dimension , 2001, Numerische Mathematik.

[20]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[21]  Zhi Li,et al.  A Computational Study of the Efficiency of Collocation Software for 1 D Parabolic PDEs with Interpolation-based Spatial Error Estimation ∗ , .

[22]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis , 2000 .

[23]  M. A. Akanbi,et al.  Numerical solution of initial value problems in differential - algebraic equations , 2005 .