Skip to content
GitLab
Commit 2388f65e authored by fbourgin's avatar fbourgin
Browse files

wip - debug

parent f1e032bf
12-evaluation-of-quantiles-prediction
1 merge request!9Draft: Resolve "evaluation of quantiles prediction"
Pipeline #63443 failed with stage
in 2 minutes and 13 seconds
Hide whitespace changes
Inline Side-by-side
Showing
with 14 additions and 14 deletions
+14 -14
include/evalhyd/detail/probabilist/intervals.hpp
+ 14
− 14
View file @ 2388f65e
......@@ -66,22 +66,22 @@ 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 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
{
xt::view(itv_bnds, xt::all(), xt::all(), i, 0, xt::all()) =
xt::view(q_prd, xt::all(), xt::all(), 0, xt::all());
xt::view(q_prd, xt::all(), xt::all(), std::min(res[0][1], res[1][1], xt::all());
xt::view(itv_bnds, xt::all(), xt::all(), i, 1, xt::all()) =
xt::view(q_prd, xt::all(), xt::all(), 1, xt::all());
//}
xt::view(q_prd, xt::all(), xt::all(), std::max(res[0][1], res[1][1], xt::all());
}
}
}
else
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment