diff --git a/include/pfasst/quadrature.hpp b/include/pfasst/quadrature.hpp index ff2b7d63..301ca72b 100644 --- a/include/pfasst/quadrature.hpp +++ b/include/pfasst/quadrature.hpp @@ -281,7 +281,9 @@ namespace pfasst if (abs(num) > 1e-32) { mat(i, j) = num / den; - } + } else { + mat(i, j) = 0.0; + } } }