Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
skip to main content
research-article

Algorithm 962: BACOLI: B-spline Adaptive Collocation Software for PDEs with Interpolation-Based Spatial Error Control

Published: 23 May 2016 Publication History

Abstract

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.

Supplementary Material

ZIP File (962.zip)
Software for BACOLI: B-spline Adaptive Collocation Software for PDEs with Interpolation-Based Spatial Error Control

References

[1]
M. Ainsworth and J. T. Oden. 2000. A Posteriori Error Estimation in Finite Element Analysis. Wiley-Interscience, New York. xx+240 pages.
[2]
T. Arsenault, T. Smith, and P. H. Muir. 2009. Superconvergent interpolants for efficient spatial error estimation in 1D PDE collocation solvers. Can. Appl. Math. Q. 17 (2009), 409--431.
[3]
T. Arsenault, T. Smith, P. H. Muir, and P. Keast. 2011. Efficient Interpolation-Based Error Estimation for 1D Time-Dependent PDE Collocation Codes. Technical Report 2011_001. Saint Mary’s University, Dept. of Mathematics and Computing Science, Halifax, Nova Scotia, Canada.
[4]
T. Arsenault, T. Smith, P. H. Muir, and J. Pew. 2012. Asymptotically correct interpolation-based spatial error estimation for 1D PDE solvers. Can. Appl. Math. Q. 20 (2012), 307--328.
[5]
U. M. Ascher, R. M. M. Mattheij, and R. D. Russell. 1995. Numerical Solution of Boundary Value Problems for Ordinary Dfferential Euations. Classics in Applied Mathematics, Vol. 13. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA.
[6]
I. Babuška, J. R. Whiteman, and T. Strouboulis. 2011. Finite Elements. Oxford University Press, Oxford. xii+323 pages.
[7]
K. E. Brenan, S. L. Campbell, and L. R. Petzold. 1996. Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations. Classics in Applied Mathematics, Vol. 14. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA. x+256 pages.
[8]
J. H. Cerutti and S. V. Parter. 1976. Collocation methods for parabolic partial differential equations in one space dimension. Numer. Math. 26, 3 (1976), 227--254.
[9]
C. de Boor. 1978. A Practical Guide to Splines. Applied Mathematical Sciences, Vol. 27. Springer-Verlag, New York. xxiv+392 pages.
[10]
J. C. Díaz, G. Fairweather, and P. Keast. 1983. Algorithm 603. COLROW and ARCECO: FORTRAN packages for solving certain almost block diagonal linear systems by modified alternate row and column elimination. ACM Trans. Math. Software 9, 3 (1983), 376--380.
[11]
J. Douglas, Jr. and T. Dupont. 1974. Collocation Methods for Parabolic Equations in a Single Space Variable. Springer-Verlag, Berlin. iv+147 pages.
[12]
D. J. Estep, M. G. Larson, and R. D. Williams. 2000. Estimating the error of numerical solutions of systems of reaction-diffusion equations. Mem. Amer. Math. Soc. 146, 696 (2000), viii+109.
[13]
W. F. Finden. 2008. An error term and uniqueness for Hermite-Birkhoff interpolation involving only function values and/or first derivative values. J. Comput. Appl. Math. 212, 1 (2008), 1--15. 10.1016/j.cam.2006.11.022
[14]
E. Hairer, S. P. Nørsett, and G. Wanner. 1993. Solving Ordinary Differential Equations. I (2nd ed.). Springer Series in Computational Mathematics, Vol. 8. Springer-Verlag, Berlin. xvi+528 pages.
[15]
E. Hairer and G. Wanner. 1996. Solving Ordinary Differential Equations. II (2nd ed.). Springer Series in Computational Mathematics, Vol. 14. Springer-Verlag, Berlin. xvi+614 pages.
[16]
Z. Li and P. H. Muir. 2013. B-spline Gaussian collocation software for two-dimensional parabolic PDEs. Adv. Appl. Math. Mech. 5 (2013), 528--547.
[17]
P. K. Moore. 2001. Interpolation error-based a posteriori error estimation for two-point boundary value problems and parabolic equations in one space dimension. Numer. Math. 90, 1 (2001), 149--177.
[18]
P. K. Moore. 2003. Applications of Lobatto polynomials to an adaptive finite element method: A posteriori error estimates for hp-adaptivity and grid-to-grid interpolation. Numer. Math. 94, 2 (2003), 367--401.
[19]
P. K. Moore. 2004a. An implicit interpolation error-based error estimation strategy for hp-adaptivity in one space dimension. J. Numer. Math. 12, 2 (2004), 143--167.
[20]
P. K. Moore. 2004b. Interpolation error-based a posteriori error estimation for hp-refinement using first and second derivative jumps. Appl. Numer. Math. 48, 1 (2004), 63--82.
[21]
P. H. Muir. 2016. Website for Paul Muir. http://cs.smu.ca/∼muir.
[22]
T. Pew, Z. Li, and P. H. Muir. 2013. A Computational Study of the Efficiency of Collocation Software for 1D Parabolic PDEs with Interpolation-Based Spatial Error Estimation. Technical Report 2013_001. Saint Mary’s University, Dept. of Mathematics and Computing Science, Halifax, Nova Scotia, Canada.
[23]
R. Wang, P. Keast, and P. H. Muir. 2004a. BACOL: B-spline adaptive collocation software for 1D parabolic PDEs. ACM Trans. Math. Software 30, 4 (2004), 454--470.
[24]
R. Wang, P. Keast, and P. H. Muir. 2004b. A comparison of adaptive software for 1D parabolic PDEs. J. Comput. Appl. Math. 169, 1 (2004), 127--150.
[25]
R. Wang, P. Keast, and P. H. Muir. 2008. Algorithm 874: BACOLR—Spatial and temporal error control software for PDEs based on high-order adaptive collocation. ACM Trans. Math. Software 34, 3 (2008), Article 15, 28.

Cited By

View all
  • (2024)Differential equation software for the computation of error-controlled continuous approximate solutionsNumerical Algorithms10.1007/s11075-024-01784-196:3(1021-1044)Online publication date: 1-Jul-2024
  • (2021)Solving Cardiac Bidomain Problems with B-spline Adaptive CollocationRecent Developments in Mathematical, Statistical and Computational Sciences10.1007/978-3-030-63591-6_28(301-311)Online publication date: 30-Aug-2021
  • (2021)bacoli_py—A Python Package for the Error Controlled Numerical Solution of 1D Time-Dependent PDEsRecent Developments in Mathematical, Statistical and Computational Sciences10.1007/978-3-030-63591-6_27(289-300)Online publication date: 30-Aug-2021
  • Show More Cited By

Index Terms

  1. Algorithm 962: BACOLI: B-spline Adaptive Collocation Software for PDEs with Interpolation-Based Spatial Error Control

    Recommendations

    Comments

    Information & Contributors

    Information

    Published In

    cover image ACM Transactions on Mathematical Software
    ACM Transactions on Mathematical Software  Volume 42, Issue 3
    June 2016
    208 pages
    ISSN:0098-3500
    EISSN:1557-7295
    DOI:10.1145/2935754
    Issue’s Table of Contents
    Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than the author(s) must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected].

    Publisher

    Association for Computing Machinery

    New York, NY, United States

    Publication History

    Published: 23 May 2016
    Accepted: 01 August 2015
    Revised: 01 March 2015
    Received: 01 February 2014
    Published in TOMS Volume 42, Issue 3

    Permissions

    Request permissions for this article.

    Check for updates

    Badges

    Author Tags

    1. 1D parabolic PDEs
    2. Numerical software
    3. collocation
    4. efficiency
    5. interpolation
    6. spatial error estimation

    Qualifiers

    • Research-article
    • Research
    • Refereed

    Funding Sources

    • Mathematics of Information Technology and Complex Systems Network
    • Saint Mary's University
    • Natural Sciences and Engineering Research Council of Canada

    Contributors

    Other Metrics

    Bibliometrics & Citations

    Bibliometrics

    Article Metrics

    • Downloads (Last 12 months)10
    • Downloads (Last 6 weeks)0
    Reflects downloads up to 28 Oct 2024

    Other Metrics

    Citations

    Cited By

    View all
    • (2024)Differential equation software for the computation of error-controlled continuous approximate solutionsNumerical Algorithms10.1007/s11075-024-01784-196:3(1021-1044)Online publication date: 1-Jul-2024
    • (2021)Solving Cardiac Bidomain Problems with B-spline Adaptive CollocationRecent Developments in Mathematical, Statistical and Computational Sciences10.1007/978-3-030-63591-6_28(301-311)Online publication date: 30-Aug-2021
    • (2021)bacoli_py—A Python Package for the Error Controlled Numerical Solution of 1D Time-Dependent PDEsRecent Developments in Mathematical, Statistical and Computational Sciences10.1007/978-3-030-63591-6_27(289-300)Online publication date: 30-Aug-2021
    • (2019)Extended BACOLIACM Transactions on Mathematical Software10.1145/330132045:1(1-19)Online publication date: 14-Mar-2019
    • (2019)Error control Gaussian collocation software for boundary value ODEs and 1D time-dependent PDEsNumerical Algorithms10.1007/s11075-019-00738-281:4(1505-1519)Online publication date: 1-Aug-2019
    • (2018)Performance analysis of error-control B-spline Gaussian collocation software for PDEsComputers & Mathematics with Applications10.1016/j.camwa.2018.11.025Online publication date: Dec-2018

    View Options

    Get Access

    Login options

    Full Access

    View options

    PDF

    View or Download as a PDF file.

    PDF

    eReader

    View online with eReader.

    eReader

    Media

    Figures

    Other

    Tables

    Share

    Share

    Share this Publication link

    Share on social media