Trisurf Monte Carlo simulator
Samo Penic
2012-07-09 443ba1603235e6cf668d0e24aa794c95941bc597
Fix in dipyramid initialization, where we had 0.25*25 as xk
3 files modified
5 ■■■■■ changed files
src/initial_distribution.c 2 ●●● patch | view | raw | blame | history
src/sh.c patch | view | raw | blame | history
src/spherical_trisurf.c 3 ●●●● patch | view | raw | blame | history
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
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;