From 93328af6a668031e77279bf336e03c9575f65640 Mon Sep 17 00:00:00 2001
From: fbourgin <francois.bourgin@inrae.fr>
Date: Thu, 24 Apr 2025 14:39:16 +0200
Subject: [PATCH] wip - debug
---
include/evalhyd/detail/probabilist/intervals.hpp | 5 +++++
1 file changed, 5 insertions(+)
diff --git a/include/evalhyd/detail/probabilist/intervals.hpp b/include/evalhyd/detail/probabilist/intervals.hpp
index e0dcd25..4e4d29a 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(
--
GitLab