Hi,
I have a question regarding the use of sqp.m in octave 3.5.0+
He example in the help works fine, but if we add inequalities
constraints it fails (or I do not know how to specify them)
The code is below
function r = g (x)
r = [ sumsq(x)-10;
x(2)*x(3)-5*x(4)*x(5);
x(1)^3+x(2)^3+1 ];
endfunction
function obj = phi (x)
obj = exp(prod(x)) - 0.5*(x(1)^3+x(2)^3+1)^2;
endfunction
function l = h(x)
l = 1 - x.^2;
endfunction
x0 = [-1.8; 1.7; 1.9; -0.8; -0.8];
[x, obj, info, iter, nf, lambda] = sqp (x0, @phi, @g, @h)
error: sqp: operator *: nonconformant arguments (op1 is 5x10, op2 is 0x0)
error: called from:
error: sqp.m at line 472, column 9
As far as I can understand, qp returns empty LAMBDA vector when the
solution is not feasible. This condition is not being checked in sqp
and LAMBDA is used when it is empty.
Any suggestion?
Thanks
--
M. Sc. Juan Pablo Carbajal
-----
PhD Student
University of Zürich
www.ailab.ch/carbajal
------------------------------------------------------------------------------
Magic Quadrant for Content-Aware Data Loss Prevention
Research study explores the data loss prevention market. Includes in-depth
analysis on the changes within the DLP market, and the criteria used to
evaluate the strengths and weaknesses of these DLP solutions.
http://www.accelacomm.com/jaw/sfnl/114/51385063/
_______________________________________________
Octave-dev mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/octave-dev
