| | |
| | | #include<stdio.h> |
| | | #include<string.h> |
| | | #include "constvol.h" |
| | | #include "globals.h" |
| | | |
| | | ts_bool single_bondflip_timestep(ts_vesicle *vesicle, ts_bond *bond, ts_double *rn){ |
| | | /*c Vertex and triangle (lm and lp) indexing for bond flip: |
| | |
| | | ts_vertex *k=bond->vtx2; |
| | | ts_uint nei,neip,neim; |
| | | ts_uint i,j; |
| | | ts_double oldenergy, delta_energy, dvol=0.0, darea=0.0; |
| | | ts_double oldenergy=0.0, delta_energy=0.0, dvol=0.0, darea=0.0; |
| | | ts_triangle *lm=NULL,*lp=NULL, *lp1=NULL, *lm2=NULL; |
| | | |
| | | ts_vertex *kp,*km; |
| | |
| | | /* end backup vertex */ |
| | | |
| | | /* Save old energy */ |
| | | oldenergy=0; |
| | | oldenergy=0.0; |
| | | oldenergy+=k->xk* k->energy; |
| | | oldenergy+=kp->xk* kp->energy; |
| | | oldenergy+=km->xk* km->energy; |
| | |
| | | |
| | | |
| | | /* Calculating the new energy */ |
| | | delta_energy=0; |
| | | delta_energy=0.0; |
| | | delta_energy+=k->xk* k->energy; |
| | | delta_energy+=kp->xk* kp->energy; |
| | | delta_energy+=km->xk* km->energy; |