diff --git a/include/evalhyd/detail/probabilist/intervals.hpp b/include/evalhyd/detail/probabilist/intervals.hpp
index 340337a4d4b373f4d908cf36a09bdfec992acf03..d1d1d537ce27cc4eb910ce69b99c8a50d7573d91 100644
--- a/include/evalhyd/detail/probabilist/intervals.hpp
+++ b/include/evalhyd/detail/probabilist/intervals.hpp
@@ -67,21 +67,21 @@ namespace evalhyd
                     for (std::size_t i = 0; i < n_itv; i++)
                     {
                     auto a = xt::broadcast(xt::view(quantiles, i), std::vector<std::size_t>({q_lvl.size(), 2}));
-                    auto b = xt::broadcast(q_lvl, std::vector<std::size_t>({2, q_lvl.size()}));
-                    auto res = xt::where(xt::equal(a, xt::transpose(b)));
-                    if (res.size() != 2) 
-                    {
-                        throw std::runtime_error(
-                                "interval-based metric requested, "
-                                "but *c_lvl* not matching *q_lvl"
-                        );
-                    } else 
-                    {
+                    //auto b = xt::broadcast(q_lvl, std::vector<std::size_t>({2, q_lvl.size()}));
+                    //auto res = xt::where(xt::equal(a, xt::transpose(b)));
+                    //if (res.size() != 2) 
+                    //{
+                    //    throw std::runtime_error(
+                    //            "interval-based metric requested, "
+                    //            "but *c_lvl* not matching *q_lvl"
+                    //    );
+                    //} else 
+                    //{
                         xt::view(itv_bnds, xt::all(), xt::all(), i, 0, xt::all()) =
                                 xt::view(q_prd, 0);
                         xt::view(itv_bnds, xt::all(), xt::all(), i, 1, xt::all()) = 
                                 xt::view(q_prd, 1);
-                    }
+                    //}
                     }
                 }
                 else