Fix in dipyramid initialization, where we had 0.25*25 as xk

Samo Penic

2012-07-09 443ba1603235e6cf668d0e24aa794c95941bc597

Fix in dipyramid initialization, where we had 0.25*25 as xk

 src/initial_distribution.c

@@ -20,7 +20,7 @@
    ts_vesicle *vesicle=init_vesicle(no_vertices,ncmax1,ncmax2,ncmax3,stepsize);

//TODO: debugging only. Please remove ASAP!
    vesicle->bending_rigidity=25.0*25.0;
    vesicle->bending_rigidity=25.0;

    vesicle->nshell=nshell;
    retval = vtx_set_global_values(vesicle);

 src/sh.c

@@ -371,7 +371,7 @@
ts_int i,j;
for(i=0;i<sph->l;i++){
    for(j=0;j<2*i+1;j++){
    sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j];
        sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j];
    }
}
    sph->N++;

 src/spherical_trisurf.c

@@ -43,7 +43,7 @@



for(i=0;i<10;i++){
for(i=0;i<10000;i++){
    centermass(vesicle);
    cell_occupation(vesicle);
    for(j=0;j<1000;j++){
@@ -89,6 +89,7 @@
        fprintf(fh,"%d,\t%d,\t%e\n", i, j-i, sph->sumUlm2[i][j]/(ts_double)sph->N);

            }
    fprintf(fh,"\n");
    }
    fclose(fh);
    return TS_SUCCESS;