alpar@209: /* -*- mode: C++; indent-tabs-mode: nil; -*-
deba@200:  *
alpar@209:  * This file is a part of LEMON, a generic C++ optimization library.
deba@200:  *
alpar@440:  * Copyright (C) 2003-2009
deba@200:  * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
deba@200:  * (Egervary Research Group on Combinatorial Optimization, EGRES).
deba@200:  *
deba@200:  * Permission to use, modify and distribute this software is granted
deba@200:  * provided that this copyright notice appears in all copies. For
deba@200:  * precise terms see the accompanying LICENSE file.
deba@200:  *
deba@200:  * This software is provided "AS IS" with no warranty of any kind,
deba@200:  * express or implied, and with no claim as to its suitability for any
deba@200:  * purpose.
deba@200:  *
deba@200:  */
deba@200: 
deba@200: #include <lemon/smart_graph.h>
deba@200: #include <lemon/list_graph.h>
kpeter@894: #include <lemon/static_graph.h>
deba@200: #include <lemon/lgf_reader.h>
deba@200: #include <lemon/error.h>
deba@200: 
deba@200: #include "test_tools.h"
deba@200: 
deba@200: using namespace std;
deba@200: using namespace lemon;
deba@200: 
kpeter@894: template <typename GR>
deba@200: void digraph_copy_test() {
deba@200:   const int nn = 10;
deba@200: 
kpeter@890:   // Build a digraph
deba@200:   SmartDigraph from;
deba@200:   SmartDigraph::NodeMap<int> fnm(from);
deba@200:   SmartDigraph::ArcMap<int> fam(from);
deba@200:   SmartDigraph::Node fn = INVALID;
deba@200:   SmartDigraph::Arc fa = INVALID;
deba@200: 
deba@200:   std::vector<SmartDigraph::Node> fnv;
deba@200:   for (int i = 0; i < nn; ++i) {
deba@200:     SmartDigraph::Node node = from.addNode();
deba@200:     fnv.push_back(node);
deba@200:     fnm[node] = i * i;
deba@200:     if (i == 0) fn = node;
deba@200:   }
deba@200: 
deba@200:   for (int i = 0; i < nn; ++i) {
deba@200:     for (int j = 0; j < nn; ++j) {
deba@200:       SmartDigraph::Arc arc = from.addArc(fnv[i], fnv[j]);
deba@200:       fam[arc] = i + j * j;
deba@200:       if (i == 0 && j == 0) fa = arc;
deba@200:     }
deba@200:   }
kpeter@894:   
kpeter@894:   // Test digraph copy
kpeter@894:   GR to;
kpeter@894:   typename GR::template NodeMap<int> tnm(to);
kpeter@894:   typename GR::template ArcMap<int> tam(to);
kpeter@894:   typename GR::Node tn;
kpeter@894:   typename GR::Arc ta;
deba@200: 
kpeter@894:   SmartDigraph::NodeMap<typename GR::Node> nr(from);
kpeter@894:   SmartDigraph::ArcMap<typename GR::Arc> er(from);
deba@200: 
kpeter@894:   typename GR::template NodeMap<SmartDigraph::Node> ncr(to);
kpeter@894:   typename GR::template ArcMap<SmartDigraph::Arc> ecr(to);
deba@200: 
kpeter@282:   digraphCopy(from, to).
kpeter@282:     nodeMap(fnm, tnm).arcMap(fam, tam).
deba@200:     nodeRef(nr).arcRef(er).
deba@200:     nodeCrossRef(ncr).arcCrossRef(ecr).
kpeter@282:     node(fn, tn).arc(fa, ta).run();
kpeter@890:   
kpeter@890:   check(countNodes(from) == countNodes(to), "Wrong copy.");
kpeter@890:   check(countArcs(from) == countArcs(to), "Wrong copy.");
deba@200: 
deba@200:   for (SmartDigraph::NodeIt it(from); it != INVALID; ++it) {
deba@200:     check(ncr[nr[it]] == it, "Wrong copy.");
deba@200:     check(fnm[it] == tnm[nr[it]], "Wrong copy.");
deba@200:   }
deba@200: 
deba@200:   for (SmartDigraph::ArcIt it(from); it != INVALID; ++it) {
deba@200:     check(ecr[er[it]] == it, "Wrong copy.");
deba@200:     check(fam[it] == tam[er[it]], "Wrong copy.");
deba@200:     check(nr[from.source(it)] == to.source(er[it]), "Wrong copy.");
deba@200:     check(nr[from.target(it)] == to.target(er[it]), "Wrong copy.");
deba@200:   }
deba@200: 
kpeter@894:   for (typename GR::NodeIt it(to); it != INVALID; ++it) {
deba@200:     check(nr[ncr[it]] == it, "Wrong copy.");
deba@200:   }
deba@200: 
kpeter@894:   for (typename GR::ArcIt it(to); it != INVALID; ++it) {
deba@200:     check(er[ecr[it]] == it, "Wrong copy.");
deba@200:   }
deba@200:   check(tn == nr[fn], "Wrong copy.");
deba@200:   check(ta == er[fa], "Wrong copy.");
kpeter@890: 
kpeter@890:   // Test repeated copy
kpeter@890:   digraphCopy(from, to).run();
kpeter@890:   
kpeter@890:   check(countNodes(from) == countNodes(to), "Wrong copy.");
kpeter@890:   check(countArcs(from) == countArcs(to), "Wrong copy.");
deba@200: }
deba@200: 
kpeter@894: template <typename GR>
deba@200: void graph_copy_test() {
deba@200:   const int nn = 10;
deba@200: 
kpeter@890:   // Build a graph
deba@200:   SmartGraph from;
deba@200:   SmartGraph::NodeMap<int> fnm(from);
deba@200:   SmartGraph::ArcMap<int> fam(from);
deba@200:   SmartGraph::EdgeMap<int> fem(from);
deba@200:   SmartGraph::Node fn = INVALID;
deba@200:   SmartGraph::Arc fa = INVALID;
deba@200:   SmartGraph::Edge fe = INVALID;
deba@200: 
deba@200:   std::vector<SmartGraph::Node> fnv;
deba@200:   for (int i = 0; i < nn; ++i) {
deba@200:     SmartGraph::Node node = from.addNode();
deba@200:     fnv.push_back(node);
deba@200:     fnm[node] = i * i;
deba@200:     if (i == 0) fn = node;
deba@200:   }
deba@200: 
deba@200:   for (int i = 0; i < nn; ++i) {
deba@200:     for (int j = 0; j < nn; ++j) {
deba@200:       SmartGraph::Edge edge = from.addEdge(fnv[i], fnv[j]);
deba@200:       fem[edge] = i * i + j * j;
deba@200:       fam[from.direct(edge, true)] = i + j * j;
deba@200:       fam[from.direct(edge, false)] = i * i + j;
deba@200:       if (i == 0 && j == 0) fa = from.direct(edge, true);
deba@200:       if (i == 0 && j == 0) fe = edge;
deba@200:     }
deba@200:   }
alpar@209: 
kpeter@890:   // Test graph copy
kpeter@894:   GR to;
kpeter@894:   typename GR::template NodeMap<int> tnm(to);
kpeter@894:   typename GR::template ArcMap<int> tam(to);
kpeter@894:   typename GR::template EdgeMap<int> tem(to);
kpeter@894:   typename GR::Node tn;
kpeter@894:   typename GR::Arc ta;
kpeter@894:   typename GR::Edge te;
deba@200: 
kpeter@894:   SmartGraph::NodeMap<typename GR::Node> nr(from);
kpeter@894:   SmartGraph::ArcMap<typename GR::Arc> ar(from);
kpeter@894:   SmartGraph::EdgeMap<typename GR::Edge> er(from);
deba@200: 
kpeter@894:   typename GR::template NodeMap<SmartGraph::Node> ncr(to);
kpeter@894:   typename GR::template ArcMap<SmartGraph::Arc> acr(to);
kpeter@894:   typename GR::template EdgeMap<SmartGraph::Edge> ecr(to);
deba@200: 
kpeter@282:   graphCopy(from, to).
kpeter@282:     nodeMap(fnm, tnm).arcMap(fam, tam).edgeMap(fem, tem).
deba@200:     nodeRef(nr).arcRef(ar).edgeRef(er).
deba@200:     nodeCrossRef(ncr).arcCrossRef(acr).edgeCrossRef(ecr).
kpeter@282:     node(fn, tn).arc(fa, ta).edge(fe, te).run();
deba@200: 
kpeter@890:   check(countNodes(from) == countNodes(to), "Wrong copy.");
kpeter@890:   check(countEdges(from) == countEdges(to), "Wrong copy.");
kpeter@890:   check(countArcs(from) == countArcs(to), "Wrong copy.");
kpeter@890: 
deba@200:   for (SmartGraph::NodeIt it(from); it != INVALID; ++it) {
deba@200:     check(ncr[nr[it]] == it, "Wrong copy.");
deba@200:     check(fnm[it] == tnm[nr[it]], "Wrong copy.");
deba@200:   }
deba@200: 
deba@200:   for (SmartGraph::ArcIt it(from); it != INVALID; ++it) {
deba@200:     check(acr[ar[it]] == it, "Wrong copy.");
deba@200:     check(fam[it] == tam[ar[it]], "Wrong copy.");
deba@200:     check(nr[from.source(it)] == to.source(ar[it]), "Wrong copy.");
deba@200:     check(nr[from.target(it)] == to.target(ar[it]), "Wrong copy.");
deba@200:   }
deba@200: 
deba@200:   for (SmartGraph::EdgeIt it(from); it != INVALID; ++it) {
deba@200:     check(ecr[er[it]] == it, "Wrong copy.");
deba@200:     check(fem[it] == tem[er[it]], "Wrong copy.");
alpar@209:     check(nr[from.u(it)] == to.u(er[it]) || nr[from.u(it)] == to.v(er[it]),
alpar@209:           "Wrong copy.");
alpar@209:     check(nr[from.v(it)] == to.u(er[it]) || nr[from.v(it)] == to.v(er[it]),
alpar@209:           "Wrong copy.");
alpar@209:     check((from.u(it) != from.v(it)) == (to.u(er[it]) != to.v(er[it])),
alpar@209:           "Wrong copy.");
deba@200:   }
deba@200: 
kpeter@894:   for (typename GR::NodeIt it(to); it != INVALID; ++it) {
deba@200:     check(nr[ncr[it]] == it, "Wrong copy.");
deba@200:   }
deba@200: 
kpeter@894:   for (typename GR::ArcIt it(to); it != INVALID; ++it) {
deba@200:     check(ar[acr[it]] == it, "Wrong copy.");
deba@200:   }
kpeter@894:   for (typename GR::EdgeIt it(to); it != INVALID; ++it) {
deba@200:     check(er[ecr[it]] == it, "Wrong copy.");
deba@200:   }
deba@200:   check(tn == nr[fn], "Wrong copy.");
deba@200:   check(ta == ar[fa], "Wrong copy.");
deba@200:   check(te == er[fe], "Wrong copy.");
kpeter@890: 
kpeter@890:   // Test repeated copy
kpeter@890:   graphCopy(from, to).run();
kpeter@890:   
kpeter@890:   check(countNodes(from) == countNodes(to), "Wrong copy.");
kpeter@890:   check(countEdges(from) == countEdges(to), "Wrong copy.");
kpeter@890:   check(countArcs(from) == countArcs(to), "Wrong copy.");
deba@200: }
deba@200: 
deba@1022: template <typename GR>
deba@1022: void bpgraph_copy_test() {
deba@1022:   const int nn = 10;
deba@1022: 
deba@1022:   // Build a graph
deba@1022:   SmartBpGraph from;
deba@1022:   SmartBpGraph::NodeMap<int> fnm(from);
deba@1026:   SmartBpGraph::RedNodeMap<int> frnm(from);
deba@1026:   SmartBpGraph::BlueNodeMap<int> fbnm(from);
deba@1022:   SmartBpGraph::ArcMap<int> fam(from);
deba@1022:   SmartBpGraph::EdgeMap<int> fem(from);
deba@1022:   SmartBpGraph::Node fn = INVALID;
deba@1025:   SmartBpGraph::RedNode frn = INVALID;
deba@1025:   SmartBpGraph::BlueNode fbn = INVALID;
deba@1022:   SmartBpGraph::Arc fa = INVALID;
deba@1022:   SmartBpGraph::Edge fe = INVALID;
deba@1022: 
deba@1025:   std::vector<SmartBpGraph::RedNode> frnv;
deba@1022:   for (int i = 0; i < nn; ++i) {
deba@1025:     SmartBpGraph::RedNode node = from.addRedNode();
deba@1022:     frnv.push_back(node);
deba@1022:     fnm[node] = i * i;
deba@1022:     frnm[node] = i + i;
deba@1025:     if (i == 0) {
deba@1025:       fn = node;
deba@1025:       frn = node;
deba@1025:     }
deba@1022:   }
deba@1022: 
deba@1025:   std::vector<SmartBpGraph::BlueNode> fbnv;
deba@1022:   for (int i = 0; i < nn; ++i) {
deba@1025:     SmartBpGraph::BlueNode node = from.addBlueNode();
deba@1022:     fbnv.push_back(node);
deba@1022:     fnm[node] = i * i;
deba@1022:     fbnm[node] = i + i;
deba@1025:     if (i == 0) fbn = node;
deba@1022:   }
deba@1022: 
deba@1022:   for (int i = 0; i < nn; ++i) {
deba@1022:     for (int j = 0; j < nn; ++j) {
deba@1022:       SmartBpGraph::Edge edge = from.addEdge(frnv[i], fbnv[j]);
deba@1022:       fem[edge] = i * i + j * j;
deba@1022:       fam[from.direct(edge, true)] = i + j * j;
deba@1022:       fam[from.direct(edge, false)] = i * i + j;
deba@1022:       if (i == 0 && j == 0) fa = from.direct(edge, true);
deba@1022:       if (i == 0 && j == 0) fe = edge;
deba@1022:     }
deba@1022:   }
deba@1022: 
deba@1022:   // Test graph copy
deba@1022:   GR to;
deba@1022:   typename GR::template NodeMap<int> tnm(to);
deba@1026:   typename GR::template RedNodeMap<int> trnm(to);
deba@1026:   typename GR::template BlueNodeMap<int> tbnm(to);
deba@1022:   typename GR::template ArcMap<int> tam(to);
deba@1022:   typename GR::template EdgeMap<int> tem(to);
deba@1022:   typename GR::Node tn;
deba@1025:   typename GR::RedNode trn;
deba@1025:   typename GR::BlueNode tbn;
deba@1022:   typename GR::Arc ta;
deba@1022:   typename GR::Edge te;
deba@1022: 
deba@1022:   SmartBpGraph::NodeMap<typename GR::Node> nr(from);
deba@1026:   SmartBpGraph::RedNodeMap<typename GR::RedNode> rnr(from);
deba@1026:   SmartBpGraph::BlueNodeMap<typename GR::BlueNode> bnr(from);
deba@1022:   SmartBpGraph::ArcMap<typename GR::Arc> ar(from);
deba@1022:   SmartBpGraph::EdgeMap<typename GR::Edge> er(from);
deba@1022: 
deba@1022:   typename GR::template NodeMap<SmartBpGraph::Node> ncr(to);
deba@1026:   typename GR::template RedNodeMap<SmartBpGraph::RedNode> rncr(to);
deba@1026:   typename GR::template BlueNodeMap<SmartBpGraph::BlueNode> bncr(to);
deba@1022:   typename GR::template ArcMap<SmartBpGraph::Arc> acr(to);
deba@1022:   typename GR::template EdgeMap<SmartBpGraph::Edge> ecr(to);
deba@1022: 
deba@1022:   bpGraphCopy(from, to).
deba@1026:     nodeMap(fnm, tnm).
deba@1026:     redNodeMap(frnm, trnm).blueNodeMap(fbnm, tbnm).
deba@1022:     arcMap(fam, tam).edgeMap(fem, tem).
deba@1022:     nodeRef(nr).redRef(rnr).blueRef(bnr).
deba@1022:     arcRef(ar).edgeRef(er).
deba@1022:     nodeCrossRef(ncr).redCrossRef(rncr).blueCrossRef(bncr).
deba@1022:     arcCrossRef(acr).edgeCrossRef(ecr).
deba@1025:     node(fn, tn).redNode(frn, trn).blueNode(fbn, tbn).
deba@1025:     arc(fa, ta).edge(fe, te).run();
deba@1022: 
deba@1022:   check(countNodes(from) == countNodes(to), "Wrong copy.");
deba@1022:   check(countRedNodes(from) == countRedNodes(to), "Wrong copy.");
deba@1022:   check(countBlueNodes(from) == countBlueNodes(to), "Wrong copy.");
deba@1022:   check(countEdges(from) == countEdges(to), "Wrong copy.");
deba@1022:   check(countArcs(from) == countArcs(to), "Wrong copy.");
deba@1022: 
deba@1022:   for (SmartBpGraph::NodeIt it(from); it != INVALID; ++it) {
deba@1022:     check(ncr[nr[it]] == it, "Wrong copy.");
deba@1022:     check(fnm[it] == tnm[nr[it]], "Wrong copy.");
deba@1025:   }
deba@1025: 
deba@1026:   for (SmartBpGraph::RedNodeIt it(from); it != INVALID; ++it) {
deba@1025:     check(ncr[nr[it]] == it, "Wrong copy.");
deba@1025:     check(fnm[it] == tnm[nr[it]], "Wrong copy.");
deba@1025:     check(rnr[it] == nr[it], "Wrong copy.");
deba@1025:     check(rncr[rnr[it]] == it, "Wrong copy.");
deba@1025:     check(frnm[it] == trnm[rnr[it]], "Wrong copy.");
deba@1025:     check(to.red(rnr[it]), "Wrong copy.");
deba@1025:   }
deba@1025: 
deba@1026:   for (SmartBpGraph::BlueNodeIt it(from); it != INVALID; ++it) {
deba@1025:     check(ncr[nr[it]] == it, "Wrong copy.");
deba@1025:     check(fnm[it] == tnm[nr[it]], "Wrong copy.");
deba@1025:     check(bnr[it] == nr[it], "Wrong copy.");
deba@1025:     check(bncr[bnr[it]] == it, "Wrong copy.");
deba@1025:     check(fbnm[it] == tbnm[bnr[it]], "Wrong copy.");
deba@1025:     check(to.blue(bnr[it]), "Wrong copy.");
deba@1022:   }
deba@1022: 
deba@1022:   for (SmartBpGraph::ArcIt it(from); it != INVALID; ++it) {
deba@1022:     check(acr[ar[it]] == it, "Wrong copy.");
deba@1022:     check(fam[it] == tam[ar[it]], "Wrong copy.");
deba@1022:     check(nr[from.source(it)] == to.source(ar[it]), "Wrong copy.");
deba@1022:     check(nr[from.target(it)] == to.target(ar[it]), "Wrong copy.");
deba@1022:   }
deba@1022: 
deba@1022:   for (SmartBpGraph::EdgeIt it(from); it != INVALID; ++it) {
deba@1022:     check(ecr[er[it]] == it, "Wrong copy.");
deba@1022:     check(fem[it] == tem[er[it]], "Wrong copy.");
deba@1022:     check(nr[from.u(it)] == to.u(er[it]) || nr[from.u(it)] == to.v(er[it]),
deba@1022:           "Wrong copy.");
deba@1022:     check(nr[from.v(it)] == to.u(er[it]) || nr[from.v(it)] == to.v(er[it]),
deba@1022:           "Wrong copy.");
deba@1022:     check((from.u(it) != from.v(it)) == (to.u(er[it]) != to.v(er[it])),
deba@1022:           "Wrong copy.");
deba@1022:   }
deba@1022: 
deba@1022:   for (typename GR::NodeIt it(to); it != INVALID; ++it) {
deba@1022:     check(nr[ncr[it]] == it, "Wrong copy.");
deba@1022:   }
deba@1026:   for (typename GR::RedNodeIt it(to); it != INVALID; ++it) {
deba@1022:     check(rncr[it] == ncr[it], "Wrong copy.");
deba@1022:     check(rnr[rncr[it]] == it, "Wrong copy.");
deba@1022:   }
deba@1026:   for (typename GR::BlueNodeIt it(to); it != INVALID; ++it) {
deba@1022:     check(bncr[it] == ncr[it], "Wrong copy.");
deba@1022:     check(bnr[bncr[it]] == it, "Wrong copy.");
deba@1022:   }
deba@1022:   for (typename GR::ArcIt it(to); it != INVALID; ++it) {
deba@1022:     check(ar[acr[it]] == it, "Wrong copy.");
deba@1022:   }
deba@1022:   for (typename GR::EdgeIt it(to); it != INVALID; ++it) {
deba@1022:     check(er[ecr[it]] == it, "Wrong copy.");
deba@1022:   }
deba@1022:   check(tn == nr[fn], "Wrong copy.");
deba@1025:   check(trn == rnr[frn], "Wrong copy.");
deba@1025:   check(tbn == bnr[fbn], "Wrong copy.");
deba@1022:   check(ta == ar[fa], "Wrong copy.");
deba@1022:   check(te == er[fe], "Wrong copy.");
deba@1022: 
deba@1022:   // Test repeated copy
deba@1022:   bpGraphCopy(from, to).run();
deba@1022:   
deba@1022:   check(countNodes(from) == countNodes(to), "Wrong copy.");
deba@1022:   check(countRedNodes(from) == countRedNodes(to), "Wrong copy.");
deba@1022:   check(countBlueNodes(from) == countBlueNodes(to), "Wrong copy.");
deba@1022:   check(countEdges(from) == countEdges(to), "Wrong copy.");
deba@1022:   check(countArcs(from) == countArcs(to), "Wrong copy.");
deba@1022: }
deba@1022: 
deba@200: 
deba@200: int main() {
kpeter@894:   digraph_copy_test<SmartDigraph>();
kpeter@894:   digraph_copy_test<ListDigraph>();
kpeter@894:   digraph_copy_test<StaticDigraph>();
kpeter@894:   graph_copy_test<SmartGraph>();
kpeter@894:   graph_copy_test<ListGraph>();
deba@1022:   bpgraph_copy_test<SmartBpGraph>();
deba@1022:   bpgraph_copy_test<ListBpGraph>();
deba@200: 
alpar@209:   return 0;
deba@200: }