histogram/test/histogram_mixed_test.cpp
Jay Gohil 10c19df918
Add fraction accumulator (for calculating efficiencies, etc) (#361)
* Added fraction accumulator
* Added binomial proportion interval computers: Wald, Wilson Score, Clopper Pearson, Jeffreys
* Width of intervals for all classes can be set via deviation or confidence_level
* Support valarray in histogram::fill
* Add fraction and count to user guide

Co-authored-by: Hans Dembinski <hans.dembinski@gmail.com>
Co-authored-by: Hans Dembinski <HDembinski@users.noreply.github.com>
2022-09-28 11:18:46 +02:00

92 lines
2.8 KiB
C++

// Copyright 2015-2017 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/axis/integer.hpp>
#include <boost/histogram/axis/ostream.hpp>
#include <boost/histogram/axis/regular.hpp>
#include <boost/histogram/histogram.hpp>
#include <boost/histogram/literals.hpp>
#include <boost/histogram/ostream.hpp>
#include <boost/histogram/storage_adaptor.hpp>
#include <boost/histogram/unlimited_storage.hpp>
#include <vector>
#include "histogram.hpp"
#include "throw_exception.hpp"
using namespace boost::histogram;
template <typename T1, typename T2>
void run_tests() {
// compare
{
auto a = make(T1{}, axis::regular<>{3, 0, 3}, axis::integer<>(0, 2));
auto b = make_s(T2{}, std::vector<unsigned>(), axis::regular<>{3, 0, 3},
axis::integer<>(0, 2));
BOOST_TEST_EQ(a, b);
auto b2 = make(T2{}, axis::integer<>{0, 3}, axis::integer<>(0, 2));
BOOST_TEST_NE(a, b2);
auto b3 = make(T2{}, axis::regular<>(3, 0, 4), axis::integer<>(0, 2));
BOOST_TEST_NE(a, b3);
}
// operators
{
auto a = make(T1{}, axis::integer<int, use_default, axis::option::none_t>{0, 2});
auto b = make(T2{}, axis::integer<int, use_default, axis::option::none_t>{0, 2});
BOOST_TEST_EQ(a, b);
a(0); // 1 0
b(1); // 0 1
a += b;
BOOST_TEST_EQ(a[0], 1);
BOOST_TEST_EQ(a[1], 1);
a *= b;
BOOST_TEST_EQ(a[0], 0);
BOOST_TEST_EQ(a[1], 1);
a -= b;
BOOST_TEST_EQ(a[0], 0);
BOOST_TEST_EQ(a[1], 0);
a[0] = 2;
a[1] = 4;
b[0] = 2;
b[1] = 2;
a /= b;
BOOST_TEST_EQ(a[0], 1);
BOOST_TEST_EQ(a[1], 2);
BOOST_TEST_THROWS(a += make(T2{}, axis::integer<>{0, 3}), std::invalid_argument);
BOOST_TEST_THROWS(a -= make(T2{}, axis::integer<>{0, 3}), std::invalid_argument);
BOOST_TEST_THROWS(a *= make(T2{}, axis::integer<>{0, 3}), std::invalid_argument);
BOOST_TEST_THROWS(a /= make(T2{}, axis::integer<>{0, 3}), std::invalid_argument);
}
// copy_assign
{
auto a = make(T1{}, axis::regular<>{3, 0, 3}, axis::integer<>{0, 2});
auto b = make_s(T2{}, std::vector<double>(), axis::regular<>{3, 0, 3},
axis::integer<>{0, 2});
a(1, 1);
BOOST_TEST_NE(a, b);
b = a;
BOOST_TEST_EQ(a, b);
}
}
int main() {
run_tests<static_tag, dynamic_tag>();
run_tests<dynamic_tag, static_tag>();
// copy assign
{
auto a = make(static_tag{}, axis::regular<>{3, 0, 3}, axis::integer<>{0, 2});
auto b = make(dynamic_tag{}, axis::regular<>{3, 0, 3}, axis::regular<>{2, 0, 2},
axis::integer<>{0, 2});
BOOST_TEST_THROWS(a = b, std::invalid_argument);
}
return boost::report_errors();
}