kwarg for stepping off bounds in projected gradient

SimPEG/Optimization.py

@@ -888,6 +888,8 @@ def __init__(self, **kwargs):
     maxIterCG = 5
     tolCG = 1e-1
 
+    stepOffBoundsFact = 0.1 # perturbation of the inactive set off the bounds
+
     lower = -np.inf
     upper = np.inf
 
@@ -998,7 +1000,8 @@ def findSearchDirection(self):
                 dm_i = max( abs( delx ) )
                 dm_a = max( abs(rhs_a) )
 
-                delx = delx + rhs_a * dm_i / dm_a /10.
+                # perturb inactive set off of bounds so that they are included in the step
+                delx = delx + self.stepOffBoundsFact * (rhs_a * dm_i / dm_a)
 
             # Only keep gradients going in the right direction on the active set
             indx = ((self.xc<=self.lower) & (delx < 0)) | ((self.xc>=self.upper) & (delx > 0))