test/maps_test.cc
author Alpar Juttner <alpar@cs.elte.hu>
Tue, 05 Jun 2012 12:06:28 +0200
changeset 1222 c40a9d94442d
parent 956 141f9c0db4a3
child 1084 d303bfa8b1ed
child 1159 7fdaa05a69a1
permissions -rw-r--r--
New features in time_measure.h (#442)
     1 /* -*- mode: C++; indent-tabs-mode: nil; -*-
     2  *
     3  * This file is a part of LEMON, a generic C++ optimization library.
     4  *
     5  * Copyright (C) 2003-2010
     6  * Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport
     7  * (Egervary Research Group on Combinatorial Optimization, EGRES).
     8  *
     9  * Permission to use, modify and distribute this software is granted
    10  * provided that this copyright notice appears in all copies. For
    11  * precise terms see the accompanying LICENSE file.
    12  *
    13  * This software is provided "AS IS" with no warranty of any kind,
    14  * express or implied, and with no claim as to its suitability for any
    15  * purpose.
    16  *
    17  */
    18 
    19 #include <deque>
    20 #include <set>
    21 
    22 #include <lemon/concept_check.h>
    23 #include <lemon/concepts/maps.h>
    24 #include <lemon/maps.h>
    25 #include <lemon/list_graph.h>
    26 #include <lemon/smart_graph.h>
    27 #include <lemon/adaptors.h>
    28 #include <lemon/dfs.h>
    29 #include <algorithm>
    30 
    31 #include "test_tools.h"
    32 
    33 using namespace lemon;
    34 using namespace lemon::concepts;
    35 
    36 struct A {};
    37 inline bool operator<(A, A) { return true; }
    38 struct B {};
    39 
    40 class C {
    41   int _x;
    42 public:
    43   C(int x) : _x(x) {}
    44   int get() const { return _x; }
    45 };
    46 inline bool operator<(C c1, C c2) { return c1.get() < c2.get(); }
    47 inline bool operator==(C c1, C c2) { return c1.get() == c2.get(); }
    48 
    49 C createC(int x) { return C(x); }
    50 
    51 template <typename T>
    52 class Less {
    53   T _t;
    54 public:
    55   Less(T t): _t(t) {}
    56   bool operator()(const T& t) const { return t < _t; }
    57 };
    58 
    59 class F {
    60 public:
    61   typedef A argument_type;
    62   typedef B result_type;
    63 
    64   B operator()(const A&) const { return B(); }
    65 private:
    66   F& operator=(const F&);
    67 };
    68 
    69 int func(A) { return 3; }
    70 
    71 int binc(int a, B) { return a+1; }
    72 
    73 template <typename T>
    74 class Sum {
    75   T& _sum;
    76 public:
    77   Sum(T& sum) : _sum(sum) {}
    78   void operator()(const T& t) { _sum += t; }
    79 };
    80 
    81 typedef ReadMap<A, double> DoubleMap;
    82 typedef ReadWriteMap<A, double> DoubleWriteMap;
    83 typedef ReferenceMap<A, double, double&, const double&> DoubleRefMap;
    84 
    85 typedef ReadMap<A, bool> BoolMap;
    86 typedef ReadWriteMap<A, bool> BoolWriteMap;
    87 typedef ReferenceMap<A, bool, bool&, const bool&> BoolRefMap;
    88 
    89 int main()
    90 {
    91   // Map concepts
    92   checkConcept<ReadMap<A,B>, ReadMap<A,B> >();
    93   checkConcept<ReadMap<A,C>, ReadMap<A,C> >();
    94   checkConcept<WriteMap<A,B>, WriteMap<A,B> >();
    95   checkConcept<WriteMap<A,C>, WriteMap<A,C> >();
    96   checkConcept<ReadWriteMap<A,B>, ReadWriteMap<A,B> >();
    97   checkConcept<ReadWriteMap<A,C>, ReadWriteMap<A,C> >();
    98   checkConcept<ReferenceMap<A,B,B&,const B&>, ReferenceMap<A,B,B&,const B&> >();
    99   checkConcept<ReferenceMap<A,C,C&,const C&>, ReferenceMap<A,C,C&,const C&> >();
   100 
   101   // NullMap
   102   {
   103     checkConcept<ReadWriteMap<A,B>, NullMap<A,B> >();
   104     NullMap<A,B> map1;
   105     NullMap<A,B> map2 = map1;
   106     map1 = nullMap<A,B>();
   107   }
   108 
   109   // ConstMap
   110   {
   111     checkConcept<ReadWriteMap<A,B>, ConstMap<A,B> >();
   112     checkConcept<ReadWriteMap<A,C>, ConstMap<A,C> >();
   113     ConstMap<A,B> map1;
   114     ConstMap<A,B> map2 = B();
   115     ConstMap<A,B> map3 = map1;
   116     map1 = constMap<A>(B());
   117     map1 = constMap<A,B>();
   118     map1.setAll(B());
   119     ConstMap<A,C> map4(C(1));
   120     ConstMap<A,C> map5 = map4;
   121     map4 = constMap<A>(C(2));
   122     map4.setAll(C(3));
   123 
   124     checkConcept<ReadWriteMap<A,int>, ConstMap<A,int> >();
   125     check(constMap<A>(10)[A()] == 10, "Something is wrong with ConstMap");
   126 
   127     checkConcept<ReadWriteMap<A,int>, ConstMap<A,Const<int,10> > >();
   128     ConstMap<A,Const<int,10> > map6;
   129     ConstMap<A,Const<int,10> > map7 = map6;
   130     map6 = constMap<A,int,10>();
   131     map7 = constMap<A,Const<int,10> >();
   132     check(map6[A()] == 10 && map7[A()] == 10,
   133           "Something is wrong with ConstMap");
   134   }
   135 
   136   // IdentityMap
   137   {
   138     checkConcept<ReadMap<A,A>, IdentityMap<A> >();
   139     IdentityMap<A> map1;
   140     IdentityMap<A> map2 = map1;
   141     map1 = identityMap<A>();
   142 
   143     checkConcept<ReadMap<double,double>, IdentityMap<double> >();
   144     check(identityMap<double>()[1.0] == 1.0 &&
   145           identityMap<double>()[3.14] == 3.14,
   146           "Something is wrong with IdentityMap");
   147   }
   148 
   149   // RangeMap
   150   {
   151     checkConcept<ReferenceMap<int,B,B&,const B&>, RangeMap<B> >();
   152     RangeMap<B> map1;
   153     RangeMap<B> map2(10);
   154     RangeMap<B> map3(10,B());
   155     RangeMap<B> map4 = map1;
   156     RangeMap<B> map5 = rangeMap<B>();
   157     RangeMap<B> map6 = rangeMap<B>(10);
   158     RangeMap<B> map7 = rangeMap(10,B());
   159 
   160     checkConcept< ReferenceMap<int, double, double&, const double&>,
   161                   RangeMap<double> >();
   162     std::vector<double> v(10, 0);
   163     v[5] = 100;
   164     RangeMap<double> map8(v);
   165     RangeMap<double> map9 = rangeMap(v);
   166     check(map9.size() == 10 && map9[2] == 0 && map9[5] == 100,
   167           "Something is wrong with RangeMap");
   168   }
   169 
   170   // SparseMap
   171   {
   172     checkConcept<ReferenceMap<A,B,B&,const B&>, SparseMap<A,B> >();
   173     SparseMap<A,B> map1;
   174     SparseMap<A,B> map2 = B();
   175     SparseMap<A,B> map3 = sparseMap<A,B>();
   176     SparseMap<A,B> map4 = sparseMap<A>(B());
   177 
   178     checkConcept< ReferenceMap<double, int, int&, const int&>,
   179                   SparseMap<double, int> >();
   180     std::map<double, int> m;
   181     SparseMap<double, int> map5(m);
   182     SparseMap<double, int> map6(m,10);
   183     SparseMap<double, int> map7 = sparseMap(m);
   184     SparseMap<double, int> map8 = sparseMap(m,10);
   185 
   186     check(map5[1.0] == 0 && map5[3.14] == 0 &&
   187           map6[1.0] == 10 && map6[3.14] == 10,
   188           "Something is wrong with SparseMap");
   189     map5[1.0] = map6[3.14] = 100;
   190     check(map5[1.0] == 100 && map5[3.14] == 0 &&
   191           map6[1.0] == 10 && map6[3.14] == 100,
   192           "Something is wrong with SparseMap");
   193   }
   194 
   195   // ComposeMap
   196   {
   197     typedef ComposeMap<DoubleMap, ReadMap<B,A> > CompMap;
   198     checkConcept<ReadMap<B,double>, CompMap>();
   199     CompMap map1 = CompMap(DoubleMap(),ReadMap<B,A>());
   200     CompMap map2 = composeMap(DoubleMap(), ReadMap<B,A>());
   201 
   202     SparseMap<double, bool> m1(false); m1[3.14] = true;
   203     RangeMap<double> m2(2); m2[0] = 3.0; m2[1] = 3.14;
   204     check(!composeMap(m1,m2)[0] && composeMap(m1,m2)[1],
   205           "Something is wrong with ComposeMap")
   206   }
   207 
   208   // CombineMap
   209   {
   210     typedef CombineMap<DoubleMap, DoubleMap, std::plus<double> > CombMap;
   211     checkConcept<ReadMap<A,double>, CombMap>();
   212     CombMap map1 = CombMap(DoubleMap(), DoubleMap());
   213     CombMap map2 = combineMap(DoubleMap(), DoubleMap(), std::plus<double>());
   214 
   215     check(combineMap(constMap<B,int,2>(), identityMap<B>(), &binc)[B()] == 3,
   216           "Something is wrong with CombineMap");
   217   }
   218 
   219   // FunctorToMap, MapToFunctor
   220   {
   221     checkConcept<ReadMap<A,B>, FunctorToMap<F,A,B> >();
   222     checkConcept<ReadMap<A,B>, FunctorToMap<F> >();
   223     FunctorToMap<F> map1;
   224     FunctorToMap<F> map2 = FunctorToMap<F>(F());
   225     B b = functorToMap(F())[A()];
   226 
   227     checkConcept<ReadMap<A,B>, MapToFunctor<ReadMap<A,B> > >();
   228     MapToFunctor<ReadMap<A,B> > map =
   229       MapToFunctor<ReadMap<A,B> >(ReadMap<A,B>());
   230 
   231     check(functorToMap(&func)[A()] == 3,
   232           "Something is wrong with FunctorToMap");
   233     check(mapToFunctor(constMap<A,int>(2))(A()) == 2,
   234           "Something is wrong with MapToFunctor");
   235     check(mapToFunctor(functorToMap(&func))(A()) == 3 &&
   236           mapToFunctor(functorToMap(&func))[A()] == 3,
   237           "Something is wrong with FunctorToMap or MapToFunctor");
   238     check(functorToMap(mapToFunctor(constMap<A,int>(2)))[A()] == 2,
   239           "Something is wrong with FunctorToMap or MapToFunctor");
   240   }
   241 
   242   // ConvertMap
   243   {
   244     checkConcept<ReadMap<double,double>,
   245       ConvertMap<ReadMap<double, int>, double> >();
   246     ConvertMap<RangeMap<bool>, int> map1(rangeMap(1, true));
   247     ConvertMap<RangeMap<bool>, int> map2 = convertMap<int>(rangeMap(2, false));
   248   }
   249 
   250   // ForkMap
   251   {
   252     checkConcept<DoubleWriteMap, ForkMap<DoubleWriteMap, DoubleWriteMap> >();
   253 
   254     typedef RangeMap<double> RM;
   255     typedef SparseMap<int, double> SM;
   256     RM m1(10, -1);
   257     SM m2(-1);
   258     checkConcept<ReadWriteMap<int, double>, ForkMap<RM, SM> >();
   259     checkConcept<ReadWriteMap<int, double>, ForkMap<SM, RM> >();
   260     ForkMap<RM, SM> map1(m1,m2);
   261     ForkMap<SM, RM> map2 = forkMap(m2,m1);
   262     map2.set(5, 10);
   263     check(m1[1] == -1 && m1[5] == 10 && m2[1] == -1 &&
   264           m2[5] == 10 && map2[1] == -1 && map2[5] == 10,
   265           "Something is wrong with ForkMap");
   266   }
   267 
   268   // Arithmetic maps:
   269   // - AddMap, SubMap, MulMap, DivMap
   270   // - ShiftMap, ShiftWriteMap, ScaleMap, ScaleWriteMap
   271   // - NegMap, NegWriteMap, AbsMap
   272   {
   273     checkConcept<DoubleMap, AddMap<DoubleMap,DoubleMap> >();
   274     checkConcept<DoubleMap, SubMap<DoubleMap,DoubleMap> >();
   275     checkConcept<DoubleMap, MulMap<DoubleMap,DoubleMap> >();
   276     checkConcept<DoubleMap, DivMap<DoubleMap,DoubleMap> >();
   277 
   278     ConstMap<int, double> c1(1.0), c2(3.14);
   279     IdentityMap<int> im;
   280     ConvertMap<IdentityMap<int>, double> id(im);
   281     check(addMap(c1,id)[0] == 1.0  && addMap(c1,id)[10] == 11.0,
   282           "Something is wrong with AddMap");
   283     check(subMap(id,c1)[0] == -1.0 && subMap(id,c1)[10] == 9.0,
   284           "Something is wrong with SubMap");
   285     check(mulMap(id,c2)[0] == 0    && mulMap(id,c2)[2]  == 6.28,
   286           "Something is wrong with MulMap");
   287     check(divMap(c2,id)[1] == 3.14 && divMap(c2,id)[2]  == 1.57,
   288           "Something is wrong with DivMap");
   289 
   290     checkConcept<DoubleMap, ShiftMap<DoubleMap> >();
   291     checkConcept<DoubleWriteMap, ShiftWriteMap<DoubleWriteMap> >();
   292     checkConcept<DoubleMap, ScaleMap<DoubleMap> >();
   293     checkConcept<DoubleWriteMap, ScaleWriteMap<DoubleWriteMap> >();
   294     checkConcept<DoubleMap, NegMap<DoubleMap> >();
   295     checkConcept<DoubleWriteMap, NegWriteMap<DoubleWriteMap> >();
   296     checkConcept<DoubleMap, AbsMap<DoubleMap> >();
   297 
   298     check(shiftMap(id, 2.0)[1] == 3.0 && shiftMap(id, 2.0)[10] == 12.0,
   299           "Something is wrong with ShiftMap");
   300     check(shiftWriteMap(id, 2.0)[1] == 3.0 &&
   301           shiftWriteMap(id, 2.0)[10] == 12.0,
   302           "Something is wrong with ShiftWriteMap");
   303     check(scaleMap(id, 2.0)[1] == 2.0 && scaleMap(id, 2.0)[10] == 20.0,
   304           "Something is wrong with ScaleMap");
   305     check(scaleWriteMap(id, 2.0)[1] == 2.0 &&
   306           scaleWriteMap(id, 2.0)[10] == 20.0,
   307           "Something is wrong with ScaleWriteMap");
   308     check(negMap(id)[1] == -1.0 && negMap(id)[-10] == 10.0,
   309           "Something is wrong with NegMap");
   310     check(negWriteMap(id)[1] == -1.0 && negWriteMap(id)[-10] == 10.0,
   311           "Something is wrong with NegWriteMap");
   312     check(absMap(id)[1] == 1.0 && absMap(id)[-10] == 10.0,
   313           "Something is wrong with AbsMap");
   314   }
   315 
   316   // Logical maps:
   317   // - TrueMap, FalseMap
   318   // - AndMap, OrMap
   319   // - NotMap, NotWriteMap
   320   // - EqualMap, LessMap
   321   {
   322     checkConcept<BoolMap, TrueMap<A> >();
   323     checkConcept<BoolMap, FalseMap<A> >();
   324     checkConcept<BoolMap, AndMap<BoolMap,BoolMap> >();
   325     checkConcept<BoolMap, OrMap<BoolMap,BoolMap> >();
   326     checkConcept<BoolMap, NotMap<BoolMap> >();
   327     checkConcept<BoolWriteMap, NotWriteMap<BoolWriteMap> >();
   328     checkConcept<BoolMap, EqualMap<DoubleMap,DoubleMap> >();
   329     checkConcept<BoolMap, LessMap<DoubleMap,DoubleMap> >();
   330 
   331     TrueMap<int> tm;
   332     FalseMap<int> fm;
   333     RangeMap<bool> rm(2);
   334     rm[0] = true; rm[1] = false;
   335     check(andMap(tm,rm)[0] && !andMap(tm,rm)[1] &&
   336           !andMap(fm,rm)[0] && !andMap(fm,rm)[1],
   337           "Something is wrong with AndMap");
   338     check(orMap(tm,rm)[0] && orMap(tm,rm)[1] &&
   339           orMap(fm,rm)[0] && !orMap(fm,rm)[1],
   340           "Something is wrong with OrMap");
   341     check(!notMap(rm)[0] && notMap(rm)[1],
   342           "Something is wrong with NotMap");
   343     check(!notWriteMap(rm)[0] && notWriteMap(rm)[1],
   344           "Something is wrong with NotWriteMap");
   345 
   346     ConstMap<int, double> cm(2.0);
   347     IdentityMap<int> im;
   348     ConvertMap<IdentityMap<int>, double> id(im);
   349     check(lessMap(id,cm)[1] && !lessMap(id,cm)[2] && !lessMap(id,cm)[3],
   350           "Something is wrong with LessMap");
   351     check(!equalMap(id,cm)[1] && equalMap(id,cm)[2] && !equalMap(id,cm)[3],
   352           "Something is wrong with EqualMap");
   353   }
   354 
   355   // LoggerBoolMap
   356   {
   357     typedef std::vector<int> vec;
   358     checkConcept<WriteMap<int, bool>, LoggerBoolMap<vec::iterator> >();
   359     checkConcept<WriteMap<int, bool>,
   360                  LoggerBoolMap<std::back_insert_iterator<vec> > >();
   361 
   362     vec v1;
   363     vec v2(10);
   364     LoggerBoolMap<std::back_insert_iterator<vec> >
   365       map1(std::back_inserter(v1));
   366     LoggerBoolMap<vec::iterator> map2(v2.begin());
   367     map1.set(10, false);
   368     map1.set(20, true);   map2.set(20, true);
   369     map1.set(30, false);  map2.set(40, false);
   370     map1.set(50, true);   map2.set(50, true);
   371     map1.set(60, true);   map2.set(60, true);
   372     check(v1.size() == 3 && v2.size() == 10 &&
   373           v1[0]==20 && v1[1]==50 && v1[2]==60 &&
   374           v2[0]==20 && v2[1]==50 && v2[2]==60,
   375           "Something is wrong with LoggerBoolMap");
   376 
   377     int i = 0;
   378     for ( LoggerBoolMap<vec::iterator>::Iterator it = map2.begin();
   379           it != map2.end(); ++it )
   380       check(v1[i++] == *it, "Something is wrong with LoggerBoolMap");
   381 
   382     typedef ListDigraph Graph;
   383     DIGRAPH_TYPEDEFS(Graph);
   384     Graph gr;
   385 
   386     Node n0 = gr.addNode();
   387     Node n1 = gr.addNode();
   388     Node n2 = gr.addNode();
   389     Node n3 = gr.addNode();
   390 
   391     gr.addArc(n3, n0);
   392     gr.addArc(n3, n2);
   393     gr.addArc(n0, n2);
   394     gr.addArc(n2, n1);
   395     gr.addArc(n0, n1);
   396 
   397     {
   398       std::vector<Node> v;
   399       dfs(gr).processedMap(loggerBoolMap(std::back_inserter(v))).run();
   400 
   401       check(v.size()==4 && v[0]==n1 && v[1]==n2 && v[2]==n0 && v[3]==n3,
   402             "Something is wrong with LoggerBoolMap");
   403     }
   404     {
   405       std::vector<Node> v(countNodes(gr));
   406       dfs(gr).processedMap(loggerBoolMap(v.begin())).run();
   407 
   408       check(v.size()==4 && v[0]==n1 && v[1]==n2 && v[2]==n0 && v[3]==n3,
   409             "Something is wrong with LoggerBoolMap");
   410     }
   411   }
   412 
   413   // IdMap, RangeIdMap
   414   {
   415     typedef ListDigraph Graph;
   416     DIGRAPH_TYPEDEFS(Graph);
   417 
   418     checkConcept<ReadMap<Node, int>, IdMap<Graph, Node> >();
   419     checkConcept<ReadMap<Arc, int>, IdMap<Graph, Arc> >();
   420     checkConcept<ReadMap<Node, int>, RangeIdMap<Graph, Node> >();
   421     checkConcept<ReadMap<Arc, int>, RangeIdMap<Graph, Arc> >();
   422 
   423     Graph gr;
   424     IdMap<Graph, Node> nmap(gr);
   425     IdMap<Graph, Arc> amap(gr);
   426     RangeIdMap<Graph, Node> nrmap(gr);
   427     RangeIdMap<Graph, Arc> armap(gr);
   428 
   429     Node n0 = gr.addNode();
   430     Node n1 = gr.addNode();
   431     Node n2 = gr.addNode();
   432 
   433     Arc a0 = gr.addArc(n0, n1);
   434     Arc a1 = gr.addArc(n0, n2);
   435     Arc a2 = gr.addArc(n2, n1);
   436     Arc a3 = gr.addArc(n2, n0);
   437 
   438     check(nmap[n0] == gr.id(n0) && nmap(gr.id(n0)) == n0, "Wrong IdMap");
   439     check(nmap[n1] == gr.id(n1) && nmap(gr.id(n1)) == n1, "Wrong IdMap");
   440     check(nmap[n2] == gr.id(n2) && nmap(gr.id(n2)) == n2, "Wrong IdMap");
   441 
   442     check(amap[a0] == gr.id(a0) && amap(gr.id(a0)) == a0, "Wrong IdMap");
   443     check(amap[a1] == gr.id(a1) && amap(gr.id(a1)) == a1, "Wrong IdMap");
   444     check(amap[a2] == gr.id(a2) && amap(gr.id(a2)) == a2, "Wrong IdMap");
   445     check(amap[a3] == gr.id(a3) && amap(gr.id(a3)) == a3, "Wrong IdMap");
   446 
   447     check(nmap.inverse()[gr.id(n0)] == n0, "Wrong IdMap::InverseMap");
   448     check(amap.inverse()[gr.id(a0)] == a0, "Wrong IdMap::InverseMap");
   449 
   450     check(nrmap.size() == 3 && armap.size() == 4,
   451           "Wrong RangeIdMap::size()");
   452 
   453     check(nrmap[n0] == 0 && nrmap(0) == n0, "Wrong RangeIdMap");
   454     check(nrmap[n1] == 1 && nrmap(1) == n1, "Wrong RangeIdMap");
   455     check(nrmap[n2] == 2 && nrmap(2) == n2, "Wrong RangeIdMap");
   456 
   457     check(armap[a0] == 0 && armap(0) == a0, "Wrong RangeIdMap");
   458     check(armap[a1] == 1 && armap(1) == a1, "Wrong RangeIdMap");
   459     check(armap[a2] == 2 && armap(2) == a2, "Wrong RangeIdMap");
   460     check(armap[a3] == 3 && armap(3) == a3, "Wrong RangeIdMap");
   461 
   462     check(nrmap.inverse()[0] == n0, "Wrong RangeIdMap::InverseMap");
   463     check(armap.inverse()[0] == a0, "Wrong RangeIdMap::InverseMap");
   464 
   465     gr.erase(n1);
   466 
   467     if (nrmap[n0] == 1) nrmap.swap(n0, n2);
   468     nrmap.swap(n2, n0);
   469     if (armap[a1] == 1) armap.swap(a1, a3);
   470     armap.swap(a3, a1);
   471 
   472     check(nrmap.size() == 2 && armap.size() == 2,
   473           "Wrong RangeIdMap::size()");
   474 
   475     check(nrmap[n0] == 1 && nrmap(1) == n0, "Wrong RangeIdMap");
   476     check(nrmap[n2] == 0 && nrmap(0) == n2, "Wrong RangeIdMap");
   477 
   478     check(armap[a1] == 1 && armap(1) == a1, "Wrong RangeIdMap");
   479     check(armap[a3] == 0 && armap(0) == a3, "Wrong RangeIdMap");
   480 
   481     check(nrmap.inverse()[0] == n2, "Wrong RangeIdMap::InverseMap");
   482     check(armap.inverse()[0] == a3, "Wrong RangeIdMap::InverseMap");
   483   }
   484 
   485   // SourceMap, TargetMap, ForwardMap, BackwardMap, InDegMap, OutDegMap
   486   {
   487     typedef ListGraph Graph;
   488     GRAPH_TYPEDEFS(Graph);
   489 
   490     checkConcept<ReadMap<Arc, Node>, SourceMap<Graph> >();
   491     checkConcept<ReadMap<Arc, Node>, TargetMap<Graph> >();
   492     checkConcept<ReadMap<Edge, Arc>, ForwardMap<Graph> >();
   493     checkConcept<ReadMap<Edge, Arc>, BackwardMap<Graph> >();
   494     checkConcept<ReadMap<Node, int>, InDegMap<Graph> >();
   495     checkConcept<ReadMap<Node, int>, OutDegMap<Graph> >();
   496 
   497     Graph gr;
   498     Node n0 = gr.addNode();
   499     Node n1 = gr.addNode();
   500     Node n2 = gr.addNode();
   501 
   502     gr.addEdge(n0,n1);
   503     gr.addEdge(n1,n2);
   504     gr.addEdge(n0,n2);
   505     gr.addEdge(n2,n1);
   506     gr.addEdge(n1,n2);
   507     gr.addEdge(n0,n1);
   508 
   509     for (EdgeIt e(gr); e != INVALID; ++e) {
   510       check(forwardMap(gr)[e] == gr.direct(e, true), "Wrong ForwardMap");
   511       check(backwardMap(gr)[e] == gr.direct(e, false), "Wrong BackwardMap");
   512     }
   513 
   514     check(mapCompare(gr,
   515           sourceMap(orienter(gr, constMap<Edge, bool>(true))),
   516           targetMap(orienter(gr, constMap<Edge, bool>(false)))),
   517           "Wrong SourceMap or TargetMap");
   518 
   519     typedef Orienter<Graph, const ConstMap<Edge, bool> > Digraph;
   520     Digraph dgr(gr, constMap<Edge, bool>(true));
   521     OutDegMap<Digraph> odm(dgr);
   522     InDegMap<Digraph> idm(dgr);
   523 
   524     check(odm[n0] == 3 && odm[n1] == 2 && odm[n2] == 1, "Wrong OutDegMap");
   525     check(idm[n0] == 0 && idm[n1] == 3 && idm[n2] == 3, "Wrong InDegMap");
   526 
   527     gr.addEdge(n2, n0);
   528 
   529     check(odm[n0] == 3 && odm[n1] == 2 && odm[n2] == 2, "Wrong OutDegMap");
   530     check(idm[n0] == 1 && idm[n1] == 3 && idm[n2] == 3, "Wrong InDegMap");
   531   }
   532 
   533   // CrossRefMap
   534   {
   535     typedef ListDigraph Graph;
   536     DIGRAPH_TYPEDEFS(Graph);
   537 
   538     checkConcept<ReadWriteMap<Node, int>,
   539                  CrossRefMap<Graph, Node, int> >();
   540     checkConcept<ReadWriteMap<Node, bool>,
   541                  CrossRefMap<Graph, Node, bool> >();
   542     checkConcept<ReadWriteMap<Node, double>,
   543                  CrossRefMap<Graph, Node, double> >();
   544 
   545     Graph gr;
   546     typedef CrossRefMap<Graph, Node, char> CRMap;
   547     CRMap map(gr);
   548 
   549     Node n0 = gr.addNode();
   550     Node n1 = gr.addNode();
   551     Node n2 = gr.addNode();
   552 
   553     map.set(n0, 'A');
   554     map.set(n1, 'B');
   555     map.set(n2, 'C');
   556 
   557     check(map[n0] == 'A' && map('A') == n0 && map.inverse()['A'] == n0,
   558           "Wrong CrossRefMap");
   559     check(map[n1] == 'B' && map('B') == n1 && map.inverse()['B'] == n1,
   560           "Wrong CrossRefMap");
   561     check(map[n2] == 'C' && map('C') == n2 && map.inverse()['C'] == n2,
   562           "Wrong CrossRefMap");
   563     check(map.count('A') == 1 && map.count('B') == 1 && map.count('C') == 1,
   564           "Wrong CrossRefMap::count()");
   565 
   566     CRMap::ValueIt it = map.beginValue();
   567     check(*it++ == 'A' && *it++ == 'B' && *it++ == 'C' &&
   568           it == map.endValue(), "Wrong value iterator");
   569 
   570     map.set(n2, 'A');
   571 
   572     check(map[n0] == 'A' && map[n1] == 'B' && map[n2] == 'A',
   573           "Wrong CrossRefMap");
   574     check(map('A') == n0 && map.inverse()['A'] == n0, "Wrong CrossRefMap");
   575     check(map('B') == n1 && map.inverse()['B'] == n1, "Wrong CrossRefMap");
   576     check(map('C') == INVALID && map.inverse()['C'] == INVALID,
   577           "Wrong CrossRefMap");
   578     check(map.count('A') == 2 && map.count('B') == 1 && map.count('C') == 0,
   579           "Wrong CrossRefMap::count()");
   580 
   581     it = map.beginValue();
   582     check(*it++ == 'A' && *it++ == 'A' && *it++ == 'B' &&
   583           it == map.endValue(), "Wrong value iterator");
   584 
   585     map.set(n0, 'C');
   586 
   587     check(map[n0] == 'C' && map[n1] == 'B' && map[n2] == 'A',
   588           "Wrong CrossRefMap");
   589     check(map('A') == n2 && map.inverse()['A'] == n2, "Wrong CrossRefMap");
   590     check(map('B') == n1 && map.inverse()['B'] == n1, "Wrong CrossRefMap");
   591     check(map('C') == n0 && map.inverse()['C'] == n0, "Wrong CrossRefMap");
   592     check(map.count('A') == 1 && map.count('B') == 1 && map.count('C') == 1,
   593           "Wrong CrossRefMap::count()");
   594 
   595     it = map.beginValue();
   596     check(*it++ == 'A' && *it++ == 'B' && *it++ == 'C' &&
   597           it == map.endValue(), "Wrong value iterator");
   598   }
   599 
   600   // CrossRefMap
   601   {
   602     typedef SmartDigraph Graph;
   603     DIGRAPH_TYPEDEFS(Graph);
   604 
   605     checkConcept<ReadWriteMap<Node, int>,
   606                  CrossRefMap<Graph, Node, int> >();
   607 
   608     Graph gr;
   609     typedef CrossRefMap<Graph, Node, char> CRMap;
   610     typedef CRMap::ValueIterator ValueIt;
   611     CRMap map(gr);
   612 
   613     Node n0 = gr.addNode();
   614     Node n1 = gr.addNode();
   615     Node n2 = gr.addNode();
   616 
   617     map.set(n0, 'A');
   618     map.set(n1, 'B');
   619     map.set(n2, 'C');
   620     map.set(n2, 'A');
   621     map.set(n0, 'C');
   622 
   623     check(map[n0] == 'C' && map[n1] == 'B' && map[n2] == 'A',
   624           "Wrong CrossRefMap");
   625     check(map('A') == n2 && map.inverse()['A'] == n2, "Wrong CrossRefMap");
   626     check(map('B') == n1 && map.inverse()['B'] == n1, "Wrong CrossRefMap");
   627     check(map('C') == n0 && map.inverse()['C'] == n0, "Wrong CrossRefMap");
   628 
   629     ValueIt it = map.beginValue();
   630     check(*it++ == 'A' && *it++ == 'B' && *it++ == 'C' &&
   631           it == map.endValue(), "Wrong value iterator");
   632   }
   633 
   634   // Iterable bool map
   635   {
   636     typedef SmartGraph Graph;
   637     typedef SmartGraph::Node Item;
   638 
   639     typedef IterableBoolMap<SmartGraph, SmartGraph::Node> Ibm;
   640     checkConcept<ReferenceMap<Item, bool, bool&, const bool&>, Ibm>();
   641 
   642     const int num = 10;
   643     Graph g;
   644     Ibm map0(g, true);
   645     std::vector<Item> items;
   646     for (int i = 0; i < num; ++i) {
   647       items.push_back(g.addNode());
   648     }
   649 
   650     Ibm map1(g, true);
   651     int n = 0;
   652     for (Ibm::TrueIt it(map1); it != INVALID; ++it) {
   653       check(map1[static_cast<Item>(it)], "Wrong TrueIt");
   654       ++n;
   655     }
   656     check(n == num, "Wrong number");
   657 
   658     n = 0;
   659     for (Ibm::ItemIt it(map1, true); it != INVALID; ++it) {
   660         check(map1[static_cast<Item>(it)], "Wrong ItemIt for true");
   661         ++n;
   662     }
   663     check(n == num, "Wrong number");
   664     check(Ibm::FalseIt(map1) == INVALID, "Wrong FalseIt");
   665     check(Ibm::ItemIt(map1, false) == INVALID, "Wrong ItemIt for false");
   666 
   667     map1[items[5]] = true;
   668 
   669     n = 0;
   670     for (Ibm::ItemIt it(map1, true); it != INVALID; ++it) {
   671         check(map1[static_cast<Item>(it)], "Wrong ItemIt for true");
   672         ++n;
   673     }
   674     check(n == num, "Wrong number");
   675 
   676     map1[items[num / 2]] = false;
   677     check(map1[items[num / 2]] == false, "Wrong map value");
   678 
   679     n = 0;
   680     for (Ibm::TrueIt it(map1); it != INVALID; ++it) {
   681         check(map1[static_cast<Item>(it)], "Wrong TrueIt for true");
   682         ++n;
   683     }
   684     check(n == num - 1, "Wrong number");
   685 
   686     n = 0;
   687     for (Ibm::FalseIt it(map1); it != INVALID; ++it) {
   688         check(!map1[static_cast<Item>(it)], "Wrong FalseIt for true");
   689         ++n;
   690     }
   691     check(n == 1, "Wrong number");
   692 
   693     map1[items[0]] = false;
   694     check(map1[items[0]] == false, "Wrong map value");
   695 
   696     map1[items[num - 1]] = false;
   697     check(map1[items[num - 1]] == false, "Wrong map value");
   698 
   699     n = 0;
   700     for (Ibm::TrueIt it(map1); it != INVALID; ++it) {
   701         check(map1[static_cast<Item>(it)], "Wrong TrueIt for true");
   702         ++n;
   703     }
   704     check(n == num - 3, "Wrong number");
   705     check(map1.trueNum() == num - 3, "Wrong number");
   706 
   707     n = 0;
   708     for (Ibm::FalseIt it(map1); it != INVALID; ++it) {
   709         check(!map1[static_cast<Item>(it)], "Wrong FalseIt for true");
   710         ++n;
   711     }
   712     check(n == 3, "Wrong number");
   713     check(map1.falseNum() == 3, "Wrong number");
   714   }
   715 
   716   // Iterable int map
   717   {
   718     typedef SmartGraph Graph;
   719     typedef SmartGraph::Node Item;
   720     typedef IterableIntMap<SmartGraph, SmartGraph::Node> Iim;
   721 
   722     checkConcept<ReferenceMap<Item, int, int&, const int&>, Iim>();
   723 
   724     const int num = 10;
   725     Graph g;
   726     Iim map0(g, 0);
   727     std::vector<Item> items;
   728     for (int i = 0; i < num; ++i) {
   729       items.push_back(g.addNode());
   730     }
   731 
   732     Iim map1(g);
   733     check(map1.size() == 0, "Wrong size");
   734 
   735     for (int i = 0; i < num; ++i) {
   736       map1[items[i]] = i;
   737     }
   738     check(map1.size() == num, "Wrong size");
   739 
   740     for (int i = 0; i < num; ++i) {
   741       Iim::ItemIt it(map1, i);
   742       check(static_cast<Item>(it) == items[i], "Wrong value");
   743       ++it;
   744       check(static_cast<Item>(it) == INVALID, "Wrong value");
   745     }
   746 
   747     for (int i = 0; i < num; ++i) {
   748       map1[items[i]] = i % 2;
   749     }
   750     check(map1.size() == 2, "Wrong size");
   751 
   752     int n = 0;
   753     for (Iim::ItemIt it(map1, 0); it != INVALID; ++it) {
   754       check(map1[static_cast<Item>(it)] == 0, "Wrong value");
   755       ++n;
   756     }
   757     check(n == (num + 1) / 2, "Wrong number");
   758 
   759     for (Iim::ItemIt it(map1, 1); it != INVALID; ++it) {
   760       check(map1[static_cast<Item>(it)] == 1, "Wrong value");
   761       ++n;
   762     }
   763     check(n == num, "Wrong number");
   764 
   765   }
   766 
   767   // Iterable value map
   768   {
   769     typedef SmartGraph Graph;
   770     typedef SmartGraph::Node Item;
   771     typedef IterableValueMap<SmartGraph, SmartGraph::Node, double> Ivm;
   772 
   773     checkConcept<ReadWriteMap<Item, double>, Ivm>();
   774 
   775     const int num = 10;
   776     Graph g;
   777     Ivm map0(g, 0.0);
   778     std::vector<Item> items;
   779     for (int i = 0; i < num; ++i) {
   780       items.push_back(g.addNode());
   781     }
   782 
   783     Ivm map1(g, 0.0);
   784     check(distance(map1.beginValue(), map1.endValue()) == 1, "Wrong size");
   785     check(*map1.beginValue() == 0.0, "Wrong value");
   786 
   787     for (int i = 0; i < num; ++i) {
   788       map1.set(items[i], static_cast<double>(i));
   789     }
   790     check(distance(map1.beginValue(), map1.endValue()) == num, "Wrong size");
   791 
   792     for (int i = 0; i < num; ++i) {
   793       Ivm::ItemIt it(map1, static_cast<double>(i));
   794       check(static_cast<Item>(it) == items[i], "Wrong value");
   795       ++it;
   796       check(static_cast<Item>(it) == INVALID, "Wrong value");
   797     }
   798 
   799     for (Ivm::ValueIt vit = map1.beginValue();
   800          vit != map1.endValue(); ++vit) {
   801       check(map1[static_cast<Item>(Ivm::ItemIt(map1, *vit))] == *vit,
   802             "Wrong ValueIt");
   803     }
   804 
   805     for (int i = 0; i < num; ++i) {
   806       map1.set(items[i], static_cast<double>(i % 2));
   807     }
   808     check(distance(map1.beginValue(), map1.endValue()) == 2, "Wrong size");
   809 
   810     int n = 0;
   811     for (Ivm::ItemIt it(map1, 0.0); it != INVALID; ++it) {
   812       check(map1[static_cast<Item>(it)] == 0.0, "Wrong value");
   813       ++n;
   814     }
   815     check(n == (num + 1) / 2, "Wrong number");
   816 
   817     for (Ivm::ItemIt it(map1, 1.0); it != INVALID; ++it) {
   818       check(map1[static_cast<Item>(it)] == 1.0, "Wrong value");
   819       ++n;
   820     }
   821     check(n == num, "Wrong number");
   822 
   823   }
   824 
   825   // Graph map utilities:
   826   // mapMin(), mapMax(), mapMinValue(), mapMaxValue()
   827   // mapFind(), mapFindIf(), mapCount(), mapCountIf()
   828   // mapCopy(), mapCompare(), mapFill()
   829   {
   830     DIGRAPH_TYPEDEFS(SmartDigraph);
   831 
   832     SmartDigraph g;
   833     Node n1 = g.addNode();
   834     Node n2 = g.addNode();
   835     Node n3 = g.addNode();
   836 
   837     SmartDigraph::NodeMap<int> map1(g);
   838     SmartDigraph::ArcMap<char> map2(g);
   839     ConstMap<Node, A> cmap1 = A();
   840     ConstMap<Arc, C> cmap2 = C(0);
   841 
   842     map1[n1] = 10;
   843     map1[n2] = 5;
   844     map1[n3] = 12;
   845 
   846     // mapMin(), mapMax(), mapMinValue(), mapMaxValue()
   847     check(mapMin(g, map1) == n2, "Wrong mapMin()");
   848     check(mapMax(g, map1) == n3, "Wrong mapMax()");
   849     check(mapMin(g, map1, std::greater<int>()) == n3, "Wrong mapMin()");
   850     check(mapMax(g, map1, std::greater<int>()) == n2, "Wrong mapMax()");
   851     check(mapMinValue(g, map1) == 5, "Wrong mapMinValue()");
   852     check(mapMaxValue(g, map1) == 12, "Wrong mapMaxValue()");
   853 
   854     check(mapMin(g, map2) == INVALID, "Wrong mapMin()");
   855     check(mapMax(g, map2) == INVALID, "Wrong mapMax()");
   856 
   857     check(mapMin(g, cmap1) != INVALID, "Wrong mapMin()");
   858     check(mapMax(g, cmap2) == INVALID, "Wrong mapMax()");
   859 
   860     Arc a1 = g.addArc(n1, n2);
   861     Arc a2 = g.addArc(n1, n3);
   862     Arc a3 = g.addArc(n2, n3);
   863     Arc a4 = g.addArc(n3, n1);
   864 
   865     map2[a1] = 'b';
   866     map2[a2] = 'a';
   867     map2[a3] = 'b';
   868     map2[a4] = 'c';
   869 
   870     // mapMin(), mapMax(), mapMinValue(), mapMaxValue()
   871     check(mapMin(g, map2) == a2, "Wrong mapMin()");
   872     check(mapMax(g, map2) == a4, "Wrong mapMax()");
   873     check(mapMin(g, map2, std::greater<int>()) == a4, "Wrong mapMin()");
   874     check(mapMax(g, map2, std::greater<int>()) == a2, "Wrong mapMax()");
   875     check(mapMinValue(g, map2, std::greater<int>()) == 'c',
   876           "Wrong mapMinValue()");
   877     check(mapMaxValue(g, map2, std::greater<int>()) == 'a',
   878           "Wrong mapMaxValue()");
   879 
   880     check(mapMin(g, cmap1) != INVALID, "Wrong mapMin()");
   881     check(mapMax(g, cmap2) != INVALID, "Wrong mapMax()");
   882     check(mapMaxValue(g, cmap2) == C(0), "Wrong mapMaxValue()");
   883 
   884     check(mapMin(g, composeMap(functorToMap(&createC), map2)) == a2,
   885           "Wrong mapMin()");
   886     check(mapMax(g, composeMap(functorToMap(&createC), map2)) == a4,
   887           "Wrong mapMax()");
   888     check(mapMinValue(g, composeMap(functorToMap(&createC), map2)) == C('a'),
   889           "Wrong mapMinValue()");
   890     check(mapMaxValue(g, composeMap(functorToMap(&createC), map2)) == C('c'),
   891           "Wrong mapMaxValue()");
   892 
   893     // mapFind(), mapFindIf()
   894     check(mapFind(g, map1, 5) == n2, "Wrong mapFind()");
   895     check(mapFind(g, map1, 6) == INVALID, "Wrong mapFind()");
   896     check(mapFind(g, map2, 'a') == a2, "Wrong mapFind()");
   897     check(mapFind(g, map2, 'e') == INVALID, "Wrong mapFind()");
   898     check(mapFind(g, cmap2, C(0)) == ArcIt(g), "Wrong mapFind()");
   899     check(mapFind(g, cmap2, C(1)) == INVALID, "Wrong mapFind()");
   900 
   901     check(mapFindIf(g, map1, Less<int>(7)) == n2,
   902           "Wrong mapFindIf()");
   903     check(mapFindIf(g, map1, Less<int>(5)) == INVALID,
   904           "Wrong mapFindIf()");
   905     check(mapFindIf(g, map2, Less<char>('d')) == ArcIt(g),
   906           "Wrong mapFindIf()");
   907     check(mapFindIf(g, map2, Less<char>('a')) == INVALID,
   908           "Wrong mapFindIf()");
   909 
   910     // mapCount(), mapCountIf()
   911     check(mapCount(g, map1, 5) == 1, "Wrong mapCount()");
   912     check(mapCount(g, map1, 6) == 0, "Wrong mapCount()");
   913     check(mapCount(g, map2, 'a') == 1, "Wrong mapCount()");
   914     check(mapCount(g, map2, 'b') == 2, "Wrong mapCount()");
   915     check(mapCount(g, map2, 'e') == 0, "Wrong mapCount()");
   916     check(mapCount(g, cmap2, C(0)) == 4, "Wrong mapCount()");
   917     check(mapCount(g, cmap2, C(1)) == 0, "Wrong mapCount()");
   918 
   919     check(mapCountIf(g, map1, Less<int>(11)) == 2,
   920           "Wrong mapCountIf()");
   921     check(mapCountIf(g, map1, Less<int>(13)) == 3,
   922           "Wrong mapCountIf()");
   923     check(mapCountIf(g, map1, Less<int>(5)) == 0,
   924           "Wrong mapCountIf()");
   925     check(mapCountIf(g, map2, Less<char>('d')) == 4,
   926           "Wrong mapCountIf()");
   927     check(mapCountIf(g, map2, Less<char>('c')) == 3,
   928           "Wrong mapCountIf()");
   929     check(mapCountIf(g, map2, Less<char>('a')) == 0,
   930           "Wrong mapCountIf()");
   931 
   932     // MapIt, ConstMapIt
   933 /*
   934 These tests can be used after applying bugfix #330
   935     typedef SmartDigraph::NodeMap<int>::MapIt MapIt;
   936     typedef SmartDigraph::NodeMap<int>::ConstMapIt ConstMapIt;
   937     check(*std::min_element(MapIt(map1), MapIt(INVALID)) == 5,
   938           "Wrong NodeMap<>::MapIt");
   939     check(*std::max_element(ConstMapIt(map1), ConstMapIt(INVALID)) == 12,
   940           "Wrong NodeMap<>::MapIt");
   941 
   942     int sum = 0;
   943     std::for_each(MapIt(map1), MapIt(INVALID), Sum<int>(sum));
   944     check(sum == 27, "Wrong NodeMap<>::MapIt");
   945     std::for_each(ConstMapIt(map1), ConstMapIt(INVALID), Sum<int>(sum));
   946     check(sum == 54, "Wrong NodeMap<>::ConstMapIt");
   947 */
   948 
   949     // mapCopy(), mapCompare(), mapFill()
   950     check(mapCompare(g, map1, map1), "Wrong mapCompare()");
   951     check(mapCompare(g, cmap2, cmap2), "Wrong mapCompare()");
   952     check(mapCompare(g, map1, shiftMap(map1, 0)), "Wrong mapCompare()");
   953     check(mapCompare(g, map2, scaleMap(map2, 1)), "Wrong mapCompare()");
   954     check(!mapCompare(g, map1, shiftMap(map1, 1)), "Wrong mapCompare()");
   955 
   956     SmartDigraph::NodeMap<int> map3(g, 0);
   957     SmartDigraph::ArcMap<char> map4(g, 'a');
   958 
   959     check(!mapCompare(g, map1, map3), "Wrong mapCompare()");
   960     check(!mapCompare(g, map2, map4), "Wrong mapCompare()");
   961 
   962     mapCopy(g, map1, map3);
   963     mapCopy(g, map2, map4);
   964 
   965     check(mapCompare(g, map1, map3), "Wrong mapCompare() or mapCopy()");
   966     check(mapCompare(g, map2, map4), "Wrong mapCompare() or mapCopy()");
   967 
   968     Undirector<SmartDigraph> ug(g);
   969     Undirector<SmartDigraph>::EdgeMap<char> umap1(ug, 'x');
   970     Undirector<SmartDigraph>::ArcMap<double> umap2(ug, 3.14);
   971 
   972     check(!mapCompare(g, map2, umap1), "Wrong mapCompare() or mapCopy()");
   973     check(!mapCompare(g, umap1, map2), "Wrong mapCompare() or mapCopy()");
   974     check(!mapCompare(ug, map2, umap1), "Wrong mapCompare() or mapCopy()");
   975     check(!mapCompare(ug, umap1, map2), "Wrong mapCompare() or mapCopy()");
   976 
   977     mapCopy(g, map2, umap1);
   978 
   979     check(mapCompare(g, map2, umap1), "Wrong mapCompare() or mapCopy()");
   980     check(mapCompare(g, umap1, map2), "Wrong mapCompare() or mapCopy()");
   981     check(mapCompare(ug, map2, umap1), "Wrong mapCompare() or mapCopy()");
   982     check(mapCompare(ug, umap1, map2), "Wrong mapCompare() or mapCopy()");
   983 
   984     mapCopy(g, map2, umap1);
   985     mapCopy(g, umap1, map2);
   986     mapCopy(ug, map2, umap1);
   987     mapCopy(ug, umap1, map2);
   988 
   989     check(!mapCompare(ug, umap1, umap2), "Wrong mapCompare() or mapCopy()");
   990     mapCopy(ug, umap1, umap2);
   991     check(mapCompare(ug, umap1, umap2), "Wrong mapCompare() or mapCopy()");
   992 
   993     check(!mapCompare(g, map1, constMap<Node>(2)), "Wrong mapCompare()");
   994     mapFill(g, map1, 2);
   995     check(mapCompare(g, constMap<Node>(2), map1), "Wrong mapFill()");
   996 
   997     check(!mapCompare(g, map2, constMap<Arc>('z')), "Wrong mapCompare()");
   998     mapCopy(g, constMap<Arc>('z'), map2);
   999     check(mapCompare(g, constMap<Arc>('z'), map2), "Wrong mapCopy()");
  1000   }
  1001 
  1002   return 0;
  1003 }