RhodeCode

/* -*- C++ -*-

 *

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

 *

 * Copyright (C) 2003-2008

 * 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.

 *

 */

#ifndef LEMON_MAPS_H

#define LEMON_MAPS_H

#include <iterator>

#include <functional>

#include <vector>

#include <lemon/bits/utility.h>

#include <lemon/bits/traits.h>

///\file

///\ingroup maps

///\brief Miscellaneous property maps

#include <map>

namespace lemon {

  /// \addtogroup maps

  /// @{

  /// Base class of maps.

  /// Base class of maps. It provides the necessary type definitions

  /// required by the map %concepts.

  template<typename K, typename V>

  class MapBase {

  public:

    /// \biref The key type of the map.

    typedef K Key;

    /// \brief The value type of the map.

    /// (The type of objects associated with the keys).

    typedef V Value;

  };

  /// Null map. (a.k.a. DoNothingMap)

  /// This map can be used if you have to provide a map only for

  /// its type definitions, or if you have to provide a writable map,

  /// but data written to it is not required (i.e. it will be sent to

  /// <tt>/dev/null</tt>).

  /// It conforms the \ref concepts::ReadWriteMap "ReadWriteMap" concept.

  ///

  /// \sa ConstMap

  template<typename K, typename V>

  class NullMap : public MapBase<K, V> {

  public:

    typedef MapBase<K, V> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Gives back a default constructed element.

    Value operator[](const Key&) const { return Value(); }

    /// Absorbs the value.

    void set(const Key&, const Value&) {}

  };

  /// Returns a \ref NullMap class

  /// This function just returns a \ref NullMap class.

  /// \relates NullMap

  template <typename K, typename V>

  NullMap<K, V> nullMap() {

    return NullMap<K, V>();

  }

  /// Constant map.

  /// This \ref concepts::ReadMap "readable map" assigns a specified

  /// value to each key.

  ///

  /// In other aspects it is equivalent to \ref NullMap.

  /// So it conforms the \ref concepts::ReadWriteMap "ReadWriteMap"

  /// concept, but it absorbs the data written to it.

  ///

  /// The simplest way of using this map is through the constMap()

  /// function.

  ///

  /// \sa NullMap

  /// \sa IdentityMap

  template<typename K, typename V>

  class ConstMap : public MapBase<K, V> {

  private:

    V _value;

  public:

    typedef MapBase<K, V> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Default constructor

    /// Default constructor.

    /// The value of the map will be default constructed.

    ConstMap() {}

    /// Constructor with specified initial value

    /// Constructor with specified initial value.

    ConstMap(const Value &v) : _value(v) {}

    /// Gives back the specified value.

    Value operator[](const Key&) const { return _value; }

    /// Absorbs the value.

    void set(const Key&, const Value&) {}

    /// Sets the value that is assigned to each key.

    void setAll(const Value &v) {

      _value = v;

    }

    template<typename V1>

    ConstMap(const ConstMap<K, V1> &, const Value &v) : _value(v) {}

  };

  /// Returns a \ref ConstMap class

  /// This function just returns a \ref ConstMap class.

  /// \relates ConstMap

  template<typename K, typename V>

  inline ConstMap<K, V> constMap(const V &v) {

    return ConstMap<K, V>(v);

  }

  template<typename T, T v>

  struct Const {};

  /// Constant map with inlined constant value.

  /// This \ref concepts::ReadMap "readable map" assigns a specified

  /// value to each key.

  ///

  /// In other aspects it is equivalent to \ref NullMap.

  /// So it conforms the \ref concepts::ReadWriteMap "ReadWriteMap"

  /// concept, but it absorbs the data written to it.

  ///

  /// The simplest way of using this map is through the constMap()

  /// function.

  ///

  /// \sa NullMap

  /// \sa IdentityMap

  template<typename K, typename V, V v>

  class ConstMap<K, Const<V, v> > : public MapBase<K, V> {

  public:

    typedef MapBase<K, V> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Constructor.

    ConstMap() {}

    /// Gives back the specified value.

    Value operator[](const Key&) const { return v; }

    /// Absorbs the value.

    void set(const Key&, const Value&) {}

  };

  /// Returns a \ref ConstMap class with inlined constant value

  /// This function just returns a \ref ConstMap class with inlined

  /// constant value.

  /// \relates ConstMap

  template<typename K, typename V, V v>

  inline ConstMap<K, Const<V, v> > constMap() {

    return ConstMap<K, Const<V, v> >();

  }

  /// Identity map.

  /// This \ref concepts::ReadMap "read-only map" gives back the given

  /// key as value without any modification.

  ///

  /// \sa ConstMap

  template <typename T>

  class IdentityMap : public MapBase<T, T> {

  public:

    typedef MapBase<T, T> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Gives back the given value without any modification.

    Value operator[](const Key &k) const {

      return k;

    }

  };

  /// Returns an \ref IdentityMap class

  /// This function just returns an \ref IdentityMap class.

  /// \relates IdentityMap

  template<typename T>

  inline IdentityMap<T> identityMap() {

    return IdentityMap<T>();

  }

  /// \brief Map for storing values for integer keys from the range

  /// <tt>[0..size-1]</tt>.

  ///

  /// This map is essentially a wrapper for \c std::vector. It assigns

  /// values to integer keys from the range <tt>[0..size-1]</tt>.

  /// It can be used with some data structures, for example

  /// \ref UnionFind, \ref BinHeap, when the used items are small

  /// integers. This map conforms the \ref concepts::ReferenceMap

  /// "ReferenceMap" concept.

  ///

  /// The simplest way of using this map is through the rangeMap()

  /// function.

  template <typename V>

  class RangeMap : public MapBase<int, V> {

    template <typename V1>

    friend class RangeMap;

  private:

    typedef std::vector<V> Vector;

    Vector _vector;

  public:

    typedef MapBase<int, V> Parent;

    /// Key type

    typedef typename Parent::Key Key;

    /// Value type

    typedef typename Parent::Value Value;

    /// Reference type

    typedef typename Vector::reference Reference;

    /// Const reference type

    typedef typename Vector::const_reference ConstReference;

    typedef True ReferenceMapTag;

  public:

    /// Constructor with specified default value.

    RangeMap(int size = 0, const Value &value = Value())

      : _vector(size, value) {}

    /// Constructs the map from an appropriate \c std::vector.

    template <typename V1>

    RangeMap(const std::vector<V1>& vector)

      : _vector(vector.begin(), vector.end()) {}

    /// Constructs the map from another \ref RangeMap.

    template <typename V1>

    RangeMap(const RangeMap<V1> &c)

      : _vector(c._vector.begin(), c._vector.end()) {}

    /// Returns the size of the map.

    int size() {

      return _vector.size();

    }

    /// Resizes the map.

    /// Resizes the underlying \c std::vector container, so changes the

    /// keyset of the map.

    /// \param size The new size of the map. The new keyset will be the

    /// range <tt>[0..size-1]</tt>.

    /// \param value The default value to assign to the new keys.

    void resize(int size, const Value &value = Value()) {

      _vector.resize(size, value);

    }

  private:

    RangeMap& operator=(const RangeMap&);

  public:

    ///\e

    Reference operator[](const Key &k) {

      return _vector[k];

    }

    ///\e

    ConstReference operator[](const Key &k) const {

      return _vector[k];

    }

    ///\e

    void set(const Key &k, const Value &v) {

      _vector[k] = v;

    }

  };

  /// Returns a \ref RangeMap class

  /// This function just returns a \ref RangeMap class.

  /// \relates RangeMap

  template<typename V>

  inline RangeMap<V> rangeMap(int size = 0, const V &value = V()) {

    return RangeMap<V>(size, value);

  }

  /// \brief Returns a \ref RangeMap class created from an appropriate

  /// \c std::vector

  /// This function just returns a \ref RangeMap class created from an

  /// appropriate \c std::vector.

  /// \relates RangeMap

  template<typename V>

  inline RangeMap<V> rangeMap(const std::vector<V> &vector) {

    return RangeMap<V>(vector);

  }

  /// Map type based on \c std::map

  /// This map is essentially a wrapper for \c std::map with addition

  /// that you can specify a default value for the keys that are not

  /// stored actually. This value can be different from the default

  /// contructed value (i.e. \c %Value()).

  /// This type conforms the \ref concepts::ReferenceMap "ReferenceMap"

    ///\e

    void setAll(const Value &v) {

      _value = v;

      _map.clear();

    }

  };

  /// Returns a \ref SparseMap class

  /// This function just returns a \ref SparseMap class with specified

  /// default value.

  /// \relates SparseMap

  template<typename K, typename V, typename Compare>

  inline SparseMap<K, V, Compare> sparseMap(const V& value = V()) {

    return SparseMap<K, V, Compare>(value);

  }

  template<typename K, typename V>

  inline SparseMap<K, V, std::less<K> > sparseMap(const V& value = V()) {

    return SparseMap<K, V, std::less<K> >(value);

  }

  /// \brief Returns a \ref SparseMap class created from an appropriate

  /// \c std::map

  /// This function just returns a \ref SparseMap class created from an

  /// appropriate \c std::map.

  /// \relates SparseMap

  template<typename K, typename V, typename Compare>

  inline SparseMap<K, V, Compare>

    sparseMap(const std::map<K, V, Compare> &map, const V& value = V())

  {

    return SparseMap<K, V, Compare>(map, value);

  }

  /// @}

  /// \addtogroup map_adaptors

  /// @{

  /// Composition of two maps

  /// This \ref concepts::ReadMap "read-only map" returns the

  /// composition of two given maps. That is to say, if \c m1 is of

  /// type \c M1 and \c m2 is of \c M2, then for

  /// \code

  ///   ComposeMap<M1, M2> cm(m1,m2);

  /// \endcode

  /// <tt>cm[x]</tt> will be equal to <tt>m1[m2[x]]</tt>.

  ///

  /// The \c Key type of the map is inherited from \c M2 and the

  /// \c Value type is from \c M1.

  /// \c M2::Value must be convertible to \c M1::Key.

  ///

  /// The simplest way of using this map is through the composeMap()

  /// function.

  ///

  /// \sa CombineMap

  ///

  /// \todo Check the requirements.

  template <typename M1, typename M2>

  class ComposeMap : public MapBase<typename M2::Key, typename M1::Value> {

    const M1 &_m1;

    const M2 &_m2;

  public:

    typedef MapBase<typename M2::Key, typename M1::Value> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Constructor

    ComposeMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}

    /// \e

    typename MapTraits<M1>::ConstReturnValue

    operator[](const Key &k) const { return _m1[_m2[k]]; }

  };

  /// Returns a \ref ComposeMap class

  /// This function just returns a \ref ComposeMap class.

  ///

  /// If \c m1 and \c m2 are maps and the \c Value type of \c m2 is

  /// convertible to the \c Key of \c m1, then <tt>composeMap(m1,m2)[x]</tt>

  /// will be equal to <tt>m1[m2[x]]</tt>.

  ///

  /// \relates ComposeMap

  template <typename M1, typename M2>

  inline ComposeMap<M1, M2> composeMap(const M1 &m1, const M2 &m2) {

    return ComposeMap<M1, M2>(m1, m2);

  }

  /// Combination of two maps using an STL (binary) functor.

  /// This \ref concepts::ReadMap "read-only map" takes two maps and a

  /// binary functor and returns the combination of the two given maps

  /// using the functor.

  /// That is to say, if \c m1 is of type \c M1 and \c m2 is of \c M2

  /// and \c f is of \c F, then for

  /// \code

  ///   CombineMap<M1,M2,F,V> cm(m1,m2,f);

  /// \endcode

  /// <tt>cm[x]</tt> will be equal to <tt>f(m1[x],m2[x])</tt>.

  ///

  /// The \c Key type of the map is inherited from \c M1 (\c M1::Key

  /// must be convertible to \c M2::Key) and the \c Value type is \c V.

  /// \c M2::Value and \c M1::Value must be convertible to the

  /// corresponding input parameter of \c F and the return type of \c F

  /// must be convertible to \c V.

  ///

  /// The simplest way of using this map is through the combineMap()

  /// function.

  ///

  /// \sa ComposeMap

  ///

  /// \todo Check the requirements.

  template<typename M1, typename M2, typename F,

	   typename V = typename F::result_type>

  class CombineMap : public MapBase<typename M1::Key, V> {

    const M1 &_m1;

    const M2 &_m2;

    F _f;

  public:

    typedef MapBase<typename M1::Key, V> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Constructor

    CombineMap(const M1 &m1, const M2 &m2, const F &f = F())

      : _m1(m1), _m2(m2), _f(f) {}

    /// \e

    Value operator[](const Key &k) const { return _f(_m1[k],_m2[k]); }

  };

  /// Returns a \ref CombineMap class

  /// This function just returns a \ref CombineMap class.

  ///

  /// For example, if \c m1 and \c m2 are both maps with \c double

  /// values, then

  /// \code

  ///   combineMap(m1,m2,std::plus<double>())

  /// \endcode

  /// is equivalent to

  /// \code

  ///   addMap(m1,m2)

  /// \endcode

  ///

  /// This function is specialized for adaptable binary function

  /// classes and C++ functions.

  ///

  /// \relates CombineMap

  template<typename M1, typename M2, typename F, typename V>

  inline CombineMap<M1, M2, F, V>

  combineMap(const M1 &m1, const M2 &m2, const F &f) {

    return CombineMap<M1, M2, F, V>(m1,m2,f);

  }

  template<typename M1, typename M2, typename F>

  inline CombineMap<M1, M2, F, typename F::result_type>

  combineMap(const M1 &m1, const M2 &m2, const F &f) {

    return combineMap<M1, M2, F, typename F::result_type>(m1,m2,f);

  }

  template<typename M1, typename M2, typename K1, typename K2, typename V>

  inline CombineMap<M1, M2, V (*)(K1, K2), V>

  combineMap(const M1 &m1, const M2 &m2, V (*f)(K1, K2)) {

    return combineMap<M1, M2, V (*)(K1, K2), V>(m1,m2,f);

  }

  /// Converts an STL style (unary) functor to a map

  /// This \ref concepts::ReadMap "read-only map" returns the value

  /// of a given functor. Actually, it just wraps the functor and

  /// provides the \c Key and \c Value typedefs.

  ///

  /// Template parameters \c K and \c V will become its \c Key and

  /// \c Value. In most cases they have to be given explicitly because

  /// a functor typically does not provide \c argument_type and

  /// \c result_type typedefs.

  /// Parameter \c F is the type of the used functor.

  ///

  /// The simplest way of using this map is through the functorToMap()

  /// function.

  ///

  /// \sa MapToFunctor

  template<typename F,

	   typename K = typename F::argument_type,

	   typename V = typename F::result_type>

  class FunctorToMap : public MapBase<K, V> {

  public:

    typedef MapBase<K, V> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Constructor

    FunctorToMap(const F &f = F()) : _f(f) {}

    /// \e

    Value operator[](const Key &k) const { return _f(k); }

  };

  /// Returns a \ref FunctorToMap class

  /// This function just returns a \ref FunctorToMap class.

  ///

  /// This function is specialized for adaptable binary function

  /// classes and C++ functions.

  ///

  /// \relates FunctorToMap

  template<typename K, typename V, typename F>

  inline FunctorToMap<F, K, V> functorToMap(const F &f) {

    return FunctorToMap<F, K, V>(f);

  }

  template <typename F>

  inline FunctorToMap<F, typename F::argument_type, typename F::result_type>

    functorToMap(const F &f)

  {

    return FunctorToMap<F, typename F::argument_type,

      typename F::result_type>(f);

  }

  template <typename K, typename V>

  inline FunctorToMap<V (*)(K), K, V> functorToMap(V (*f)(K)) {

    return FunctorToMap<V (*)(K), K, V>(f);

  }

  /// Converts a map to an STL style (unary) functor

  /// This class converts a map to an STL style (unary) functor.

  /// That is it provides an <tt>operator()</tt> to read its values.

  ///

  /// For the sake of convenience it also works as a usual

  /// \ref concepts::ReadMap "readable map", i.e. <tt>operator[]</tt>

  /// and the \c Key and \c Value typedefs also exist.

  ///

  /// The simplest way of using this map is through the mapToFunctor()

  /// function.

  ///

  ///\sa FunctorToMap

  template <typename M>

  class MapToFunctor : public MapBase<typename M::Key, typename M::Value> {

    const M &_m;

  public:

    typedef MapBase<typename M::Key, typename M::Value> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    typedef typename Parent::Key argument_type;

    typedef typename Parent::Value result_type;

    /// Constructor

    MapToFunctor(const M &m) : _m(m) {}

    /// \e

    Value operator()(const Key &k) const { return _m[k]; }

    /// \e

    Value operator[](const Key &k) const { return _m[k]; }

  };

  /// Returns a \ref MapToFunctor class

  /// This function just returns a \ref MapToFunctor class.

  /// \relates MapToFunctor

  template<typename M>

  inline MapToFunctor<M> mapToFunctor(const M &m) {

    return MapToFunctor<M>(m);

  }

  /// \brief Map adaptor to convert the \c Value type of a map to

  /// another type using the default conversion.

  /// Map adaptor to convert the \c Value type of a \ref concepts::ReadMap

  /// "readable map" to another type using the default conversion.

  /// The \c Key type of it is inherited from \c M and the \c Value

  /// type is \c V.

  /// This type conforms the \ref concepts::ReadMap "ReadMap" concept.

  ///

  /// The simplest way of using this map is through the convertMap()

  /// function.

  template <typename M, typename V>

  class ConvertMap : public MapBase<typename M::Key, V> {

    const M &_m;

  public:

    typedef MapBase<typename M::Key, V> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Constructor

    /// Constructor.

    /// \param m The underlying map.

    ConvertMap(const M &m) : _m(m) {}

    /// \e

    Value operator[](const Key &k) const { return _m[k]; }

  };

  /// Returns a \ref ConvertMap class

  /// This function just returns a \ref ConvertMap class.

  /// \relates ConvertMap

  template<typename V, typename M>

  inline ConvertMap<M, V> convertMap(const M &map) {

    return ConvertMap<M, V>(map);

  }

  /// Applies all map setting operations to two maps

  /// This map has two \ref concepts::WriteMap "writable map" parameters

  /// and each write request will be passed to both of them.

  /// If \c M1 is also \ref concepts::ReadMap "readable", then the read

  /// operations will return the corresponding values of \c M1.

  ///

  /// The \c Key and \c Value types are inherited from \c M1.

  /// The \c Key and \c Value of \c M2 must be convertible from those

  /// of \c M1.

  ///

  /// The simplest way of using this map is through the forkMap()

  /// function.

  template<typename  M1, typename M2>

  class ForkMap : public MapBase<typename M1::Key, typename M1::Value> {

    M1 &_m1;

    M2 &_m2;

  public:

    typedef MapBase<typename M1::Key, typename M1::Value> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Constructor

    ForkMap(M1 &m1, M2 &m2) : _m1(m1), _m2(m2) {}

    /// Returns the value associated with the given key in the first map.

    Value operator[](const Key &k) const { return _m1[k]; }

    /// Sets the value associated with the given key in both maps.

    void set(const Key &k, const Value &v) { _m1.set(k,v); _m2.set(k,v); }

  };

  /// Returns a \ref ForkMap class

  /// This function just returns a \ref ForkMap class.

  /// \relates ForkMap

  template <typename M1, typename M2>

  inline ForkMap<M1,M2> forkMap(M1 &m1, M2 &m2) {

    return ForkMap<M1,M2>(m1,m2);

  }

  /// Sum of two maps

  /// This \ref concepts::ReadMap "read-only map" returns the sum

  /// of the values of the two given maps.

  /// Its \c Key and \c Value types are inherited from \c M1.

  /// The \c Key and \c Value of \c M2 must be convertible to those of

  /// \c M1.

  ///

  /// If \c m1 is of type \c M1 and \c m2 is of \c M2, then for

  /// \code

  ///   AddMap<M1,M2> am(m1,m2);

  /// \endcode

  /// <tt>am[x]</tt> will be equal to <tt>m1[x]+m2[x]</tt>.

  ///

  /// The simplest way of using this map is through the addMap()

  /// function.

  ///

  /// \sa SubMap, MulMap, DivMap

  /// \sa ShiftMap, ShiftWriteMap

  template<typename M1, typename M2>

  class AddMap : public MapBase<typename M1::Key, typename M1::Value> {

    const M1 &_m1;

    const M2 &_m2;

  public:

    typedef MapBase<typename M1::Key, typename M1::Value> Parent;

    typedef typename Parent::Key Key;

    typedef typename Parent::Value Value;

    /// Constructor

    AddMap(const M1 &m1, const M2 &m2) : _m1(m1), _m2(m2) {}

    /// \e

    Value operator[](const Key &k) const { return _m1[k]+_m2[k]; }

  };

/* -*- C++ -*-

 *

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

 *

 * Copyright (C) 2003-2008

 * 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 <deque>

#include <set>

#include <lemon/concept_check.h>

#include <lemon/concepts/maps.h>

#include <lemon/maps.h>

#include "test_tools.h"

using namespace lemon;

using namespace lemon::concepts;

struct A {};

inline bool operator<(A, A) { return true; }

struct B {};

class C {

  int x;

public:

  C(int _x) : x(_x) {}

};

class F {

public:

  typedef A argument_type;

  typedef B result_type;

  B operator()(const A&) const { return B(); }

private:

  F& operator=(const F&);

};

int func(A) { return 3; }

int binc(int a, B) { return a+1; }

typedef ReadMap<A, double> DoubleMap;

typedef ReadWriteMap<A, double> DoubleWriteMap;

typedef ReferenceMap<A, double, double&, const double&> DoubleRefMap;

typedef ReadMap<A, bool> BoolMap;

typedef ReadWriteMap<A, bool> BoolWriteMap;

typedef ReferenceMap<A, bool, bool&, const bool&> BoolRefMap;

int main()

{

  // Map concepts

  checkConcept<ReadMap<A,B>, ReadMap<A,B> >();

  checkConcept<ReadMap<A,C>, ReadMap<A,C> >();

  checkConcept<WriteMap<A,B>, WriteMap<A,B> >();

  checkConcept<WriteMap<A,C>, WriteMap<A,C> >();

  checkConcept<ReadWriteMap<A,B>, ReadWriteMap<A,B> >();

  checkConcept<ReadWriteMap<A,C>, ReadWriteMap<A,C> >();

  checkConcept<ReferenceMap<A,B,B&,const B&>, ReferenceMap<A,B,B&,const B&> >();

  checkConcept<ReferenceMap<A,C,C&,const C&>, ReferenceMap<A,C,C&,const C&> >();

  // NullMap

  {

    checkConcept<ReadWriteMap<A,B>, NullMap<A,B> >();

    NullMap<A,B> map1;

    NullMap<A,B> map2 = map1;

    map1 = nullMap<A,B>();

  }

  // ConstMap

  {

    checkConcept<ReadWriteMap<A,B>, ConstMap<A,B> >();

    ConstMap<A,B> map1;

    ConstMap<A,B> map2(B());

    ConstMap<A,B> map3 = map1;

    map1 = constMap<A>(B());

    map1.setAll(B());

    checkConcept<ReadWriteMap<A,int>, ConstMap<A,int> >();

    check(constMap<A>(10)[A()] == 10, "Something is wrong with ConstMap");

    checkConcept<ReadWriteMap<A,int>, ConstMap<A,Const<int,10> > >();

  }

  // IdentityMap

  {

    checkConcept<ReadMap<A,A>, IdentityMap<A> >();

    IdentityMap<A> map1;

    IdentityMap<A> map2 = map1;

    map1 = identityMap<A>();

    checkConcept<ReadMap<double,double>, IdentityMap<double> >();

    check(identityMap<double>()[1.0] == 1.0 && identityMap<double>()[3.14] == 3.14,

          "Something is wrong with IdentityMap");

  }

  // RangeMap

  {

    checkConcept<ReferenceMap<int,B,B&,const B&>, RangeMap<B> >();

    RangeMap<B> map1;

    RangeMap<B> map2(10);

    RangeMap<B> map3(10,B());

    RangeMap<B> map4 = map1;

    RangeMap<B> map5 = rangeMap<B>();

    RangeMap<B> map6 = rangeMap<B>(10);

    RangeMap<B> map7 = rangeMap(10,B());

    checkConcept< ReferenceMap<int, double, double&, const double&>,

                  RangeMap<double> >();

    std::vector<double> v(10, 0);

    v[5] = 100;

    RangeMap<double> map8(v);

    RangeMap<double> map9 = rangeMap(v);

    check(map9.size() == 10 && map9[2] == 0 && map9[5] == 100,

          "Something is wrong with RangeMap");

  }

  // SparseMap

  {

    checkConcept<ReferenceMap<A,B,B&,const B&>, SparseMap<A,B> >();

    SparseMap<A,B> map1;

    SparseMap<A,B> map2(B());

    SparseMap<A,B> map3 = sparseMap<A,B>();

    SparseMap<A,B> map4 = sparseMap<A>(B());

    checkConcept< ReferenceMap<double, int, int&, const int&>,

                  SparseMap<double, int> >();

    std::map<double, int> m;

    SparseMap<double, int> map5(m);

    SparseMap<double, int> map6(m,10);

    SparseMap<double, int> map7 = sparseMap(m);

    SparseMap<double, int> map8 = sparseMap(m,10);

    check(map5[1.0] == 0 && map5[3.14] == 0 && map6[1.0] == 10 && map6[3.14] == 10,

          "Something is wrong with SparseMap");

    map5[1.0] = map6[3.14] = 100;

    check(map5[1.0] == 100 && map5[3.14] == 0 && map6[1.0] == 10 && map6[3.14] == 100,

          "Something is wrong with SparseMap");

  }

  // ComposeMap

  {

    typedef ComposeMap<DoubleMap, ReadMap<B,A> > CompMap;

    checkConcept<ReadMap<B,double>, CompMap>();

    CompMap map1(DoubleMap(),ReadMap<B,A>());

    CompMap map2 = composeMap(DoubleMap(), ReadMap<B,A>());

    SparseMap<double, bool> m1(false); m1[3.14] = true;

    RangeMap<double> m2(2); m2[0] = 3.0; m2[1] = 3.14;

    check(!composeMap(m1,m2)[0] && composeMap(m1,m2)[1], "Something is wrong with ComposeMap")

  }

  // CombineMap

  {

    typedef CombineMap<DoubleMap, DoubleMap, std::plus<double> > CombMap;

    checkConcept<ReadMap<A,double>, CombMap>();

    CombMap map1(DoubleMap(), DoubleMap());

    CombMap map2 = combineMap(DoubleMap(), DoubleMap(), std::plus<double>());

    check(combineMap(constMap<B,int,2>(), identityMap<B>(), &binc)[B()] == 3,

          "Something is wrong with CombineMap");

  }

  // FunctorToMap, MapToFunctor

  {

    checkConcept<ReadMap<A,B>, FunctorToMap<F,A,B> >();

    checkConcept<ReadMap<A,B>, FunctorToMap<F> >();

    FunctorToMap<F> map1;

    FunctorToMap<F> map2(F());

    B b = functorToMap(F())[A()];

    checkConcept<ReadMap<A,B>, MapToFunctor<ReadMap<A,B> > >();

    MapToFunctor<ReadMap<A,B> > map(ReadMap<A,B>());

    check(functorToMap(&func)[A()] == 3, "Something is wrong with FunctorToMap");

    check(mapToFunctor(constMap<A,int>(2))(A()) == 2, "Something is wrong with MapToFunctor");

    check(mapToFunctor(functorToMap(&func))(A()) == 3 && mapToFunctor(functorToMap(&func))[A()] == 3,

          "Something is wrong with FunctorToMap or MapToFunctor");

    check(functorToMap(mapToFunctor(constMap<A,int>(2)))[A()] == 2,

          "Something is wrong with FunctorToMap or MapToFunctor");

  }

  // ConvertMap

  {

    checkConcept<ReadMap<double,double>, ConvertMap<ReadMap<double, int>, double> >();

    ConvertMap<RangeMap<bool>, int> map1(rangeMap(1, true));

    ConvertMap<RangeMap<bool>, int> map2 = convertMap<int>(rangeMap(2, false));

  }

  // ForkMap

  {

    checkConcept<DoubleWriteMap, ForkMap<DoubleWriteMap, DoubleWriteMap> >();

    typedef RangeMap<double> RM;

    typedef SparseMap<int, double> SM;

    RM m1(10, -1);

    SM m2(-1);

    checkConcept<ReadWriteMap<int, double>, ForkMap<RM, SM> >();

    checkConcept<ReadWriteMap<int, double>, ForkMap<SM, RM> >();

    ForkMap<RM, SM> map1(m1,m2);

    ForkMap<SM, RM> map2 = forkMap(m2,m1);

    map2.set(5, 10);

    check(m1[1] == -1 && m1[5] == 10 && m2[1] == -1 && m2[5] == 10 && map2[1] == -1 && map2[5] == 10,

          "Something is wrong with ForkMap");

  }

  // Arithmetic maps:

  // - AddMap, SubMap, MulMap, DivMap

  // - ShiftMap, ShiftWriteMap, ScaleMap, ScaleWriteMap

  // - NegMap, NegWriteMap, AbsMap

  {

    checkConcept<DoubleMap, AddMap<DoubleMap,DoubleMap> >();

    checkConcept<DoubleMap, SubMap<DoubleMap,DoubleMap> >();

    checkConcept<DoubleMap, MulMap<DoubleMap,DoubleMap> >();

    checkConcept<DoubleMap, DivMap<DoubleMap,DoubleMap> >();

    ConstMap<int, double> c1(1.0), c2(3.14);

    IdentityMap<int> im;

    ConvertMap<IdentityMap<int>, double> id(im);

    check(addMap(c1,id)[0] == 1.0  && addMap(c1,id)[10] == 11.0, "Something is wrong with AddMap");

    check(subMap(id,c1)[0] == -1.0 && subMap(id,c1)[10] == 9.0,  "Something is wrong with SubMap");

    check(mulMap(id,c2)[0] == 0    && mulMap(id,c2)[2]  == 6.28, "Something is wrong with MulMap");

    check(divMap(c2,id)[1] == 3.14 && divMap(c2,id)[2]  == 1.57, "Something is wrong with DivMap");

    checkConcept<DoubleMap, ShiftMap<DoubleMap> >();

    checkConcept<DoubleWriteMap, ShiftWriteMap<DoubleWriteMap> >();

    checkConcept<DoubleMap, ScaleMap<DoubleMap> >();

    checkConcept<DoubleWriteMap, ScaleWriteMap<DoubleWriteMap> >();

    checkConcept<DoubleMap, NegMap<DoubleMap> >();

    checkConcept<DoubleWriteMap, NegWriteMap<DoubleWriteMap> >();

    checkConcept<DoubleMap, AbsMap<DoubleMap> >();

    check(shiftMap(id, 2.0)[1] == 3.0 && shiftMap(id, 2.0)[10] == 12.0,

          "Something is wrong with ShiftMap");

    check(shiftWriteMap(id, 2.0)[1] == 3.0 && shiftWriteMap(id, 2.0)[10] == 12.0,

          "Something is wrong with ShiftWriteMap");

    check(scaleMap(id, 2.0)[1] == 2.0 && scaleMap(id, 2.0)[10] == 20.0,

          "Something is wrong with ScaleMap");

    check(scaleWriteMap(id, 2.0)[1] == 2.0 && scaleWriteMap(id, 2.0)[10] == 20.0,

          "Something is wrong with ScaleWriteMap");

    check(negMap(id)[1] == -1.0 && negMap(id)[-10] == 10.0,

          "Something is wrong with NegMap");

    check(negWriteMap(id)[1] == -1.0 && negWriteMap(id)[-10] == 10.0,

          "Something is wrong with NegWriteMap");

    check(absMap(id)[1] == 1.0 && absMap(id)[-10] == 10.0,

          "Something is wrong with AbsMap");

  }

  // Logical maps:

  // - TrueMap, FalseMap

  // - AndMap, OrMap

  // - NotMap, NotWriteMap

  // - EqualMap, LessMap

  {

    checkConcept<BoolMap, TrueMap<A> >();

    checkConcept<BoolMap, FalseMap<A> >();

    checkConcept<BoolMap, AndMap<BoolMap,BoolMap> >();

    checkConcept<BoolMap, OrMap<BoolMap,BoolMap> >();

    checkConcept<BoolMap, NotMap<BoolMap> >();

    checkConcept<BoolWriteMap, NotWriteMap<BoolWriteMap> >();

    checkConcept<BoolMap, EqualMap<DoubleMap,DoubleMap> >();

    checkConcept<BoolMap, LessMap<DoubleMap,DoubleMap> >();

    TrueMap<int> tm;

    FalseMap<int> fm;

    RangeMap<bool> rm(2);

    rm[0] = true; rm[1] = false;

    check(andMap(tm,rm)[0] && !andMap(tm,rm)[1] && !andMap(fm,rm)[0] && !andMap(fm,rm)[1],

          "Something is wrong with AndMap");

    check(orMap(tm,rm)[0] && orMap(tm,rm)[1] && orMap(fm,rm)[0] && !orMap(fm,rm)[1],

          "Something is wrong with OrMap");

    check(!notMap(rm)[0] && notMap(rm)[1], "Something is wrong with NotMap");

    check(!notWriteMap(rm)[0] && notWriteMap(rm)[1], "Something is wrong with NotWriteMap");