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 ++++--
 1 files changed, 4 insertions(+), 2 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);

--
Gitblit v1.9.3