From 8e58b474db332bf54c3e91684d8ac8d3f14ef11b Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Mon, 09 Jul 2012 17:15:37 +0000
Subject: [PATCH] Maybe found a bug... Lash sh was not initialized to 0 after each calculation
---
src/co_test.c | 6 ++++--
src/sh.c | 3 ++-
2 files changed, 6 insertions(+), 3 deletions(-)
diff --git a/src/co_test.c b/src/co_test.c
index 54e18b7..61a6a3a 100644
--- a/src/co_test.c
+++ b/src/co_test.c
@@ -40,8 +40,10 @@
preparationSh(vesicle,r0);
calculateYlmi(vesicle);
+/*
ts_coord *coord=(ts_coord *)malloc(sizeof(ts_coord));
ts_double fi, theta;
+
for(i=0;i<vesicle->vlist->n;i++){
cart2sph(coord,vesicle->vlist->vtx[i]->x, vesicle->vlist->vtx[i]->y, vesicle->vlist->vtx[i]->z);
@@ -54,15 +56,15 @@
}
free(coord);
+*/
calculateUlm(vesicle);
-
for(i=0;i<vesicle->sphHarmonics->l;i++){
for(j=0;j<2*i+1;j++){
fprintf(stderr,"ulm(%d,%d)=%e\n",i,j+1,vesicle->sphHarmonics->ulm[i][j]);
}
}
-
+storeUlm2(vesicle);
sph_free(vesicle->sphHarmonics);
vesicle_free(vesicle);
diff --git a/src/sh.c b/src/sh.c
index baeb8f9..0c4b06d 100644
--- a/src/sh.c
+++ b/src/sh.c
@@ -343,7 +343,7 @@
ts_uint i,j,k;
ts_vertex *cvtx;
for(i=0;i<vesicle->sphHarmonics->l;i++){
- for(j=0;j<2*i;j++) vesicle->sphHarmonics->ulm[i][j]=0.0;
+ for(j=0;j<2*i+1;j++) vesicle->sphHarmonics->ulm[i][j]=0.0;
}
//TODO: call calculateYlmi !!!
@@ -372,6 +372,7 @@
ts_int i,j;
for(i=0;i<sph->l;i++){
for(j=0;j<2*i+1;j++){
+ /* DEBUG fprintf(stderr,"sph->sumUlm2[%d][%d]=%e\n",i,j,sph->ulm[i][j]* sph->ulm[i][j]); */
sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j];
}
}
--
Gitblit v1.9.3