|
From: | Ben Abbott |
Subject: | Re: normalized ALF (Assotiated Legendre Function) |
Date: | Tue, 12 Feb 2008 20:42:35 -0500 |
The method I used detects an underflow for orders greater than 151. octave:104> result_octave1 = legendre (152, [-1:0.1:1]); warning: legendre: results may be unstable for high orders. I've attached the script. Ben
legendre.m
Description: Binary data
On Feb 12, 2008, at 1:38 PM, Ben Abbott wrote:
On Tuesday, February 12, 2008, at 10:02AM, "Marco Caliari" <address@hidden > wrote:Hi.The normalized Lagrange functions should (almost) never give NaNs of Infs.The enclosed script fixes a problem in the previous: octave:1> legendreold(151,-0.9,"norm")(end-1:end) ans = -3.3248e-53 Inf octave:2> legendre(151,-0.9,"norm")(end-1:end) ans = -3.3248e-53 9.2660e-55 MarcoI've tried to combine the various implementations. I still need to look into adding a warning in then event of an overflow.We should also add more tests to include the current improvements. Ben <legendre.m>_______________________________________________ Help-octave mailing list address@hidden https://www.cae.wisc.edu/mailman/listinfo/help-octave
[Prev in Thread] | Current Thread | [Next in Thread] |