Confinement planes, not debugged

Samo Penic

2018-09-23 267cf14b1ced6fce480e4292131a987a1802e28a

 src/vertexmove.c

@@ -99,7 +99,13 @@
        return TS_FAIL;
    } 
   
 
// plane confinement test
   if(vesicle->tape->plane_confinement_switch){
      if(vtx->z >= vesicle->confinement_plane.z_max || vtx->z <= vesicle->confinement_plane.z_min){
            vtx=memcpy((void *)vtx,(void *)&backupvtx[0],sizeof(ts_vertex));
            return TS_FAIL;
      }
   }
//if all the tests are successful, then energy for vtx and neighbours is calculated
   for(i=0;i<vtx->neigh_no;i++){
   memcpy((void *)&backupvtx[i+1],(void *)vtx->neigh[i],sizeof(ts_vertex));
@@ -201,6 +207,16 @@
//   fprintf(stderr, "DE=%f\n",delta_energy);
    //MONTE CARLOOOOOOOO
//   if(vtx->c!=0.0) printf("DE=%f\n",delta_energy);
// plane confinement
   if(vesicle->tape->plane_confinement_switch){
      //if planes are not close enough, then repusion force is on
      if(vesicle->confinement_plane.z_max-vesicle->confinement_plane.z_min > vesicle->tape->plane_d){	
         delta_energy+=vesicle->tape->plane_F * 1.0/( (backupvtx->z-vesicle->confinement_plane.z_min) + (backupvtx->z-vesicle->confinement_plane.z_max) );
         delta_energy+=-(vesicle->tape->plane_F * 1.0/( (vtx->z-vesicle->confinement_plane.z_min) + (vtx->z-vesicle->confinement_plane.z_max) ) );

      }
   }
// end plane confinement
    if(delta_energy>=0){
#ifdef TS_DOUBLE_DOUBLE
        if(exp(-delta_energy)< drand48())