Confidence intervals (spectrums) now support 2D

treden · Oct 9, 2023 · 62f4326 · 62f4326
1 parent 0f4140b
commit 62f4326
Show file tree

Hide file tree

Showing 2 changed files with 212 additions and 164 deletions.
diff --git a/utide/_solve.py b/utide/_solve.py
@@ -609,26 +609,25 @@ def _solv2(tin, uin, vin, lat, **opts):
             coef.theta = np.hstack((coef.theta, theta))
             coef.g = np.hstack((coef.g, g))
 
-    # TODO: 2D OLS Confidence intervals
-    # if opt["conf_int"]:
-    #     coef = _confidence(
-    #         coef,
-    #         cnstit,
-    #         opt,
-    #         t,
-    #         e,
-    #         tin,
-    #         elor,
-    #         xraw,
-    #         xmod,
-    #         W,
-    #         m,
-    #         B,
-    #         Xu,
-    #         Yu,
-    #         Xv,
-    #         Yv,
-    #     )
+    if opt["conf_int"]:
+        coef = _confidence(
+            coef,
+            cnstit,
+            opt,
+            t,
+            e,
+            tin,
+            elor,
+            xraw,
+            xmod,
+            W,
+            m,
+            B,
+            Xu,
+            Yu,
+            Xv,
+            Yv,
+        )
 
     # Diagnostics.
     # TODO: 2D Diagnostics