Merge pull request #307 from MiraGeoscience/GEOPY-2569

GEOPY-2569: Wrong normalization for MT

Author: domfournier

simpeg_drivers/components/data.py

@@ -300,7 +300,7 @@ def get_normalizations(self):
                 if comp in ["gz", "bz", "gxz", "gyz", "bxz", "byz"]:
                     normalizations[chan][comp] = -1 * np.ones(self.mask.sum())
                 elif self.params.inversion_type in ["magnetotellurics"]:
-                    normalizations[chan][comp] = -1 * np.ones(self.mask.sum())
+                    normalizations[chan][comp] = np.ones(self.mask.sum())
                 elif self.params.inversion_type in ["tipper"]:
                     if "imag" in comp:
                         normalizations[chan][comp] = -1 * np.ones(self.mask.sum())

tests/run_tests/driver_mt_test.py

@@ -41,7 +41,7 @@
 # To test the full run and validate the inversion.
 # Move this file out of the test directory and run.
 
-target_run = {"data_norm": 0.032649770, "phi_d": 6.68, "phi_m": 263}
+target_run = {"data_norm": 0.032649770, "phi_d": 7.13, "phi_m": 282}
 
 
 def setup_data(workspace, survey):
@@ -187,7 +187,7 @@ def test_magnetotellurics_run(tmp_path: Path, max_iterations=1, pytest=True):
         output["data"] = orig_zyy_real_1
         assert not run_ws.get_entity("Iteration_0_sensitivities")[0]
         if pytest:
-            check_target(output, target_run, tolerance=0.2)
+            check_target(output, target_run, tolerance=0.1)
             nan_ind = np.isnan(run_ws.get_entity("Iteration_0_model")[0].values)
             inactive_ind = run_ws.get_entity("active_cells")[0].values == 0
             assert np.all(nan_ind == inactive_ind)