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
