diff --git a/include/evalhyd/detail/probabilist/intervals.hpp b/include/evalhyd/detail/probabilist/intervals.hpp
index e0dcd25c9a2976e13f6b19b9b599909e414b554d..4e4d29aff4528104ee132715f4a45042c2eda024 100644
--- a/include/evalhyd/detail/probabilist/intervals.hpp
+++ b/include/evalhyd/detail/probabilist/intervals.hpp
@@ -69,11 +69,16 @@ namespace evalhyd
                     auto a = xt::broadcast(xt::view(quantiles, i), std::vector<std::size_t>({q_lvl.size(), 2}));
                     auto b = xt::broadcast(q_lvl / 100., std::vector<std::size_t>({2, q_lvl.size()}));
                     auto res = xt::where(xt::equal(a, xt::transpose(b)));
+                    
                     std::cout << "res: " << res.size() << std::endl;
                     std::cout << "res_00: " << res[0][0] << std::endl;
                     std::cout << "res_01: " << res[0][1] << std::endl;
                     std::cout << "res_10: " << res[1][0] << std::endl;
                     std::cout << "res_11: " << res[1][1] << std::endl;
+
+                    std::cout << "res_min: " <<  std::min(res[0][0], res[0][1]) << std::endl;
+                    std::cout << "res_max: " << std::max(res[0][0], res[0][1]) << std::endl;
+
                     if (res.size() != 2) 
                     {
                         throw std::runtime_error(