aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorFranz Zotter <fzotter@users.sourceforge.net>2009-01-12 19:14:05 +0000
committerFranz Zotter <fzotter@users.sourceforge.net>2009-01-12 19:14:05 +0000
commitca0d6cb8ef5e7ad9ada495fb38892e2d6562ec81 (patch)
tree35241d566124517e802972961b2b6199ee4e11ea /src
parent55253e4e15b093a0826fc10b376f7724cea2ba7e (diff)
corrected some misbehavior of the spherical harmoncis mathematical functions
svn path=/trunk/externals/iem/iemmatrix/; revision=10518
Diffstat (limited to 'src')
-rw-r--r--src/mtx_sh/legendre_a.c11
-rw-r--r--src/mtx_sh/sharmonics.c10
2 files changed, 11 insertions, 10 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;
diff --git a/src/mtx_sh/sharmonics.c b/src/mtx_sh/sharmonics.c
index 57f7367..8593182 100644
--- a/src/mtx_sh/sharmonics.c
+++ b/src/mtx_sh/sharmonics.c
@@ -56,7 +56,15 @@ static void sharmonics_multcheby12(SHWorkSpace *ws) {
const int yincr=(ws->nmax+1)*(ws->nmax+1);
for (n=0,ny0=0; n<=ws->nmax; n++) {
- for (m=0; m<=n; m++) {
+ m=0;
+ ly0=0;
+ lt0=nt0;
+ for (l=0; l<ws->l; l++) {
+ ws->y[ly0+ny0+m]*= ws->wc->t[lt0+m];
+ ly0+=yincr;
+ lt0+=tincr;
+ }
+ for (m=1; m<=n; m++) {
ly0=0;
lt0=nt0;
for (l=0; l<ws->l; l++) {