Before adding constant volume subrutines
| | |
| | | #include<stdio.h> |
| | | #include<math.h> |
| | | #include<stdlib.h> |
| | | #include "general.h" |
| | | #include "vertex.h" |
| | | #include "bond.h" |
| | |
| | | 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 */ |
| | |
| | | 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); |
| | |
| | | if(retval==TS_SUCCESS) cell_remove_vertex(backupvtx[0].cell,vtx); |
| | | |
| | | } |
| | | |
| | | if(vesicle->tape->constvolswitch == 1); |
| | | // if(oldcellidx); |
| | | //END MONTE CARLOOOOOOO |
| | | return TS_SUCCESS; |