Skip to content

Commit 5748fc3

Browse files
committed
Add small constant. Fix issue with lower_bound
1 parent 0c7fa28 commit 5748fc3

1 file changed

Lines changed: 6 additions & 5 deletions

File tree

simpeg_drivers/components/models.py

Lines changed: 6 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -233,6 +233,7 @@ def reference_model(self) -> np.ndarray | None:
233233
self.reference_inclination,
234234
self.reference_declination,
235235
)
236+
ref_model += 1e-8 # To avoid zeroing the angles
236237
ref_model = (field_vecs.T * ref_model).flatten()
237238

238239
return ref_model
@@ -271,11 +272,11 @@ def lower_bound(self) -> np.ndarray | None:
271272
else:
272273
bound_model = self._lower_bound.model
273274

274-
if (
275-
self.driver.params.inversion_type == "magnetic vector"
276-
and self._upper_bound.model is not None
277-
):
278-
bound_model = -self._upper_bound.model
275+
if self.driver.params.inversion_type == "magnetic vector":
276+
bound_model = None
277+
278+
if self._upper_bound.model is not None:
279+
bound_model = -self._upper_bound.model
279280

280281
if bound_model is None:
281282
lbound = np.full(self.n_active, -np.inf)

0 commit comments

Comments
 (0)