X-Git-Url: https://fleuret.org/cgi-bin/gitweb/gitweb.cgi?a=blobdiff_plain;f=mtp_graph.h;h=2bc59336cad9bb7cf1aa25e063ebfef6fb6889d3;hb=7f847d18f82e92c8af842f50e0ea16fa66bf9584;hp=b07e4618c62e44dd9c120ad261667e282898c71d;hpb=90a6bc97162d199e3b840bf178f49854e5df91fa;p=mtp.git

diff --git a/mtp_graph.h b/mtp_graph.h
index b07e461..2bc5933 100644
--- a/mtp_graph.h
+++ b/mtp_graph.h
@@ -37,22 +37,39 @@ class Vertex;
 class Edge;
 
 class MTPGraph {
+  // Uses the estimated vertex distances to the source to make all the
+  // edge lengths positive, resulting in an identical added value to
+  // all the paths from the same initial node to the same final node
+  // (in particular from source to sink)
   void update_positivized_lengths();
+
   // It may happen that numerical errors in update_positivized_lengths
   // make the resulting lengths negative, albeit very small. The
-  // following method force such length to zero, and print the
-  // correction when compiled in VERBOSE mode.
+  // following method forces all negative lengths to zero, and prints
+  // the total correction when compiled in VERBOSE mode.
   void force_positivized_lengths();
-  // Set the edge pred_edge_toward_source correspondingly to the path
-  // of shortest length. The current implementation is not Dijkstra's!
+
+  void decrease_distance_in_heap(Vertex *v);
+  void increase_distance_in_heap(Vertex *v);
+
+  // Visit the vertices according to _dp_order and simply update their
+  // distance to the source
+  void dp_distance_propagation();
+
+  // Set in every vertex pred_edge_toward_source correspondingly to
+  // the path of shortest length. The current implementation is
+  // Dijkstra with a Binary Heap (and not with Fibonnaci heap (yet))
   void find_shortest_path();
-  // Follows the path starting on edge e and returns its length. If
-  // nodes is non-null, stores in it the nodes met along the path.
+
+  // Follows the path starting on edge e and returns the number of
+  // nodes to reach the sink. If path is non-null, stores in it the
+  // nodes met along the path, and computes path->length properly.
   int retrieve_one_path(Edge *e, Path *path);
-  // Returns if the graph is a DAG
-  int is_dag();
 
-  Vertex **_front, **_new_front;
+  // Returns if the graph is a DAG, and set the last_change field of
+  // each vertex to the maximum number of iterations required to reach
+  // it
+  int check_DAG_and_set_last_change();
 
   int _nb_vertices, _nb_edges;
   Vertex *_source, *_sink;
@@ -60,6 +77,13 @@ class MTPGraph {
   Edge *_edges;
   Vertex *_vertices;
 
+  // For the shortest path search
+
+  Vertex **_heap;
+  int _heap_size;
+
+  Vertex **_dp_order;
+
 public:
 
   // These variables are filled when retrieve_disjoint_paths is called
@@ -71,7 +95,12 @@ public:
 
   ~MTPGraph();
 
+  // Compute the family of path with minimum total length, set the
+  // edge occupied fields accordingly.
   void find_best_paths(scalar_t *lengths);
+
+  // Retrieve the path corresponding to the occupied edges, and save
+  // the result in the nb_paths and paths fields.
   void retrieve_disjoint_paths();
 
   void print(ostream *os);