histogram/test/storage_adaptor_test.cpp

// 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 <array>
#include <boost/core/lightweight_test.hpp>
#include <boost/histogram/accumulators/mean.hpp>
#include <boost/histogram/accumulators/weight.hpp>
#include <boost/histogram/adaptive_storage.hpp>
#include <boost/histogram/histogram_fwd.hpp>
#include <boost/histogram/storage_adaptor.hpp>
#include <deque>
#include <limits>
#include <vector>
#include "is_close.hpp"
#include "utility_allocator.hpp"

using namespace boost::histogram;

template <typename T>
void tests() {
  // ctor, copy, move
  {
    storage_adaptor<T> a;
    a.reset(2);
    storage_adaptor<T> b(a);
    storage_adaptor<T> c;
    c = a;
    BOOST_TEST_EQ(a.size(), 2);
    BOOST_TEST_EQ(b.size(), 2);
    BOOST_TEST_EQ(c.size(), 2);

    storage_adaptor<T> d(std::move(a));
    BOOST_TEST_EQ(d.size(), 2);
    storage_adaptor<T> e;
    e = std::move(d);
    BOOST_TEST_EQ(e.size(), 2);

    T t;
    storage_adaptor<T> g(t); // tests converting ctor
  }

  // increment, add, set, reset
  {
    storage_adaptor<T> a;
    a.reset(1);
    a(0);
    a(0);
    BOOST_TEST_EQ(a[0], 2);
    a.reset(2);
    BOOST_TEST_EQ(a.size(), 2);
    a(0);
    a.add(0, 2);
    a.add(1, 5);
    BOOST_TEST_EQ(a[0], 3);
    BOOST_TEST_EQ(a[1], 5);
    a[1] = 9;
    BOOST_TEST_EQ(a[0], 3);
    BOOST_TEST_EQ(a[1], 9);
    a.reset(0);
    BOOST_TEST_EQ(a.size(), 0);
  }

  // multiply
  {
    storage_adaptor<T> a;
    a.reset(2);
    a(0);
    a *= 3;
    BOOST_TEST_EQ(a[0], 3);
    BOOST_TEST_EQ(a[1], 0);
    a.add(1, 2);
    BOOST_TEST_EQ(a[0], 3);
    BOOST_TEST_EQ(a[1], 2);
    a *= 3;
    BOOST_TEST_EQ(a[0], 9);
    BOOST_TEST_EQ(a[1], 6);
  }

  // copy
  {
    storage_adaptor<T> a;
    a.reset(1);
    a(0);
    decltype(a) b;
    b.reset(2);
    BOOST_TEST(!(a == b));
    b = a;
    BOOST_TEST(a == b);
    BOOST_TEST_EQ(b.size(), 1);
    BOOST_TEST_EQ(b[0], 1);

    decltype(a) c(a);
    BOOST_TEST(a == c);
    BOOST_TEST_EQ(c.size(), 1);
    BOOST_TEST_EQ(c[0], 1);
  }

  // move
  {
    storage_adaptor<T> a;
    a.reset(1);
    a(0);
    decltype(a) b;
    BOOST_TEST(!(a == b));
    b = std::move(a);
    BOOST_TEST_EQ(b.size(), 1);
    BOOST_TEST_EQ(b[0], 1);
    decltype(a) c(std::move(b));
    BOOST_TEST_EQ(c.size(), 1);
    BOOST_TEST_EQ(c[0], 1);
  }

  // add
  {
    storage_adaptor<T> a;
    a.reset(2);
    a(1);
    auto b = a;
    b += a;
    BOOST_TEST_EQ(b[0], 0);
    BOOST_TEST_EQ(b[1], 2);
    a += a;
    // also test self-add
    BOOST_TEST_EQ(a[0], 0);
    BOOST_TEST_EQ(a[1], 2);
  }

  // multiply
  {
    storage_adaptor<T> a;
    a.reset(2);
    a(1);
    a *= 2;
    BOOST_TEST_EQ(a[0], 0);
    BOOST_TEST_EQ(a[1], 2);
  }
}

template <typename A, typename B>
void mixed_tests() {
  // comparison
  {
    A a, b;
    a.reset(1);
    b.reset(1);
    B c, d;
    c.reset(1);
    d.reset(2);
    a(0);
    b(0);
    c(0);
    c(0);
    d(0);
    d.add(1, 5);
    d.add(0, 2);
    BOOST_TEST_EQ(a[0], 1);
    BOOST_TEST_EQ(b[0], 1);
    BOOST_TEST_EQ(c[0], 2);
    BOOST_TEST_EQ(d[0], 3);
    BOOST_TEST_EQ(d[1], 5);
    BOOST_TEST(a == a);
    BOOST_TEST(a == b);
    BOOST_TEST(!(a == c));
    BOOST_TEST(!(a == d));
  }

  // ctor, copy, move
  {
    A a;
    a.reset(2);
    a(1);
    B b(a);
    B c;
    c = a;
    BOOST_TEST_EQ(c[0], 0);
    BOOST_TEST_EQ(c[1], 1);
    B d(std::move(a));
    B e;
    e = std::move(d);
    BOOST_TEST_EQ(e[0], 0);
    BOOST_TEST_EQ(e[1], 1);
  }
}

int main() {
  tests<std::vector<unsigned>>();
  tests<std::array<unsigned, 100>>();
  tests<std::deque<int>>();
  tests<std::map<std::size_t, unsigned>>();
  tests<std::unordered_map<std::size_t, unsigned>>();

  mixed_tests<storage_adaptor<std::vector<unsigned>>,
              storage_adaptor<std::array<double, 100>>>();
  mixed_tests<adaptive_storage<>, storage_adaptor<std::vector<unsigned>>>();
  mixed_tests<storage_adaptor<std::vector<unsigned>>, adaptive_storage<>>();

  // with accumulators::weight
  {
    auto a = storage_adaptor<std::vector<accumulators::weight<double>>>();
    a.reset(1);
    a(0);
    a.add(0, 1);
    a.add(0, accumulators::weight<double>(1, 0));
    BOOST_TEST_EQ(a[0].value(), 3);
    BOOST_TEST_EQ(a[0].variance(), 2);
    a(0, weight(2));
    BOOST_TEST_EQ(a[0].value(), 5);
    BOOST_TEST_EQ(a[0].variance(), 6);
  }

  // with accumulators::mean
  {
    auto a = storage_adaptor<std::vector<accumulators::mean<double>>>();
    a.reset(1);
    a(0, 1);
    a(0, weight(2), 2);
    a.add(0, accumulators::mean<double>(1, 0, 0));
    BOOST_TEST_EQ(a[0].sum(), 4);
    BOOST_TEST_IS_CLOSE(a[0].value(), 1.25, 1e-3);
    BOOST_TEST_IS_CLOSE(a[0].variance(), 0.916, 1e-3);
  }

  // exceeding array capacity
  {
    auto a = storage_adaptor<std::array<int, 10>>();
    a.reset(10); // should not throw
    BOOST_TEST_THROWS(a.reset(11), std::runtime_error);
  }

  // test sparsity of map backend
  {
    tracing_allocator_db db;
    tracing_allocator<char> alloc(db);
    using map = std::map<std::size_t, double, std::less<std::size_t>,
                         tracing_allocator<std::pair<const std::size_t, double>>>;
    using A = storage_adaptor<map>;
    auto a = A(map(alloc));
    // MSVC implementation allocates some structures for debugging
    const auto baseline = db.sum.first;
    a.reset(10);
    BOOST_TEST_EQ(db.sum.first, baseline); // nothing allocated yet
    // query on const reference does not allocate
    BOOST_TEST_EQ(static_cast<const A&>(a)[0], 0);
    // query on writeable reference allocates
    BOOST_TEST_EQ(a[9], 0);
    BOOST_TEST_EQ(db.sum.first, baseline + 1);

    a(5);
    BOOST_TEST_EQ(a[0], 0);
    BOOST_TEST_EQ(a[5], 1);
    BOOST_TEST_EQ(a[9], 0);
    BOOST_TEST_EQ(db.sum.first, baseline + 3);
    a *= 2;
    BOOST_TEST_EQ(a[5], 2);
    BOOST_TEST_EQ(db.sum.first, baseline + 3);

    auto b = storage_adaptor<std::vector<int>>();
    b.reset(5);
    b(2);
    a = b;
    // only one new allocation for non-zero value
    BOOST_TEST_EQ(db.sum.first, baseline + 4);
  }

  return boost::report_errors();
}