Trisurf Monte Carlo simulator
Samo Penic
2014-04-28 2b14da2e8f012d7b1bf669af085ab75d87f378f1
Before adding constant volume subrutines
2 files modified
7 ■■■■ changed files
src/main.c 3 ●●●● patch | view | raw | blame | history
src/vertexmove.c 4 ●●● patch | view | raw | blame | history
src/main.c
@@ -1,5 +1,6 @@
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include "general.h"
#include "vertex.h"
#include "bond.h"
@@ -51,7 +52,7 @@
        vesicle->pressure=tape->pressure;
        vesicle->dmax=tape->dmax*tape->dmax;
        poly_assign_filament_xi(vesicle,tape);
        tape_free(vesicle->tape);
        free(vesicle->tape);
        vesicle->tape=tape;
        vesicle->clist->dmin_interspecies = tape->dmin_interspecies*tape->dmin_interspecies;
        /* spherical harmonics */
src/vertexmove.c
@@ -145,7 +145,7 @@
    return TS_FAIL; 
    }
}
    //accepted
//    oldcellidx=vertex_self_avoidance(vesicle, &backupvtx[0]);
    if(vtx->cell!=vesicle->clist->cell[cellidx]){
        retval=cell_add_vertex(vesicle->clist->cell[cellidx],vtx);
@@ -153,6 +153,8 @@
        if(retval==TS_SUCCESS) cell_remove_vertex(backupvtx[0].cell,vtx);
        
    }
    if(vesicle->tape->constvolswitch == 1);
//    if(oldcellidx);
    //END MONTE CARLOOOOOOO
    return TS_SUCCESS;