chiark / gitweb /
do not actually compute vertex areas
authorIan Jackson <ian@davenant.greenend.org.uk>
Sun, 6 Jan 2008 21:50:07 +0000 (21:50 +0000)
committerIan Jackson <ian@davenant.greenend.org.uk>
Sun, 6 Jan 2008 21:50:07 +0000 (21:50 +0000)
energy.c

index 4723e0c..e9b93fd 100644 (file)
--- a/energy.c
+++ b/energy.c
@@ -72,7 +72,8 @@ void compute_edge_lengths(const Vertices vertices) {
 }
 
 void compute_vertex_areas(const Vertices vertices) {
-  int v0,v1,v2, e1,e2, k;
+  int v0,v1,v2, e1,e2;
+//  int k;
 
   FOR_VERTEX(v0) {
     double total= 0.0, edges_total=0;
@@ -85,13 +86,13 @@ void compute_vertex_areas(const Vertices vertices) {
 
       edges_total += edge_lengths[v0][e1];
 
-      double e1v[D3], e2v[D3], av[D3];
-      K {
-       e1v[k]= vertices[v1][k] - vertices[v0][k];
-       e2v[k]= vertices[v2][k] - vertices[v0][k];
-      }
-      xprod(av, e1v, e2v);
-      total += magnD(av);
+//      double e1v[D3], e2v[D3], av[D3];
+//      K {
+//     e1v[k]= vertices[v1][k] - vertices[v0][k];
+//     e2v[k]= vertices[v2][k] - vertices[v0][k];
+//      }
+//      xprod(av, e1v, e2v);
+//      total += magnD(av);
       
       count++;
     }