Fix Hipparchus-Math#316 EmpiricalDistribution#density() shouldn't ret…

…urn NaN see discussion in Hipparchus-Math#315
MaximeJo · Mar 6, 2024 · 57f8f3f · 57f8f3f
1 parent a8fe2cf
commit 57f8f3f
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/hipparchus-stat/src/main/java/org/hipparchus/stat/fitting/EmpiricalDistribution.java b/hipparchus-stat/src/main/java/org/hipparchus/stat/fitting/EmpiricalDistribution.java
@@ -578,7 +578,11 @@ public double density(double x) {
         }
         final int binIndex = findBin(x);
         final RealDistribution kernel = getKernel(binStats.get(binIndex));
-        return kernel.density(x) * pB(binIndex) / kB(binIndex);
+        double kernelDensity = kernel.density(x);
+        if (kernelDensity == 0d) {
+            return 0d;
+        }
+        return kernelDensity * pB(binIndex) / kB(binIndex);
     }
 
     /**