diff --git a/include/evalhyd/detail/probabilist/brier.hpp b/include/evalhyd/detail/probabilist/brier.hpp index 013395e24d8fb9a42f795e1a8c2a0a63e2447f04..5e7e73db9f65b635c57b906d1280f94e7baa9d9b 100644 --- a/include/evalhyd/detail/probabilist/brier.hpp +++ b/include/evalhyd/detail/probabilist/brier.hpp @@ -508,7 +508,8 @@ namespace evalhyd for (std::size_t e = 0; e < n_exp; e++) { // retrieve length of period - auto l = xt::view(t_counts, xt::all(), xt::all(), m, e); + auto l = xt::view(t_counts, xt::all(), xt::all(), + m, xt::newaxis(), e); // retrieve range of forecast values $y_i$ auto y_i = xt::eval( @@ -670,7 +671,8 @@ namespace evalhyd xt::all(), b_exp[e]); // retrieve length of period - auto l = xt::view(t_counts, xt::all(), xt::all(), m, e); + auto l = xt::view(t_counts, xt::all(), xt::all(), + m, xt::newaxis(), e); // compute mask to subsample time steps belonging to same observation bin // (where bins are defined as the range of forecast values) diff --git a/include/evalhyd/detail/probabilist/evaluator.hpp b/include/evalhyd/detail/probabilist/evaluator.hpp index c3dc9669eb758d20e7a7f5cce341f85cdcab4e02..88eb1361203384311048e7534e5b1d369d4a67c9 100644 --- a/include/evalhyd/detail/probabilist/evaluator.hpp +++ b/include/evalhyd/detail/probabilist/evaluator.hpp @@ -390,7 +390,7 @@ namespace evalhyd ); } return qs.value(); - };; + }; xt::xtensor<double, 3> get_crps() {