Trisurf Monte Carlo simulator
Samo Penic
2014-03-06 86f5e72cbecb70197e77fb5e411874d6b015ebdd
Dump and restore complete and debugged. Still, use with caution
3 files modified
21 ■■■■■ changed files
src/io.c 6 ●●●●● patch | view | raw | blame | history
src/main.c 13 ●●●● patch | view | raw | blame | history
src/tape 2 ●●● patch | view | raw | blame | history
src/io.c
@@ -22,7 +22,7 @@
    FILE *fh=fopen("dump.bin","wb");
    /* dump vesicle */
    fwrite(vesicle, sizeof(vesicle),1,fh);
    fwrite(vesicle, sizeof(ts_vesicle),1,fh);
    /* dump vertex list */
    fwrite(vesicle->vlist, sizeof(ts_vertex_list),1,fh);
    /* dump bond list */
@@ -144,7 +144,9 @@
/* we restore all the data from the dump */
    /* restore vesicle */
    ts_vesicle *vesicle=(ts_vesicle *)calloc(1,sizeof(ts_vesicle));
    retval=fread(vesicle, sizeof(vesicle),1,fh);
    retval=fread(vesicle, sizeof(ts_vesicle),1,fh);
//    fprintf(stderr,"was here! %e\n",vesicle->dmax);
    /* restore vertex list */
    vesicle->vlist=(ts_vertex_list *)malloc(sizeof(ts_vertex_list));
    retval=fread(vesicle->vlist, sizeof(ts_vertex_list),1,fh);
src/main.c
@@ -72,18 +72,11 @@
//poly_list_free(vesicle->poly_list);
/*End testing*/
//write_vertex_xml_file(vesicle,0);
//vesicle_free(vesicle);
dump_state(vesicle1);
vesicle=restore_state();
//write_vertex_xml_file(vesicle,1);
//write_master_xml_file("test.pvd");
//vesicle_free(vesicle1);
//vesicle_free(vesicle);
vesicle_free(vesicle1);
fprintf(stderr,"Starting simulation...");
run_simulation(vesicle, mcsweeps, inititer, iterations);
vesicle=restore_state();
//vesicle_free(vesicle1);
run_simulation(vesicle1, mcsweeps, inititer, iterations);
write_master_xml_file("test.pvd");
write_dout_fcompat_file(vesicle,"dout");
vesicle_free(vesicle);
src/tape
@@ -28,7 +28,7 @@
#how many initial mcsweeps*inititer MC sweeps before recording to disk?
inititer=1
#how many records do you want on the disk iteration are there in a run?
iterations=10
iterations=10000
#shut up if we are using cluster!!!