Core/src/cpp/base.cpp at 37380b751202c28fc15e409e948a63a4a86e387f · StatisKit/Core · GitHub

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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
#include <string>
#include <iostream>
#include <limits>
#include <memory>

#include "base.h"
#include "sample_space.h"

namespace statiskit
{
    namespace __impl
    {
        double reldiff(const double& prev, const double& curr)
        {
            return fabs(curr-prev) / fabs(prev);
        }

        double reldiff(const Eigen::VectorXd& prev, const Eigen::VectorXd& curr)
        {
            double norm = 0;
            for (Index i = 0, max_index = prev.size(); i < max_index; ++i) {
                if (prev[i] == 0) {
                    norm +=  pow(curr[i], 2);
                } else {
                    norm +=  pow((prev[i]-curr[i])/prev[i], 2);
                }
            }
            return sqrt(norm);
        }

        double reldiff(const Eigen::MatrixXd& prev, const Eigen::MatrixXd& curr)
        {
            double norm = 0;
            for (Index i = 0, max_index = prev.rows(); i < max_index; ++i) {
                for (Index j = 0, max_index = prev.cols(); j < max_index; ++j) {
                    if (prev(i,j) == 0) {
                        norm +=  pow(curr(i,j), 2);
                    } else {
                        norm +=  pow((prev(i,j)-curr(i,j))/prev(i,j), 2);
                    }
                }
            }
            return sqrt(norm);
        }

        boost::mt19937 _random_generator = boost::mt19937(0);

        boost::mt19937& get_random_generator()
        {
            return _random_generator;
        }
    }

    void set_seed()
    {
        __impl::_random_generator.seed();
    }

    void set_seed(const Index& seed)
    {
        __impl::_random_generator.seed(seed);
    }

    not_implemented_error::not_implemented_error(const std::string& function, const std::string& file, const unsigned int& line) : std::runtime_error("'" + function + "' in file '" + file + "' at line " + __impl::to_string(line) + " is not implemented")
    {}

    parameter_error::parameter_error(const std::string& parameter, const std::string& error) : std::runtime_error("'" + parameter + "' parameter: " + error)
    {}

    member_error::member_error(const std::string& member, const std::string& error) : std::runtime_error("'" + member + "' member: " + error)
    {}

    size_error::size_error(const std::string& parameter, const Index& self, const Index& other, const size_type& size) : parameter_error(parameter, "size of " + __impl::to_string(self) + " instead of " + __impl::to_string(other))
    {}

    nullptr_error::nullptr_error(const std::string& parameter) : parameter_error(parameter, "cannot be set to nullptr")
    {}
}