From f8e6baf4cc4dc4529eeda076556096e63da2dec8 Mon Sep 17 00:00:00 2001
From: Miha <mihaf@hapi.(none)>
Date: Wed, 12 Feb 2014 21:20:57 +0000
Subject: [PATCH] Just added some fprint info during run.
---
src/energy.c | 12 ++++++++++--
1 files changed, 10 insertions(+), 2 deletions(-)
diff --git a/src/energy.c b/src/energy.c
index 4d93753..1e37e08 100644
--- a/src/energy.c
+++ b/src/energy.c
@@ -20,6 +20,12 @@
}
+inline ts_bool bond_energy(ts_bond *bond,ts_poly *poly){
+ bond->energy=poly->k*pow(bond->bond_length-1,2);
+ return TS_SUCCESS;
+};
+
+
inline ts_bool energy_vertex(ts_vertex *vtx){
// ts_vertex *vtx=&vlist->vertex[n]-1; // Caution! 0 Indexed value!
// ts_triangle *tristar=vtx->tristar-1;
@@ -28,7 +34,7 @@
ts_uint jjp,jjm;
ts_vertex *j,*jp, *jm;
ts_triangle *jt;
- ts_double s=0,xh=0,yh=0,zh=0,txn=0,tyn=0,tzn=0;
+ ts_double s=0.0,xh=0.0,yh=0.0,zh=0.0,txn=0.0,tyn=0.0,tzn=0.0;
ts_double x1,x2,x3,ctp,ctm,tot,xlen;
ts_double h,ht;
for(jj=1; jj<=vtx->neigh_no;jj++){
@@ -72,7 +78,9 @@
#endif
tot=ctp+ctm;
tot=0.5*tot;
+
xlen=vtx_distance_sq(j,vtx);
+/*
#ifdef TS_DOUBLE_DOUBLE
vtx->bond[jj-1]->bond_length=sqrt(xlen);
#endif
@@ -84,7 +92,7 @@
#endif
vtx->bond[jj-1]->bond_length_dual=tot*vtx->bond[jj-1]->bond_length;
-
+*/
s+=tot*xlen;
xh+=tot*(j->x - vtx->x);
yh+=tot*(j->y - vtx->y);
--
Gitblit v1.9.3