sa_colours_read_graph.cc

Go to the documentation of this file.

#include <iostream>
#include <lemon/smart_graph.h>
#include <utils/graph_utils.hh>
#include <minlca/optimisation-sa-colours.hh>
#include <utils/timer.hh>
#include <minlca/profiler.hh>
#include <minlca/greedy.hh>
#include <minlca/optimisation-profiler.hh>
using namespace std;
using namespace lemon;
using namespace minlca;
using namespace minlca::utils;


int main(int argc, char **argv)
{
  GRAPH_TYPEDEFS(SmartGraph);
  SmartGraph g;
  ifstream graphFile(argv[1]);
  GraphFormat format = RMF;

  if (string(argv[2]) == "GRA") {
    format = GRA;
  }

  Timer timer;
  readGraph(g, format, graphFile);

  graphFile.close();

  cout << "Graph reading time: " << timer.elapsed() << endl;
  cout << endl;

  int n = countNodes(g);
  cout << "Total nodes: " << n << endl;
  cout << "Total edges: " << countEdges(g) << endl;

  cout << endl;
  cout << "Using greedy nearest:" << endl;

  MinLCAProfiler<MinLCAGreedyNearestBfs<SmartGraph>> greedyNearest(g);
  cout << "Creation elapsed time: " << greedyNearest.creationTime() << endl;
  greedyNearest.init();
  cout << "Initialisation elapsed time: " << greedyNearest.initTime() << endl;

  if (greedyNearest.run() == SOLUTION_FOUND) {
    cout << "Algorithm running time: " << greedyNearest.runTime() << endl;
    cout << "Max colours: " << greedyNearest.maxK() << endl;
    cout << "Colours used: " << greedyNearest.totalColours() << endl;
    cout << "LCA value: " << greedyNearest.lcaValue() << endl;
  }
  else {
    cout << "Could not find a solution!" << endl;
  }

  cout << endl;
  cout << "Using simulated annealing:" << endl;

  MinLCAOptProfiler<MinLCAColoursSA<SmartGraph>>
      coloursSA(g, 0, atoi(argv[3]));
  cout << "Creation elapsed time: " << coloursSA.creationTime() << endl;
  coloursSA.init(atoi(argv[4]), atoi(argv[5]), atoi(argv[6]), atof(argv[7]));
  cout << "Initialisation elapsed time: " << coloursSA.initTime() << endl;

  coloursSA.initialSolution(greedyNearest.colouring());
  cout << "Initial solution elapsed time: " << coloursSA.initialSolutionTime() <<
       endl;

  if (coloursSA.run() == SOLUTION_FOUND) {
    cout << "Algorithm running time: " << coloursSA.runTime() << endl;
    cout << "Max colours: " << coloursSA.maxK() << endl;
    cout << "Colours used: " << coloursSA.totalColours() << endl;
    cout << "LCA value: " << coloursSA.lcaValue() << endl;
  }
  else {
    cout << "Could not find a solution!" << endl;
  }
}