Show simple item record

dc.contributor.authorHolmes, S.
dc.contributor.authorFeatherstone, Will
dc.identifier.citationHolmes, S.A. and Featherstone, W.E.. 2002. A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalised associated Legendre functions. Journal of Geodesy 76 (5): 279-299.

Spherical harmonic expansions form partial sums of fully normalised associated Legendre functions (ALFs). However, when evaluated increasingly close to the poles, the ultra-high degree and order (e.g. 2700) ALFs range over thousands of orders of magnitude. This causes existing recursion techniques for computing values of individual ALFs and their derivatives to fail. A common solution in geodesy is to evaluate these expansions using Clenshaw's method, which does not compute individual ALFs or their derivatives. Straightforward numerical principles govern the stability of this technique. Elementary algebra is employed to illustrate how these principles are implemented in Clenshaw's method. It is also demonstrated how existing recursion algorithms for computing ALFs and their first derivatives are easily modified to incorporate these same numerical principles. These modified recursions yield scaled ALFs and first derivatives, which can then be combined using Horner's scheme to compute partial sums, complete to degree and order 2700, for all latitudes (except at the poles for first derivatives). This exceeds any previously published result. Numerical tests suggest that this new approach is at least as precise and efficient as Clenshaw's method. However, the principal strength of the new techniques lies in their simplicity of formulation and implementation, since this quality should simplify the task of extending the approach to other uses, such as spherical harmonic analysis.

dc.subjectSpherical Harmonic Expansions - Fully Normalised Associated Legendre Functions - Clenshaw Summation - Recursion - Horner's Scheme
dc.titleA unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalised associated Legendre functions
dc.typeJournal Article
dcterms.source.titleJournal of Geodesy

Originally published in Journal of Geodesy 2002 76(5) pp.279-299.


Copyright Springer-Verlag


The original article is available at

curtin.accessStatusOpen access
curtin.facultyDivision of Resources and Environment
curtin.facultyDepartment of Spatial Sciences

Files in this item


This item appears in the following Collection(s)

Show simple item record