0
17
0
1
1
1
1
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
///\file |
| 20 | 20 |
///\brief Some basic non inline function and static global data. |
| 21 | 21 |
|
| 22 | 22 |
#include<lemon/tolerance.h> |
| 23 | 23 |
#include<lemon/bits/invalid.h> |
| 24 | 24 |
namespace lemon {
|
| 25 | 25 |
|
| 26 | 26 |
float Tolerance<float>::def_epsilon = 1e-4; |
| 27 | 27 |
double Tolerance<double>::def_epsilon = 1e-10; |
| 28 | 28 |
long double Tolerance<long double>::def_epsilon = 1e-14; |
| 29 | 29 |
|
| 30 | 30 |
#ifndef LEMON_ONLY_TEMPLATES |
| 31 | 31 |
const Invalid INVALID = Invalid(); |
| 32 | 32 |
#endif |
| 33 | 33 |
|
| 34 | 34 |
} //namespace lemon |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#ifndef LEMON_BITS_INVALID_H |
| 20 | 20 |
#define LEMON_BITS_INVALID_H |
| 21 | 21 |
|
| 22 | 22 |
///\file |
| 23 | 23 |
///\brief Definition of INVALID. |
| 24 | 24 |
|
| 25 | 25 |
namespace lemon {
|
| 26 | 26 |
|
| 27 | 27 |
/// \brief Dummy type to make it easier to create invalid iterators. |
| 28 | 28 |
/// |
| 29 | 29 |
/// See \ref INVALID for the usage. |
| 30 | 30 |
struct Invalid {
|
| 31 | 31 |
public: |
| 32 | 32 |
bool operator==(Invalid) { return true; }
|
| 33 | 33 |
bool operator!=(Invalid) { return false; }
|
| 34 | 34 |
bool operator< (Invalid) { return false; }
|
| 35 | 35 |
}; |
| 36 | 36 |
|
| 37 | 37 |
/// \brief Invalid iterators. |
| 38 | 38 |
/// |
| 39 | 39 |
/// \ref Invalid is a global type that converts to each iterator |
| 40 | 40 |
/// in such a way that the value of the target iterator will be invalid. |
| 41 | 41 |
|
| 42 | 42 |
//Some people didn't like this: |
| 43 | 43 |
//const Invalid &INVALID = *(Invalid *)0; |
| 44 | 44 |
|
| 45 | 45 |
#ifdef LEMON_ONLY_TEMPLATES |
| 46 | 46 |
const Invalid INVALID = Invalid(); |
| 47 | 47 |
#else |
| 48 | 48 |
extern const Invalid INVALID; |
| 49 | 49 |
#endif |
| 50 | 50 |
|
| 51 | 51 |
} //namespace lemon |
| 52 | 52 |
|
| 53 | 53 |
#endif |
| 54 | 54 |
|
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
// This file contains a modified version of the enable_if library from BOOST. |
| 20 | 20 |
// See the appropriate copyright notice below. |
| 21 | 21 |
|
| 22 | 22 |
// Boost enable_if library |
| 23 | 23 |
|
| 24 | 24 |
// Copyright 2003 � The Trustees of Indiana University. |
| 25 | 25 |
|
| 26 | 26 |
// Use, modification, and distribution is subject to the Boost Software |
| 27 | 27 |
// License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at |
| 28 | 28 |
// http://www.boost.org/LICENSE_1_0.txt) |
| 29 | 29 |
|
| 30 | 30 |
// Authors: Jaakko J�rvi (jajarvi at osl.iu.edu) |
| 31 | 31 |
// Jeremiah Willcock (jewillco at osl.iu.edu) |
| 32 | 32 |
// Andrew Lumsdaine (lums at osl.iu.edu) |
| 33 | 33 |
|
| 34 | 34 |
|
| 35 | 35 |
#ifndef LEMON_BITS_UTILITY_H |
| 36 | 36 |
#define LEMON_BITS_UTILITY_H |
| 37 | 37 |
|
| 38 | 38 |
///\file |
| 39 | 39 |
///\brief Miscellaneous basic utilities |
| 40 | 40 |
/// |
| 41 | 41 |
///\todo Please rethink the organisation of the basic files like this. |
| 42 | 42 |
///E.g. this file might be merged with invalid.h. |
| 43 | 43 |
|
| 44 | 44 |
|
| 45 | 45 |
namespace lemon |
| 46 | 46 |
{
|
| 47 | 47 |
|
| 48 | 48 |
/// Basic type for defining "tags". A "YES" condition for \c enable_if. |
| 49 | 49 |
|
| 50 | 50 |
/// Basic type for defining "tags". A "YES" condition for \c enable_if. |
| 51 | 51 |
/// |
| 52 | 52 |
///\sa False |
| 53 | 53 |
/// |
| 54 | 54 |
/// \todo This should go to a separate "basic_types.h" (or something) |
| 55 | 55 |
/// file. |
| 56 | 56 |
struct True {
|
| 57 | 57 |
///\e |
| 58 | 58 |
static const bool value = true; |
| 59 | 59 |
}; |
| 60 | 60 |
|
| 61 | 61 |
/// Basic type for defining "tags". A "NO" condition for \c enable_if. |
| 62 | 62 |
|
| 63 | 63 |
/// Basic type for defining "tags". A "NO" condition for \c enable_if. |
| 64 | 64 |
/// |
| 65 | 65 |
///\sa True |
| 66 | 66 |
struct False {
|
| 67 | 67 |
///\e |
| 68 | 68 |
static const bool value = false; |
| 69 | 69 |
}; |
| 70 | 70 |
|
| 71 | 71 |
|
| 72 | 72 |
struct InvalidType {
|
| 73 | 73 |
}; |
| 74 | 74 |
|
| 75 | 75 |
template <typename T> |
| 76 | 76 |
struct Wrap {
|
| 77 | 77 |
const T &value; |
| 78 | 78 |
Wrap(const T &t) : value(t) {}
|
| 79 | 79 |
}; |
| 80 | 80 |
|
| 81 | 81 |
/**************** dummy class to avoid ambiguity ****************/ |
| 82 | 82 |
|
| 83 | 83 |
template<int T> struct dummy { dummy(int) {} };
|
| 84 | 84 |
|
| 85 | 85 |
/**************** enable_if from BOOST ****************/ |
| 86 | 86 |
|
| 87 | 87 |
template <typename Type, typename T = void> |
| 88 | 88 |
struct exists {
|
| 89 | 89 |
typedef T type; |
| 90 | 90 |
}; |
| 91 | 91 |
|
| 92 | 92 |
|
| 93 | 93 |
template <bool B, class T = void> |
| 94 | 94 |
struct enable_if_c {
|
| 95 | 95 |
typedef T type; |
| 96 | 96 |
}; |
| 97 | 97 |
|
| 98 | 98 |
template <class T> |
| 99 | 99 |
struct enable_if_c<false, T> {};
|
| 100 | 100 |
|
| 101 | 101 |
template <class Cond, class T = void> |
| 102 | 102 |
struct enable_if : public enable_if_c<Cond::value, T> {};
|
| 103 | 103 |
|
| 104 | 104 |
template <bool B, class T> |
| 105 | 105 |
struct lazy_enable_if_c {
|
| 106 | 106 |
typedef typename T::type type; |
| 107 | 107 |
}; |
| 108 | 108 |
|
| 109 | 109 |
template <class T> |
| 110 | 110 |
struct lazy_enable_if_c<false, T> {};
|
| 111 | 111 |
|
| 112 | 112 |
template <class Cond, class T> |
| 113 | 113 |
struct lazy_enable_if : public lazy_enable_if_c<Cond::value, T> {};
|
| 114 | 114 |
|
| 115 | 115 |
|
| 116 | 116 |
template <bool B, class T = void> |
| 117 | 117 |
struct disable_if_c {
|
| 118 | 118 |
typedef T type; |
| 119 | 119 |
}; |
| 120 | 120 |
|
| 121 | 121 |
template <class T> |
| 122 | 122 |
struct disable_if_c<true, T> {};
|
| 123 | 123 |
|
| 124 | 124 |
template <class Cond, class T = void> |
| 125 | 125 |
struct disable_if : public disable_if_c<Cond::value, T> {};
|
| 126 | 126 |
|
| 127 | 127 |
template <bool B, class T> |
| 128 | 128 |
struct lazy_disable_if_c {
|
| 129 | 129 |
typedef typename T::type type; |
| 130 | 130 |
}; |
| 131 | 131 |
|
| 132 | 132 |
template <class T> |
| 133 | 133 |
struct lazy_disable_if_c<true, T> {};
|
| 134 | 134 |
|
| 135 | 135 |
template <class Cond, class T> |
| 136 | 136 |
struct lazy_disable_if : public lazy_disable_if_c<Cond::value, T> {};
|
| 137 | 137 |
|
| 138 | 138 |
} // namespace lemon |
| 139 | 139 |
|
| 140 | 140 |
#endif |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
// This file contains a modified version of the concept checking |
| 20 | 20 |
// utility from BOOST. |
| 21 | 21 |
// See the appropriate copyright notice below. |
| 22 | 22 |
|
| 23 | 23 |
// (C) Copyright Jeremy Siek 2000. |
| 24 | 24 |
// Distributed under the Boost Software License, Version 1.0. (See |
| 25 | 25 |
// accompanying file LICENSE_1_0.txt or copy at |
| 26 | 26 |
// http://www.boost.org/LICENSE_1_0.txt) |
| 27 | 27 |
// |
| 28 | 28 |
// Revision History: |
| 29 | 29 |
// 05 May 2001: Workarounds for HP aCC from Thomas Matelich. (Jeremy Siek) |
| 30 | 30 |
// 02 April 2001: Removed limits header altogether. (Jeremy Siek) |
| 31 | 31 |
// 01 April 2001: Modified to use new <boost/limits.hpp> header. (JMaddock) |
| 32 | 32 |
// |
| 33 | 33 |
|
| 34 | 34 |
// See http://www.boost.org/libs/concept_check for documentation. |
| 35 | 35 |
|
| 36 | 36 |
#ifndef LEMON_CONCEPT_CHECKS_H |
| 37 | 37 |
#define LEMON_CONCEPT_CHECKS_H |
| 38 | 38 |
|
| 39 | 39 |
namespace lemon {
|
| 40 | 40 |
|
| 41 | 41 |
/* |
| 42 | 42 |
"inline" is used for ignore_unused_variable_warning() |
| 43 | 43 |
and function_requires() to make sure there is no |
| 44 | 44 |
overtarget with g++. |
| 45 | 45 |
*/ |
| 46 | 46 |
|
| 47 | 47 |
template <class T> inline void ignore_unused_variable_warning(const T&) { }
|
| 48 | 48 |
|
| 49 | 49 |
template <class Concept> |
| 50 | 50 |
inline void function_requires() |
| 51 | 51 |
{
|
| 52 | 52 |
#if !defined(NDEBUG) |
| 53 | 53 |
void (Concept::*x)() = & Concept::constraints; |
| 54 | 54 |
ignore_unused_variable_warning(x); |
| 55 | 55 |
#endif |
| 56 | 56 |
} |
| 57 | 57 |
|
| 58 | 58 |
template <typename Concept, typename Type> |
| 59 | 59 |
inline void checkConcept() {
|
| 60 | 60 |
#if !defined(NDEBUG) |
| 61 | 61 |
typedef typename Concept::template Constraints<Type> ConceptCheck; |
| 62 | 62 |
void (ConceptCheck::*x)() = & ConceptCheck::constraints; |
| 63 | 63 |
ignore_unused_variable_warning(x); |
| 64 | 64 |
#endif |
| 65 | 65 |
} |
| 66 | 66 |
|
| 67 | 67 |
#define BOOST_CLASS_REQUIRE(type_var, ns, concept) \ |
| 68 | 68 |
typedef void (ns::concept <type_var>::* func##type_var##concept)(); \ |
| 69 | 69 |
template <func##type_var##concept Tp1_> \ |
| 70 | 70 |
struct concept_checking_##type_var##concept { }; \
|
| 71 | 71 |
typedef concept_checking_##type_var##concept< \ |
| 72 | 72 |
BOOST_FPTR ns::concept<type_var>::constraints> \ |
| 73 | 73 |
concept_checking_typedef_##type_var##concept |
| 74 | 74 |
|
| 75 | 75 |
#define BOOST_CLASS_REQUIRE2(type_var1, type_var2, ns, concept) \ |
| 76 | 76 |
typedef void (ns::concept <type_var1,type_var2>::* \ |
| 77 | 77 |
func##type_var1##type_var2##concept)(); \ |
| 78 | 78 |
template <func##type_var1##type_var2##concept Tp1_> \ |
| 79 | 79 |
struct concept_checking_##type_var1##type_var2##concept { }; \
|
| 80 | 80 |
typedef concept_checking_##type_var1##type_var2##concept< \ |
| 81 | 81 |
BOOST_FPTR ns::concept<type_var1,type_var2>::constraints> \ |
| 82 | 82 |
concept_checking_typedef_##type_var1##type_var2##concept |
| 83 | 83 |
|
| 84 | 84 |
#define BOOST_CLASS_REQUIRE3(tv1, tv2, tv3, ns, concept) \ |
| 85 | 85 |
typedef void (ns::concept <tv1,tv2,tv3>::* \ |
| 86 | 86 |
func##tv1##tv2##tv3##concept)(); \ |
| 87 | 87 |
template <func##tv1##tv2##tv3##concept Tp1_> \ |
| 88 | 88 |
struct concept_checking_##tv1##tv2##tv3##concept { }; \
|
| 89 | 89 |
typedef concept_checking_##tv1##tv2##tv3##concept< \ |
| 90 | 90 |
BOOST_FPTR ns::concept<tv1,tv2,tv3>::constraints> \ |
| 91 | 91 |
concept_checking_typedef_##tv1##tv2##tv3##concept |
| 92 | 92 |
|
| 93 | 93 |
#define BOOST_CLASS_REQUIRE4(tv1, tv2, tv3, tv4, ns, concept) \ |
| 94 | 94 |
typedef void (ns::concept <tv1,tv2,tv3,tv4>::* \ |
| 95 | 95 |
func##tv1##tv2##tv3##tv4##concept)(); \ |
| 96 | 96 |
template <func##tv1##tv2##tv3##tv4##concept Tp1_> \ |
| 97 | 97 |
struct concept_checking_##tv1##tv2##tv3##tv4##concept { }; \
|
| 98 | 98 |
typedef concept_checking_##tv1##tv2##tv3##tv4##concept< \ |
| 99 | 99 |
BOOST_FPTR ns::concept<tv1,tv2,tv3,tv4>::constraints> \ |
| 100 | 100 |
concept_checking_typedef_##tv1##tv2##tv3##tv4##concept |
| 101 | 101 |
|
| 102 | 102 |
|
| 103 | 103 |
} // namespace lemon |
| 104 | 104 |
|
| 105 | 105 |
#endif // LEMON_CONCEPT_CHECKS_H |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#ifndef LEMON_CONCEPT_MAPS_H |
| 20 | 20 |
#define LEMON_CONCEPT_MAPS_H |
| 21 | 21 |
|
| 22 | 22 |
#include <lemon/bits/utility.h> |
| 23 | 23 |
#include <lemon/concept_check.h> |
| 24 | 24 |
|
| 25 | 25 |
///\ingroup concept |
| 26 | 26 |
///\file |
| 27 | 27 |
///\brief Map concepts checking classes for testing and documenting. |
| 28 | 28 |
|
| 29 | 29 |
namespace lemon {
|
| 30 | 30 |
|
| 31 | 31 |
namespace concepts {
|
| 32 | 32 |
|
| 33 | 33 |
/// \addtogroup concept |
| 34 | 34 |
/// @{
|
| 35 | 35 |
|
| 36 | 36 |
/// Readable map concept |
| 37 | 37 |
|
| 38 | 38 |
/// Readable map concept. |
| 39 | 39 |
/// |
| 40 | 40 |
template<typename K, typename T> |
| 41 | 41 |
class ReadMap |
| 42 | 42 |
{
|
| 43 | 43 |
public: |
| 44 | 44 |
/// The key type of the map. |
| 45 | 45 |
typedef K Key; |
| 46 | 46 |
/// The value type of the map. (The type of objects associated with the keys). |
| 47 | 47 |
typedef T Value; |
| 48 | 48 |
|
| 49 | 49 |
/// Returns the value associated with a key. |
| 50 | 50 |
|
| 51 | 51 |
/// \bug Value shouldn't need to be default constructible. |
| 52 | 52 |
/// |
| 53 | 53 |
Value operator[](const Key &) const {return Value();}
|
| 54 | 54 |
|
| 55 | 55 |
template<typename _ReadMap> |
| 56 | 56 |
struct Constraints {
|
| 57 | 57 |
|
| 58 | 58 |
void constraints() {
|
| 59 | 59 |
Value val = m[key]; |
| 60 | 60 |
val = m[key]; |
| 61 | 61 |
typename _ReadMap::Value own_val = m[own_key]; |
| 62 | 62 |
own_val = m[own_key]; |
| 63 | 63 |
|
| 64 | 64 |
ignore_unused_variable_warning(val); |
| 65 | 65 |
ignore_unused_variable_warning(own_val); |
| 66 | 66 |
ignore_unused_variable_warning(key); |
| 67 | 67 |
} |
| 68 | 68 |
Key& key; |
| 69 | 69 |
typename _ReadMap::Key& own_key; |
| 70 | 70 |
_ReadMap& m; |
| 71 | 71 |
}; |
| 72 | 72 |
|
| 73 | 73 |
}; |
| 74 | 74 |
|
| 75 | 75 |
|
| 76 | 76 |
/// Writable map concept |
| 77 | 77 |
|
| 78 | 78 |
/// Writable map concept. |
| 79 | 79 |
/// |
| 80 | 80 |
template<typename K, typename T> |
| 81 | 81 |
class WriteMap |
| 82 | 82 |
{
|
| 83 | 83 |
public: |
| 84 | 84 |
/// The key type of the map. |
| 85 | 85 |
typedef K Key; |
| 86 | 86 |
/// The value type of the map. (The type of objects associated with the keys). |
| 87 | 87 |
typedef T Value; |
| 88 | 88 |
|
| 89 | 89 |
/// Sets the value associated with a key. |
| 90 | 90 |
void set(const Key &,const Value &) {}
|
| 91 | 91 |
|
| 92 | 92 |
///Default constructor |
| 93 | 93 |
WriteMap() {}
|
| 94 | 94 |
|
| 95 | 95 |
template <typename _WriteMap> |
| 96 | 96 |
struct Constraints {
|
| 97 | 97 |
void constraints() {
|
| 98 | 98 |
// No constraints for constructor. |
| 99 | 99 |
m.set(key, val); |
| 100 | 100 |
m.set(own_key, own_val); |
| 101 | 101 |
ignore_unused_variable_warning(key); |
| 102 | 102 |
ignore_unused_variable_warning(val); |
| 103 | 103 |
ignore_unused_variable_warning(own_key); |
| 104 | 104 |
ignore_unused_variable_warning(own_val); |
| 105 | 105 |
} |
| 106 | 106 |
|
| 107 | 107 |
Value& val; |
| 108 | 108 |
typename _WriteMap::Value own_val; |
| 109 | 109 |
Key& key; |
| 110 | 110 |
typename _WriteMap::Key& own_key; |
| 111 | 111 |
_WriteMap& m; |
| 112 | 112 |
|
| 113 | 113 |
}; |
| 114 | 114 |
}; |
| 115 | 115 |
|
| 116 | 116 |
/// Read/Writable map concept |
| 117 | 117 |
|
| 118 | 118 |
/// Read/writable map concept. |
| 119 | 119 |
/// |
| 120 | 120 |
template<typename K, typename T> |
| 121 | 121 |
class ReadWriteMap : public ReadMap<K,T>, |
| 122 | 122 |
public WriteMap<K,T> |
| 123 | 123 |
{
|
| 124 | 124 |
public: |
| 125 | 125 |
/// The key type of the map. |
| 126 | 126 |
typedef K Key; |
| 127 | 127 |
/// The value type of the map. (The type of objects associated with the keys). |
| 128 | 128 |
typedef T Value; |
| 129 | 129 |
|
| 130 | 130 |
/// Returns the value associated with a key. |
| 131 | 131 |
Value operator[](const Key &) const {return Value();}
|
| 132 | 132 |
/// Sets the value associated with a key. |
| 133 | 133 |
void set(const Key & ,const Value &) {}
|
| 134 | 134 |
|
| 135 | 135 |
template<typename _ReadWriteMap> |
| 136 | 136 |
struct Constraints {
|
| 137 | 137 |
void constraints() {
|
| 138 | 138 |
checkConcept<ReadMap<K, T>, _ReadWriteMap >(); |
| 139 | 139 |
checkConcept<WriteMap<K, T>, _ReadWriteMap >(); |
| 140 | 140 |
} |
| 141 | 141 |
}; |
| 142 | 142 |
}; |
| 143 | 143 |
|
| 144 | 144 |
|
| 145 | 145 |
/// Dereferable map concept |
| 146 | 146 |
|
| 147 | 147 |
/// Dereferable map concept. |
| 148 | 148 |
/// |
| 149 | 149 |
template<typename K, typename T, typename R, typename CR> |
| 150 | 150 |
class ReferenceMap : public ReadWriteMap<K,T> |
| 151 | 151 |
{
|
| 152 | 152 |
public: |
| 153 | 153 |
/// Tag for reference maps. |
| 154 | 154 |
typedef True ReferenceMapTag; |
| 155 | 155 |
/// The key type of the map. |
| 156 | 156 |
typedef K Key; |
| 157 | 157 |
/// The value type of the map. (The type of objects associated with the keys). |
| 158 | 158 |
typedef T Value; |
| 159 | 159 |
/// The reference type of the map. |
| 160 | 160 |
typedef R Reference; |
| 161 | 161 |
/// The const reference type of the map. |
| 162 | 162 |
typedef CR ConstReference; |
| 163 | 163 |
|
| 164 | 164 |
protected: |
| 165 | 165 |
Value tmp; |
| 166 | 166 |
public: |
| 167 | 167 |
|
| 168 | 168 |
///Returns a reference to the value associated to a key. |
| 169 | 169 |
Reference operator[](const Key &) { return tmp; }
|
| 170 | 170 |
///Returns a const reference to the value associated to a key. |
| 171 | 171 |
ConstReference operator[](const Key &) const { return tmp; }
|
| 172 | 172 |
/// Sets the value associated with a key. |
| 173 | 173 |
void set(const Key &k,const Value &t) { operator[](k)=t; }
|
| 174 | 174 |
|
| 175 | 175 |
/// \todo Rethink this concept. |
| 176 | 176 |
template<typename _ReferenceMap> |
| 177 | 177 |
struct ReferenceMapConcept {
|
| 178 | 178 |
|
| 179 | 179 |
void constraints() {
|
| 180 | 180 |
checkConcept<ReadWriteMap, _ReferenceMap >(); |
| 181 | 181 |
m[key] = val; |
| 182 | 182 |
val = m[key]; |
| 183 | 183 |
m[key] = ref; |
| 184 | 184 |
ref = m[key]; |
| 185 | 185 |
m[own_key] = own_val; |
| 186 | 186 |
own_val = m[own_key]; |
| 187 | 187 |
m[own_key] = own_ref; |
| 188 | 188 |
own_ref = m[own_key]; |
| 189 | 189 |
} |
| 190 | 190 |
|
| 191 | 191 |
typename _ReferenceMap::Key& own_key; |
| 192 | 192 |
typename _ReferenceMap::Value& own_val; |
| 193 | 193 |
typename _ReferenceMap::Reference& own_ref; |
| 194 | 194 |
Key& key; |
| 195 | 195 |
Value& val; |
| 196 | 196 |
Reference& ref; |
| 197 | 197 |
_ReferenceMap& m; |
| 198 | 198 |
}; |
| 199 | 199 |
}; |
| 200 | 200 |
|
| 201 | 201 |
// @} |
| 202 | 202 |
|
| 203 | 203 |
} //namespace concepts |
| 204 | 204 |
|
| 205 | 205 |
} //namespace lemon |
| 206 | 206 |
|
| 207 | 207 |
#endif // LEMON_CONCEPT_MAPS_H |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#ifndef LEMON_DIM2_H |
| 20 | 20 |
#define LEMON_DIM2_H |
| 21 | 21 |
|
| 22 | 22 |
#include <iostream> |
| 23 | 23 |
#include <lemon/bits/utility.h> |
| 24 | 24 |
|
| 25 | 25 |
///\ingroup misc |
| 26 | 26 |
///\file |
| 27 | 27 |
///\brief A simple two dimensional vector and a bounding box implementation |
| 28 | 28 |
/// |
| 29 | 29 |
/// The class \ref lemon::dim2::Point "dim2::Point" implements |
| 30 | 30 |
///a two dimensional vector with the usual |
| 31 | 31 |
/// operations. |
| 32 | 32 |
/// |
| 33 | 33 |
/// The class \ref lemon::dim2::BoundingBox "dim2::BoundingBox" |
| 34 | 34 |
/// can be used to determine |
| 35 | 35 |
/// the rectangular bounding box of a set of |
| 36 | 36 |
/// \ref lemon::dim2::Point "dim2::Point"'s. |
| 37 | 37 |
|
| 38 | 38 |
namespace lemon {
|
| 39 | 39 |
|
| 40 | 40 |
///Tools for handling two dimensional coordinates |
| 41 | 41 |
|
| 42 | 42 |
///This namespace is a storage of several |
| 43 | 43 |
///tools for handling two dimensional coordinates |
| 44 | 44 |
namespace dim2 {
|
| 45 | 45 |
|
| 46 | 46 |
/// \addtogroup misc |
| 47 | 47 |
/// @{
|
| 48 | 48 |
|
| 49 | 49 |
/// A simple two dimensional vector (plainvector) implementation |
| 50 | 50 |
|
| 51 | 51 |
/// A simple two dimensional vector (plainvector) implementation |
| 52 | 52 |
///with the usual vector |
| 53 | 53 |
/// operators. |
| 54 | 54 |
/// |
| 55 | 55 |
template<typename T> |
| 56 | 56 |
class Point {
|
| 57 | 57 |
|
| 58 | 58 |
public: |
| 59 | 59 |
|
| 60 | 60 |
typedef T Value; |
| 61 | 61 |
|
| 62 | 62 |
///First coordinate |
| 63 | 63 |
T x; |
| 64 | 64 |
///Second coordinate |
| 65 | 65 |
T y; |
| 66 | 66 |
|
| 67 | 67 |
///Default constructor |
| 68 | 68 |
Point() {}
|
| 69 | 69 |
|
| 70 | 70 |
///Construct an instance from coordinates |
| 71 | 71 |
Point(T a, T b) : x(a), y(b) { }
|
| 72 | 72 |
|
| 73 | 73 |
///The dimension of the vector. |
| 74 | 74 |
|
| 75 | 75 |
///The dimension of the vector. |
| 76 | 76 |
///This function always returns 2. |
| 77 | 77 |
int size() const { return 2; }
|
| 78 | 78 |
|
| 79 | 79 |
///Subscripting operator |
| 80 | 80 |
|
| 81 | 81 |
///\c p[0] is \c p.x and \c p[1] is \c p.y |
| 82 | 82 |
/// |
| 83 | 83 |
T& operator[](int idx) { return idx == 0 ? x : y; }
|
| 84 | 84 |
|
| 85 | 85 |
///Const subscripting operator |
| 86 | 86 |
|
| 87 | 87 |
///\c p[0] is \c p.x and \c p[1] is \c p.y |
| 88 | 88 |
/// |
| 89 | 89 |
const T& operator[](int idx) const { return idx == 0 ? x : y; }
|
| 90 | 90 |
|
| 91 | 91 |
///Conversion constructor |
| 92 | 92 |
template<class TT> Point(const Point<TT> &p) : x(p.x), y(p.y) {}
|
| 93 | 93 |
|
| 94 | 94 |
///Give back the square of the norm of the vector |
| 95 | 95 |
T normSquare() const {
|
| 96 | 96 |
return x*x+y*y; |
| 97 | 97 |
} |
| 98 | 98 |
|
| 99 | 99 |
///Increment the left hand side by u |
| 100 | 100 |
Point<T>& operator +=(const Point<T>& u) {
|
| 101 | 101 |
x += u.x; |
| 102 | 102 |
y += u.y; |
| 103 | 103 |
return *this; |
| 104 | 104 |
} |
| 105 | 105 |
|
| 106 | 106 |
///Decrement the left hand side by u |
| 107 | 107 |
Point<T>& operator -=(const Point<T>& u) {
|
| 108 | 108 |
x -= u.x; |
| 109 | 109 |
y -= u.y; |
| 110 | 110 |
return *this; |
| 111 | 111 |
} |
| 112 | 112 |
|
| 113 | 113 |
///Multiply the left hand side with a scalar |
| 114 | 114 |
Point<T>& operator *=(const T &u) {
|
| 115 | 115 |
x *= u; |
| 116 | 116 |
y *= u; |
| 117 | 117 |
return *this; |
| 118 | 118 |
} |
| 119 | 119 |
|
| 120 | 120 |
///Divide the left hand side by a scalar |
| 121 | 121 |
Point<T>& operator /=(const T &u) {
|
| 122 | 122 |
x /= u; |
| 123 | 123 |
y /= u; |
| 124 | 124 |
return *this; |
| 125 | 125 |
} |
| 126 | 126 |
|
| 127 | 127 |
///Return the scalar product of two vectors |
| 128 | 128 |
T operator *(const Point<T>& u) const {
|
| 129 | 129 |
return x*u.x+y*u.y; |
| 130 | 130 |
} |
| 131 | 131 |
|
| 132 | 132 |
///Return the sum of two vectors |
| 133 | 133 |
Point<T> operator+(const Point<T> &u) const {
|
| 134 | 134 |
Point<T> b=*this; |
| 135 | 135 |
return b+=u; |
| 136 | 136 |
} |
| 137 | 137 |
|
| 138 | 138 |
///Return the negative of the vector |
| 139 | 139 |
Point<T> operator-() const {
|
| 140 | 140 |
Point<T> b=*this; |
| 141 | 141 |
b.x=-b.x; b.y=-b.y; |
| 142 | 142 |
return b; |
| 143 | 143 |
} |
| 144 | 144 |
|
| 145 | 145 |
///Return the difference of two vectors |
| 146 | 146 |
Point<T> operator-(const Point<T> &u) const {
|
| 147 | 147 |
Point<T> b=*this; |
| 148 | 148 |
return b-=u; |
| 149 | 149 |
} |
| 150 | 150 |
|
| 151 | 151 |
///Return a vector multiplied by a scalar |
| 152 | 152 |
Point<T> operator*(const T &u) const {
|
| 153 | 153 |
Point<T> b=*this; |
| 154 | 154 |
return b*=u; |
| 155 | 155 |
} |
| 156 | 156 |
|
| 157 | 157 |
///Return a vector divided by a scalar |
| 158 | 158 |
Point<T> operator/(const T &u) const {
|
| 159 | 159 |
Point<T> b=*this; |
| 160 | 160 |
return b/=u; |
| 161 | 161 |
} |
| 162 | 162 |
|
| 163 | 163 |
///Test equality |
| 164 | 164 |
bool operator==(const Point<T> &u) const {
|
| 165 | 165 |
return (x==u.x) && (y==u.y); |
| 166 | 166 |
} |
| 167 | 167 |
|
| 168 | 168 |
///Test inequality |
| 169 | 169 |
bool operator!=(Point u) const {
|
| 170 | 170 |
return (x!=u.x) || (y!=u.y); |
| 171 | 171 |
} |
| 172 | 172 |
|
| 173 | 173 |
}; |
| 174 | 174 |
|
| 175 | 175 |
///Return a Point |
| 176 | 176 |
|
| 177 | 177 |
///Return a Point. |
| 178 | 178 |
///\relates Point |
| 179 | 179 |
template <typename T> |
| 180 | 180 |
inline Point<T> makePoint(const T& x, const T& y) {
|
| 181 | 181 |
return Point<T>(x, y); |
| 182 | 182 |
} |
| 183 | 183 |
|
| 184 | 184 |
///Return a vector multiplied by a scalar |
| 185 | 185 |
|
| 186 | 186 |
///Return a vector multiplied by a scalar. |
| 187 | 187 |
///\relates Point |
| 188 | 188 |
template<typename T> Point<T> operator*(const T &u,const Point<T> &x) {
|
| 189 | 189 |
return x*u; |
| 190 | 190 |
} |
| 191 | 191 |
|
| 192 | 192 |
///Read a plainvector from a stream |
| 193 | 193 |
|
| 194 | 194 |
///Read a plainvector from a stream. |
| 195 | 195 |
///\relates Point |
| 196 | 196 |
/// |
| 197 | 197 |
template<typename T> |
| 198 | 198 |
inline std::istream& operator>>(std::istream &is, Point<T> &z) {
|
| 199 | 199 |
char c; |
| 200 | 200 |
if (is >> c) {
|
| 201 | 201 |
if (c != '(') is.putback(c);
|
| 202 | 202 |
} else {
|
| 203 | 203 |
is.clear(); |
| 204 | 204 |
} |
| 205 | 205 |
if (!(is >> z.x)) return is; |
| 206 | 206 |
if (is >> c) {
|
| 207 | 207 |
if (c != ',') is.putback(c); |
| 208 | 208 |
} else {
|
| 209 | 209 |
is.clear(); |
| 210 | 210 |
} |
| 211 | 211 |
if (!(is >> z.y)) return is; |
| 212 | 212 |
if (is >> c) {
|
| 213 | 213 |
if (c != ')') is.putback(c); |
| 214 | 214 |
} else {
|
| 215 | 215 |
is.clear(); |
| 216 | 216 |
} |
| 217 | 217 |
return is; |
| 218 | 218 |
} |
| 219 | 219 |
|
| 220 | 220 |
///Write a plainvector to a stream |
| 221 | 221 |
|
| 222 | 222 |
///Write a plainvector to a stream. |
| 223 | 223 |
///\relates Point |
| 224 | 224 |
/// |
| 225 | 225 |
template<typename T> |
| 226 | 226 |
inline std::ostream& operator<<(std::ostream &os, const Point<T>& z) |
| 227 | 227 |
{
|
| 228 | 228 |
os << "(" << z.x << ", " << z.y << ")";
|
| 229 | 229 |
return os; |
| 230 | 230 |
} |
| 231 | 231 |
|
| 232 | 232 |
///Rotate by 90 degrees |
| 233 | 233 |
|
| 234 | 234 |
///Returns the parameter rotated by 90 degrees in positive direction. |
| 235 | 235 |
///\relates Point |
| 236 | 236 |
/// |
| 237 | 237 |
template<typename T> |
| 238 | 238 |
inline Point<T> rot90(const Point<T> &z) |
| 239 | 239 |
{
|
| 240 | 240 |
return Point<T>(-z.y,z.x); |
| 241 | 241 |
} |
| 242 | 242 |
|
| 243 | 243 |
///Rotate by 180 degrees |
| 244 | 244 |
|
| 245 | 245 |
///Returns the parameter rotated by 180 degrees. |
| 246 | 246 |
///\relates Point |
| 247 | 247 |
/// |
| 248 | 248 |
template<typename T> |
| 249 | 249 |
inline Point<T> rot180(const Point<T> &z) |
| 250 | 250 |
{
|
| 251 | 251 |
return Point<T>(-z.x,-z.y); |
| 252 | 252 |
} |
| 253 | 253 |
|
| 254 | 254 |
///Rotate by 270 degrees |
| 255 | 255 |
|
| 256 | 256 |
///Returns the parameter rotated by 90 degrees in negative direction. |
| 257 | 257 |
///\relates Point |
| 258 | 258 |
/// |
| 259 | 259 |
template<typename T> |
| 260 | 260 |
inline Point<T> rot270(const Point<T> &z) |
| 261 | 261 |
{
|
| 262 | 262 |
return Point<T>(z.y,-z.x); |
| 263 | 263 |
} |
| 264 | 264 |
|
| 265 | 265 |
|
| 266 | 266 |
|
| 267 | 267 |
/// A class to calculate or store the bounding box of plainvectors. |
| 268 | 268 |
|
| 269 | 269 |
/// A class to calculate or store the bounding box of plainvectors. |
| 270 | 270 |
/// |
| 271 | 271 |
template<typename T> |
| 272 | 272 |
class BoundingBox {
|
| 273 | 273 |
Point<T> bottom_left, top_right; |
| 274 | 274 |
bool _empty; |
| 275 | 275 |
public: |
| 276 | 276 |
|
| 277 | 277 |
///Default constructor: creates an empty bounding box |
| 278 | 278 |
BoundingBox() { _empty = true; }
|
| 279 | 279 |
|
| 280 | 280 |
///Construct an instance from one point |
| 281 | 281 |
BoundingBox(Point<T> a) { bottom_left=top_right=a; _empty = false; }
|
| 282 | 282 |
|
| 283 | 283 |
///Construct an instance from two points |
| 284 | 284 |
|
| 285 | 285 |
///Construct an instance from two points. |
| 286 | 286 |
///\param a The bottom left corner. |
| 287 | 287 |
///\param b The top right corner. |
| 288 | 288 |
///\warning The coordinates of the bottom left corner must be no more |
| 289 | 289 |
///than those of the top right one. |
| 290 | 290 |
BoundingBox(Point<T> a,Point<T> b) |
| 291 | 291 |
{
|
| 292 | 292 |
bottom_left=a; |
| 293 | 293 |
top_right=b; |
| 294 | 294 |
_empty = false; |
| 295 | 295 |
} |
| 296 | 296 |
|
| 297 | 297 |
///Construct an instance from four numbers |
| 298 | 298 |
|
| 299 | 299 |
///Construct an instance from four numbers. |
| 300 | 300 |
///\param l The left side of the box. |
| 301 | 301 |
///\param b The bottom of the box. |
| 302 | 302 |
///\param r The right side of the box. |
| 303 | 303 |
///\param t The top of the box. |
| 304 | 304 |
///\warning The left side must be no more than the right side and |
| 305 | 305 |
///bottom must be no more than the top. |
| 306 | 306 |
BoundingBox(T l,T b,T r,T t) |
| 307 | 307 |
{
|
| 308 | 308 |
bottom_left=Point<T>(l,b); |
| 309 | 309 |
top_right=Point<T>(r,t); |
| 310 | 310 |
_empty = false; |
| 311 | 311 |
} |
| 312 | 312 |
|
| 313 | 313 |
///Return \c true if the bounding box is empty. |
| 314 | 314 |
|
| 315 | 315 |
///Return \c true if the bounding box is empty (i.e. return \c false |
| 316 | 316 |
///if at least one point was added to the box or the coordinates of |
| 317 | 317 |
///the box were set). |
| 318 | 318 |
///The coordinates of an empty bounding box are not defined. |
| 319 | 319 |
bool empty() const {
|
| 320 | 320 |
return _empty; |
| 321 | 321 |
} |
| 322 | 322 |
|
| 323 | 323 |
///Make the BoundingBox empty |
| 324 | 324 |
void clear() {
|
| 325 | 325 |
_empty=1; |
| 326 | 326 |
} |
| 327 | 327 |
|
| 328 | 328 |
///Give back the bottom left corner |
| 329 | 329 |
|
| 330 | 330 |
///Give back the bottom left corner. |
| 331 | 331 |
///If the bounding box is empty, then the return value is not defined. |
| 332 | 332 |
Point<T> bottomLeft() const {
|
| 333 | 333 |
return bottom_left; |
| 334 | 334 |
} |
| 335 | 335 |
|
| 336 | 336 |
///Set the bottom left corner |
| 337 | 337 |
|
| 338 | 338 |
///Set the bottom left corner. |
| 339 | 339 |
///It should only be used for non-empty box. |
| 340 | 340 |
void bottomLeft(Point<T> p) {
|
| 341 | 341 |
bottom_left = p; |
| 342 | 342 |
} |
| 343 | 343 |
|
| 344 | 344 |
///Give back the top right corner |
| 345 | 345 |
|
| 346 | 346 |
///Give back the top right corner. |
| 347 | 347 |
///If the bounding box is empty, then the return value is not defined. |
| 348 | 348 |
Point<T> topRight() const {
|
| 349 | 349 |
return top_right; |
| 350 | 350 |
} |
| 351 | 351 |
|
| 352 | 352 |
///Set the top right corner |
| 353 | 353 |
|
| 354 | 354 |
///Set the top right corner. |
| 355 | 355 |
///It should only be used for non-empty box. |
| 356 | 356 |
void topRight(Point<T> p) {
|
| 357 | 357 |
top_right = p; |
| 358 | 358 |
} |
| 359 | 359 |
|
| 360 | 360 |
///Give back the bottom right corner |
| 361 | 361 |
|
| 362 | 362 |
///Give back the bottom right corner. |
| 363 | 363 |
///If the bounding box is empty, then the return value is not defined. |
| 364 | 364 |
Point<T> bottomRight() const {
|
| 365 | 365 |
return Point<T>(top_right.x,bottom_left.y); |
| 366 | 366 |
} |
| 367 | 367 |
|
| 368 | 368 |
///Set the bottom right corner |
| 369 | 369 |
|
| 370 | 370 |
///Set the bottom right corner. |
| 371 | 371 |
///It should only be used for non-empty box. |
| 372 | 372 |
void bottomRight(Point<T> p) {
|
| 373 | 373 |
top_right.x = p.x; |
| 374 | 374 |
bottom_left.y = p.y; |
| 375 | 375 |
} |
| 376 | 376 |
|
| 377 | 377 |
///Give back the top left corner |
| 378 | 378 |
|
| 379 | 379 |
///Give back the top left corner. |
| 380 | 380 |
///If the bounding box is empty, then the return value is not defined. |
| 381 | 381 |
Point<T> topLeft() const {
|
| 382 | 382 |
return Point<T>(bottom_left.x,top_right.y); |
| 383 | 383 |
} |
| 384 | 384 |
|
| 385 | 385 |
///Set the top left corner |
| 386 | 386 |
|
| 387 | 387 |
///Set the top left corner. |
| 388 | 388 |
///It should only be used for non-empty box. |
| 389 | 389 |
void topLeft(Point<T> p) {
|
| 1 |
/* -*- C++ -*- |
|
| 2 |
* |
|
| 3 |
* This file is a part of LEMON, a generic C++ optimization library |
|
| 4 |
* |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 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 |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#ifndef LEMON_MAPS_H |
| 20 | 20 |
#define LEMON_MAPS_H |
| 21 | 21 |
|
| 22 | 22 |
#include <iterator> |
| 23 | 23 |
#include <functional> |
| 24 | 24 |
#include <vector> |
| 25 | 25 |
|
| 26 | 26 |
#include <lemon/bits/utility.h> |
| 27 | 27 |
// #include <lemon/bits/traits.h> |
| 28 | 28 |
|
| 29 | 29 |
///\file |
| 30 | 30 |
///\ingroup maps |
| 31 | 31 |
///\brief Miscellaneous property maps |
| 32 | 32 |
/// |
| 33 | 33 |
#include <map> |
| 34 | 34 |
|
| 35 | 35 |
namespace lemon {
|
| 36 | 36 |
|
| 37 | 37 |
/// \addtogroup maps |
| 38 | 38 |
/// @{
|
| 39 | 39 |
|
| 40 | 40 |
/// Base class of maps. |
| 41 | 41 |
|
| 42 | 42 |
/// Base class of maps. |
| 43 | 43 |
/// It provides the necessary <tt>typedef</tt>s required by the map concept. |
| 44 | 44 |
template<typename K, typename T> |
| 45 | 45 |
class MapBase {
|
| 46 | 46 |
public: |
| 47 | 47 |
/// The key type of the map. |
| 48 | 48 |
typedef K Key; |
| 49 | 49 |
/// The value type of the map. (The type of objects associated with the keys). |
| 50 | 50 |
typedef T Value; |
| 51 | 51 |
}; |
| 52 | 52 |
|
| 53 | 53 |
/// Null map. (a.k.a. DoNothingMap) |
| 54 | 54 |
|
| 55 | 55 |
/// This map can be used if you have to provide a map only for |
| 56 | 56 |
/// its type definitions, or if you have to provide a writable map, |
| 57 | 57 |
/// but data written to it is not required (i.e. it will be sent to |
| 58 | 58 |
/// <tt>/dev/null</tt>). |
| 59 | 59 |
template<typename K, typename T> |
| 60 | 60 |
class NullMap : public MapBase<K, T> {
|
| 61 | 61 |
public: |
| 62 | 62 |
typedef MapBase<K, T> Parent; |
| 63 | 63 |
typedef typename Parent::Key Key; |
| 64 | 64 |
typedef typename Parent::Value Value; |
| 65 | 65 |
|
| 66 | 66 |
/// Gives back a default constructed element. |
| 67 | 67 |
T operator[](const K&) const { return T(); }
|
| 68 | 68 |
/// Absorbs the value. |
| 69 | 69 |
void set(const K&, const T&) {}
|
| 70 | 70 |
}; |
| 71 | 71 |
|
| 72 | 72 |
///Returns a \c NullMap class |
| 73 | 73 |
|
| 74 | 74 |
///This function just returns a \c NullMap class. |
| 75 | 75 |
///\relates NullMap |
| 76 | 76 |
template <typename K, typename V> |
| 77 | 77 |
NullMap<K, V> nullMap() {
|
| 78 | 78 |
return NullMap<K, V>(); |
| 79 | 79 |
} |
| 80 | 80 |
|
| 81 | 81 |
|
| 82 | 82 |
/// Constant map. |
| 83 | 83 |
|
| 84 | 84 |
/// This is a readable map which assigns a specified value to each key. |
| 85 | 85 |
/// In other aspects it is equivalent to the \c NullMap. |
| 86 | 86 |
template<typename K, typename T> |
| 87 | 87 |
class ConstMap : public MapBase<K, T> {
|
| 88 | 88 |
private: |
| 89 | 89 |
T v; |
| 90 | 90 |
public: |
| 91 | 91 |
|
| 92 | 92 |
typedef MapBase<K, T> Parent; |
| 93 | 93 |
typedef typename Parent::Key Key; |
| 94 | 94 |
typedef typename Parent::Value Value; |
| 95 | 95 |
|
| 96 | 96 |
/// Default constructor |
| 97 | 97 |
|
| 98 | 98 |
/// Default constructor. |
| 99 | 99 |
/// The value of the map will be uninitialized. |
| 100 | 100 |
/// (More exactly it will be default constructed.) |
| 101 | 101 |
ConstMap() {}
|
| 102 | 102 |
|
| 103 | 103 |
/// Constructor with specified initial value |
| 104 | 104 |
|
| 105 | 105 |
/// Constructor with specified initial value. |
| 106 | 106 |
/// \param _v is the initial value of the map. |
| 107 | 107 |
ConstMap(const T &_v) : v(_v) {}
|
| 108 | 108 |
|
| 109 | 109 |
///\e |
| 110 | 110 |
T operator[](const K&) const { return v; }
|
| 111 | 111 |
|
| 112 | 112 |
///\e |
| 113 | 113 |
void setAll(const T &t) {
|
| 114 | 114 |
v = t; |
| 115 | 115 |
} |
| 116 | 116 |
|
| 117 | 117 |
template<typename T1> |
| 118 | 118 |
struct rebind {
|
| 119 | 119 |
typedef ConstMap<K, T1> other; |
| 120 | 120 |
}; |
| 121 | 121 |
|
| 122 | 122 |
template<typename T1> |
| 123 | 123 |
ConstMap(const ConstMap<K, T1> &, const T &_v) : v(_v) {}
|
| 124 | 124 |
}; |
| 125 | 125 |
|
| 126 | 126 |
///Returns a \c ConstMap class |
| 127 | 127 |
|
| 128 | 128 |
///This function just returns a \c ConstMap class. |
| 129 | 129 |
///\relates ConstMap |
| 130 | 130 |
template<typename K, typename V> |
| 131 | 131 |
inline ConstMap<K, V> constMap(const V &v) {
|
| 132 | 132 |
return ConstMap<K, V>(v); |
| 133 | 133 |
} |
| 134 | 134 |
|
| 135 | 135 |
|
| 136 | 136 |
template<typename T, T v> |
| 137 | 137 |
struct Const { };
|
| 138 | 138 |
|
| 139 | 139 |
/// Constant map with inlined constant value. |
| 140 | 140 |
|
| 141 | 141 |
/// This is a readable map which assigns a specified value to each key. |
| 142 | 142 |
/// In other aspects it is equivalent to the \c NullMap. |
| 143 | 143 |
template<typename K, typename V, V v> |
| 144 | 144 |
class ConstMap<K, Const<V, v> > : public MapBase<K, V> {
|
| 145 | 145 |
public: |
| 146 | 146 |
typedef MapBase<K, V> Parent; |
| 147 | 147 |
typedef typename Parent::Key Key; |
| 148 | 148 |
typedef typename Parent::Value Value; |
| 149 | 149 |
|
| 150 | 150 |
ConstMap() { }
|
| 151 | 151 |
///\e |
| 152 | 152 |
V operator[](const K&) const { return v; }
|
| 153 | 153 |
///\e |
| 154 | 154 |
void set(const K&, const V&) { }
|
| 155 | 155 |
}; |
| 156 | 156 |
|
| 157 | 157 |
///Returns a \c ConstMap class |
| 158 | 158 |
|
| 159 | 159 |
///This function just returns a \c ConstMap class with inlined value. |
| 160 | 160 |
///\relates ConstMap |
| 161 | 161 |
template<typename K, typename V, V v> |
| 162 | 162 |
inline ConstMap<K, Const<V, v> > constMap() {
|
| 163 | 163 |
return ConstMap<K, Const<V, v> >(); |
| 164 | 164 |
} |
| 165 | 165 |
|
| 166 | 166 |
///Map based on std::map |
| 167 | 167 |
|
| 168 | 168 |
///This is essentially a wrapper for \c std::map with addition that |
| 169 | 169 |
///you can specify a default value different from \c Value(). |
| 170 | 170 |
template <typename K, typename T, typename Compare = std::less<K> > |
| 171 | 171 |
class StdMap {
|
| 172 | 172 |
template <typename K1, typename T1, typename C1> |
| 173 | 173 |
friend class StdMap; |
| 174 | 174 |
public: |
| 175 | 175 |
|
| 176 | 176 |
typedef True ReferenceMapTag; |
| 177 | 177 |
///\e |
| 178 | 178 |
typedef K Key; |
| 179 | 179 |
///\e |
| 180 | 180 |
typedef T Value; |
| 181 | 181 |
///\e |
| 182 | 182 |
typedef T& Reference; |
| 183 | 183 |
///\e |
| 184 | 184 |
typedef const T& ConstReference; |
| 185 | 185 |
|
| 186 | 186 |
private: |
| 187 | 187 |
|
| 188 | 188 |
typedef std::map<K, T, Compare> Map; |
| 189 | 189 |
Value _value; |
| 190 | 190 |
Map _map; |
| 191 | 191 |
|
| 192 | 192 |
public: |
| 193 | 193 |
|
| 194 | 194 |
/// Constructor with specified default value |
| 195 | 195 |
StdMap(const T& value = T()) : _value(value) {}
|
| 196 | 196 |
/// \brief Constructs the map from an appropriate std::map, and explicitly |
| 197 | 197 |
/// specifies a default value. |
| 198 | 198 |
template <typename T1, typename Comp1> |
| 199 | 199 |
StdMap(const std::map<Key, T1, Comp1> &map, const T& value = T()) |
| 200 | 200 |
: _map(map.begin(), map.end()), _value(value) {}
|
| 201 | 201 |
|
| 202 | 202 |
/// \brief Constructs a map from an other StdMap. |
| 203 | 203 |
template<typename T1, typename Comp1> |
| 204 | 204 |
StdMap(const StdMap<Key, T1, Comp1> &c) |
| 205 | 205 |
: _map(c._map.begin(), c._map.end()), _value(c._value) {}
|
| 206 | 206 |
|
| 207 | 207 |
private: |
| 208 | 208 |
|
| 209 | 209 |
StdMap& operator=(const StdMap&); |
| 210 | 210 |
|
| 211 | 211 |
public: |
| 212 | 212 |
|
| 213 | 213 |
///\e |
| 214 | 214 |
Reference operator[](const Key &k) {
|
| 215 | 215 |
typename Map::iterator it = _map.lower_bound(k); |
| 216 | 216 |
if (it != _map.end() && !_map.key_comp()(k, it->first)) |
| 217 | 217 |
return it->second; |
| 218 | 218 |
else |
| 219 | 219 |
return _map.insert(it, std::make_pair(k, _value))->second; |
| 220 | 220 |
} |
| 221 | 221 |
|
| 222 | 222 |
/// \e |
| 223 | 223 |
ConstReference operator[](const Key &k) const {
|
| 224 | 224 |
typename Map::const_iterator it = _map.find(k); |
| 225 | 225 |
if (it != _map.end()) |
| 226 | 226 |
return it->second; |
| 227 | 227 |
else |
| 228 | 228 |
return _value; |
| 229 | 229 |
} |
| 230 | 230 |
|
| 231 | 231 |
/// \e |
| 232 | 232 |
void set(const Key &k, const T &t) {
|
| 233 | 233 |
typename Map::iterator it = _map.lower_bound(k); |
| 234 | 234 |
if (it != _map.end() && !_map.key_comp()(k, it->first)) |
| 235 | 235 |
it->second = t; |
| 236 | 236 |
else |
| 237 | 237 |
_map.insert(it, std::make_pair(k, t)); |
| 238 | 238 |
} |
| 239 | 239 |
|
| 240 | 240 |
/// \e |
| 241 | 241 |
void setAll(const T &t) {
|
| 242 | 242 |
_value = t; |
| 243 | 243 |
_map.clear(); |
| 244 | 244 |
} |
| 245 | 245 |
|
| 246 | 246 |
template <typename T1, typename C1 = std::less<T1> > |
| 247 | 247 |
struct rebind {
|
| 248 | 248 |
typedef StdMap<Key, T1, C1> other; |
| 249 | 249 |
}; |
| 250 | 250 |
}; |
| 251 | 251 |
|
| 252 | 252 |
/// \brief Map for storing values for keys from the range <tt>[0..size-1]</tt> |
| 253 | 253 |
/// |
| 254 | 254 |
/// The current map has the <tt>[0..size-1]</tt> keyset and the values |
| 255 | 255 |
/// are stored in a \c std::vector<T> container. It can be used with |
| 256 | 256 |
/// some data structures, for example \c UnionFind, \c BinHeap, when |
| 257 | 257 |
/// the used items are small integer numbers. |
| 258 | 258 |
/// |
| 259 | 259 |
/// \todo Revise its name |
| 260 | 260 |
template <typename T> |
| 261 | 261 |
class IntegerMap {
|
| 262 | 262 |
|
| 263 | 263 |
template <typename T1> |
| 264 | 264 |
friend class IntegerMap; |
| 265 | 265 |
|
| 266 | 266 |
public: |
| 267 | 267 |
|
| 268 | 268 |
typedef True ReferenceMapTag; |
| 269 | 269 |
///\e |
| 270 | 270 |
typedef int Key; |
| 271 | 271 |
///\e |
| 272 | 272 |
typedef T Value; |
| 273 | 273 |
///\e |
| 274 | 274 |
typedef T& Reference; |
| 275 | 275 |
///\e |
| 276 | 276 |
typedef const T& ConstReference; |
| 277 | 277 |
|
| 278 | 278 |
private: |
| 279 | 279 |
|
| 280 | 280 |
typedef std::vector<T> Vector; |
| 281 | 281 |
Vector _vector; |
| 282 | 282 |
|
| 283 | 283 |
public: |
| 284 | 284 |
|
| 285 | 285 |
/// Constructor with specified default value |
| 286 | 286 |
IntegerMap(int size = 0, const T& value = T()) : _vector(size, value) {}
|
| 287 | 287 |
|
| 288 | 288 |
/// \brief Constructs the map from an appropriate std::vector. |
| 289 | 289 |
template <typename T1> |
| 290 | 290 |
IntegerMap(const std::vector<T1>& vector) |
| 291 | 291 |
: _vector(vector.begin(), vector.end()) {}
|
| 292 | 292 |
|
| 293 | 293 |
/// \brief Constructs a map from an other IntegerMap. |
| 294 | 294 |
template <typename T1> |
| 295 | 295 |
IntegerMap(const IntegerMap<T1> &c) |
| 296 | 296 |
: _vector(c._vector.begin(), c._vector.end()) {}
|
| 297 | 297 |
|
| 298 | 298 |
/// \brief Resize the container |
| 299 | 299 |
void resize(int size, const T& value = T()) {
|
| 300 | 300 |
_vector.resize(size, value); |
| 301 | 301 |
} |
| 302 | 302 |
|
| 303 | 303 |
private: |
| 304 | 304 |
|
| 305 | 305 |
IntegerMap& operator=(const IntegerMap&); |
| 306 | 306 |
|
| 307 | 307 |
public: |
| 308 | 308 |
|
| 309 | 309 |
///\e |
| 310 | 310 |
Reference operator[](Key k) {
|
| 311 | 311 |
return _vector[k]; |
| 312 | 312 |
} |
| 313 | 313 |
|
| 314 | 314 |
/// \e |
| 315 | 315 |
ConstReference operator[](Key k) const {
|
| 316 | 316 |
return _vector[k]; |
| 317 | 317 |
} |
| 318 | 318 |
|
| 319 | 319 |
/// \e |
| 320 | 320 |
void set(const Key &k, const T& t) {
|
| 321 | 321 |
_vector[k] = t; |
| 322 | 322 |
} |
| 323 | 323 |
|
| 324 | 324 |
}; |
| 325 | 325 |
|
| 326 | 326 |
/// @} |
| 327 | 327 |
|
| 328 | 328 |
/// \addtogroup map_adaptors |
| 329 | 329 |
/// @{
|
| 330 | 330 |
|
| 331 | 331 |
/// \brief Identity map. |
| 332 | 332 |
/// |
| 333 | 333 |
/// This map gives back the given key as value without any |
| 334 | 334 |
/// modification. |
| 335 | 335 |
template <typename T> |
| 336 | 336 |
class IdentityMap : public MapBase<T, T> {
|
| 337 | 337 |
public: |
| 338 | 338 |
typedef MapBase<T, T> Parent; |
| 339 | 339 |
typedef typename Parent::Key Key; |
| 340 | 340 |
typedef typename Parent::Value Value; |
| 341 | 341 |
|
| 342 | 342 |
/// \e |
| 343 | 343 |
const T& operator[](const T& t) const {
|
| 344 | 344 |
return t; |
| 345 | 345 |
} |
| 346 | 346 |
}; |
| 347 | 347 |
|
| 348 | 348 |
///Returns an \c IdentityMap class |
| 349 | 349 |
|
| 350 | 350 |
///This function just returns an \c IdentityMap class. |
| 351 | 351 |
///\relates IdentityMap |
| 352 | 352 |
template<typename T> |
| 353 | 353 |
inline IdentityMap<T> identityMap() {
|
| 354 | 354 |
return IdentityMap<T>(); |
| 355 | 355 |
} |
| 356 | 356 |
|
| 357 | 357 |
|
| 358 | 358 |
///\brief Convert the \c Value of a map to another type using |
| 359 | 359 |
///the default conversion. |
| 360 | 360 |
/// |
| 361 | 361 |
///This \c concepts::ReadMap "read only map" |
| 362 | 362 |
///converts the \c Value of a map to type \c T. |
| 363 | 363 |
///Its \c Key is inherited from \c M. |
| 364 | 364 |
template <typename M, typename T> |
| 365 | 365 |
class ConvertMap : public MapBase<typename M::Key, T> {
|
| 366 | 366 |
const M& m; |
| 367 | 367 |
public: |
| 368 | 368 |
typedef MapBase<typename M::Key, T> Parent; |
| 369 | 369 |
typedef typename Parent::Key Key; |
| 370 | 370 |
typedef typename Parent::Value Value; |
| 371 | 371 |
|
| 372 | 372 |
///Constructor |
| 373 | 373 |
|
| 374 | 374 |
///Constructor. |
| 375 | 375 |
///\param _m is the underlying map. |
| 376 | 376 |
ConvertMap(const M &_m) : m(_m) {};
|
| 377 | 377 |
|
| 378 | 378 |
/// \brief The subscript operator. |
| 379 | 379 |
/// |
| 380 | 380 |
/// The subscript operator. |
| 381 | 381 |
Value operator[](const Key& k) const {return m[k];}
|
| 382 | 382 |
}; |
| 383 | 383 |
|
| 384 | 384 |
///Returns a \c ConvertMap class |
| 385 | 385 |
|
| 386 | 386 |
///This function just returns a \c ConvertMap class. |
| 387 | 387 |
///\relates ConvertMap |
| 388 | 388 |
template<typename T, typename M> |
| 389 | 389 |
inline ConvertMap<M, T> convertMap(const M &m) {
|
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
///\file |
| 20 | 20 |
///\brief Instantiation of the Random class. |
| 21 | 21 |
|
| 22 | 22 |
#include <lemon/random.h> |
| 23 | 23 |
|
| 24 | 24 |
namespace lemon {
|
| 25 | 25 |
/// \brief Global random number generator instance |
| 26 | 26 |
/// |
| 27 | 27 |
/// A global Mersenne Twister random number generator instance. |
| 28 | 28 |
Random rnd; |
| 29 | 29 |
} |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
/* |
| 20 | 20 |
* This file contains the reimplemented version of the Mersenne Twister |
| 21 | 21 |
* Generator of Matsumoto and Nishimura. |
| 22 | 22 |
* |
| 23 | 23 |
* See the appropriate copyright notice below. |
| 24 | 24 |
* |
| 25 | 25 |
* Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, |
| 26 | 26 |
* All rights reserved. |
| 27 | 27 |
* |
| 28 | 28 |
* Redistribution and use in source and binary forms, with or without |
| 29 | 29 |
* modification, are permitted provided that the following conditions |
| 30 | 30 |
* are met: |
| 31 | 31 |
* |
| 32 | 32 |
* 1. Redistributions of source code must retain the above copyright |
| 33 | 33 |
* notice, this list of conditions and the following disclaimer. |
| 34 | 34 |
* |
| 35 | 35 |
* 2. Redistributions in binary form must reproduce the above copyright |
| 36 | 36 |
* notice, this list of conditions and the following disclaimer in the |
| 37 | 37 |
* documentation and/or other materials provided with the distribution. |
| 38 | 38 |
* |
| 39 | 39 |
* 3. The names of its contributors may not be used to endorse or promote |
| 40 | 40 |
* products derived from this software without specific prior written |
| 41 | 41 |
* permission. |
| 42 | 42 |
* |
| 43 | 43 |
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS |
| 44 | 44 |
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT |
| 45 | 45 |
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS |
| 46 | 46 |
* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE |
| 47 | 47 |
* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, |
| 48 | 48 |
* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES |
| 49 | 49 |
* (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR |
| 50 | 50 |
* SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) |
| 51 | 51 |
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, |
| 52 | 52 |
* STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) |
| 53 | 53 |
* ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED |
| 54 | 54 |
* OF THE POSSIBILITY OF SUCH DAMAGE. |
| 55 | 55 |
* |
| 56 | 56 |
* |
| 57 | 57 |
* Any feedback is very welcome. |
| 58 | 58 |
* http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html |
| 59 | 59 |
* email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space) |
| 60 | 60 |
*/ |
| 61 | 61 |
|
| 62 | 62 |
#ifndef LEMON_RANDOM_H |
| 63 | 63 |
#define LEMON_RANDOM_H |
| 64 | 64 |
|
| 65 | 65 |
#include <algorithm> |
| 66 | 66 |
#include <iterator> |
| 67 | 67 |
#include <vector> |
| 68 | 68 |
|
| 69 | 69 |
#include <ctime> |
| 70 | 70 |
#include <cmath> |
| 71 | 71 |
|
| 72 | 72 |
#include <lemon/dim2.h> |
| 73 | 73 |
///\ingroup misc |
| 74 | 74 |
///\file |
| 75 | 75 |
///\brief Mersenne Twister random number generator |
| 76 | 76 |
|
| 77 | 77 |
namespace lemon {
|
| 78 | 78 |
|
| 79 | 79 |
namespace _random_bits {
|
| 80 | 80 |
|
| 81 | 81 |
template <typename _Word, int _bits = std::numeric_limits<_Word>::digits> |
| 82 | 82 |
struct RandomTraits {};
|
| 83 | 83 |
|
| 84 | 84 |
template <typename _Word> |
| 85 | 85 |
struct RandomTraits<_Word, 32> {
|
| 86 | 86 |
|
| 87 | 87 |
typedef _Word Word; |
| 88 | 88 |
static const int bits = 32; |
| 89 | 89 |
|
| 90 | 90 |
static const int length = 624; |
| 91 | 91 |
static const int shift = 397; |
| 92 | 92 |
|
| 93 | 93 |
static const Word mul = 0x6c078965u; |
| 94 | 94 |
static const Word arrayInit = 0x012BD6AAu; |
| 95 | 95 |
static const Word arrayMul1 = 0x0019660Du; |
| 96 | 96 |
static const Word arrayMul2 = 0x5D588B65u; |
| 97 | 97 |
|
| 98 | 98 |
static const Word mask = 0x9908B0DFu; |
| 99 | 99 |
static const Word loMask = (1u << 31) - 1; |
| 100 | 100 |
static const Word hiMask = ~loMask; |
| 101 | 101 |
|
| 102 | 102 |
|
| 103 | 103 |
static Word tempering(Word rnd) {
|
| 104 | 104 |
rnd ^= (rnd >> 11); |
| 105 | 105 |
rnd ^= (rnd << 7) & 0x9D2C5680u; |
| 106 | 106 |
rnd ^= (rnd << 15) & 0xEFC60000u; |
| 107 | 107 |
rnd ^= (rnd >> 18); |
| 108 | 108 |
return rnd; |
| 109 | 109 |
} |
| 110 | 110 |
|
| 111 | 111 |
}; |
| 112 | 112 |
|
| 113 | 113 |
template <typename _Word> |
| 114 | 114 |
struct RandomTraits<_Word, 64> {
|
| 115 | 115 |
|
| 116 | 116 |
typedef _Word Word; |
| 117 | 117 |
static const int bits = 64; |
| 118 | 118 |
|
| 119 | 119 |
static const int length = 312; |
| 120 | 120 |
static const int shift = 156; |
| 121 | 121 |
|
| 122 | 122 |
static const Word mul = Word(0x5851F42Du) << 32 | Word(0x4C957F2Du); |
| 123 | 123 |
static const Word arrayInit = Word(0x00000000u) << 32 |Word(0x012BD6AAu); |
| 124 | 124 |
static const Word arrayMul1 = Word(0x369DEA0Fu) << 32 |Word(0x31A53F85u); |
| 125 | 125 |
static const Word arrayMul2 = Word(0x27BB2EE6u) << 32 |Word(0x87B0B0FDu); |
| 126 | 126 |
|
| 127 | 127 |
static const Word mask = Word(0xB5026F5Au) << 32 | Word(0xA96619E9u); |
| 128 | 128 |
static const Word loMask = (Word(1u) << 31) - 1; |
| 129 | 129 |
static const Word hiMask = ~loMask; |
| 130 | 130 |
|
| 131 | 131 |
static Word tempering(Word rnd) {
|
| 132 | 132 |
rnd ^= (rnd >> 29) & (Word(0x55555555u) << 32 | Word(0x55555555u)); |
| 133 | 133 |
rnd ^= (rnd << 17) & (Word(0x71D67FFFu) << 32 | Word(0xEDA60000u)); |
| 134 | 134 |
rnd ^= (rnd << 37) & (Word(0xFFF7EEE0u) << 32 | Word(0x00000000u)); |
| 135 | 135 |
rnd ^= (rnd >> 43); |
| 136 | 136 |
return rnd; |
| 137 | 137 |
} |
| 138 | 138 |
|
| 139 | 139 |
}; |
| 140 | 140 |
|
| 141 | 141 |
template <typename _Word> |
| 142 | 142 |
class RandomCore {
|
| 143 | 143 |
public: |
| 144 | 144 |
|
| 145 | 145 |
typedef _Word Word; |
| 146 | 146 |
|
| 147 | 147 |
private: |
| 148 | 148 |
|
| 149 | 149 |
static const int bits = RandomTraits<Word>::bits; |
| 150 | 150 |
|
| 151 | 151 |
static const int length = RandomTraits<Word>::length; |
| 152 | 152 |
static const int shift = RandomTraits<Word>::shift; |
| 153 | 153 |
|
| 154 | 154 |
public: |
| 155 | 155 |
|
| 156 | 156 |
void initState() {
|
| 157 | 157 |
static const Word seedArray[4] = {
|
| 158 | 158 |
0x12345u, 0x23456u, 0x34567u, 0x45678u |
| 159 | 159 |
}; |
| 160 | 160 |
|
| 161 | 161 |
initState(seedArray, seedArray + 4); |
| 162 | 162 |
} |
| 163 | 163 |
|
| 164 | 164 |
void initState(Word seed) {
|
| 165 | 165 |
|
| 166 | 166 |
static const Word mul = RandomTraits<Word>::mul; |
| 167 | 167 |
|
| 168 | 168 |
current = state; |
| 169 | 169 |
|
| 170 | 170 |
Word *curr = state + length - 1; |
| 171 | 171 |
curr[0] = seed; --curr; |
| 172 | 172 |
for (int i = 1; i < length; ++i) {
|
| 173 | 173 |
curr[0] = (mul * ( curr[1] ^ (curr[1] >> (bits - 2)) ) + i); |
| 174 | 174 |
--curr; |
| 175 | 175 |
} |
| 176 | 176 |
} |
| 177 | 177 |
|
| 178 | 178 |
template <typename Iterator> |
| 179 | 179 |
void initState(Iterator begin, Iterator end) {
|
| 180 | 180 |
|
| 181 | 181 |
static const Word init = RandomTraits<Word>::arrayInit; |
| 182 | 182 |
static const Word mul1 = RandomTraits<Word>::arrayMul1; |
| 183 | 183 |
static const Word mul2 = RandomTraits<Word>::arrayMul2; |
| 184 | 184 |
|
| 185 | 185 |
|
| 186 | 186 |
Word *curr = state + length - 1; --curr; |
| 187 | 187 |
Iterator it = begin; int cnt = 0; |
| 188 | 188 |
int num; |
| 189 | 189 |
|
| 190 | 190 |
initState(init); |
| 191 | 191 |
|
| 192 | 192 |
num = length > end - begin ? length : end - begin; |
| 193 | 193 |
while (num--) {
|
| 194 | 194 |
curr[0] = (curr[0] ^ ((curr[1] ^ (curr[1] >> (bits - 2))) * mul1)) |
| 195 | 195 |
+ *it + cnt; |
| 196 | 196 |
++it; ++cnt; |
| 197 | 197 |
if (it == end) {
|
| 198 | 198 |
it = begin; cnt = 0; |
| 199 | 199 |
} |
| 200 | 200 |
if (curr == state) {
|
| 201 | 201 |
curr = state + length - 1; curr[0] = state[0]; |
| 202 | 202 |
} |
| 203 | 203 |
--curr; |
| 204 | 204 |
} |
| 205 | 205 |
|
| 206 | 206 |
num = length - 1; cnt = length - (curr - state) - 1; |
| 207 | 207 |
while (num--) {
|
| 208 | 208 |
curr[0] = (curr[0] ^ ((curr[1] ^ (curr[1] >> (bits - 2))) * mul2)) |
| 209 | 209 |
- cnt; |
| 210 | 210 |
--curr; ++cnt; |
| 211 | 211 |
if (curr == state) {
|
| 212 | 212 |
curr = state + length - 1; curr[0] = state[0]; --curr; |
| 213 | 213 |
cnt = 1; |
| 214 | 214 |
} |
| 215 | 215 |
} |
| 216 | 216 |
|
| 217 | 217 |
state[length - 1] = Word(1) << (bits - 1); |
| 218 | 218 |
} |
| 219 | 219 |
|
| 220 | 220 |
void copyState(const RandomCore& other) {
|
| 221 | 221 |
std::copy(other.state, other.state + length, state); |
| 222 | 222 |
current = state + (other.current - other.state); |
| 223 | 223 |
} |
| 224 | 224 |
|
| 225 | 225 |
Word operator()() {
|
| 226 | 226 |
if (current == state) fillState(); |
| 227 | 227 |
--current; |
| 228 | 228 |
Word rnd = *current; |
| 229 | 229 |
return RandomTraits<Word>::tempering(rnd); |
| 230 | 230 |
} |
| 231 | 231 |
|
| 232 | 232 |
private: |
| 233 | 233 |
|
| 234 | 234 |
|
| 235 | 235 |
void fillState() {
|
| 236 | 236 |
static const Word mask[2] = { 0x0ul, RandomTraits<Word>::mask };
|
| 237 | 237 |
static const Word loMask = RandomTraits<Word>::loMask; |
| 238 | 238 |
static const Word hiMask = RandomTraits<Word>::hiMask; |
| 239 | 239 |
|
| 240 | 240 |
current = state + length; |
| 241 | 241 |
|
| 242 | 242 |
register Word *curr = state + length - 1; |
| 243 | 243 |
register long num; |
| 244 | 244 |
|
| 245 | 245 |
num = length - shift; |
| 246 | 246 |
while (num--) {
|
| 247 | 247 |
curr[0] = (((curr[0] & hiMask) | (curr[-1] & loMask)) >> 1) ^ |
| 248 | 248 |
curr[- shift] ^ mask[curr[-1] & 1ul]; |
| 249 | 249 |
--curr; |
| 250 | 250 |
} |
| 251 | 251 |
num = shift - 1; |
| 252 | 252 |
while (num--) {
|
| 253 | 253 |
curr[0] = (((curr[0] & hiMask) | (curr[-1] & loMask)) >> 1) ^ |
| 254 | 254 |
curr[length - shift] ^ mask[curr[-1] & 1ul]; |
| 255 | 255 |
--curr; |
| 256 | 256 |
} |
| 257 | 257 |
curr[0] = (((curr[0] & hiMask) | (curr[length - 1] & loMask)) >> 1) ^ |
| 258 | 258 |
curr[length - shift] ^ mask[curr[length - 1] & 1ul]; |
| 259 | 259 |
|
| 260 | 260 |
} |
| 261 | 261 |
|
| 262 | 262 |
|
| 263 | 263 |
Word *current; |
| 264 | 264 |
Word state[length]; |
| 265 | 265 |
|
| 266 | 266 |
}; |
| 267 | 267 |
|
| 268 | 268 |
|
| 269 | 269 |
template <typename Result, |
| 270 | 270 |
int shift = (std::numeric_limits<Result>::digits + 1) / 2> |
| 271 | 271 |
struct Masker {
|
| 272 | 272 |
static Result mask(const Result& result) {
|
| 273 | 273 |
return Masker<Result, (shift + 1) / 2>:: |
| 274 | 274 |
mask(static_cast<Result>(result | (result >> shift))); |
| 275 | 275 |
} |
| 276 | 276 |
}; |
| 277 | 277 |
|
| 278 | 278 |
template <typename Result> |
| 279 | 279 |
struct Masker<Result, 1> {
|
| 280 | 280 |
static Result mask(const Result& result) {
|
| 281 | 281 |
return static_cast<Result>(result | (result >> 1)); |
| 282 | 282 |
} |
| 283 | 283 |
}; |
| 284 | 284 |
|
| 285 | 285 |
template <typename Result, typename Word, |
| 286 | 286 |
int rest = std::numeric_limits<Result>::digits, int shift = 0, |
| 287 | 287 |
bool last = rest <= std::numeric_limits<Word>::digits> |
| 288 | 288 |
struct IntConversion {
|
| 289 | 289 |
static const int bits = std::numeric_limits<Word>::digits; |
| 290 | 290 |
|
| 291 | 291 |
static Result convert(RandomCore<Word>& rnd) {
|
| 292 | 292 |
return static_cast<Result>(rnd() >> (bits - rest)) << shift; |
| 293 | 293 |
} |
| 294 | 294 |
|
| 295 | 295 |
}; |
| 296 | 296 |
|
| 297 | 297 |
template <typename Result, typename Word, int rest, int shift> |
| 298 | 298 |
struct IntConversion<Result, Word, rest, shift, false> {
|
| 299 | 299 |
static const int bits = std::numeric_limits<Word>::digits; |
| 300 | 300 |
|
| 301 | 301 |
static Result convert(RandomCore<Word>& rnd) {
|
| 302 | 302 |
return (static_cast<Result>(rnd()) << shift) | |
| 303 | 303 |
IntConversion<Result, Word, rest - bits, shift + bits>::convert(rnd); |
| 304 | 304 |
} |
| 305 | 305 |
}; |
| 306 | 306 |
|
| 307 | 307 |
|
| 308 | 308 |
template <typename Result, typename Word, |
| 309 | 309 |
bool one_word = (std::numeric_limits<Word>::digits < |
| 310 | 310 |
std::numeric_limits<Result>::digits) > |
| 311 | 311 |
struct Mapping {
|
| 312 | 312 |
static Result map(RandomCore<Word>& rnd, const Result& bound) {
|
| 313 | 313 |
Word max = Word(bound - 1); |
| 314 | 314 |
Result mask = Masker<Result>::mask(bound - 1); |
| 315 | 315 |
Result num; |
| 316 | 316 |
do {
|
| 317 | 317 |
num = IntConversion<Result, Word>::convert(rnd) & mask; |
| 318 | 318 |
} while (num > max); |
| 319 | 319 |
return num; |
| 320 | 320 |
} |
| 321 | 321 |
}; |
| 322 | 322 |
|
| 323 | 323 |
template <typename Result, typename Word> |
| 324 | 324 |
struct Mapping<Result, Word, false> {
|
| 325 | 325 |
static Result map(RandomCore<Word>& rnd, const Result& bound) {
|
| 326 | 326 |
Word max = Word(bound - 1); |
| 327 | 327 |
Word mask = Masker<Word, (std::numeric_limits<Result>::digits + 1) / 2> |
| 328 | 328 |
::mask(max); |
| 329 | 329 |
Word num; |
| 330 | 330 |
do {
|
| 331 | 331 |
num = rnd() & mask; |
| 332 | 332 |
} while (num > max); |
| 333 | 333 |
return num; |
| 334 | 334 |
} |
| 335 | 335 |
}; |
| 336 | 336 |
|
| 337 | 337 |
template <typename Result, int exp, bool pos = (exp >= 0)> |
| 338 | 338 |
struct ShiftMultiplier {
|
| 339 | 339 |
static const Result multiplier() {
|
| 340 | 340 |
Result res = ShiftMultiplier<Result, exp / 2>::multiplier(); |
| 341 | 341 |
res *= res; |
| 342 | 342 |
if ((exp & 1) == 1) res *= static_cast<Result>(2.0); |
| 343 | 343 |
return res; |
| 344 | 344 |
} |
| 345 | 345 |
}; |
| 346 | 346 |
|
| 347 | 347 |
template <typename Result, int exp> |
| 348 | 348 |
struct ShiftMultiplier<Result, exp, false> {
|
| 349 | 349 |
static const Result multiplier() {
|
| 350 | 350 |
Result res = ShiftMultiplier<Result, exp / 2>::multiplier(); |
| 351 | 351 |
res *= res; |
| 352 | 352 |
if ((exp & 1) == 1) res *= static_cast<Result>(0.5); |
| 353 | 353 |
return res; |
| 354 | 354 |
} |
| 355 | 355 |
}; |
| 356 | 356 |
|
| 357 | 357 |
template <typename Result> |
| 358 | 358 |
struct ShiftMultiplier<Result, 0, true> {
|
| 359 | 359 |
static const Result multiplier() {
|
| 360 | 360 |
return static_cast<Result>(1.0); |
| 361 | 361 |
} |
| 362 | 362 |
}; |
| 363 | 363 |
|
| 364 | 364 |
template <typename Result> |
| 365 | 365 |
struct ShiftMultiplier<Result, -20, true> {
|
| 366 | 366 |
static const Result multiplier() {
|
| 367 | 367 |
return static_cast<Result>(1.0/1048576.0); |
| 368 | 368 |
} |
| 369 | 369 |
}; |
| 370 | 370 |
|
| 371 | 371 |
template <typename Result> |
| 372 | 372 |
struct ShiftMultiplier<Result, -32, true> {
|
| 373 | 373 |
static const Result multiplier() {
|
| 374 | 374 |
return static_cast<Result>(1.0/424967296.0); |
| 375 | 375 |
} |
| 376 | 376 |
}; |
| 377 | 377 |
|
| 378 | 378 |
template <typename Result> |
| 379 | 379 |
struct ShiftMultiplier<Result, -53, true> {
|
| 380 | 380 |
static const Result multiplier() {
|
| 381 | 381 |
return static_cast<Result>(1.0/9007199254740992.0); |
| 382 | 382 |
} |
| 383 | 383 |
}; |
| 384 | 384 |
|
| 385 | 385 |
template <typename Result> |
| 386 | 386 |
struct ShiftMultiplier<Result, -64, true> {
|
| 387 | 387 |
static const Result multiplier() {
|
| 388 | 388 |
return static_cast<Result>(1.0/18446744073709551616.0); |
| 389 | 389 |
} |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#ifndef LEMON_TOLERANCE_H |
| 20 | 20 |
#define LEMON_TOLERANCE_H |
| 21 | 21 |
|
| 22 | 22 |
///\ingroup misc |
| 23 | 23 |
///\file |
| 24 | 24 |
///\brief A basic tool to handle the anomalies of calculation with |
| 25 | 25 |
///floating point numbers. |
| 26 | 26 |
/// |
| 27 | 27 |
///\todo It should be in a module like "Basic tools" |
| 28 | 28 |
|
| 29 | 29 |
|
| 30 | 30 |
namespace lemon {
|
| 31 | 31 |
|
| 32 | 32 |
/// \addtogroup misc |
| 33 | 33 |
/// @{
|
| 34 | 34 |
|
| 35 | 35 |
///\brief A class to provide a basic way to |
| 36 | 36 |
///handle the comparison of numbers that are obtained |
| 37 | 37 |
///as a result of a probably inexact computation. |
| 38 | 38 |
/// |
| 39 | 39 |
///Tolerance is a class to provide a basic way to |
| 40 | 40 |
///handle the comparison of numbers that are obtained |
| 41 | 41 |
///as a result of a probably inexact computation. |
| 42 | 42 |
/// |
| 43 | 43 |
///This is an abstract class, it should be specialized for all numerical |
| 44 | 44 |
///data types. These specialized classes like \ref Tolerance\<double\> |
| 45 | 45 |
///may offer additional tuning parameters. |
| 46 | 46 |
/// |
| 47 | 47 |
///\sa Tolerance<float> |
| 48 | 48 |
///\sa Tolerance<double> |
| 49 | 49 |
///\sa Tolerance<long double> |
| 50 | 50 |
///\sa Tolerance<int> |
| 51 | 51 |
#if defined __GNUC__ && !defined __STRICT_ANSI__ |
| 52 | 52 |
///\sa Tolerance<long long int> |
| 53 | 53 |
#endif |
| 54 | 54 |
///\sa Tolerance<unsigned int> |
| 55 | 55 |
#if defined __GNUC__ && !defined __STRICT_ANSI__ |
| 56 | 56 |
///\sa Tolerance<unsigned long long int> |
| 57 | 57 |
#endif |
| 58 | 58 |
|
| 59 | 59 |
template<class T> |
| 60 | 60 |
class Tolerance |
| 61 | 61 |
{
|
| 62 | 62 |
public: |
| 63 | 63 |
typedef T Value; |
| 64 | 64 |
|
| 65 | 65 |
///\name Comparisons |
| 66 | 66 |
///The concept is that these bool functions return with \c true only if |
| 67 | 67 |
///the related comparisons hold even if some numerical error appeared |
| 68 | 68 |
///during the computations. |
| 69 | 69 |
|
| 70 | 70 |
///@{
|
| 71 | 71 |
|
| 72 | 72 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 73 | 73 |
static bool less(Value a,Value b) {return false;}
|
| 74 | 74 |
///Returns \c true if \c a is \e surely different from \c b |
| 75 | 75 |
static bool different(Value a,Value b) {return false;}
|
| 76 | 76 |
///Returns \c true if \c a is \e surely positive |
| 77 | 77 |
static bool positive(Value a) {return false;}
|
| 78 | 78 |
///Returns \c true if \c a is \e surely negative |
| 79 | 79 |
static bool negative(Value a) {return false;}
|
| 80 | 80 |
///Returns \c true if \c a is \e surely non-zero |
| 81 | 81 |
static bool nonZero(Value a) {return false;}
|
| 82 | 82 |
|
| 83 | 83 |
///@} |
| 84 | 84 |
|
| 85 | 85 |
///Returns the zero value. |
| 86 | 86 |
static Value zero() {return T();}
|
| 87 | 87 |
|
| 88 | 88 |
// static bool finite(Value a) {}
|
| 89 | 89 |
// static Value big() {}
|
| 90 | 90 |
// static Value negativeBig() {}
|
| 91 | 91 |
}; |
| 92 | 92 |
|
| 93 | 93 |
|
| 94 | 94 |
///Float specialization of \ref Tolerance. |
| 95 | 95 |
|
| 96 | 96 |
///Float specialization of \ref Tolerance. |
| 97 | 97 |
///\sa Tolerance |
| 98 | 98 |
///\relates Tolerance |
| 99 | 99 |
template<> |
| 100 | 100 |
class Tolerance<float> |
| 101 | 101 |
{
|
| 102 | 102 |
static float def_epsilon; |
| 103 | 103 |
float _epsilon; |
| 104 | 104 |
public: |
| 105 | 105 |
///\e |
| 106 | 106 |
typedef float Value; |
| 107 | 107 |
|
| 108 | 108 |
///Constructor setting the epsilon tolerance to the default value. |
| 109 | 109 |
Tolerance() : _epsilon(def_epsilon) {}
|
| 110 | 110 |
///Constructor setting the epsilon tolerance. |
| 111 | 111 |
Tolerance(float e) : _epsilon(e) {}
|
| 112 | 112 |
|
| 113 | 113 |
///Return the epsilon value. |
| 114 | 114 |
Value epsilon() const {return _epsilon;}
|
| 115 | 115 |
///Set the epsilon value. |
| 116 | 116 |
void epsilon(Value e) {_epsilon=e;}
|
| 117 | 117 |
|
| 118 | 118 |
///Return the default epsilon value. |
| 119 | 119 |
static Value defaultEpsilon() {return def_epsilon;}
|
| 120 | 120 |
///Set the default epsilon value. |
| 121 | 121 |
static void defaultEpsilon(Value e) {def_epsilon=e;}
|
| 122 | 122 |
|
| 123 | 123 |
///\name Comparisons |
| 124 | 124 |
///See class Tolerance for more details. |
| 125 | 125 |
|
| 126 | 126 |
///@{
|
| 127 | 127 |
|
| 128 | 128 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 129 | 129 |
bool less(Value a,Value b) const {return a+_epsilon<b;}
|
| 130 | 130 |
///Returns \c true if \c a is \e surely different from \c b |
| 131 | 131 |
bool different(Value a,Value b) const { return less(a,b)||less(b,a); }
|
| 132 | 132 |
///Returns \c true if \c a is \e surely positive |
| 133 | 133 |
bool positive(Value a) const { return _epsilon<a; }
|
| 134 | 134 |
///Returns \c true if \c a is \e surely negative |
| 135 | 135 |
bool negative(Value a) const { return -_epsilon>a; }
|
| 136 | 136 |
///Returns \c true if \c a is \e surely non-zero |
| 137 | 137 |
bool nonZero(Value a) const { return positive(a)||negative(a); }
|
| 138 | 138 |
|
| 139 | 139 |
///@} |
| 140 | 140 |
|
| 141 | 141 |
///Returns zero |
| 142 | 142 |
static Value zero() {return 0;}
|
| 143 | 143 |
}; |
| 144 | 144 |
|
| 145 | 145 |
///Double specialization of \ref Tolerance. |
| 146 | 146 |
|
| 147 | 147 |
///Double specialization of \ref Tolerance. |
| 148 | 148 |
///\sa Tolerance |
| 149 | 149 |
///\relates Tolerance |
| 150 | 150 |
template<> |
| 151 | 151 |
class Tolerance<double> |
| 152 | 152 |
{
|
| 153 | 153 |
static double def_epsilon; |
| 154 | 154 |
double _epsilon; |
| 155 | 155 |
public: |
| 156 | 156 |
///\e |
| 157 | 157 |
typedef double Value; |
| 158 | 158 |
|
| 159 | 159 |
///Constructor setting the epsilon tolerance to the default value. |
| 160 | 160 |
Tolerance() : _epsilon(def_epsilon) {}
|
| 161 | 161 |
///Constructor setting the epsilon tolerance. |
| 162 | 162 |
Tolerance(double e) : _epsilon(e) {}
|
| 163 | 163 |
|
| 164 | 164 |
///Return the epsilon value. |
| 165 | 165 |
Value epsilon() const {return _epsilon;}
|
| 166 | 166 |
///Set the epsilon value. |
| 167 | 167 |
void epsilon(Value e) {_epsilon=e;}
|
| 168 | 168 |
|
| 169 | 169 |
///Return the default epsilon value. |
| 170 | 170 |
static Value defaultEpsilon() {return def_epsilon;}
|
| 171 | 171 |
///Set the default epsilon value. |
| 172 | 172 |
static void defaultEpsilon(Value e) {def_epsilon=e;}
|
| 173 | 173 |
|
| 174 | 174 |
///\name Comparisons |
| 175 | 175 |
///See class Tolerance for more details. |
| 176 | 176 |
|
| 177 | 177 |
///@{
|
| 178 | 178 |
|
| 179 | 179 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 180 | 180 |
bool less(Value a,Value b) const {return a+_epsilon<b;}
|
| 181 | 181 |
///Returns \c true if \c a is \e surely different from \c b |
| 182 | 182 |
bool different(Value a,Value b) const { return less(a,b)||less(b,a); }
|
| 183 | 183 |
///Returns \c true if \c a is \e surely positive |
| 184 | 184 |
bool positive(Value a) const { return _epsilon<a; }
|
| 185 | 185 |
///Returns \c true if \c a is \e surely negative |
| 186 | 186 |
bool negative(Value a) const { return -_epsilon>a; }
|
| 187 | 187 |
///Returns \c true if \c a is \e surely non-zero |
| 188 | 188 |
bool nonZero(Value a) const { return positive(a)||negative(a); }
|
| 189 | 189 |
|
| 190 | 190 |
///@} |
| 191 | 191 |
|
| 192 | 192 |
///Returns zero |
| 193 | 193 |
static Value zero() {return 0;}
|
| 194 | 194 |
}; |
| 195 | 195 |
|
| 196 | 196 |
///Long double specialization of \ref Tolerance. |
| 197 | 197 |
|
| 198 | 198 |
///Long double specialization of \ref Tolerance. |
| 199 | 199 |
///\sa Tolerance |
| 200 | 200 |
///\relates Tolerance |
| 201 | 201 |
template<> |
| 202 | 202 |
class Tolerance<long double> |
| 203 | 203 |
{
|
| 204 | 204 |
static long double def_epsilon; |
| 205 | 205 |
long double _epsilon; |
| 206 | 206 |
public: |
| 207 | 207 |
///\e |
| 208 | 208 |
typedef long double Value; |
| 209 | 209 |
|
| 210 | 210 |
///Constructor setting the epsilon tolerance to the default value. |
| 211 | 211 |
Tolerance() : _epsilon(def_epsilon) {}
|
| 212 | 212 |
///Constructor setting the epsilon tolerance. |
| 213 | 213 |
Tolerance(long double e) : _epsilon(e) {}
|
| 214 | 214 |
|
| 215 | 215 |
///Return the epsilon value. |
| 216 | 216 |
Value epsilon() const {return _epsilon;}
|
| 217 | 217 |
///Set the epsilon value. |
| 218 | 218 |
void epsilon(Value e) {_epsilon=e;}
|
| 219 | 219 |
|
| 220 | 220 |
///Return the default epsilon value. |
| 221 | 221 |
static Value defaultEpsilon() {return def_epsilon;}
|
| 222 | 222 |
///Set the default epsilon value. |
| 223 | 223 |
static void defaultEpsilon(Value e) {def_epsilon=e;}
|
| 224 | 224 |
|
| 225 | 225 |
///\name Comparisons |
| 226 | 226 |
///See class Tolerance for more details. |
| 227 | 227 |
|
| 228 | 228 |
///@{
|
| 229 | 229 |
|
| 230 | 230 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 231 | 231 |
bool less(Value a,Value b) const {return a+_epsilon<b;}
|
| 232 | 232 |
///Returns \c true if \c a is \e surely different from \c b |
| 233 | 233 |
bool different(Value a,Value b) const { return less(a,b)||less(b,a); }
|
| 234 | 234 |
///Returns \c true if \c a is \e surely positive |
| 235 | 235 |
bool positive(Value a) const { return _epsilon<a; }
|
| 236 | 236 |
///Returns \c true if \c a is \e surely negative |
| 237 | 237 |
bool negative(Value a) const { return -_epsilon>a; }
|
| 238 | 238 |
///Returns \c true if \c a is \e surely non-zero |
| 239 | 239 |
bool nonZero(Value a) const { return positive(a)||negative(a); }
|
| 240 | 240 |
|
| 241 | 241 |
///@} |
| 242 | 242 |
|
| 243 | 243 |
///Returns zero |
| 244 | 244 |
static Value zero() {return 0;}
|
| 245 | 245 |
}; |
| 246 | 246 |
|
| 247 | 247 |
///Integer specialization of \ref Tolerance. |
| 248 | 248 |
|
| 249 | 249 |
///Integer specialization of \ref Tolerance. |
| 250 | 250 |
///\sa Tolerance |
| 251 | 251 |
template<> |
| 252 | 252 |
class Tolerance<int> |
| 253 | 253 |
{
|
| 254 | 254 |
public: |
| 255 | 255 |
///\e |
| 256 | 256 |
typedef int Value; |
| 257 | 257 |
|
| 258 | 258 |
///\name Comparisons |
| 259 | 259 |
///See \ref Tolerance for more details. |
| 260 | 260 |
|
| 261 | 261 |
///@{
|
| 262 | 262 |
|
| 263 | 263 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 264 | 264 |
static bool less(Value a,Value b) { return a<b;}
|
| 265 | 265 |
///Returns \c true if \c a is \e surely different from \c b |
| 266 | 266 |
static bool different(Value a,Value b) { return a!=b; }
|
| 267 | 267 |
///Returns \c true if \c a is \e surely positive |
| 268 | 268 |
static bool positive(Value a) { return 0<a; }
|
| 269 | 269 |
///Returns \c true if \c a is \e surely negative |
| 270 | 270 |
static bool negative(Value a) { return 0>a; }
|
| 271 | 271 |
///Returns \c true if \c a is \e surely non-zero |
| 272 | 272 |
static bool nonZero(Value a) { return a!=0; }
|
| 273 | 273 |
|
| 274 | 274 |
///@} |
| 275 | 275 |
|
| 276 | 276 |
///Returns zero |
| 277 | 277 |
static Value zero() {return 0;}
|
| 278 | 278 |
}; |
| 279 | 279 |
|
| 280 | 280 |
///Unsigned integer specialization of \ref Tolerance. |
| 281 | 281 |
|
| 282 | 282 |
///Unsigned integer specialization of \ref Tolerance. |
| 283 | 283 |
///\sa Tolerance |
| 284 | 284 |
template<> |
| 285 | 285 |
class Tolerance<unsigned int> |
| 286 | 286 |
{
|
| 287 | 287 |
public: |
| 288 | 288 |
///\e |
| 289 | 289 |
typedef unsigned int Value; |
| 290 | 290 |
|
| 291 | 291 |
///\name Comparisons |
| 292 | 292 |
///See \ref Tolerance for more details. |
| 293 | 293 |
|
| 294 | 294 |
///@{
|
| 295 | 295 |
|
| 296 | 296 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 297 | 297 |
static bool less(Value a,Value b) { return a<b;}
|
| 298 | 298 |
///Returns \c true if \c a is \e surely different from \c b |
| 299 | 299 |
static bool different(Value a,Value b) { return a!=b; }
|
| 300 | 300 |
///Returns \c true if \c a is \e surely positive |
| 301 | 301 |
static bool positive(Value a) { return 0<a; }
|
| 302 | 302 |
///Returns \c true if \c a is \e surely negative |
| 303 | 303 |
static bool negative(Value) { return false; }
|
| 304 | 304 |
///Returns \c true if \c a is \e surely non-zero |
| 305 | 305 |
static bool nonZero(Value a) { return a!=0; }
|
| 306 | 306 |
|
| 307 | 307 |
///@} |
| 308 | 308 |
|
| 309 | 309 |
///Returns zero |
| 310 | 310 |
static Value zero() {return 0;}
|
| 311 | 311 |
}; |
| 312 | 312 |
|
| 313 | 313 |
|
| 314 | 314 |
///Long integer specialization of \ref Tolerance. |
| 315 | 315 |
|
| 316 | 316 |
///Long integer specialization of \ref Tolerance. |
| 317 | 317 |
///\sa Tolerance |
| 318 | 318 |
template<> |
| 319 | 319 |
class Tolerance<long int> |
| 320 | 320 |
{
|
| 321 | 321 |
public: |
| 322 | 322 |
///\e |
| 323 | 323 |
typedef long int Value; |
| 324 | 324 |
|
| 325 | 325 |
///\name Comparisons |
| 326 | 326 |
///See \ref Tolerance for more details. |
| 327 | 327 |
|
| 328 | 328 |
///@{
|
| 329 | 329 |
|
| 330 | 330 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 331 | 331 |
static bool less(Value a,Value b) { return a<b;}
|
| 332 | 332 |
///Returns \c true if \c a is \e surely different from \c b |
| 333 | 333 |
static bool different(Value a,Value b) { return a!=b; }
|
| 334 | 334 |
///Returns \c true if \c a is \e surely positive |
| 335 | 335 |
static bool positive(Value a) { return 0<a; }
|
| 336 | 336 |
///Returns \c true if \c a is \e surely negative |
| 337 | 337 |
static bool negative(Value a) { return 0>a; }
|
| 338 | 338 |
///Returns \c true if \c a is \e surely non-zero |
| 339 | 339 |
static bool nonZero(Value a) { return a!=0;}
|
| 340 | 340 |
|
| 341 | 341 |
///@} |
| 342 | 342 |
|
| 343 | 343 |
///Returns zero |
| 344 | 344 |
static Value zero() {return 0;}
|
| 345 | 345 |
}; |
| 346 | 346 |
|
| 347 | 347 |
///Unsigned long integer specialization of \ref Tolerance. |
| 348 | 348 |
|
| 349 | 349 |
///Unsigned long integer specialization of \ref Tolerance. |
| 350 | 350 |
///\sa Tolerance |
| 351 | 351 |
template<> |
| 352 | 352 |
class Tolerance<unsigned long int> |
| 353 | 353 |
{
|
| 354 | 354 |
public: |
| 355 | 355 |
///\e |
| 356 | 356 |
typedef unsigned long int Value; |
| 357 | 357 |
|
| 358 | 358 |
///\name Comparisons |
| 359 | 359 |
///See \ref Tolerance for more details. |
| 360 | 360 |
|
| 361 | 361 |
///@{
|
| 362 | 362 |
|
| 363 | 363 |
///Returns \c true if \c a is \e surely strictly less than \c b |
| 364 | 364 |
static bool less(Value a,Value b) { return a<b;}
|
| 365 | 365 |
///Returns \c true if \c a is \e surely different from \c b |
| 366 | 366 |
static bool different(Value a,Value b) { return a!=b; }
|
| 367 | 367 |
///Returns \c true if \c a is \e surely positive |
| 368 | 368 |
static bool positive(Value a) { return 0<a; }
|
| 369 | 369 |
///Returns \c true if \c a is \e surely negative |
| 370 | 370 |
static bool negative(Value) { return false; }
|
| 371 | 371 |
///Returns \c true if \c a is \e surely non-zero |
| 372 | 372 |
static bool nonZero(Value a) { return a!=0;}
|
| 373 | 373 |
|
| 374 | 374 |
///@} |
| 375 | 375 |
|
| 376 | 376 |
///Returns zero |
| 377 | 377 |
static Value zero() {return 0;}
|
| 378 | 378 |
}; |
| 379 | 379 |
|
| 380 | 380 |
#if defined __GNUC__ && !defined __STRICT_ANSI__ |
| 381 | 381 |
|
| 382 | 382 |
///Long long integer specialization of \ref Tolerance. |
| 383 | 383 |
|
| 384 | 384 |
///Long long integer specialization of \ref Tolerance. |
| 385 | 385 |
///\warning This class (more exactly, type <tt>long long</tt>) |
| 386 | 386 |
///is not ansi compatible. |
| 387 | 387 |
///\sa Tolerance |
| 388 | 388 |
template<> |
| 389 | 389 |
class Tolerance<long long int> |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#include <lemon/dim2.h> |
| 20 | 20 |
#include <iostream> |
| 21 | 21 |
#include "test_tools.h" |
| 22 | 22 |
|
| 23 | 23 |
using namespace std; |
| 24 | 24 |
using namespace lemon; |
| 25 | 25 |
|
| 26 | 26 |
int main() |
| 27 | 27 |
{
|
| 28 | 28 |
cout << "Testing classes 'dim2::Point' and 'dim2::BoundingBox'." << endl; |
| 29 | 29 |
|
| 30 | 30 |
typedef dim2::Point<int> Point; |
| 31 | 31 |
|
| 32 | 32 |
Point p; |
| 33 | 33 |
check(p.size()==2, "Wrong vector initialization."); |
| 34 | 34 |
|
| 35 | 35 |
Point a(1,2); |
| 36 | 36 |
Point b(3,4); |
| 37 | 37 |
check(a[0]==1 && a[1]==2, "Wrong vector initialization."); |
| 38 | 38 |
|
| 39 | 39 |
p = a+b; |
| 40 | 40 |
check(p.x==4 && p.y==6, "Wrong vector addition."); |
| 41 | 41 |
|
| 42 | 42 |
p = a-b; |
| 43 | 43 |
check(p.x==-2 && p.y==-2, "Wrong vector subtraction."); |
| 44 | 44 |
|
| 45 | 45 |
check(a.normSquare()==5,"Wrong vector norm calculation."); |
| 46 | 46 |
check(a*b==11, "Wrong vector scalar product."); |
| 47 | 47 |
|
| 48 | 48 |
int l=2; |
| 49 | 49 |
p = a*l; |
| 50 | 50 |
check(p.x==2 && p.y==4, "Wrong vector multiplication by a scalar."); |
| 51 | 51 |
|
| 52 | 52 |
p = b/l; |
| 53 | 53 |
check(p.x==1 && p.y==2, "Wrong vector division by a scalar."); |
| 54 | 54 |
|
| 55 | 55 |
typedef dim2::BoundingBox<int> BB; |
| 56 | 56 |
BB box1; |
| 57 | 57 |
check(box1.empty(), "It should be empty."); |
| 58 | 58 |
|
| 59 | 59 |
box1.add(a); |
| 60 | 60 |
check(!box1.empty(), "It should not be empty."); |
| 61 | 61 |
box1.add(b); |
| 62 | 62 |
|
| 63 | 63 |
check(box1.bottomLeft().x==1 && |
| 64 | 64 |
box1.bottomLeft().y==2 && |
| 65 | 65 |
box1.topRight().x==3 && |
| 66 | 66 |
box1.topRight().y==4, |
| 67 | 67 |
"Wrong addition of points to box."); |
| 68 | 68 |
|
| 69 | 69 |
p.x=2; p.y=3; |
| 70 | 70 |
check(box1.inside(p), "It should be inside."); |
| 71 | 71 |
|
| 72 | 72 |
p.x=1; p.y=3; |
| 73 | 73 |
check(box1.inside(p), "It should be inside."); |
| 74 | 74 |
|
| 75 | 75 |
p.x=0; p.y=3; |
| 76 | 76 |
check(!box1.inside(p), "It should not be inside."); |
| 77 | 77 |
|
| 78 | 78 |
BB box2(p); |
| 79 | 79 |
check(!box2.empty(), |
| 80 | 80 |
"It should not be empty. Constructed from 1 point."); |
| 81 | 81 |
|
| 82 | 82 |
box2.add(box1); |
| 83 | 83 |
check(box2.inside(p), |
| 84 | 84 |
"It should be inside. Incremented a box with another one."); |
| 85 | 85 |
|
| 86 | 86 |
return 0; |
| 87 | 87 |
} |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#include <deque> |
| 20 | 20 |
#include <set> |
| 21 | 21 |
|
| 22 | 22 |
#include <lemon/concept_check.h> |
| 23 | 23 |
#include <lemon/concepts/maps.h> |
| 24 | 24 |
#include <lemon/maps.h> |
| 25 | 25 |
|
| 26 | 26 |
#include "test_tools.h" |
| 27 | 27 |
|
| 28 | 28 |
using namespace lemon; |
| 29 | 29 |
using namespace lemon::concepts; |
| 30 | 30 |
|
| 31 | 31 |
struct A {};
|
| 32 | 32 |
inline bool operator<(A, A) { return true; }
|
| 33 | 33 |
struct B {};
|
| 34 | 34 |
|
| 35 | 35 |
class F {
|
| 36 | 36 |
public: |
| 37 | 37 |
typedef A argument_type; |
| 38 | 38 |
typedef B result_type; |
| 39 | 39 |
|
| 40 | 40 |
B operator()(const A &) const {return B();}
|
| 41 | 41 |
}; |
| 42 | 42 |
|
| 43 | 43 |
int func(A) {return 3;}
|
| 44 | 44 |
|
| 45 | 45 |
int binc(int, B) {return 4;}
|
| 46 | 46 |
|
| 47 | 47 |
typedef ReadMap<A,double> DoubleMap; |
| 48 | 48 |
typedef ReadWriteMap<A, double> WriteDoubleMap; |
| 49 | 49 |
|
| 50 | 50 |
typedef ReadMap<A,bool> BoolMap; |
| 51 | 51 |
typedef ReadWriteMap<A, bool> BoolWriteMap; |
| 52 | 52 |
|
| 53 | 53 |
int main() |
| 54 | 54 |
{ // checking graph components
|
| 55 | 55 |
|
| 56 | 56 |
checkConcept<ReadMap<A,B>, ReadMap<A,B> >(); |
| 57 | 57 |
checkConcept<WriteMap<A,B>, WriteMap<A,B> >(); |
| 58 | 58 |
checkConcept<ReadWriteMap<A,B>, ReadWriteMap<A,B> >(); |
| 59 | 59 |
checkConcept<ReferenceMap<A,B,B&,const B&>, ReferenceMap<A,B,B&,const B&> >(); |
| 60 | 60 |
|
| 61 | 61 |
checkConcept<ReadMap<A,double>, AddMap<DoubleMap,DoubleMap> >(); |
| 62 | 62 |
checkConcept<ReadMap<A,double>, SubMap<DoubleMap,DoubleMap> >(); |
| 63 | 63 |
checkConcept<ReadMap<A,double>, MulMap<DoubleMap,DoubleMap> >(); |
| 64 | 64 |
checkConcept<ReadMap<A,double>, DivMap<DoubleMap,DoubleMap> >(); |
| 65 | 65 |
checkConcept<ReadMap<A,double>, NegMap<DoubleMap> >(); |
| 66 | 66 |
checkConcept<ReadWriteMap<A,double>, NegWriteMap<WriteDoubleMap> >(); |
| 67 | 67 |
checkConcept<ReadMap<A,double>, AbsMap<DoubleMap> >(); |
| 68 | 68 |
checkConcept<ReadMap<A,double>, ShiftMap<DoubleMap> >(); |
| 69 | 69 |
checkConcept<ReadWriteMap<A,double>, ShiftWriteMap<WriteDoubleMap> >(); |
| 70 | 70 |
checkConcept<ReadMap<A,double>, ScaleMap<DoubleMap> >(); |
| 71 | 71 |
checkConcept<ReadWriteMap<A,double>, ScaleWriteMap<WriteDoubleMap> >(); |
| 72 | 72 |
checkConcept<ReadMap<A,double>, ForkMap<DoubleMap, DoubleMap> >(); |
| 73 | 73 |
checkConcept<ReadWriteMap<A,double>, |
| 74 | 74 |
ForkWriteMap<WriteDoubleMap, WriteDoubleMap> >(); |
| 75 | 75 |
|
| 76 | 76 |
checkConcept<ReadMap<B,double>, ComposeMap<DoubleMap,ReadMap<B,A> > >(); |
| 77 | 77 |
|
| 78 | 78 |
checkConcept<ReadMap<A,B>, FunctorMap<F, A, B> >(); |
| 79 | 79 |
|
| 80 | 80 |
checkConcept<ReadMap<A, bool>, NotMap<BoolMap> >(); |
| 81 | 81 |
checkConcept<ReadWriteMap<A, bool>, NotWriteMap<BoolWriteMap> >(); |
| 82 | 82 |
|
| 83 | 83 |
checkConcept<WriteMap<A, bool>, StoreBoolMap<A*> >(); |
| 84 | 84 |
checkConcept<WriteMap<A, bool>, BackInserterBoolMap<std::deque<A> > >(); |
| 85 | 85 |
checkConcept<WriteMap<A, bool>, FrontInserterBoolMap<std::deque<A> > >(); |
| 86 | 86 |
checkConcept<WriteMap<A, bool>, InserterBoolMap<std::set<A> > >(); |
| 87 | 87 |
checkConcept<WriteMap<A, bool>, FillBoolMap<WriteMap<A, B> > >(); |
| 88 | 88 |
checkConcept<WriteMap<A, bool>, SettingOrderBoolMap<WriteMap<A, int> > >(); |
| 89 | 89 |
|
| 90 | 90 |
int a; |
| 91 | 91 |
|
| 92 | 92 |
a=mapFunctor(constMap<A,int>(2))(A()); |
| 93 | 93 |
check(a==2,"Something is wrong with mapFunctor"); |
| 94 | 94 |
|
| 95 | 95 |
B b; |
| 96 | 96 |
b=functorMap(F())[A()]; |
| 97 | 97 |
|
| 98 | 98 |
a=functorMap(&func)[A()]; |
| 99 | 99 |
check(a==3,"Something is wrong with functorMap"); |
| 100 | 100 |
|
| 101 | 101 |
a=combineMap(constMap<B, int, 1>(), identityMap<B>(), &binc)[B()]; |
| 102 | 102 |
check(a==4,"Something is wrong with combineMap"); |
| 103 | 103 |
|
| 104 | 104 |
|
| 105 | 105 |
std::cout << __FILE__ ": All tests passed.\n"; |
| 106 | 106 |
|
| 107 | 107 |
return 0; |
| 108 | 108 |
} |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#include <lemon/random.h> |
| 20 | 20 |
#include "test_tools.h" |
| 21 | 21 |
|
| 22 | 22 |
///\file \brief Test cases for random.h |
| 23 | 23 |
/// |
| 24 | 24 |
///\todo To be extended |
| 25 | 25 |
/// |
| 26 | 26 |
|
| 27 | 27 |
int main() |
| 28 | 28 |
{
|
| 29 | 29 |
double a=lemon::rnd(); |
| 30 | 30 |
check(a<1.0&&a>0.0,"This should be in [0,1)"); |
| 31 | 31 |
a=lemon::rnd.gauss(); |
| 32 | 32 |
a=lemon::rnd.gamma(3.45,0); |
| 33 | 33 |
a=lemon::rnd.gamma(4); |
| 34 | 34 |
//Does gamma work with integer k? |
| 35 | 35 |
a=lemon::rnd.gamma(4.0,0); |
| 36 | 36 |
} |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#ifndef LEMON_TEST_TEST_TOOLS_H |
| 20 | 20 |
#define LEMON_TEST_TEST_TOOLS_H |
| 21 | 21 |
|
| 22 | 22 |
#include <iostream> |
| 23 | 23 |
|
| 24 | 24 |
//! \ingroup misc |
| 25 | 25 |
//! \file |
| 26 | 26 |
//! \brief Some utilities to write test programs. |
| 27 | 27 |
|
| 28 | 28 |
|
| 29 | 29 |
///If \c rc is fail, writes an error message end exit. |
| 30 | 30 |
|
| 31 | 31 |
///If \c rc is fail, writes an error message end exit. |
| 32 | 32 |
///The error message contains the file name and the line number of the |
| 33 | 33 |
///source code in a standard from, which makes it possible to go there |
| 34 | 34 |
///using good source browsers like e.g. \c emacs. |
| 35 | 35 |
/// |
| 36 | 36 |
///For example |
| 37 | 37 |
///\code check(0==1,"This is obviously false.");\endcode will |
| 38 | 38 |
///print this (and then exits). |
| 39 | 39 |
///\verbatim graph_test.cc:123: error: This is obviously false. \endverbatim |
| 40 | 40 |
/// |
| 41 | 41 |
///\todo It should be in \c error.h |
| 42 | 42 |
#define check(rc, msg) \ |
| 43 | 43 |
if(!(rc)) { \
|
| 44 | 44 |
std::cerr << __FILE__ ":" << __LINE__ << ": error: " << msg << std::endl; \ |
| 45 | 45 |
abort(); \ |
| 46 | 46 |
} else { } \
|
| 47 | 47 |
|
| 48 | 48 |
#endif |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#include "test_tools.h" |
| 20 | 20 |
|
| 21 | 21 |
int main() |
| 22 | 22 |
{
|
| 23 | 23 |
check(false, "Don't panic. Failing is the right behaviour here."); |
| 24 | 24 |
return 0; |
| 25 | 25 |
} |
| 1 | 1 |
/* -*- C++ -*- |
| 2 | 2 |
* |
| 3 | 3 |
* This file is a part of LEMON, a generic C++ optimization library |
| 4 | 4 |
* |
| 5 |
* Copyright (C) 2003- |
|
| 5 |
* Copyright (C) 2003-2008 |
|
| 6 | 6 |
* Egervary Jeno Kombinatorikus Optimalizalasi Kutatocsoport |
| 7 | 7 |
* (Egervary Research Group on Combinatorial Optimization, EGRES). |
| 8 | 8 |
* |
| 9 | 9 |
* Permission to use, modify and distribute this software is granted |
| 10 | 10 |
* provided that this copyright notice appears in all copies. For |
| 11 | 11 |
* precise terms see the accompanying LICENSE file. |
| 12 | 12 |
* |
| 13 | 13 |
* This software is provided "AS IS" with no warranty of any kind, |
| 14 | 14 |
* express or implied, and with no claim as to its suitability for any |
| 15 | 15 |
* purpose. |
| 16 | 16 |
* |
| 17 | 17 |
*/ |
| 18 | 18 |
|
| 19 | 19 |
#include "test_tools.h" |
| 20 | 20 |
|
| 21 | 21 |
int main() |
| 22 | 22 |
{
|
| 23 | 23 |
check(true, "It should pass."); |
| 24 | 24 |
return 0; |
| 25 | 25 |
} |
0 comments (0 inline)