fix mistake in RANK_DIAG computation

afa26b47 · Thibault Hallouin · 3743f8e3 · afa26b47
Commit afa26b47 authored 2 years ago by Thibault Hallouin
Hide whitespace changes
Inline Side-by-side

Showing

with 4 additions and 6 deletions
+4 -6
--- a/include/evalhyd/detail/probabilist/ranks.hpp
+++ b/include/evalhyd/detail/probabilist/ranks.hpp
@@ -240,12 +240,10 @@ namespace evalhyd
                        // compute the rank diagram
                        xt::view(REL_DIAG, xt::all(), xt::all(), m, e, xt::all()) =
-                                xt::sum(
+                                xt::view(o_j, xt::all(), xt::all(),
-                                        xt::view(o_j, xt::all(), xt::all(),
+                                         m, e, xt::all())
-                                                 m, e, xt::all())
+                                / l
-                                        / l,
+                        ;
-                                        -1
-                                );
                    }
                }