We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent eecff65 commit 9a9c6e8Copy full SHA for 9a9c6e8
include/albatross/src/stats/incomplete_gamma.hpp
@@ -35,7 +35,12 @@ constexpr double INCOMPLETE_GAMMA_EQUALITY_TRESHOLD = 1e-12;
35
36
inline double incomplete_gamma_quadrature_inp_vals(double lb, double ub,
37
std::size_t counter) {
38
- ALBATROSS_ASSERT(counter < gauss_legendre_50_points.size());
+ if (counter >= gauss_legendre_50_points.size()) {
39
+ std::cerr << "ERROR: Counter value " << counter
40
+ << " is higher than gauss_legendre_50_points size "
41
+ << gauss_legendre_50_points.size() << std::endl;
42
+ ALBATROSS_ASSERT(false);
43
+ }
44
return (ub - lb) * 0.5 * gauss_legendre_50_points[counter] + 0.5 * (ub + lb);
45
}
46
0 commit comments