From 48bb92b75443c452e2a873d01be92bfe419e26bc Mon Sep 17 00:00:00 2001
From: mihaf <miha.fosnaric@gmail.com>
Date: Thu, 06 Feb 2014 12:04:07 +0000
Subject: [PATCH] Adding polymers MC steps (not done yet)...
---
src/poly.c | 13 +++++++++++++
1 files changed, 13 insertions(+), 0 deletions(-)
diff --git a/src/poly.c b/src/poly.c
index 05566d0..edab18f 100644
--- a/src/poly.c
+++ b/src/poly.c
@@ -5,6 +5,15 @@
#include"bond.h"
#include<math.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 +28,10 @@
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));
+ }
+
return poly;
}
--
Gitblit v1.9.3