diff --git a/include/evalhyd/detail/probabilist/intervals.hpp b/include/evalhyd/detail/probabilist/intervals.hpp
index f5ed84657d3f6a1d1da2af728ac02d5f2bcd5f39..a66396930cab9da0ae8974ee4073a9ac807027bf 100644
--- a/include/evalhyd/detail/probabilist/intervals.hpp
+++ b/include/evalhyd/detail/probabilist/intervals.hpp
@@ -68,7 +68,7 @@ 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, std::vector<std::size_t>({2, q_lvl.size()}));
-                    auto res = xt::where(xt::equal(a, xt::transpose(b));
+                    auto res = xt::where(xt::equal(a, xt::transpose(b)));
                     if (res.size() != 2) 
                     {
                         throw std::runtime_error(