diff --git a/include/evalhyd/detail/probabilist/contingency.hpp b/include/evalhyd/detail/probabilist/contingency.hpp index 771f7681c80c3edf25f92301ee94ece5e9d524b0..533e4ba5f3c15273a3c88a1563fbbfd79cbaa823 100644 --- a/include/evalhyd/detail/probabilist/contingency.hpp +++ b/include/evalhyd/detail/probabilist/contingency.hpp @@ -56,9 +56,10 @@ namespace evalhyd auto alert_lvl = xt::arange<double>(double(n_mbr + 1)); // determine whether forecast yield alert - return sum_f_k >= - xt::view(alert_lvl, xt::newaxis(), xt::newaxis(), - xt::all(), xt::newaxis(), xt::newaxis()); + return xt::view(sum_f_k, xt::all(), xt::all(), xt::newaxis(), + xt::all(), xt::all()) + >= xt::view(alert_lvl, xt::newaxis(), xt::newaxis(), + xt::all(), xt::newaxis(), xt::newaxis()); } /// Determine hits ('a' in contingency table).