| | |
| | | |
| | | delta_energy=0; |
| | | |
| | | |
| | | // vesicle_volume(vesicle); |
| | | // fprintf(stderr,"Volume in the beginning=%1.16e\n", vesicle->volume); |
| | | |
| | |
| | | pow(sqrt(vtx_distance_sq(&backupvtx[0], vtx->grafted_poly->vlist->vtx[0])-1),2)) *vtx->grafted_poly->k; |
| | | } |
| | | */ |
| | | |
| | | // plane confinement energy due to compressing force |
| | | if(vesicle->tape->plane_confinement_switch){ |
| | | if(vesicle->confinement_plane.force_switch){ |
| | | //substract old energy |
| | | if(abs(vesicle->tape->plane_d/2.0-vesicle->confinement_plane.z_max)>1e-10) { |
| | | delta_energy-=vesicle->tape->plane_F / pow(vesicle->confinement_plane.z_max-backupvtx[0].z,2); |
| | | delta_energy+=vesicle->tape->plane_F / pow(vesicle->confinement_plane.z_max-vtx->z,2); |
| | | } |
| | | if(abs(-vesicle->tape->plane_d/2.0-vesicle->confinement_plane.z_min)>1e-10) { |
| | | delta_energy-=vesicle->tape->plane_F / pow(vesicle->confinement_plane.z_min-backupvtx[0].z,2); |
| | | delta_energy+=vesicle->tape->plane_F / pow(vesicle->confinement_plane.z_min-vtx->z,2); |
| | | } |
| | | } |
| | | } |
| | | |
| | | // fprintf(stderr, "DE=%f\n",delta_energy); |
| | | //MONTE CARLOOOOOOOO |
| | | // if(vtx->c!=0.0) printf("DE=%f\n",delta_energy); |