qp() in Octave 3.0.0 returns result egregiously violating input constraints
Gabriele Pannocchia
g.pannocchia at ing.unipi.it
Sun Apr 6 15:41:04 CDT 2008
> It appears to work for me.
>
> octave:44> INFO
> INFO =
> {
> solveiter = 200
> info = 3
> }
>
> octave:45> OBJ
> OBJ = 2.5372e-06
> octave:46>
Do you also know if the starting point computed by glpk (if you do not
supply the starting point I sent) is feasible or not?
> I also tried to determine which rows were unique, by using the code
> below.
>
> [jnk, n] = unique ([A;Q';H;LB';UB';A_IN]', 'rows');
>
> When I solved the problem using only the unique rows (there are 119
> unique rows) ...
>
> n = sort(n);
> [X0, H, Q, A, B, LB, UB, A_LB, A_IN, A_UB] = mod_qp(n);
> [X, OBJ, INFO, LAMBDA] = qp (X0, H, Q, A, B, LB, UB, A_LB, A_IN, A_UB);
>
> The result is
>
> OBJ = 2.5114e-05
> octave:50> INFO
> INFO =
> {
> solveiter = 56
> info = 0
> }
>
> Assuming I haven't got something wrong ... Is this sort of behavior to
> be expected?
It is hard to tell... You removed 40 variables from the problem, so that
this is not the same QP.
Gabriele
More information about the Bug-octave
mailing list