Updated rules.

adolfomunoz · Jun 5, 2024 · 6d9d2f5 · 6d9d2f5
1 parent 7806ea3
commit 6d9d2f5
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/newton-cotes/rules.h b/src/newton-cotes/rules.h
@@ -149,7 +149,9 @@ struct Simpson {
  template<typename Float, typename T,typename Norm = NormDefault> 
  constexpr Float pdf(Float t, const std::array<T,samples>& p, Float a, Float b, const Norm& norm = Norm()) const {
  Float num = norm(at(t,p));
- return (num<1.e-10)?Float(0):(num/pdf_integral_subrange(a,b,p,norm));
+ Float den = pdf_integral_subrange(a,b,p,norm);
+ if (den<1.e-10) return 1.0/(b-a); 
+ else return num/den;
  }
 
  //My functions are below this line