wip - debug

include/evalhyd/detail/probabilist/intervals.hpp

@@ -67,11 +67,11 @@ 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 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: " << res[0][1] << std::endl;
-                    std::cout << "res: " << res[1][1] << std::endl;
+                    std::cout << "res: " << res[0] << std::endl;
+                    std::cout << "res: " << res[1] << std::endl;
                     if (res.size() != 2) 
                     {
                         throw std::runtime_error(