From 317dd1daf7f06932abe2ba90a291b084c3502437 Mon Sep 17 00:00:00 2001 From: dominiquef Date: Wed, 21 Feb 2024 10:02:22 -0800 Subject: [PATCH 1/2] Fix issue with duplicate log on reference --- geoapps/inversion/components/models.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/geoapps/inversion/components/models.py b/geoapps/inversion/components/models.py index ef6f4e74a..3cf0b40bd 100644 --- a/geoapps/inversion/components/models.py +++ b/geoapps/inversion/components/models.py @@ -133,8 +133,8 @@ def reference(self) -> np.ndarray | None: return mref if mref is None or (self.is_sigma and all(mref == 0)): - mref = self.starting self.driver.params.alpha_s = 0.0 + return self.starting.copy() ref_model = mref.copy() ref_model = np.log(ref_model) if self.is_sigma else ref_model From afb886cfef8404eaef83f760d0bca896e727dd78 Mon Sep 17 00:00:00 2001 From: dominiquef Date: Wed, 21 Feb 2024 10:11:54 -0800 Subject: [PATCH 2/2] Add Expmap on joint --- geoapps/inversion/joint/joint_surveys/driver.py | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/geoapps/inversion/joint/joint_surveys/driver.py b/geoapps/inversion/joint/joint_surveys/driver.py index 46aab255b..6e4265377 100644 --- a/geoapps/inversion/joint/joint_surveys/driver.py +++ b/geoapps/inversion/joint/joint_surveys/driver.py @@ -98,9 +98,15 @@ def directives(self): ) self._directives = DirectivesFactory(self) + global_model_save = self._directives.save_iteration_model_directive + if self.models.is_sigma: + global_model_save.transforms += [ + maps.ExpMap(self.inversion_mesh.mesh) + ] + self._directives.directive_list = ( self._directives.inversion_directives - + [self._directives.save_iteration_model_directive] + + [global_model_save] + directives_list ) return self._directives