Maybe found a bug... Lash sh was not initialized to 0 after each calculation

Samo Penic

2012-07-09 8e58b474db332bf54c3e91684d8ac8d3f14ef11b

Maybe found a bug... Lash sh was not initialized to 0 after each calculation

 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);

 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];
    }
}