From: Brian Gough <bjg@gnu.org>
To: Patrick Alken <patrick.alken@colorado.edu>
Cc: gsl-discuss@sourceware.org
Subject: Re: New extension for Associated Legendre Functions (ALFs)
Date: Tue, 17 Mar 2009 11:37:00 -0000 [thread overview]
Message-ID: <m3wsaojsln.wl%bjg@network-theory.co.uk> (raw)
In-Reply-To: <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
next prev parent reply other threads:[~2009-03-17 11:37 UTC|newest]
Thread overview: 3+ messages / expand[flat|nested] mbox.gz Atom feed top
2009-03-16 21:29 Patrick Alken
2009-03-17 11:37 ` Brian Gough [this message]
2009-03-17 16:23 ` Patrick Alken
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=m3wsaojsln.wl%bjg@network-theory.co.uk \
--to=bjg@gnu.org \
--cc=gsl-discuss@sourceware.org \
--cc=patrick.alken@colorado.edu \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).