diff options
| author | Keith Packard <keithp@keithp.com> | 2017-10-12 00:30:23 -0700 |
|---|---|---|
| committer | Keith Packard <keithp@keithp.com> | 2017-10-12 00:30:23 -0700 |
| commit | 77d1aee917306ad59492c4c8352fe2125b430d0c (patch) | |
| tree | a07f2aece79df2fd6cc267d4ba76399db6ef6939 | |
| parent | e04679ba52761d5531037117a21ab1f1896358b0 (diff) | |
altoslib: Fix time series filter window computation
Small floating point rounding errors could lead to NaNs.
Signed-off-by: Keith Packard <keithp@keithp.com>
| -rw-r--r-- | altoslib/AltosTimeSeries.java | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/altoslib/AltosTimeSeries.java b/altoslib/AltosTimeSeries.java index 7208c176..c6a780a3 100644 --- a/altoslib/AltosTimeSeries.java +++ b/altoslib/AltosTimeSeries.java @@ -294,7 +294,7 @@ public class AltosTimeSeries implements Iterable<AltosTimeValue>, Comparable<Alt double alpha = m / 2; double t = (n - alpha) / alpha; - if (t > 1) + if (t > 1 || t < -1) t = 1; double k = i0 (beta * Math.sqrt (1 - t*t)) / i0(beta); return k; |