Changeset 100:f1de2ab64e1c in lemon-0.x for src/work/edmonds_karp.hh

Timestamp:

02/18/04 18:27:13 (20 years ago)

Author:

marci

Branch:

default

Phase:

public

Convert:

svn:c9d7d8f5-90d6-0310-b91f-818b3a526b0e/lemon/trunk@129

Message:

.

File:

• src/work/edmonds_karp.hh (modified) (2 diffs)

src/work/edmonds_karp.hh

@@ -7 +7 @@
 
 #include <bfs_iterator.hh>
+#include <time_measure.h>
 
 namespace marci {
@@ -444 +445 @@
       return _augment;
     }
+    bool augmentWithBlockingFlow() {
+      BfsIterator4< Graph, OutEdgeIt, typename Graph::NodeMap<bool> > bfs(G);
+      bfs.pushAndSetReached(s);
+      typename Graph::NodeMap<int> dist(G); //filled up with 0's
+      while ( !bfs.finished() ) { 
+        OutEdgeIt e=OutEdgeIt(bfs);
+        if (e.valid() && bfs.isBNodeNewlyReached()) {
+          dist.set(G.head(e), dist.get(G.tail(e))+1);
+          //NodeIt v=res_graph.tail(e);
+          //NodeIt w=res_graph.head(e);
+          //pred.set(w, e);
+          //if (pred.get(v).valid()) {
+          //  free.set(w, std::min(free.get(v), e.free()));
+          //} else {
+          //  free.set(w, e.free()); 
+          //}
+          //if (res_graph.head(e)==t) { _augment=true; break; }
+        }
+        
+        ++bfs;
+      } //end of searching augmenting path
+
+      double pre_time_copy=currTime();
+      typedef Graph MutableGraph;
+      MutableGraph F;
+      typename Graph::NodeMap<NodeIt> G_to_F(G);
+      for(typename Graph::EachNodeIt n=G.template first<typename Graph::EachNodeIt>(); n.valid(); ++n) {
+        G_to_F.set(n, F.addNode());
+      }
+      for(typename Graph::EachEdgeIt e=G.template first<typename Graph::EachEdgeIt>(); e.valid(); ++e) {
+        if (dist.get(G.head(e))==dist.get(G.tail(e))+1) {
+          F.addEdge(G_to_F.get(G.tail(e)), G_to_F.get(G.head(e)));
+        }
+      }
+      double post_time_copy=currTime();
+      std::cout << "copy time: " << post_time_copy-pre_time_copy << " sec"<< std::endl; 
+
+      return 0;
+    }
     void run() {
-      while (augment()) { } 
+      //int num_of_augmentations=0;
+      while (augment()) { 
+        //std::cout << ++num_of_augmentations << " ";
+        //std::cout<<std::endl;
+      } 
     }
     Number flowValue() {