A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalised associated Legendre functions

By S. Holmes, Will Featherstone. 2002

221 citations recorded by the source corpus.

How researchers see this paper

The researchers whose trusted networks rank it highest, strongest advocate first. Every researcher gets a different order.

Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies (2012, methodology): “Wenzel single normalization (1998), or Horner’s schemes (in cos φ) on partial sums such as in Holmes and Featherstone (2002),”
A computationally compact representation of Magnetic-Apex and Quasi-Dipole coordinates with smooth base vectors (2009, methodology): “Holmes and Featherstone [2002] noted that the standard forward column method, applied at double precision, will produce numerical underflow for m > 1900; we therefore expect that our algorithm will produce reliable output for orders considerably higher than m = 72.”
Gravity and Topography of the Terrestrial Planets (2015, methodology): “If starting values used in the recursion are appropriately scaled, as is summarized by Holmes and Featherstone (2002), these can be computed to high accuracy up to a maximum spherical harmonic degree of about 2700.”