// Patch for bug /. Detecting this sort of issue earlier
// (via abnormalities in the indicators) would be better.
bounds = np.array(bounds) // again, this should have been the standard form
lb = bounds[:, 0]
ub = bounds[:, 1]
lb[np.equal(lb, None)] = -np.inf
ub[np.equal(ub, None)] = np.inf
return x, fun, slack, con, lb, ub
After Change
if i in no_adjust:
continue
lbi = bi[0]
ubi = bi[1]
if lbi == -np.inf and ubi == np.inf:
n_unbounded += 1
x[i] = x[i] - x[n_x + n_unbounded - 1]
else: