From 34d3de02dc183d38bf661789e98e43d9c1b4c0bd Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Thu, 23 Feb 2012 15:13:28 +0000
Subject: [PATCH] cell data removed
---
src/energy.c | 18 ++++++++++--------
1 files changed, 10 insertions(+), 8 deletions(-)
diff --git a/src/energy.c b/src/energy.c
index e30295b..7cd96f5 100644
--- a/src/energy.c
+++ b/src/energy.c
@@ -74,24 +74,24 @@
tot=0.5*tot;
xlen=vtx_distance_sq(j,vtx);
#ifdef TS_DOUBLE_DOUBLE
- data->bond[jj-1]->data->bond_length=sqrt(xlen);
+ data->bond[jj-1]->bond_length=sqrt(xlen);
#endif
#ifdef TS_DOUBLE_FLOAT
- data->bond[jj-1]->data->bond_length=sqrtf(xlen);
+ data->bond[jj-1]->bond_length=sqrtf(xlen);
#endif
#ifdef TS_DOUBLE_LONGDOUBLE
- data->bond[jj-1]->data->bond_length=sqrtl(xlen);
+ data->bond[jj-1]->bond_length=sqrtl(xlen);
#endif
- data->bond[jj-1]->data->bond_length_dual=tot*data->bond[jj-1]->data->bond_length;
+ data->bond[jj-1]->bond_length_dual=tot*data->bond[jj-1]->bond_length;
s+=tot*xlen;
xh+=tot*(j->data->x - data->x);
yh+=tot*(j->data->y - data->y);
zh+=tot*(j->data->z - data->z);
- txn+=jt->data->xnorm;
- tyn+=jt->data->ynorm;
- tzn+=jt->data->znorm;
+ txn+=jt->xnorm;
+ tyn+=jt->ynorm;
+ tzn+=jt->znorm;
}
h=xh*xh+yh*yh+zh*zh;
@@ -118,7 +118,9 @@
data->curvature=-sqrtl(h);
}
#endif
-//TODO: MAJOR!!!! What is vtx->data->c?????????????? Here it is 0!
+// What is vtx->data->c?????????????? Here it is 0!
+// c is forced curvature energy for each vertex. Should be set to zero for
+// norman circumstances.
data->energy=0.5*s*(data->curvature/s-data->c)*(data->curvature/s-data->c);
return TS_SUCCESS;
--
Gitblit v1.9.3