From ca0d6cb8ef5e7ad9ada495fb38892e2d6562ec81 Mon Sep 17 00:00:00 2001
From: Franz Zotter <fzotter@users.sourceforge.net>
Date: Mon, 12 Jan 2009 19:14:05 +0000
Subject: corrected some misbehavior of the spherical harmoncis mathematical
 functions

svn path=/trunk/externals/iem/iemmatrix/; revision=10518
---
 src/mtx_sh/sharmonics.c | 10 +++++++++-
 1 file changed, 9 insertions(+), 1 deletion(-)

(limited to 'src/mtx_sh/sharmonics.c')

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++) {
-- 
cgit v1.2.1