commit 9f17f6de64d5f8376bb0eda601d8c2d730e70d5e
author Gilles Sadowski <gilleseran@gmail.com> Wed May 27 01:52:14 2020 +0200
committer Gilles Sadowski <gilleseran@gmail.com> Wed May 27 01:52:14 2020 +0200
tree 44d59bf094c21e1116e1666a5413b9bf8b9d6e69
parent 52b374d85e37235805cc0820c8f49cb139432d02

MATH-1534: Use value that probably conforms to the original intent.

Thanks to Christoph Läubrich.

src/main/java/org/apache/commons/math4/fitting/GaussianCurveFitter.java

diff --git a/src/main/java/org/apache/commons/math4/fitting/GaussianCurveFitter.java b/src/main/java/org/apache/commons/math4/fitting/GaussianCurveFitter.java
index 711fc1d..4af1200 100644
--- a/src/main/java/org/apache/commons/math4/fitting/GaussianCurveFitter.java
+++ b/src/main/java/org/apache/commons/math4/fitting/GaussianCurveFitter.java
@@ -294,11 +294,10 @@
         private double[] basicGuess(WeightedObservedPoint[] points) {
             final int maxYIdx = findMaxY(points);
             final double n = points[maxYIdx].getY();
-            final double m = points[maxYIdx].getX();
 
             double fwhmApprox;
             try {
-                final double halfY = n + ((m - n) / 2);
+                final double halfY = 0.5 * n;
                 final double fwhmX1 = interpolateXAtY(points, maxYIdx, -1, halfY);
                 final double fwhmX2 = interpolateXAtY(points, maxYIdx, 1, halfY);
                 fwhmApprox = fwhmX2 - fwhmX1;
@@ -308,7 +307,7 @@
             }
             final double s = fwhmApprox / (2 * FastMath.sqrt(2 * FastMath.log(2)));
 
-            return new double[] { n, m, s };
+            return new double[] { n, points[maxYIdx].getX(), s };
         }
 
         /**