Fix in dipyramid initialization, where we had 0.25*25 as xk
| | |
| | | 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); |
| | |
| | | |
| | | |
| | | |
| | | for(i=0;i<10;i++){ |
| | | for(i=0;i<10000;i++){ |
| | | centermass(vesicle); |
| | | cell_occupation(vesicle); |
| | | for(j=0;j<1000;j++){ |
| | |
| | | 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; |