catches error in the computation of confidence intervals for null tim…

…e series
wesleybowman · Jun 21, 2018 · f00d018 · f00d018 · wesleybowman · Jun 22, 2018
1 parent d5548e4
commit f00d018
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/utide/confidence.py b/utide/confidence.py
@@ -252,7 +252,9 @@ def _confidence(coef, cnstit, opt, t, e, tin, elor, xraw, xmod, W, m, B,
             varcov_mCw[c, :2, :2] = Duu
 
             if not opt.white:
-                Duu = Puu[c] * Duu / np.trace(Duu)
+                Duu = Puu[c] * Duu 
+                if np.trace(Duu) != 0. :
+                    Duu = Duu / np.trace(Duu)
                 varcov_mCc[c, :2, :2] = Duu
 
             if not opt.twodim: