Constant volume and constant area fix. Removed spherical harmonics calculat...

Samo Penic

2019-09-25 d5d78a49c6b51fb1f6f0661d063034eaf65e02f7

 src/bondflip.c

@@ -176,7 +176,6 @@
/*    vesicle_volume(vesicle);
    fprintf(stderr,"Volume in the beginning=%1.16e\n", vesicle->volume);
*/

/* fix data structure for flipped bond */
    ts_flip_bond(k,it,km,kp, bond,lm, lp, lm2, lp1, vesicle->tape->w);

@@ -189,6 +188,12 @@
  delta_energy+=it->xk* it->energy;
  delta_energy+=bond->energy; /* attraction with neighboring vertices, that have spontaneous curvature */
  //Neigbours of k, it, km, kp don't change its energy.
   if(vesicle->tape->stretchswitch==1){
      oldenergy+=lm->energy+lp->energy;
      stretchenergy(vesicle,lm);
      stretchenergy(vesicle,lp);
      delta_energy+=lm->energy+lp->energy;
   }

    delta_energy-=oldenergy;
   if(vesicle->pswitch == 1 || vesicle->tape->constvolswitch>0){
@@ -334,6 +339,10 @@
//      fprintf(stderr,"Restoration complete!!!\n");
//    vesicle_volume(vesicle);
//    fprintf(stderr,"Volume after fail=%1.16e\n", vesicle->volume);
   if(vesicle->tape->stretchswitch==1){
      stretchenergy(vesicle,lm);
      stretchenergy(vesicle,lp);
   }

      return TS_FAIL;
        }