Create copy of models

geoapps/inversion/components/models.py

@@ -115,28 +115,31 @@ def active_cells(self, active_cells):
 
     @property
     def starting(self):
-        mstart = self._starting.model
+        mstart = self._starting.model.copy()
 
         if mstart is not None and self.is_sigma:
+            mstart = mstart.copy()
             mstart = np.log(mstart)
 
         return mstart
 
     @property
     def reference(self):
-        mref = self._reference.model
+        mref = self._reference.model.copy()
 
         if self.driver.params.forward_only:
             return mref
 
         if mref is None:
             mref = self.starting
             self.driver.params.alpha_s = 0.0
-        elif self.is_sigma & (all(mref == 0)):
-            mref = self.starting
-            self.driver.params.alpha_s = 0.0
         else:
-            mref = np.log(mref) if self.is_sigma else mref
+            mref = mref.copy()
+            if self.is_sigma & (all(mref == 0)):
+                mref = self.starting
+                self.driver.params.alpha_s = 0.0
+            else:
+                mref = np.log(mref) if self.is_sigma else mref
         return mref
 
     @property
@@ -146,6 +149,8 @@ def lower_bound(self):
         if lbound is None:
             return -np.inf
 
+        lbound = lbound.copy()
+
         if self.is_sigma:
             is_finite = np.isfinite(lbound)
             lbound[is_finite] = np.log(lbound[is_finite])
@@ -158,6 +163,7 @@ def upper_bound(self):
         if ubound is None:
             return np.inf
 
+        ubound = ubound.copy()
         if self.is_sigma:
             is_finite = np.isfinite(ubound)
             ubound[is_finite] = np.log(ubound[is_finite])
@@ -169,6 +175,7 @@ def conductivity(self):
         mstart = self._conductivity.model
 
         if mstart is not None and self.is_sigma:
+            mstart = mstart.copy()
             mstart = np.log(mstart)
 
         return mstart