Polymers thermicaly move now! HArmonic potential added.

Miha

2014-02-06 fedf2b217113800a15f367e111e2dbcad4a3c92c

 src/poly.c

@@ -4,7 +4,17 @@
#include"vertex.h"
#include"bond.h"
#include<math.h>
#include"energy.h"

ts_bool poly_assign_spring_const(ts_vesicle *vesicle){
   ts_uint i;

   for(i=0;i<vesicle->poly_list->n;i++){
    vesicle->poly_list->poly[i]->k = vesicle->spring_constant;
       }
	
   return TS_SUCCESS;
}

ts_poly   *init_poly(ts_uint n, ts_vertex *grafted_vtx){
   ts_poly   *poly=(ts_poly *)calloc(1,sizeof(ts_poly));
@@ -19,6 +29,11 @@
      vtx_add_neighbour(poly->vlist->vtx[i+1], poly->vlist->vtx[i]);
   }

   for(i=0;i<poly->blist->n;i++){
   poly->blist->bond[i]->bond_length=sqrt(vtx_distance_sq(poly->blist->bond[i]->vtx1,poly->blist->bond[i]->vtx2));
   bond_energy(poly->blist->bond[i],poly);
   }

   return poly;
}