It works

Samo Penic

2013-11-30 bc4e1ee61725fdf7ec81a9583b7fa545f847e07b

It works

 src/tape

@@ -21,7 +21,7 @@
#how many initial mcsweeps*inititer MC sweeps before recording to disk?
inititer=100
#how many records do you want on the disk iteration are there in a run?
iterations=1000
iterations=10


#shut up if we are using cluster!!!

 src/timestep.c

@@ -30,7 +30,7 @@
ts_bool single_timestep(ts_vesicle *vesicle){
    ts_bool retval;
    ts_double rnvec[3];
    ts_uint i, b;
    ts_uint i; //b;
    for(i=0;i<vesicle->vlist->n;i++){
        rnvec[0]=drand48();
        rnvec[1]=drand48();
@@ -45,10 +45,10 @@
        rnvec[1]=drand48();
        rnvec[2]=drand48();
*/ 
    b=rand() % vesicle->blist->n;
//    b=rand() % vesicle->blist->n;
        //find a bond and return a pointer to a bond...
        //call single_bondflip_timestep...
        retval=single_bondflip_timestep(vesicle,vesicle->blist->bond[b],rnvec);
   //     retval=single_bondflip_timestep(vesicle,vesicle->blist->bond[b],rnvec);
//    if(retval==TS_SUCCESS) cnt++;        
    } 
//    printf("Bondflip success rate in one sweep: %d/%d=%e\n", cnt,vesicle->blist->n,(double)cnt/(double)vesicle->blist->n);

 src/vertex.c

@@ -99,7 +99,7 @@
    return TS_SUCCESS;
}


// Add neighbour just in one direction (add vtx2 as a neigh of vtx1 and not another way around!
ts_bool vtx_add_neighbour(ts_vertex *vtx1, ts_vertex *vtx2){
    ts_uint i;
    if(vtx1==NULL || vtx2==NULL) return TS_FAIL;