[SciPy-User] fmin_slsqp constraint problem
Oleksandr Huziy
guziy.sasha at gmail.com
Fri Oct 26 11:04:03 EDT 2012
What is your obj_fnc, I know it is naive, bu still, is it possible that it
returns a list?
Cheers
--
Oleksandr (Sasha) Huziy
2012/10/26 Neal Becker <ndbecker2 at gmail.com>
> I have a ineq constraint:
>
> ## constrain poles to be inside unit circle
> def c(coef, len_z, len_p, dz, dp):
> p = compose ((coef/opt.scale)[len_z:-1], dp)
> return np.abs(p) - 1
>
> So this will return a 1D array where each value should satisfy the
> constraint.
> fmin_slsqp will not accept this directly:
>
> e = fmin_slsqp (obj_fnc, coef*opt.scale, ieqcons=[lambda coef, len_z,
> len_p,
> dz, dp: -c(coef, len_z, len_p, dz, dp)], args=(len(lz), len(lp), dz, dp),
> eqcons=[lambda coef, len_z, len_p, dz, dp: h(coef, len_z, len_p, dz, dp)],
> full_output=True)
>
> Traceback (most recent call last):
> File "./optimize_pll5.3.2.py", line 519, in <module>
> run_line (sys.argv)
> File "./optimize_pll5.3.2.py", line 498, in run_line
> e = fmin_slsqp (obj_fnc, coef*opt.scale, ieqcons=[lambda coef, len_z,
> len_p,
> dz, dp: -c(coef, len_z, len_p, dz, dp)], args=(len(lz), len(lp), dz, dp),
> eqcons=[lambda coef, len_z, len_p, dz, dp: h(coef, len_z, len_p, dz, dp)],
> full_output=True)
> File "/usr/lib64/python2.7/site-packages/scipy/optimize/slsqp.py", line
> 334,
> in fmin_slsqp
> a_ieq[i] = ieqcons_prime[i](x)
> File "/usr/lib64/python2.7/site-packages/scipy/optimize/optimize.py",
> line
> 176, in function_wrapper
> return function(x, *args)
> File "/usr/lib64/python2.7/site-packages/scipy/optimize/optimize.py",
> line
> 398, in approx_fprime
> grad[k] = (f(*((xk+ei,)+args)) - f0)/epsilon
> ValueError: setting an array element with a sequence.
>
> Any ideas on this?
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20121026/377377b0/attachment.html>
More information about the SciPy-User
mailing list