From fa638251d7d080fb9a0ba1b4fefb1a58c24d039e Mon Sep 17 00:00:00 2001
From: "E. Joshua Rigler" <erigler@usgs.gov>
Date: Mon, 1 Oct 2018 13:16:05 -0600
Subject: [PATCH] Small fix to SqDistAlgorithm initiization
Use NumPy's `isnan(var)` function instead of `var is np.nan`, since
the latter doesn't work if var is a 64bit float NaN.
---
geomagio/algorithm/SqDistAlgorithm.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/geomagio/algorithm/SqDistAlgorithm.py b/geomagio/algorithm/SqDistAlgorithm.py
index a3e65e911..acc66392a 100644
--- a/geomagio/algorithm/SqDistAlgorithm.py
+++ b/geomagio/algorithm/SqDistAlgorithm.py
@@ -332,7 +332,7 @@ class SqDistAlgorithm(Algorithm):
# set some default values
if l0 is None:
l = np.nanmean(yobs[0:int(m)])
- if l is np.nan:
+ if np.isnan(l):
l = 0.
else:
l = l0
--
GitLab