[CIG-MC] CitcomS Legendre polynomial normalization
Magali Billen
mibillen at ucdavis.edu
Wed Mar 3 17:26:21 PST 2010
Hello Shijie,
Yes, its the factors in the recursion relationship that we were
surprised to find, as the rest the
recursion relationship steps are the same as that in Numerical
recipes...I'm glad to know that this
is stable for large L as we'll be going up to about L=360. Any chance
that there's a reference
for the recursion relationship (or maybe point us in the right
direction, we've looked up several others
and haven't found the same thing)? No hurry, but it seems like this is
a good thing to have reference for general
CitcomS background information too.
Thanks for answering our e-mail so quickly,
Magali
On Mar 3, 2010, at 5:13 PM, Shijie Zhong wrote:
> You are right that the recursion relation for the Legendre
> polynomials used in
> CitcomS is different from the standard in Numerical Recipes by
> Press. The
> standard one in Numerical Recipes blows up for spherical harmonic
> degree l >
> 20 or something, even with double precision computation. I re-
> derived a
> recursion relation in 1997 that is stable for nearly all l's and is
> used in CitcomS.
-----------------------------
Associate Professor, U.C. Davis
Department of Geology/KeckCAVEs
Earth & Physical Sciences Bldg, rm 2129
Davis, CA 95616
-----------------
mibillen at ucdavis.edu
(530) 754-5696
--------------------------
** Note new e-mail, building, office
information as of Sept. 2009 **
-----------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-mc/attachments/20100303/a4553430/attachment.htm
More information about the CIG-MC
mailing list