This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: New extension for Associated Legendre Functions (ALFs)
- From: Brian Gough <bjg at gnu dot org>
- To: Patrick Alken <patrick dot alken at colorado dot edu>
- Cc: gsl-discuss at sourceware dot org
- Date: Tue, 17 Mar 2009 11:28:52 +0000
- Subject: Re: New extension for Associated Legendre Functions (ALFs)
- References: <20090316213137.GA4277@hippogriff.homeunix.org>
At Mon, 16 Mar 2009 15:31:37 -0600,
Patrick Alken wrote:
> Most current implementations of associated Legendre functions
> (ALFs), including the current GSL version fail due to underflow
> above degree and order (L,M) 1900. A 2002 paper by Holmes and
> Featherstone shows that with a very simple scaling modification
> this limit can be extended to 2700 in double precision. I have
> recently needed to compute high degree ALFs so I have implemented
> their algorithm and made a new GSL extension if others are interested.
Great work!
Incidentally, I found a copy of their paper online at
http://www.cage.curtin.edu.au/~geogrp/research.html
http://www.cage.curtin.edu.au/~will/20760279.pdf
--
Brian Gough