From: Ian Jackson <ian@davenant.greenend.org.uk>
Date: Sun, 6 Jan 2008 21:50:07 +0000 (+0000)
Subject: do not actually compute vertex areas
X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?p=moebius2.git;a=commitdiff_plain;h=93490dcad45ee8fbcf8e2e5c5fcb59ad6f69f01a

do not actually compute vertex areas
---

diff --git a/energy.c b/energy.c
index 4723e0c..e9b93fd 100644
--- 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++;
     }