Check if standard deviation is zero with np.isclose

sblunt · Sep 24, 2024 · c409f49 · c409f49
1 parent ed80f58
commit c409f49
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/orbitize/plot.py b/orbitize/plot.py
@@ -112,7 +112,7 @@ def plot_corner(results, param_list=None, **corner_kwargs):
         index_num = results.param_idx[label_key]
 
         # only plot non-fixed parameters
-        if np.std(results.post[:, index_num]) > 0:
+        if not np.isclose(0.0, np.std(results.post[:, index_num])):
             param_indices.append(index_num)
             if (
                 label_key.startswith("aop")