Set lower bound for A0

kandread · Oct 15, 2019 · 3951313 · 3951313
1 parent b4418f0
commit 3951313
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/confluence/integrator.py b/confluence/integrator.py
@@ -85,6 +85,7 @@ def integrate(self):
         ubnds = [0.0] * (self.nreaches * 2)
         lbnds[1::2] = [0.01] * self.nreaches
         ubnds[1::2] = [0.09] * self.nreaches
+        lbnds[::2] = np.abs(np.amin(self.data[-1], axis=0))
         ubnds[::2] = [1e4] * self.nreaches
         bnds = list(zip(lbnds, ubnds))
         x0 = [0.0] * (self.nreaches * 2)