| | |
| | | ts_bool single_timestep(ts_vesicle *vesicle,ts_double *vmsr, ts_double *bfsr){ |
| | | ts_bool retval; |
| | | ts_double rnvec[3]; |
| | | ts_uint i,j,b; |
| | | ts_uint i,j;//b; |
| | | ts_uint vmsrcnt=0; |
| | | for(i=0;i<vesicle->vlist->n;i++){ |
| | | rnvec[0]=drand48(); |
| | |
| | | |
| | | ts_int bfsrcnt=0; |
| | | for(i=0;i<3*vesicle->vlist->n;i++){ |
| | | 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); |
| | | retval=TS_FAIL; |
| | | if(retval==TS_SUCCESS) bfsrcnt++; |
| | | } |
| | | |