/* -*- mode: C++; indent-tabs-mode: nil; -*-

  *

  * This file is a part of LEMON, a generic C++ optimization library.

  *

  * Copyright (C) 2003-2009

  * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport

  * (Egervary Research Group on Combinatorial Optimization, EGRES).

  *

  * Permission to use, modify and distribute this software is granted

  * provided that this copyright notice appears in all copies. For

  * precise terms see the accompanying LICENSE file.

  *

  * This software is provided "AS IS" with no warranty of any kind,

  * express or implied, and with no claim as to its suitability for any

  * purpose.

  *

  */

 #include <lemon/concepts/graph.h>

 #include <lemon/list_graph.h>

 #include <lemon/smart_graph.h>

 #include <lemon/full_graph.h>

 #include <lemon/grid_graph.h>

 #include <lemon/hypercube_graph.h>

 #include "test_tools.h"

 #include "graph_test.h"

 using namespace lemon;

 using namespace lemon::concepts;

 template <class Graph>

 void checkGraphBuild() {

   TEMPLATE_GRAPH_TYPEDEFS(Graph);

   Graph G;

   checkGraphNodeList(G, 0);

   checkGraphEdgeList(G, 0);

   checkGraphArcList(G, 0);

   G.reserveNode(3);

   G.reserveEdge(3);

   Node

     n1 = G.addNode(),

     n2 = G.addNode(),

     n3 = G.addNode();

   checkGraphNodeList(G, 3);

   checkGraphEdgeList(G, 0);

   checkGraphArcList(G, 0);

   Edge e1 = G.addEdge(n1, n2);

   check((G.u(e1) == n1 && G.v(e1) == n2) || (G.u(e1) == n2 && G.v(e1) == n1),

         "Wrong edge");

   checkGraphNodeList(G, 3);

   checkGraphEdgeList(G, 1);

   checkGraphArcList(G, 2);

   checkGraphIncEdgeArcLists(G, n1, 1);

   checkGraphIncEdgeArcLists(G, n2, 1);

   checkGraphIncEdgeArcLists(G, n3, 0);

   checkGraphConEdgeList(G, 1);

   checkGraphConArcList(G, 2);

   Edge e2 = G.addEdge(n2, n1),

        e3 = G.addEdge(n2, n3);

   checkGraphNodeList(G, 3);

   checkGraphEdgeList(G, 3);

   checkGraphArcList(G, 6);

   checkGraphIncEdgeArcLists(G, n1, 2);

   checkGraphIncEdgeArcLists(G, n2, 3);

   checkGraphIncEdgeArcLists(G, n3, 1);

   checkGraphConEdgeList(G, 3);

   checkGraphConArcList(G, 6);

   checkArcDirections(G);

   checkNodeIds(G);

   checkArcIds(G);

   checkEdgeIds(G);

   checkGraphNodeMap(G);

   checkGraphArcMap(G);

   checkGraphEdgeMap(G);

 }

 template <class Graph>

 void checkGraphAlter() {

   TEMPLATE_GRAPH_TYPEDEFS(Graph);

   Graph G;

   Node n1 = G.addNode(), n2 = G.addNode(),

        n3 = G.addNode(), n4 = G.addNode();

   Edge e1 = G.addEdge(n1, n2), e2 = G.addEdge(n2, n1),

        e3 = G.addEdge(n2, n3), e4 = G.addEdge(n1, n4),

        e5 = G.addEdge(n4, n3);

   checkGraphNodeList(G, 4);

   checkGraphEdgeList(G, 5);

   checkGraphArcList(G, 10);

   // Check changeU() and changeV()

   if (G.u(e2) == n2) {

     G.changeU(e2, n3);

   } else {

     G.changeV(e2, n3);

   }

   checkGraphNodeList(G, 4);

   checkGraphEdgeList(G, 5);

   checkGraphArcList(G, 10);

   checkGraphIncEdgeArcLists(G, n1, 3);

   checkGraphIncEdgeArcLists(G, n2, 2);

   checkGraphIncEdgeArcLists(G, n3, 3);

   checkGraphIncEdgeArcLists(G, n4, 2);

   checkGraphConEdgeList(G, 5);

   checkGraphConArcList(G, 10);

   if (G.u(e2) == n1) {

     G.changeU(e2, n2);

   } else {

     G.changeV(e2, n2);

   }

   checkGraphNodeList(G, 4);

   checkGraphEdgeList(G, 5);

   checkGraphArcList(G, 10);

   checkGraphIncEdgeArcLists(G, n1, 2);

   checkGraphIncEdgeArcLists(G, n2, 3);

   checkGraphIncEdgeArcLists(G, n3, 3);

   checkGraphIncEdgeArcLists(G, n4, 2);

   checkGraphConEdgeList(G, 5);

   checkGraphConArcList(G, 10);

   // Check contract()

   G.contract(n1, n4, false);

   checkGraphNodeList(G, 3);

   checkGraphEdgeList(G, 5);

   checkGraphArcList(G, 10);

   checkGraphIncEdgeArcLists(G, n1, 4);

   checkGraphIncEdgeArcLists(G, n2, 3);

   checkGraphIncEdgeArcLists(G, n3, 3);

   checkGraphConEdgeList(G, 5);

   checkGraphConArcList(G, 10);

   G.contract(n2, n3);

   checkGraphNodeList(G, 2);

   checkGraphEdgeList(G, 3);

   checkGraphArcList(G, 6);

   checkGraphIncEdgeArcLists(G, n1, 4);

   checkGraphIncEdgeArcLists(G, n2, 2);

   checkGraphConEdgeList(G, 3);

   checkGraphConArcList(G, 6);

 }

 template <class Graph>

 void checkGraphErase() {

   TEMPLATE_GRAPH_TYPEDEFS(Graph);

   Graph G;

   Node n1 = G.addNode(), n2 = G.addNode(),

        n3 = G.addNode(), n4 = G.addNode();

   Edge e1 = G.addEdge(n1, n2), e2 = G.addEdge(n2, n1),

        e3 = G.addEdge(n2, n3), e4 = G.addEdge(n1, n4),

        e5 = G.addEdge(n4, n3);

   // Check edge deletion

   G.erase(e2);

   checkGraphNodeList(G, 4);

   checkGraphEdgeList(G, 4);

   checkGraphArcList(G, 8);

   checkGraphIncEdgeArcLists(G, n1, 2);

   checkGraphIncEdgeArcLists(G, n2, 2);

   checkGraphIncEdgeArcLists(G, n3, 2);

   checkGraphIncEdgeArcLists(G, n4, 2);

   checkGraphConEdgeList(G, 4);

   checkGraphConArcList(G, 8);

   // Check node deletion

   G.erase(n3);

   checkGraphNodeList(G, 3);

   checkGraphEdgeList(G, 2);

   checkGraphArcList(G, 4);

   checkGraphIncEdgeArcLists(G, n1, 2);

   checkGraphIncEdgeArcLists(G, n2, 1);

   checkGraphIncEdgeArcLists(G, n4, 1);

   checkGraphConEdgeList(G, 2);

   checkGraphConArcList(G, 4);

 }

 template <class Graph>

 void checkGraphSnapshot() {

   TEMPLATE_GRAPH_TYPEDEFS(Graph);

   Graph G;

   Node n1 = G.addNode(), n2 = G.addNode(), n3 = G.addNode();

   Edge e1 = G.addEdge(n1, n2), e2 = G.addEdge(n2, n1),

        e3 = G.addEdge(n2, n3);

   checkGraphNodeList(G, 3);

   checkGraphEdgeList(G, 3);

   checkGraphArcList(G, 6);

   typename Graph::Snapshot snapshot(G);

   Node n = G.addNode();

   G.addEdge(n3, n);

   G.addEdge(n, n3);

   G.addEdge(n3, n2);

   checkGraphNodeList(G, 4);

   checkGraphEdgeList(G, 6);

   checkGraphArcList(G, 12);

   snapshot.restore();

   checkGraphNodeList(G, 3);

   checkGraphEdgeList(G, 3);

   checkGraphArcList(G, 6);

   checkGraphIncEdgeArcLists(G, n1, 2);

   checkGraphIncEdgeArcLists(G, n2, 3);

   checkGraphIncEdgeArcLists(G, n3, 1);

   checkGraphConEdgeList(G, 3);

   checkGraphConArcList(G, 6);

   checkNodeIds(G);

   checkEdgeIds(G);

   checkArcIds(G);

   checkGraphNodeMap(G);

   checkGraphEdgeMap(G);

   checkGraphArcMap(G);

   G.addNode();

   snapshot.save(G);

   G.addEdge(G.addNode(), G.addNode());

   snapshot.restore();

   snapshot.save(G);

   checkGraphNodeList(G, 4);

   checkGraphEdgeList(G, 3);

   checkGraphArcList(G, 6);

   G.addEdge(G.addNode(), G.addNode());

   snapshot.restore();

   checkGraphNodeList(G, 4);

   checkGraphEdgeList(G, 3);

   checkGraphArcList(G, 6);

 }

 void checkFullGraph(int num) {

   typedef FullGraph Graph;

   GRAPH_TYPEDEFS(Graph);

   Graph G(num);

   check(G.nodeNum() == num && G.edgeNum() == num * (num - 1) / 2,

         "Wrong size");

   G.resize(num);

   check(G.nodeNum() == num && G.edgeNum() == num * (num - 1) / 2,

         "Wrong size");

   checkGraphNodeList(G, num);

   checkGraphEdgeList(G, num * (num - 1) / 2);

   for (NodeIt n(G); n != INVALID; ++n) {

     checkGraphOutArcList(G, n, num - 1);

     checkGraphInArcList(G, n, num - 1);

     checkGraphIncEdgeList(G, n, num - 1);

   }

   checkGraphConArcList(G, num * (num - 1));

   checkGraphConEdgeList(G, num * (num - 1) / 2);

   checkArcDirections(G);

   checkNodeIds(G);

   checkArcIds(G);

   checkEdgeIds(G);

   checkGraphNodeMap(G);

   checkGraphArcMap(G);

   checkGraphEdgeMap(G);

   for (int i = 0; i < G.nodeNum(); ++i) {

     check(G.index(G(i)) == i, "Wrong index");

   }

   for (NodeIt u(G); u != INVALID; ++u) {

     for (NodeIt v(G); v != INVALID; ++v) {

       Edge e = G.edge(u, v);

       Arc a = G.arc(u, v);

       if (u == v) {

         check(e == INVALID, "Wrong edge lookup");

         check(a == INVALID, "Wrong arc lookup");

       } else {

         check((G.u(e) == u && G.v(e) == v) ||

               (G.u(e) == v && G.v(e) == u), "Wrong edge lookup");

         check(G.source(a) == u && G.target(a) == v, "Wrong arc lookup");

       }

     }

   }

 }

 void checkConcepts() {

   { // Checking graph components

     checkConcept<BaseGraphComponent, BaseGraphComponent >();

     checkConcept<IDableGraphComponent<>,

       IDableGraphComponent<> >();

     checkConcept<IterableGraphComponent<>,

       IterableGraphComponent<> >();

     checkConcept<MappableGraphComponent<>,

       MappableGraphComponent<> >();

   }

   { // Checking skeleton graph

     checkConcept<Graph, Graph>();

   }

   { // Checking ListGraph

     checkConcept<Graph, ListGraph>();

     checkConcept<AlterableGraphComponent<>, ListGraph>();

     checkConcept<ExtendableGraphComponent<>, ListGraph>();

     checkConcept<ClearableGraphComponent<>, ListGraph>();

     checkConcept<ErasableGraphComponent<>, ListGraph>();

   }

   { // Checking SmartGraph

     checkConcept<Graph, SmartGraph>();

     checkConcept<AlterableGraphComponent<>, SmartGraph>();

     checkConcept<ExtendableGraphComponent<>, SmartGraph>();

     checkConcept<ClearableGraphComponent<>, SmartGraph>();

   }

   { // Checking FullGraph

     checkConcept<Graph, FullGraph>();

   }

   { // Checking GridGraph

     checkConcept<Graph, GridGraph>();

   }

   { // Checking HypercubeGraph

     checkConcept<Graph, HypercubeGraph>();

   }

 }

 template <typename Graph>

 void checkGraphValidity() {

   TEMPLATE_GRAPH_TYPEDEFS(Graph);

   Graph g;

   Node

     n1 = g.addNode(),

     n2 = g.addNode(),

     n3 = g.addNode();

   Edge

     e1 = g.addEdge(n1, n2),

     e2 = g.addEdge(n2, n3);

   check(g.valid(n1), "Wrong validity check");

   check(g.valid(e1), "Wrong validity check");

   check(g.valid(g.direct(e1, true)), "Wrong validity check");

   check(!g.valid(g.nodeFromId(-1)), "Wrong validity check");

   check(!g.valid(g.edgeFromId(-1)), "Wrong validity check");

   check(!g.valid(g.arcFromId(-1)), "Wrong validity check");

 }

 template <typename Graph>

 void checkGraphValidityErase() {

   TEMPLATE_GRAPH_TYPEDEFS(Graph);

   Graph g;

   Node

     n1 = g.addNode(),

     n2 = g.addNode(),

     n3 = g.addNode();

   Edge

     e1 = g.addEdge(n1, n2),

     e2 = g.addEdge(n2, n3);

   check(g.valid(n1), "Wrong validity check");

   check(g.valid(e1), "Wrong validity check");

   check(g.valid(g.direct(e1, true)), "Wrong validity check");

   g.erase(n1);

   check(!g.valid(n1), "Wrong validity check");

   check(g.valid(n2), "Wrong validity check");

   check(g.valid(n3), "Wrong validity check");

   check(!g.valid(e1), "Wrong validity check");

   check(g.valid(e2), "Wrong validity check");

   check(!g.valid(g.nodeFromId(-1)), "Wrong validity check");

   check(!g.valid(g.edgeFromId(-1)), "Wrong validity check");

   check(!g.valid(g.arcFromId(-1)), "Wrong validity check");

 }

 void checkGridGraph(int width, int height) {

   typedef GridGraph Graph;

   GRAPH_TYPEDEFS(Graph);

   Graph G(width, height);

   check(G.width() == width, "Wrong column number");

   check(G.height() == height, "Wrong row number");

   G.resize(width, height);

   check(G.width() == width, "Wrong column number");

   check(G.height() == height, "Wrong row number");

   for (int i = 0; i < width; ++i) {

     for (int j = 0; j < height; ++j) {

       check(G.col(G(i, j)) == i, "Wrong column");

       check(G.row(G(i, j)) == j, "Wrong row");

       check(G.pos(G(i, j)).x == i, "Wrong column");

       check(G.pos(G(i, j)).y == j, "Wrong row");

     }

   }

   for (int j = 0; j < height; ++j) {

     for (int i = 0; i < width - 1; ++i) {

       check(G.source(G.right(G(i, j))) == G(i, j), "Wrong right");

       check(G.target(G.right(G(i, j))) == G(i + 1, j), "Wrong right");

     }

     check(G.right(G(width - 1, j)) == INVALID, "Wrong right");

   }

   for (int j = 0; j < height; ++j) {

     for (int i = 1; i < width; ++i) {

       check(G.source(G.left(G(i, j))) == G(i, j), "Wrong left");

       check(G.target(G.left(G(i, j))) == G(i - 1, j), "Wrong left");

     }

     check(G.left(G(0, j)) == INVALID, "Wrong left");

   }

   for (int i = 0; i < width; ++i) {

     for (int j = 0; j < height - 1; ++j) {

       check(G.source(G.up(G(i, j))) == G(i, j), "Wrong up");

       check(G.target(G.up(G(i, j))) == G(i, j + 1), "Wrong up");

     }

     check(G.up(G(i, height - 1)) == INVALID, "Wrong up");

   }

   for (int i = 0; i < width; ++i) {

     for (int j = 1; j < height; ++j) {

       check(G.source(G.down(G(i, j))) == G(i, j), "Wrong down");

       check(G.target(G.down(G(i, j))) == G(i, j - 1), "Wrong down");

     }

     check(G.down(G(i, 0)) == INVALID, "Wrong down");

   }

   checkGraphNodeList(G, width * height);

   checkGraphEdgeList(G, width * (height - 1) + (width - 1) * height);

   checkGraphArcList(G, 2 * (width * (height - 1) + (width - 1) * height));

   for (NodeIt n(G); n != INVALID; ++n) {

     int nb = 4;

     if (G.col(n) == 0) --nb;

     if (G.col(n) == width - 1) --nb;

     if (G.row(n) == 0) --nb;

     if (G.row(n) == height - 1) --nb;

     checkGraphOutArcList(G, n, nb);

     checkGraphInArcList(G, n, nb);

     checkGraphIncEdgeList(G, n, nb);

   }

   checkArcDirections(G);

   checkGraphConArcList(G, 2 * (width * (height - 1) + (width - 1) * height));

   checkGraphConEdgeList(G, width * (height - 1) + (width - 1) * height);

   checkNodeIds(G);

   checkArcIds(G);

   checkEdgeIds(G);

   checkGraphNodeMap(G);

   checkGraphArcMap(G);

   checkGraphEdgeMap(G);

 }

 void checkHypercubeGraph(int dim) {

   GRAPH_TYPEDEFS(HypercubeGraph);

   HypercubeGraph G(dim);

   check(G.dimension() == dim, "Wrong dimension");

   G.resize(dim);

   check(G.dimension() == dim, "Wrong dimension");

   checkGraphNodeList(G, 1 << dim);

   checkGraphEdgeList(G, dim * (1 << (dim-1)));

   checkGraphArcList(G, dim * (1 << dim));

   Node n = G.nodeFromId(dim);

   for (NodeIt n(G); n != INVALID; ++n) {

     checkGraphIncEdgeList(G, n, dim);

     for (IncEdgeIt e(G, n); e != INVALID; ++e) {

       check( (G.u(e) == n &&

               G.id(G.v(e)) == (G.id(n) ^ (1 << G.dimension(e)))) ||

              (G.v(e) == n &&

               G.id(G.u(e)) == (G.id(n) ^ (1 << G.dimension(e)))),

              "Wrong edge or wrong dimension");

     }

     checkGraphOutArcList(G, n, dim);

     for (OutArcIt a(G, n); a != INVALID; ++a) {

       check(G.source(a) == n &&

             G.id(G.target(a)) == (G.id(n) ^ (1 << G.dimension(a))),

             "Wrong arc or wrong dimension");

     }

     checkGraphInArcList(G, n, dim);

     for (InArcIt a(G, n); a != INVALID; ++a) {

       check(G.target(a) == n &&

             G.id(G.source(a)) == (G.id(n) ^ (1 << G.dimension(a))),

             "Wrong arc or wrong dimension");

     }

   }

   checkGraphConArcList(G, (1 << dim) * dim);

   checkGraphConEdgeList(G, dim * (1 << (dim-1)));

   checkArcDirections(G);

   checkNodeIds(G);

   checkArcIds(G);

   checkEdgeIds(G);

   checkGraphNodeMap(G);

   checkGraphArcMap(G);

   checkGraphEdgeMap(G);

 }

 void checkGraphs() {

   { // Checking ListGraph

     checkGraphBuild<ListGraph>();

     checkGraphAlter<ListGraph>();

     checkGraphErase<ListGraph>();

     checkGraphSnapshot<ListGraph>();

     checkGraphValidityErase<ListGraph>();

   }

   { // Checking SmartGraph

     checkGraphBuild<SmartGraph>();

     checkGraphSnapshot<SmartGraph>();

     checkGraphValidity<SmartGraph>();

   }

   { // Checking FullGraph

     checkFullGraph(7);

     checkFullGraph(8);

   }

   { // Checking GridGraph

     checkGridGraph(5, 8);

     checkGridGraph(8, 5);

     checkGridGraph(5, 5);

     checkGridGraph(0, 0);

     checkGridGraph(1, 1);

   }

   { // Checking HypercubeGraph

     checkHypercubeGraph(1);

     checkHypercubeGraph(2);

     checkHypercubeGraph(3);

     checkHypercubeGraph(4);

   }

 }

 int main() {

   checkConcepts();

   checkGraphs();

   return 0;

 }