Trisurf Monte Carlo simulator
Samo Penic
2016-06-21 e5858ff1dcf5cecfcfe81602f52d255de9d508f4
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;
}