[PD-cvs] SF.net SVN: pure-data:[10518] trunk/externals/iem/iemmatrix/src/mtx_sh

[fzotter at users.sourceforge.net]

Revision: 10518
          http://pure-data.svn.sourceforge.net/pure-data/?rev=10518&view=rev
Author:   fzotter
Date:     2009-01-12 19:14:05 +0000 (Mon, 12 Jan 2009)

Log Message:
-----------
corrected some misbehavior of the spherical harmoncis mathematical functions

Modified Paths:
--------------
    trunk/externals/iem/iemmatrix/src/mtx_sh/legendre_a.c
    trunk/externals/iem/iemmatrix/src/mtx_sh/sharmonics.c

Modified: trunk/externals/iem/iemmatrix/src/mtx_sh/legendre_a.c
===================================================================
--- trunk/externals/iem/iemmatrix/src/mtx_sh/legendre_a.c	2009-01-12 18:19:54 UTC (rev 10517)
+++ trunk/externals/iem/iemmatrix/src/mtx_sh/legendre_a.c	2009-01-12 19:14:05 UTC (rev 10518)
@@ -25,16 +25,9 @@
    // 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];
-	  }
+      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];
-	  }
-      }
       nmo0=n0;
       n0+=n+1;
    }

Modified: trunk/externals/iem/iemmatrix/src/mtx_sh/sharmonics.c
===================================================================
--- trunk/externals/iem/iemmatrix/src/mtx_sh/sharmonics.c	2009-01-12 18:19:54 UTC (rev 10517)
+++ trunk/externals/iem/iemmatrix/src/mtx_sh/sharmonics.c	2009-01-12 19:14:05 UTC (rev 10518)
@@ -56,7 +56,15 @@
    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++) {


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.