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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
// Copyright 2015-2018 Hans Dembinski
//
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt
// or copy at http://www.boost.org/LICENSE_1_0.txt)
#include <boost/core/lightweight_test.hpp>
#include <boost/histogram/accumulators/mean.hpp>
#include <boost/histogram/accumulators/ostream.hpp>
#include <boost/histogram/detail/accumulator_traits.hpp>
#include <boost/histogram/weight.hpp>
#include <sstream>
#include <tuple>
#include <type_traits>
#include "is_close.hpp"
#include "str.hpp"
#include "throw_exception.hpp"
using namespace boost::histogram;
using namespace std::literals;
int main() {
using m_t = accumulators::mean<double>;
using traits = detail::accumulator_traits<m_t>;
static_assert(traits::weight_support, "");
static_assert(std::is_same<traits::args, std::tuple<const double&>>::value, "");
// basic interface, string conversion
{
m_t a;
BOOST_TEST_EQ(a.count(), 0);
BOOST_TEST_EQ(a, m_t{});
a(4);
a(7);
a(13);
a(16);
BOOST_TEST_EQ(a.count(), 4);
BOOST_TEST_EQ(a.value(), 10);
BOOST_TEST_EQ(a.variance(), 30);
BOOST_TEST_EQ(str(a), "mean(4, 10, 30)"s);
BOOST_TEST_EQ(str(a, 20, false), " mean(4, 10, 30)"s);
BOOST_TEST_EQ(str(a, 20, true), "mean(4, 10, 30) "s);
}
// small variation on large number
{
m_t a;
a(1e8 + 4);
a(1e8 + 7);
a(1e8 + 13);
a(1e8 + 16);
BOOST_TEST_EQ(a.count(), 4);
BOOST_TEST_EQ(a.value(), 1e8 + 10);
BOOST_TEST_EQ(a.variance(), 30);
}
// addition of zero element
{
BOOST_TEST_EQ(m_t() += m_t(), m_t());
BOOST_TEST_EQ(m_t(1, 2, 3) += m_t(), m_t(1, 2, 3));
BOOST_TEST_EQ(m_t() += m_t(1, 2, 3), m_t(1, 2, 3));
}
// addition
{
m_t a, b, c;
a(1);
a(2);
a(3);
BOOST_TEST_EQ(a.count(), 3);
BOOST_TEST_EQ(a.value(), 2);
BOOST_TEST_EQ(a.variance(), 1);
b(4);
b(6);
BOOST_TEST_EQ(b.count(), 2);
BOOST_TEST_EQ(b.value(), 5);
BOOST_TEST_EQ(b.variance(), 2);
c(1);
c(2);
c(3);
c(4);
c(6);
auto d = a;
d += b;
BOOST_TEST_EQ(c.count(), d.count());
BOOST_TEST_EQ(c.value(), d.value());
BOOST_TEST_IS_CLOSE(c.variance(), d.variance(), 1e-3);
}
// using weight=2 must be same as adding all samples twice
{
m_t a, b;
for (int i = 0; i < 2; ++i) {
a(4);
a(7);
a(13);
a(16);
}
b(weight(2), 4);
b(weight(2), 7);
b(weight(2), 13);
b(weight(2), 16);
BOOST_TEST_EQ(a.count(), b.count());
BOOST_TEST_EQ(a.value(), b.value());
BOOST_TEST_IS_CLOSE(a.variance(), b.variance(), 1e-3);
}
return boost::report_errors();
}
|
