×

Fast computation of incomplete elliptic integral of first kind by half argument transformation. (English) Zbl 1201.65035

Summary: We develop a new method to calculate the incomplete elliptic integral of the first kind, \({F(\varphi|m)}\), by using the half argument formulas of Jacobian elliptic functions. The method reduces the magnitude of \({\varphi}\) by repeated usage of the formulas while fixing \(m\). The method is sufficiently precise in the sense that the maximum relative error is 3–5 machine epsilons at most. Thanks to the simplicity of the half argument formulas, the new procedure is significantly faster than the existing procedures. For example, it runs 20–60% faster than Bulirsch’ function, el1, and 1.9–2.2 times faster than the method using Carlson’s function, \(R _{F}\).

MSC:

65D20 Computation of special functions and constants, construction of tables
33E05 Elliptic functions and integrals
33F05 Numerical approximation and evaluation of special functions
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Abramowitz, M. Stegun, I.A. (eds.): Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Chapter 17. National Bureau of Standards, Washington (1964) · Zbl 0171.38503
[2] Bulirsch R.: Numerical computation of elliptic integrals and elliptic functions. Numer. Math. 7, 78–90 (1965a) · Zbl 0133.08702 · doi:10.1007/BF01397975
[3] Bulirsch R.: Numerical computation of elliptic integrals and elliptic functions II. Numer. Math. 7, 353–354 (1965b) · Zbl 0128.37204 · doi:10.1007/BF01436529
[4] Bulirsch R.: An extension of the Bartky-transformation to incomplete elliptic integrals of the third kind. Numer. Math. 13, 266–284 (1969a) · Zbl 0175.46001 · doi:10.1007/BF02167558
[5] Bulirsch R.: Numerical computation of elliptic integrals and elliptic functions III. Numer. Math. 13, 305–315 (1969b) · Zbl 0181.17502 · doi:10.1007/BF02165405
[6] Byrd P.F., Friedman M.D.: Handbook on Elliptic Integrals for Engineers and Physicists, 2nd edn. Springer, Berlin (1971) · Zbl 0213.16602
[7] Carlson B.C.: On computing elliptic integrals and functions. J. Math. Phys. 44, 332–345 (1965) · Zbl 0163.39601
[8] Carlson B.C.: Elliptic integrals of the first kind. SIAM J. Math. Anal. 8, 231–242 (1977) · Zbl 0356.33002 · doi:10.1137/0508016
[9] Carlson B.C.: Short proofs of three theorems on elliptic integrals. SIAM J. Math. Anal. 9, 524–528 (1978) · Zbl 0401.33001 · doi:10.1137/0509033
[10] Carlson B.C.: Computing elliptic integrals by duplication. Numer. Math. 33, 1–16 (1979) · Zbl 0438.65029 · doi:10.1007/BF01396491
[11] Carlson B.C., Notis E.M.: Algorithm 577. Algorithms for incomplete elliptic integrals. ACM Trans. Math. Software 7, 398–403 (1981) · Zbl 0464.65008 · doi:10.1145/355958.355970
[12] Cody W.J.: Chebyshev approximations for the complete elliptic integrals K and E. Math. Comp. 19, 105–112 (1965a) · Zbl 0137.33502
[13] Cody W.J.: Chebyshev polynomial expansions of complete elliptic integrals K and E. Math. Comp. 19, 249–259 (1965b) · Zbl 0134.33301 · doi:10.1090/S0025-5718-1965-0178563-0
[14] Cody W.J.: Corrigenda: Chebyshev approximations for the complete elliptic integrals K and E. Math. Comp. 20, 207 (1966) · doi:10.1090/S0025-5718-66-99936-4
[15] Didonato A.R., Hershey A.V.: New formulas for computing incomplete elliptic integrals of the first and second kind. J. Assoc. Comput. Mach. 6, 515–526 (1959) · Zbl 0117.10804
[16] Fukushima T.: Gaussian element formulation of short-axis-mode rotation of a rigid body. Astron. J. 136, 649–653 (2008a) · doi:10.1088/0004-6256/136/2/649
[17] Fukushima T.: Canonical and universal elements of rotational motion of triaxial rigid body. Astron. J. 136, 1728–1735 (2008b) · doi:10.1088/0004-6256/136/4/1728
[18] Fukushima T.: Fast computation of Jacobian elliptic functions and incomplete elliptic integrals for constant values of elliptic parameter and elliptic characteristic. Celest. Mech. Dyn. Astron. 105, 245–260 (2009a) · Zbl 1223.70004 · doi:10.1007/s10569-008-9177-y
[19] Fukushima T.: Efficient solution of initial-value problem of torque-free rotation. Astron. J. 137, 210–218 (2009b) · doi:10.1088/0004-6256/138/1/210
[20] Fukushima T.: Fast computation of complete elliptic integrals and Jacobian elliptic functions. Celest. Mech. Dyn. Astron. 105, 305–328 (2009c) · Zbl 1223.70005 · doi:10.1007/s10569-009-9228-z
[21] Fukushima T., Ishizaki H.: Numerical computation of incomplete elliptic integrals of a general form. Celest. Mech. Dyn. Astron. 59, 237–251 (1994) · Zbl 0818.33013 · doi:10.1007/BF00692874
[22] Glynn, E.F.: efg’s Computer Lab and Reference Library. http://www.efg2.com/Lab/Library/mathematics.htm (2009)
[23] Hastings C. Jr: Approximations for Digital Computers. Princeton University Press, Princeton (1955) · Zbl 0066.10704
[24] Hofsommer D.J., van de Riet R.P.: On the Numerical Calculation of Elliptic Integrals of the First and Second Kind and Elliptic Functions of Jacobi. Numer. Math. 5, 291–302 (1963) · Zbl 0123.13204 · doi:10.1007/BF01385899
[25] Morris, A.H. Jr.: NSWC Library of Mathematics Subroutines, Tech. Rep. NSWCDD/TR-92/425, 107-110. Naval Surface Warfare Center, Dahlgren (1993)
[26] Press W.H., Flannery B.P., Teukolsky S.A., Vetterling W.T.: Numerical Recipes: the Art of Scientific Computing. Cambridge University Press, Cambridge (1986) · Zbl 0587.65003
[27] Press W.H., Teukolsky S.A., Vetterling W.T., Flannery B.P.: Numerical Recipes: the Art of Scientific Computing, 3rd edn. Cambridge University Press, Cambridge (2007) · Zbl 1132.65001
[28] Vande Vel H.: On the series expansion method for computing incomplete elliptic integrals of the first and second kinds. Math. Comp. 23, 61–69 (1969) · Zbl 0187.10201 · doi:10.2307/2005054
[29] Wolfram S.: The Mathematica Book, 5th edn. Wolfram Research Inc./Cambridge University Press, Cambridge (2003)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.