restoring vertices with spontaneous curvatures seems to work, however, the ...

Samo Penic

2016-07-03 def8b5262e183297ccee42ebbf9ebdb89c10f44c

 src/energy.c

@@ -185,3 +185,25 @@

    return TS_SUCCESS;
}



ts_bool sweep_attraction_bond_energy(ts_vesicle *vesicle){
   int i;
   for(i=0;i<vesicle->blist->n;i++){
      attraction_bond_energy(vesicle->blist->bond[i], vesicle->tape->w);
   }
   return TS_SUCCESS;
}


inline ts_bool attraction_bond_energy(ts_bond *bond, ts_double w){

   if(fabs(bond->vtx1->c)>1e-16 && fabs(bond->vtx2->c)>1e-16){
      bond->energy=-w;
   }
   else {
      bond->energy=0.0;
   }
   return TS_SUCCESS;
}