diff options
| author | Franz Zotter <fzotter@users.sourceforge.net> | 2009-01-12 19:14:05 +0000 |
|---|---|---|
| committer | Franz Zotter <fzotter@users.sourceforge.net> | 2009-01-12 19:14:05 +0000 |
| commit | ca0d6cb8ef5e7ad9ada495fb38892e2d6562ec81 (patch) | |
| tree | 35241d566124517e802972961b2b6199ee4e11ea /src/mtx_sh/legendre_a.c | |
| parent | 55253e4e15b093a0826fc10b376f7724cea2ba7e (diff) | |
corrected some misbehavior of the spherical harmoncis mathematical functions
svn path=/trunk/externals/iem/iemmatrix/; revision=10518
Diffstat (limited to 'src/mtx_sh/legendre_a.c')
| -rw-r--r-- | src/mtx_sh/legendre_a.c | 11 |
1 files changed, 2 insertions, 9 deletions
diff --git a/src/mtx_sh/legendre_a.c b/src/mtx_sh/legendre_a.c index 687cb41..d4e53f4 100644 --- a/src/mtx_sh/legendre_a.c +++ b/src/mtx_sh/legendre_a.c @@ -25,15 +25,8 @@ static void legendre_first_recurrence (double *sintheta, LegendreWorkSpace *wl) // computes the legendre functions P_n^m(costheta) for m=n // from P_0^0 for (n=1; n<=wl->nmax; n++) { - if ((n%2)==1) { - for (l=0,l0=0; l<wl->l; l++,l0+=incr) { - wl->p[l0+n0+n] = -(2*n-1) * wl->p[l0+nmo0+n-1] * sintheta[l]; - } - } - else { - for (l=0,l0=0; l<wl->l; l++,l0+=incr) { - wl->p[l0+n0+n] = -(2*n-1) * wl->p[l0+nmo0+n-1]; - } + for (l=0,l0=0; l<wl->l; l++,l0+=incr) { + wl->p[l0+n0+n] = -(2*n-1) * wl->p[l0+nmo0+n-1] * sintheta[l]; } nmo0=n0; n0+=n+1; |