Solved bug.

adolfomunoz · May 10, 2024 · f3cc24f · f3cc24f
1 parent fc9e1cd
commit f3cc24f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/newton-cotes/region.h b/src/newton-cotes/region.h
@@ -267,8 +267,8 @@ class Region {
  const std::array<Float,DIMSUB>& b,const Norm& norm = Norm()) const {
  static_assert(DIM>=DIMSUB,"Cannot calculate the subrange pdf for that many dimensions, as the region has less dimensions");
  Float den = this->pdf_integral_subrange(a,b,norm);
- if (den<1.e-10) den=1;
- return Range(a,b).volume()*norm(this->approximation_at(pos))/den;
+ if (den<1.e-10) return Float(0);
+ else return Range(a,b).volume()*norm(this->approximation_at(pos))/den;
  }
 
  template<std::size_t DIMSUB,typename Norm = NormDefault>