source: lemon-benchmark/tests/circulation.cc @ 9:7768d68909e8

#include <lemon/circulation.h>
#include <lemon/smart_graph.h>
#include <lemon/lgf_reader.h>
#include <lemon/dimacs.h>
#include <lemon/time_measure.h>
#include <lemon/error.h>

const char test_name[] = "circulation";

using namespace lemon;

int main(int argc, char **argv)
{
  if(argc!=2) exit(1);

  typedef int Value;
  
  SmartDigraph g;
  
  SmartDigraph::ArcMap<Value> lo_cap(g);
  SmartDigraph::ArcMap<Value> up_cap(g);
  SmartDigraph::ArcMap<Value> cost(g);
  SmartDigraph::NodeMap<Value> sup(g);
  
  Timer ti;
  try {
    std::ifstream input;
    input.open(argv[1]);
    readDimacsMin(input, g, lo_cap, up_cap, cost, sup);
  } catch (FormatError& error) {
    std::cerr << error.what() << std::endl;
    return 1;
  }
  std::cerr << "Read the file: " << ti << '\n';
  ti.restart();
  
  Circulation<SmartDigraph,
    SmartDigraph::ArcMap<Value>,SmartDigraph::ArcMap<Value>,
    SmartDigraph::NodeMap<Value> > circ(g,lo_cap,up_cap,sup);
  std::cerr << "Setup Circulation class: " << ti << '\n';
  ti.restart();
  bool res = circ.run();
  if(res)
    {
      std::cerr << "A feasible circulation is found: " << ti << "\n";
      ti.restart();
      bool res2 = circ.checkFlow();
      std::cerr << "Checked in time " << ti << "\n";
      if(res2)
        std::cerr << "Success!\nn";
      else
        std::cerr << "Oops!!!!\n\n";
    }
  else
    {
      std::cerr << "A dual solution is found: " << ti << "\n";
      ti.restart();
      bool res2 = circ.checkBarrier();
      std::cerr << "Checked in time " << ti << "\n";
      if(res2)
        std::cerr << "Success!\nn";
      else
        std::cerr << "Dual-Oops!!!!\n\n";

    }
}