math.hpp 3.15 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36


//| This file is a part of the sferes2 framework.
//| Copyright 2009, ISIR / Universite Pierre et Marie Curie (UPMC)
//| Main contributor(s): Jean-Baptiste Mouret, mouret@isir.fr
//|
//| This software is a computer program whose purpose is to facilitate
//| experiments in evolutionary computation and evolutionary robotics.
//|
//| This software is governed by the CeCILL license under French law
//| and abiding by the rules of distribution of free software.  You
//| can use, modify and/ or redistribute the software under the terms
//| of the CeCILL license as circulated by CEA, CNRS and INRIA at the
//| following URL "http://www.cecill.info".
//|
//| As a counterpart to the access to the source code and rights to
//| copy, modify and redistribute granted by the license, users are
//| provided only with a limited warranty and the software's author,
//| the holder of the economic rights, and the successive licensors
//| have only limited liability.
//|
//| In this respect, the user's attention is drawn to the risks
//| associated with loading, using, modifying and/or developing or
//| reproducing the software by the user in light of its specific
//| status of free software, that may mean that it is complicated to
//| manipulate, and that also therefore means that it is reserved for
//| developers and experienced professionals having in-depth computer
//| knowledge. Users are therefore encouraged to load and test the
//| software's suitability as regards their requirements in conditions
//| enabling the security of their systems and/or data to be ensured
//| and, more generally, to use and operate it in the same conditions
//| as regards security.
//|
//| The fact that you are presently reading this means that you have
//| had knowledge of the CeCILL license and that you accept its terms.

37
38
#ifndef LIMBO_TOOLS_MATH_HPP
#define LIMBO_TOOLS_MATH_HPP
39
40
41

#include <cstdlib>
#include <cmath>
42
#include <ctime>
43
44
45
#include <list>
#include <stdlib.h>
#include <random>
46
#include <utility>
47

48
namespace limbo {
49
    namespace tools {
50
51
52
53
54
55
56
        // usage :
        // auto rgen = make_rgen();
        // double r = rand_double(rgen);
        using rand_gen_t = std::mt19937;
        using uniform_dist_t = std::uniform_real_distribution<double>;
        using rdist_t = std::pair<rand_gen_t, uniform_dist_t>;
        inline rdist_t make_rgen(double min, double max)
57
        {
58
59
          std::random_device rd;
          return std::make_pair(std::mt19937(rd()), std::uniform_real_distribution<double>(min, max));
60
        }
61
        inline double rand_double(rdist_t& rgen)
62
        {
63
64
65
          rand_gen_t& gen = std::get<0>(rgen);
          uniform_dist_t& dist = std::get<1>(rgen);
          return dist(gen);
66
        }
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85

        // get sign of number
        template <typename T>
        inline constexpr int signum(T x, std::false_type is_signed)
        {
            return T(0) < x;
        }

        template <typename T>
        inline constexpr int signum(T x, std::true_type is_signed)
        {
            return (T(0) < x) - (x < T(0));
        }

        template <typename T>
        inline constexpr int signum(T x)
        {
            return signum(x, std::is_signed<T>());
        }
86
87
    }
}
88

89
#endif