Affiliation:
1. University of Bergen, Bergen, Norway
Abstract
We discuss how to modify a recently published Matlab code, coteglob, so that the excellent performance this code demonstrates for low and intermediate accuracy requests is retained while the performance is improved for high accuracy requests. coteglob is a globally adaptive code using a 5 and 9 point pair of Newton-Cotes rules. Combining an extended sequence of rules using 5, 9, 17 and 33 points with a doubly adaptive bisection strategy is the main focus of the paper. We also discuss local versus global adaptivity and conclude that globally adaptive codes are to be preferred. Based on this we develop several new globally adaptive codes that all compare favorably both with coteglob, with Matlab's best currently available quadrature software quadl and the general purpose QUADPACK codes dqk15 and dqk21. We include the results from extensive testing using both a Lyness-Kaganove testing technique and a battery test.
Publisher
Association for Computing Machinery (ACM)
Subject
Applied Mathematics,Software
Reference23 articles.
1. Berntsen J. 1986. A test of some well known quadrature routines. Reports in Informatics 20 Dept. of Informatics University of Bergen. Berntsen J. 1986. A test of some well known quadrature routines. Reports in Informatics 20 Dept. of Informatics University of Bergen.
2. Error estimation in automatic quadrature routines
3. An adaptive algorithm for the approximate calculation of multiple integrals
4. de Boor C. 1971. On writing an automatic integration algorithm. In Mathematical Software J. R. Rice Ed. Academic Press New York NY. de Boor C. 1971. On writing an automatic integration algorithm. In Mathematical Software J. R. Rice Ed. Academic Press New York NY.
Cited by
11 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献