/usr/include/boost/histogram/detail
Edit: /usr/include/boost/histogram/detail/index_translator.hpp (3437B)
// Copyright 2019 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)
#ifndef BOOST_HISTOGRAM_DETAIL_INDEX_TRANSLATOR_HPP
#define BOOST_HISTOGRAM_DETAIL_INDEX_TRANSLATOR_HPP
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
namespace boost {
namespace histogram {
namespace detail {
template
struct index_translator {
using index_type = axis::index_type;
using multi_index_type = multi_index::value>;
using cref = const A&;
cref dst, src;
bool pass_through[buffer_size::value];
index_translator(cref d, cref s) : dst{d}, src{s} { init(dst, src); }
template
void init(const T& a, const T& b) {
std::transform(a.begin(), a.end(), b.begin(), pass_through,
[](const auto& a, const auto& b) {
return axis::visit(
[&](const auto& a) {
using U = std::decay_t;
return relaxed_equal{}(a, axis::get(b));
},
a);
});
}
template
void init(const std::tuple& a, const std::tuple& b) {
using Seq = mp11::mp_iota_c;
mp11::mp_for_each([&](auto I) {
pass_through[I] = relaxed_equal{}(std::get(a), std::get(b));
});
}
template
static index_type translate(const T& dst, const T& src, index_type i) noexcept {
assert(axis::traits::is_continuous::value == false); // LCOV_EXCL_LINE: unreachable
return dst.index(src.value(i));
}
template
void impl(const std::tuple& a, const std::tuple& b, It i,
index_type* j) const noexcept {
using Seq = mp11::mp_iota_c;
mp11::mp_for_each([&](auto I) {
if (pass_through[I])
*(j + I) = *(i + I);
else
*(j + I) = this->translate(std::get(a), std::get(b), *(i + I));
});
}
template
void impl(const T& a, const T& b, It i, index_type* j) const noexcept {
const bool* p = pass_through;
for (unsigned k = 0; k < a.size(); ++k, ++i, ++j, ++p) {
if (*p)
*j = *i;
else {
const auto& bk = b[k];
axis::visit(
[&](const auto& ak) {
using U = std::decay_t;
*j = this->translate(ak, axis::get(bk), *i);
},
a[k]);
}
}
}
template
auto operator()(const Indices& seq) const noexcept {
auto mi = multi_index_type::create(seq.size());
impl(dst, src, seq.begin(), mi.begin());
return mi;
}
};
template
auto make_index_translator(const Axes& dst, const Axes& src) noexcept {
return index_translator{dst, src};
}
} // namespace detail
} // namespace histogram
} // namespace boost
#endif