test/digraph_test.cc
author Peter Kovacs <kpeter@inf.elte.hu>
Fri, 03 Apr 2009 18:59:15 +0200
changeset 655 6ac5d9ae1d3d
parent 388 2d87dbd7f8c8
child 783 2e20aad15754
child 821 f4b5c2d5449d
child 1157 761fe0846f49
permissions -rw-r--r--
Support real types + numerical stability fix in NS (#254)

- Real types are supported by appropriate inicialization.
- A feature of the XTI spanning tree structure is removed to ensure
numerical stability (could cause problems using integer types).
The node potentials are updated always on the lower subtree,
in order to prevent overflow problems.
The former method isn't notably faster during to our tests.
alpar@209
     1
/* -*- mode: C++; indent-tabs-mode: nil; -*-
deba@57
     2
 *
alpar@209
     3
 * This file is a part of LEMON, a generic C++ optimization library.
deba@57
     4
 *
alpar@463
     5
 * Copyright (C) 2003-2009
deba@57
     6
 * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
deba@57
     7
 * (Egervary Research Group on Combinatorial Optimization, EGRES).
deba@57
     8
 *
deba@57
     9
 * Permission to use, modify and distribute this software is granted
deba@57
    10
 * provided that this copyright notice appears in all copies. For
deba@57
    11
 * precise terms see the accompanying LICENSE file.
deba@57
    12
 *
deba@57
    13
 * This software is provided "AS IS" with no warranty of any kind,
deba@57
    14
 * express or implied, and with no claim as to its suitability for any
deba@57
    15
 * purpose.
deba@57
    16
 *
deba@57
    17
 */
deba@57
    18
deba@57
    19
#include <lemon/concepts/digraph.h>
deba@57
    20
#include <lemon/list_graph.h>
kpeter@171
    21
#include <lemon/smart_graph.h>
deba@365
    22
#include <lemon/full_graph.h>
deba@57
    23
deba@57
    24
#include "test_tools.h"
kpeter@171
    25
#include "graph_test.h"
deba@57
    26
deba@57
    27
using namespace lemon;
deba@57
    28
using namespace lemon::concepts;
deba@57
    29
deba@228
    30
template <class Digraph>
kpeter@387
    31
void checkDigraphBuild() {
deba@228
    32
  TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
deba@228
    33
  Digraph G;
deba@228
    34
deba@228
    35
  checkGraphNodeList(G, 0);
deba@228
    36
  checkGraphArcList(G, 0);
deba@228
    37
deba@228
    38
  Node
deba@228
    39
    n1 = G.addNode(),
deba@228
    40
    n2 = G.addNode(),
deba@228
    41
    n3 = G.addNode();
deba@228
    42
  checkGraphNodeList(G, 3);
deba@228
    43
  checkGraphArcList(G, 0);
deba@228
    44
deba@228
    45
  Arc a1 = G.addArc(n1, n2);
deba@228
    46
  check(G.source(a1) == n1 && G.target(a1) == n2, "Wrong arc");
deba@228
    47
  checkGraphNodeList(G, 3);
deba@228
    48
  checkGraphArcList(G, 1);
deba@228
    49
deba@228
    50
  checkGraphOutArcList(G, n1, 1);
deba@228
    51
  checkGraphOutArcList(G, n2, 0);
deba@228
    52
  checkGraphOutArcList(G, n3, 0);
deba@228
    53
deba@228
    54
  checkGraphInArcList(G, n1, 0);
deba@228
    55
  checkGraphInArcList(G, n2, 1);
deba@228
    56
  checkGraphInArcList(G, n3, 0);
deba@228
    57
deba@228
    58
  checkGraphConArcList(G, 1);
deba@228
    59
kpeter@387
    60
  Arc a2 = G.addArc(n2, n1),
kpeter@387
    61
      a3 = G.addArc(n2, n3),
kpeter@387
    62
      a4 = G.addArc(n2, n3);
kpeter@387
    63
kpeter@387
    64
  checkGraphNodeList(G, 3);
kpeter@387
    65
  checkGraphArcList(G, 4);
kpeter@387
    66
kpeter@387
    67
  checkGraphOutArcList(G, n1, 1);
kpeter@387
    68
  checkGraphOutArcList(G, n2, 3);
kpeter@387
    69
  checkGraphOutArcList(G, n3, 0);
kpeter@387
    70
kpeter@387
    71
  checkGraphInArcList(G, n1, 1);
kpeter@387
    72
  checkGraphInArcList(G, n2, 1);
kpeter@387
    73
  checkGraphInArcList(G, n3, 2);
kpeter@387
    74
kpeter@387
    75
  checkGraphConArcList(G, 4);
kpeter@387
    76
kpeter@387
    77
  checkNodeIds(G);
kpeter@387
    78
  checkArcIds(G);
kpeter@387
    79
  checkGraphNodeMap(G);
kpeter@387
    80
  checkGraphArcMap(G);
kpeter@387
    81
}
kpeter@387
    82
kpeter@387
    83
template <class Digraph>
kpeter@387
    84
void checkDigraphSplit() {
kpeter@387
    85
  TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
kpeter@387
    86
kpeter@387
    87
  Digraph G;
kpeter@387
    88
  Node n1 = G.addNode(), n2 = G.addNode(), n3 = G.addNode();
kpeter@387
    89
  Arc a1 = G.addArc(n1, n2), a2 = G.addArc(n2, n1),
kpeter@387
    90
      a3 = G.addArc(n2, n3), a4 = G.addArc(n2, n3);
kpeter@387
    91
kpeter@387
    92
  Node n4 = G.split(n2);
kpeter@387
    93
kpeter@387
    94
  check(G.target(OutArcIt(G, n2)) == n4 &&
kpeter@387
    95
        G.source(InArcIt(G, n4)) == n2,
kpeter@387
    96
        "Wrong split.");
kpeter@387
    97
kpeter@387
    98
  checkGraphNodeList(G, 4);
kpeter@387
    99
  checkGraphArcList(G, 5);
kpeter@387
   100
kpeter@387
   101
  checkGraphOutArcList(G, n1, 1);
kpeter@387
   102
  checkGraphOutArcList(G, n2, 1);
kpeter@387
   103
  checkGraphOutArcList(G, n3, 0);
kpeter@387
   104
  checkGraphOutArcList(G, n4, 3);
kpeter@387
   105
kpeter@387
   106
  checkGraphInArcList(G, n1, 1);
kpeter@387
   107
  checkGraphInArcList(G, n2, 1);
kpeter@387
   108
  checkGraphInArcList(G, n3, 2);
kpeter@387
   109
  checkGraphInArcList(G, n4, 1);
kpeter@387
   110
kpeter@387
   111
  checkGraphConArcList(G, 5);
kpeter@387
   112
}
kpeter@387
   113
kpeter@387
   114
template <class Digraph>
kpeter@387
   115
void checkDigraphAlter() {
kpeter@387
   116
  TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
kpeter@387
   117
kpeter@387
   118
  Digraph G;
kpeter@387
   119
  Node n1 = G.addNode(), n2 = G.addNode(),
kpeter@387
   120
       n3 = G.addNode(), n4 = G.addNode();
kpeter@387
   121
  Arc a1 = G.addArc(n1, n2), a2 = G.addArc(n4, n1),
kpeter@387
   122
      a3 = G.addArc(n4, n3), a4 = G.addArc(n4, n3),
kpeter@387
   123
      a5 = G.addArc(n2, n4);
kpeter@387
   124
kpeter@387
   125
  checkGraphNodeList(G, 4);
kpeter@387
   126
  checkGraphArcList(G, 5);
kpeter@387
   127
kpeter@387
   128
  // Check changeSource() and changeTarget()
kpeter@387
   129
  G.changeTarget(a4, n1);
kpeter@387
   130
kpeter@387
   131
  checkGraphNodeList(G, 4);
kpeter@387
   132
  checkGraphArcList(G, 5);
kpeter@387
   133
kpeter@387
   134
  checkGraphOutArcList(G, n1, 1);
kpeter@387
   135
  checkGraphOutArcList(G, n2, 1);
kpeter@387
   136
  checkGraphOutArcList(G, n3, 0);
kpeter@387
   137
  checkGraphOutArcList(G, n4, 3);
kpeter@387
   138
kpeter@387
   139
  checkGraphInArcList(G, n1, 2);
kpeter@387
   140
  checkGraphInArcList(G, n2, 1);
kpeter@387
   141
  checkGraphInArcList(G, n3, 1);
kpeter@387
   142
  checkGraphInArcList(G, n4, 1);
kpeter@387
   143
kpeter@387
   144
  checkGraphConArcList(G, 5);
kpeter@387
   145
kpeter@387
   146
  G.changeSource(a4, n3);
kpeter@387
   147
kpeter@387
   148
  checkGraphNodeList(G, 4);
kpeter@387
   149
  checkGraphArcList(G, 5);
kpeter@387
   150
kpeter@387
   151
  checkGraphOutArcList(G, n1, 1);
kpeter@387
   152
  checkGraphOutArcList(G, n2, 1);
kpeter@387
   153
  checkGraphOutArcList(G, n3, 1);
kpeter@387
   154
  checkGraphOutArcList(G, n4, 2);
kpeter@387
   155
kpeter@387
   156
  checkGraphInArcList(G, n1, 2);
kpeter@387
   157
  checkGraphInArcList(G, n2, 1);
kpeter@387
   158
  checkGraphInArcList(G, n3, 1);
kpeter@387
   159
  checkGraphInArcList(G, n4, 1);
kpeter@387
   160
kpeter@387
   161
  checkGraphConArcList(G, 5);
kpeter@387
   162
kpeter@387
   163
  // Check contract()
kpeter@387
   164
  G.contract(n2, n4, false);
kpeter@387
   165
kpeter@387
   166
  checkGraphNodeList(G, 3);
kpeter@387
   167
  checkGraphArcList(G, 5);
kpeter@387
   168
kpeter@387
   169
  checkGraphOutArcList(G, n1, 1);
kpeter@387
   170
  checkGraphOutArcList(G, n2, 3);
kpeter@387
   171
  checkGraphOutArcList(G, n3, 1);
kpeter@387
   172
kpeter@387
   173
  checkGraphInArcList(G, n1, 2);
kpeter@387
   174
  checkGraphInArcList(G, n2, 2);
kpeter@387
   175
  checkGraphInArcList(G, n3, 1);
kpeter@387
   176
kpeter@387
   177
  checkGraphConArcList(G, 5);
kpeter@387
   178
kpeter@387
   179
  G.contract(n2, n1);
kpeter@387
   180
kpeter@387
   181
  checkGraphNodeList(G, 2);
kpeter@387
   182
  checkGraphArcList(G, 3);
kpeter@387
   183
kpeter@387
   184
  checkGraphOutArcList(G, n2, 2);
kpeter@387
   185
  checkGraphOutArcList(G, n3, 1);
kpeter@387
   186
kpeter@387
   187
  checkGraphInArcList(G, n2, 2);
kpeter@387
   188
  checkGraphInArcList(G, n3, 1);
kpeter@387
   189
kpeter@387
   190
  checkGraphConArcList(G, 3);
kpeter@387
   191
}
kpeter@387
   192
kpeter@387
   193
template <class Digraph>
kpeter@387
   194
void checkDigraphErase() {
kpeter@387
   195
  TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
kpeter@387
   196
kpeter@387
   197
  Digraph G;
kpeter@387
   198
  Node n1 = G.addNode(), n2 = G.addNode(),
kpeter@387
   199
       n3 = G.addNode(), n4 = G.addNode();
kpeter@387
   200
  Arc a1 = G.addArc(n1, n2), a2 = G.addArc(n4, n1),
kpeter@387
   201
      a3 = G.addArc(n4, n3), a4 = G.addArc(n3, n1),
kpeter@387
   202
      a5 = G.addArc(n2, n4);
kpeter@387
   203
kpeter@387
   204
  // Check arc deletion
kpeter@387
   205
  G.erase(a1);
kpeter@387
   206
kpeter@387
   207
  checkGraphNodeList(G, 4);
kpeter@387
   208
  checkGraphArcList(G, 4);
kpeter@387
   209
kpeter@387
   210
  checkGraphOutArcList(G, n1, 0);
kpeter@387
   211
  checkGraphOutArcList(G, n2, 1);
kpeter@387
   212
  checkGraphOutArcList(G, n3, 1);
kpeter@387
   213
  checkGraphOutArcList(G, n4, 2);
kpeter@387
   214
kpeter@387
   215
  checkGraphInArcList(G, n1, 2);
kpeter@387
   216
  checkGraphInArcList(G, n2, 0);
kpeter@387
   217
  checkGraphInArcList(G, n3, 1);
kpeter@387
   218
  checkGraphInArcList(G, n4, 1);
kpeter@387
   219
kpeter@387
   220
  checkGraphConArcList(G, 4);
kpeter@387
   221
kpeter@387
   222
  // Check node deletion
kpeter@387
   223
  G.erase(n4);
kpeter@387
   224
kpeter@387
   225
  checkGraphNodeList(G, 3);
kpeter@387
   226
  checkGraphArcList(G, 1);
kpeter@387
   227
kpeter@387
   228
  checkGraphOutArcList(G, n1, 0);
kpeter@387
   229
  checkGraphOutArcList(G, n2, 0);
kpeter@387
   230
  checkGraphOutArcList(G, n3, 1);
kpeter@387
   231
  checkGraphOutArcList(G, n4, 0);
kpeter@387
   232
kpeter@387
   233
  checkGraphInArcList(G, n1, 1);
kpeter@387
   234
  checkGraphInArcList(G, n2, 0);
kpeter@387
   235
  checkGraphInArcList(G, n3, 0);
kpeter@387
   236
  checkGraphInArcList(G, n4, 0);
kpeter@387
   237
kpeter@387
   238
  checkGraphConArcList(G, 1);
kpeter@387
   239
}
kpeter@387
   240
kpeter@387
   241
kpeter@387
   242
template <class Digraph>
kpeter@387
   243
void checkDigraphSnapshot() {
kpeter@387
   244
  TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
kpeter@387
   245
kpeter@387
   246
  Digraph G;
kpeter@387
   247
  Node n1 = G.addNode(), n2 = G.addNode(), n3 = G.addNode();
kpeter@387
   248
  Arc a1 = G.addArc(n1, n2), a2 = G.addArc(n2, n1),
kpeter@387
   249
      a3 = G.addArc(n2, n3), a4 = G.addArc(n2, n3);
kpeter@387
   250
kpeter@387
   251
  typename Digraph::Snapshot snapshot(G);
kpeter@387
   252
kpeter@387
   253
  Node n = G.addNode();
kpeter@387
   254
  G.addArc(n3, n);
kpeter@387
   255
  G.addArc(n, n3);
kpeter@387
   256
kpeter@387
   257
  checkGraphNodeList(G, 4);
kpeter@387
   258
  checkGraphArcList(G, 6);
kpeter@387
   259
kpeter@387
   260
  snapshot.restore();
kpeter@387
   261
deba@228
   262
  checkGraphNodeList(G, 3);
deba@228
   263
  checkGraphArcList(G, 4);
deba@228
   264
deba@228
   265
  checkGraphOutArcList(G, n1, 1);
deba@228
   266
  checkGraphOutArcList(G, n2, 3);
deba@228
   267
  checkGraphOutArcList(G, n3, 0);
deba@228
   268
deba@228
   269
  checkGraphInArcList(G, n1, 1);
deba@228
   270
  checkGraphInArcList(G, n2, 1);
deba@228
   271
  checkGraphInArcList(G, n3, 2);
deba@228
   272
deba@228
   273
  checkGraphConArcList(G, 4);
deba@228
   274
deba@228
   275
  checkNodeIds(G);
deba@228
   276
  checkArcIds(G);
deba@228
   277
  checkGraphNodeMap(G);
deba@228
   278
  checkGraphArcMap(G);
deba@228
   279
kpeter@387
   280
  G.addNode();
kpeter@387
   281
  snapshot.save(G);
kpeter@387
   282
kpeter@387
   283
  G.addArc(G.addNode(), G.addNode());
kpeter@387
   284
kpeter@387
   285
  snapshot.restore();
kpeter@387
   286
kpeter@387
   287
  checkGraphNodeList(G, 4);
kpeter@387
   288
  checkGraphArcList(G, 4);
deba@228
   289
}
deba@228
   290
deba@228
   291
void checkConcepts() {
kpeter@171
   292
  { // Checking digraph components
deba@57
   293
    checkConcept<BaseDigraphComponent, BaseDigraphComponent >();
deba@57
   294
alpar@209
   295
    checkConcept<IDableDigraphComponent<>,
deba@57
   296
      IDableDigraphComponent<> >();
deba@57
   297
alpar@209
   298
    checkConcept<IterableDigraphComponent<>,
deba@57
   299
      IterableDigraphComponent<> >();
deba@57
   300
alpar@209
   301
    checkConcept<MappableDigraphComponent<>,
deba@57
   302
      MappableDigraphComponent<> >();
deba@57
   303
  }
kpeter@171
   304
  { // Checking skeleton digraph
deba@57
   305
    checkConcept<Digraph, Digraph>();
deba@57
   306
  }
kpeter@171
   307
  { // Checking ListDigraph
kpeter@171
   308
    checkConcept<Digraph, ListDigraph>();
deba@57
   309
    checkConcept<AlterableDigraphComponent<>, ListDigraph>();
deba@57
   310
    checkConcept<ExtendableDigraphComponent<>, ListDigraph>();
deba@57
   311
    checkConcept<ClearableDigraphComponent<>, ListDigraph>();
deba@57
   312
    checkConcept<ErasableDigraphComponent<>, ListDigraph>();
kpeter@171
   313
  }
kpeter@171
   314
  { // Checking SmartDigraph
kpeter@171
   315
    checkConcept<Digraph, SmartDigraph>();
kpeter@171
   316
    checkConcept<AlterableDigraphComponent<>, SmartDigraph>();
kpeter@171
   317
    checkConcept<ExtendableDigraphComponent<>, SmartDigraph>();
kpeter@171
   318
    checkConcept<ClearableDigraphComponent<>, SmartDigraph>();
kpeter@171
   319
  }
kpeter@366
   320
  { // Checking FullDigraph
kpeter@366
   321
    checkConcept<Digraph, FullDigraph>();
kpeter@366
   322
  }
kpeter@171
   323
}
deba@57
   324
kpeter@171
   325
template <typename Digraph>
deba@228
   326
void checkDigraphValidity() {
kpeter@171
   327
  TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
kpeter@171
   328
  Digraph g;
kpeter@171
   329
kpeter@171
   330
  Node
kpeter@171
   331
    n1 = g.addNode(),
kpeter@171
   332
    n2 = g.addNode(),
kpeter@171
   333
    n3 = g.addNode();
kpeter@171
   334
kpeter@171
   335
  Arc
kpeter@171
   336
    e1 = g.addArc(n1, n2),
kpeter@171
   337
    e2 = g.addArc(n2, n3);
kpeter@171
   338
kpeter@171
   339
  check(g.valid(n1), "Wrong validity check");
kpeter@171
   340
  check(g.valid(e1), "Wrong validity check");
kpeter@171
   341
kpeter@171
   342
  check(!g.valid(g.nodeFromId(-1)), "Wrong validity check");
kpeter@171
   343
  check(!g.valid(g.arcFromId(-1)), "Wrong validity check");
kpeter@171
   344
}
kpeter@171
   345
kpeter@171
   346
template <typename Digraph>
deba@228
   347
void checkDigraphValidityErase() {
kpeter@171
   348
  TEMPLATE_DIGRAPH_TYPEDEFS(Digraph);
kpeter@171
   349
  Digraph g;
kpeter@171
   350
kpeter@171
   351
  Node
kpeter@171
   352
    n1 = g.addNode(),
kpeter@171
   353
    n2 = g.addNode(),
kpeter@171
   354
    n3 = g.addNode();
kpeter@171
   355
kpeter@171
   356
  Arc
kpeter@171
   357
    e1 = g.addArc(n1, n2),
kpeter@171
   358
    e2 = g.addArc(n2, n3);
kpeter@171
   359
kpeter@171
   360
  check(g.valid(n1), "Wrong validity check");
kpeter@171
   361
  check(g.valid(e1), "Wrong validity check");
kpeter@171
   362
kpeter@171
   363
  g.erase(n1);
kpeter@171
   364
kpeter@171
   365
  check(!g.valid(n1), "Wrong validity check");
kpeter@171
   366
  check(g.valid(n2), "Wrong validity check");
kpeter@171
   367
  check(g.valid(n3), "Wrong validity check");
kpeter@171
   368
  check(!g.valid(e1), "Wrong validity check");
kpeter@171
   369
  check(g.valid(e2), "Wrong validity check");
kpeter@171
   370
kpeter@171
   371
  check(!g.valid(g.nodeFromId(-1)), "Wrong validity check");
kpeter@171
   372
  check(!g.valid(g.arcFromId(-1)), "Wrong validity check");
kpeter@171
   373
}
kpeter@171
   374
alpar@388
   375
void checkFullDigraph(int num) {
alpar@388
   376
  typedef FullDigraph Digraph;
alpar@388
   377
  DIGRAPH_TYPEDEFS(Digraph);
alpar@388
   378
  Digraph G(num);
alpar@388
   379
alpar@388
   380
  checkGraphNodeList(G, num);
alpar@388
   381
  checkGraphArcList(G, num * num);
alpar@388
   382
alpar@388
   383
  for (NodeIt n(G); n != INVALID; ++n) {
alpar@388
   384
    checkGraphOutArcList(G, n, num);
alpar@388
   385
    checkGraphInArcList(G, n, num);
alpar@388
   386
  }
alpar@388
   387
alpar@388
   388
  checkGraphConArcList(G, num * num);
alpar@388
   389
alpar@388
   390
  checkNodeIds(G);
alpar@388
   391
  checkArcIds(G);
alpar@388
   392
  checkGraphNodeMap(G);
alpar@388
   393
  checkGraphArcMap(G);
alpar@388
   394
alpar@388
   395
  for (int i = 0; i < G.nodeNum(); ++i) {
alpar@388
   396
    check(G.index(G(i)) == i, "Wrong index");
alpar@388
   397
  }
alpar@388
   398
alpar@388
   399
  for (NodeIt s(G); s != INVALID; ++s) {
alpar@388
   400
    for (NodeIt t(G); t != INVALID; ++t) {
alpar@388
   401
      Arc a = G.arc(s, t);
alpar@388
   402
      check(G.source(a) == s && G.target(a) == t, "Wrong arc lookup");
alpar@388
   403
    }
alpar@388
   404
  }
alpar@388
   405
}
alpar@388
   406
deba@228
   407
void checkDigraphs() {
kpeter@171
   408
  { // Checking ListDigraph
kpeter@387
   409
    checkDigraphBuild<ListDigraph>();
kpeter@387
   410
    checkDigraphSplit<ListDigraph>();
kpeter@387
   411
    checkDigraphAlter<ListDigraph>();
kpeter@387
   412
    checkDigraphErase<ListDigraph>();
kpeter@387
   413
    checkDigraphSnapshot<ListDigraph>();
deba@228
   414
    checkDigraphValidityErase<ListDigraph>();
deba@57
   415
  }
kpeter@171
   416
  { // Checking SmartDigraph
kpeter@387
   417
    checkDigraphBuild<SmartDigraph>();
kpeter@387
   418
    checkDigraphSplit<SmartDigraph>();
kpeter@387
   419
    checkDigraphSnapshot<SmartDigraph>();
deba@228
   420
    checkDigraphValidity<SmartDigraph>();
kpeter@171
   421
  }
deba@365
   422
  { // Checking FullDigraph
deba@365
   423
    checkFullDigraph(8);
deba@365
   424
  }
kpeter@171
   425
}
deba@57
   426
kpeter@171
   427
int main() {
deba@228
   428
  checkDigraphs();
deba@228
   429
  checkConcepts();
deba@57
   430
  return 0;
deba@57
   431
}