experiment with other scalings for shortest path

[moebius2.git] / bgl.cpp

diff --git a/bgl.cpp b/bgl.cpp
index cbe45b52d798865b2ad1ce0b7ad09218cd689ed7..5d1e0ad2540d63f8248e2ce28813dc5a8b082a00 100644 (file)
--- a/bgl.cpp
+++ b/bgl.cpp
@@ -190,7 +190,7 @@ double graph_layout_cost(const Vertices v, const double vertex_areas[N]) {
    */
   static const double d2_epsilon= 1e-6;
   
-  double edge_weights[N*V6], vertex_distances[N], total_cost=0;
+  double edge_weights[V6<<ESHIFT], vertex_distances[N], total_cost=0;
   int v1,v2,e,f;
 
   FOR_VERTEX(v1)
@@ -200,21 +200,22 @@ double graph_layout_cost(const Vertices v, const double vertex_areas[N]) {
       edge_weights[f]= hypotD(v[v1], v[v2]);
 
   FOR_VERTEX(v1) {
-    double a1= vertex_areas[v1];
+    //double a1= vertex_areas[v1];
     single_source_shortest_paths(v1, edge_weights, vertex_distances);
     FOR_VERTEX(v2) {
       if (v1 == v2) continue;
-      double a2= vertex_areas[v2];
+      //double a2= vertex_areas[v2];
       double d2= hypotD2plus(v[v1],v[v2], d2_epsilon);
       double s= vertex_distances[v2];
       double s2= s*s + d2_epsilon;
       double sd2= s2 / d2;
-      double cost_contrib= a1*a2 * (sd2 - 1) / (d2*d2);
-      if (cost_contrib < -1e-4) {
-       printf("layout %03x..%03x (a=%g,%g) s=%g s2=%g d2=%g sd2=%g"
-              " cost+=%g\n", v1,v2, a1,a2, s,s2,d2,sd2, cost_contrib);
-       abort();
-      }
+      //double cost_contrib= a1*a2 * (sd2 - 1) / (d2*d2);
+      double cost_contrib= sd2;
+      //if (cost_contrib < -1e-4) {
+      //       printf("layout %03x..%03x (a=%g,%g) s=%g s2=%g d2=%g sd2=%g"
+      //              " cost+=%g\n", v1,v2, a1,a2, s,s2,d2,sd2, cost_contrib);
+      //       abort();
+      //      }
       total_cost += cost_contrib;
     }
   }