add tests for integral solutions

Author: MuellerSeb

tests/test_anaflow.py

@@ -221,6 +221,74 @@ def test_tpl_3d(self):
             self.assertTrue(dec(time_arr, self.delta))
         self.assertAlmostEqual(np.abs(np.max(steady - transient[-1, :])), 0.0, places=2)
 
+    def test_int(self):
+        steady = af.ext_thiem_int(
+            rad=self.rad,
+            r_ref=self.r_ref,
+            cond_gmean=self.t_gmean,
+            var=self.var,
+            len_scale=self.len_scale,
+            lat_ext=self.lat_ext,
+            roughness=self.hurst,
+            dim=self.frac_dim,
+            rate=self.rate,
+            h_ref=self.h_ref,
+        )
+        transient = af.ext_theis_int(
+            time=self.time,
+            rad=self.rad,
+            storage=self.storage,
+            cond_gmean=self.t_gmean,
+            var=self.var,
+            len_scale=self.len_scale,
+            lat_ext=self.lat_ext,
+            roughness=self.hurst,
+            dim=self.frac_dim,
+            rate=self.rate,
+            r_bound=self.r_ref,
+            h_bound=self.h_ref,
+        )
+        self.assertTrue(inc(steady, self.delta))
+        for rad_arr in transient:
+            self.assertTrue(inc(rad_arr, self.delta))
+        for time_arr in transient.T:
+            self.assertTrue(dec(time_arr, self.delta))
+        self.assertAlmostEqual(np.abs(np.max(steady - transient[-1, :])), 0.0, places=3)
+
+    def test_int_3d(self):
+        steady = af.ext_thiem_int_3d(
+            rad=self.rad,
+            r_ref=self.r_ref,
+            cond_gmean=self.t_gmean,
+            var=self.var,
+            len_scale=self.len_scale,
+            lat_ext=self.lat_ext,
+            roughness=self.hurst,
+            anis=self.anis,
+            rate=self.rate,
+            h_ref=self.h_ref,
+        )
+        transient = af.ext_theis_int_3d(
+            time=self.time,
+            rad=self.rad,
+            storage=self.storage,
+            cond_gmean=self.t_gmean,
+            var=self.var,
+            len_scale=self.len_scale,
+            lat_ext=self.lat_ext,
+            roughness=self.hurst,
+            anis=self.anis,
+            rate=self.rate,
+            r_bound=self.r_ref,
+            h_bound=self.h_ref,
+        )
+        self.assertTrue(inc(steady, self.delta))
+        for rad_arr in transient:
+            self.assertTrue(inc(rad_arr, self.delta))
+        for time_arr in transient.T:
+            self.assertTrue(dec(time_arr, self.delta))
+        self.assertAlmostEqual(np.abs(np.max(steady - transient[-1, :])), 0.0, places=2)
+
         # plt.plot(self.rad, steady)
         # for rad_arr in transient:
         #     plt.plot(self.rad, rad_arr, dashes=[1, 1])