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