Only reaction on successful acceptance of MC step needs to be programmed.

Samo Penic

2014-04-28 958e0e46e248cc7fd408fc2c4f7eb119ed93b1e5

Only reaction on successful acceptance of MC step needs to be programmed.

 src/vertexmove.c

@@ -23,7 +23,7 @@
    ts_double delta_energy,oenergy,dvol=0.0;
    ts_double costheta,sintheta,phi,r;
    //This will hold all the information of vtx and its neighbours
    ts_vertex backupvtx[20], *constvol_vtx_moved, *constvol_vtx_backup;
    ts_vertex backupvtx[20], *constvol_vtx_moved=NULL, *constvol_vtx_backup=NULL;
    memcpy((void *)&backupvtx[0],(void *)vtx,sizeof(ts_vertex));

    //Some stupid tests for debugging cell occupation!
@@ -157,7 +157,7 @@
   for(i=0;i<vtx->tristar_no;i++) triangle_normal_vector(vtx->tristar[i]);

    if(vesicle->tape->constvolswitch == 1){
        ts_bool constvolumerestore(constvol_vtx_backup);
        constvolumerestore(constvol_vtx_moved,constvol_vtx_backup);
    }

    return TS_FAIL; 
@@ -173,7 +173,7 @@
    }

    if(vesicle->tape->constvolswitch == 1){
        ts_bool constvolumeaccept(constvol_vtx_backup);
        constvolumeaccept(constvol_vtx_moved,constvol_vtx_backup);
    }
//    if(oldcellidx);
    //END MONTE CARLOOOOOOO