Mercurial > octave-nkf
diff src/DLD-FUNCTIONS/__qp__.cc @ 11879:57c9b3f2b822 release-3-0-x
__qp__.cc fix
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 31 Oct 2008 07:52:55 +0100 |
parents | 8e6dc8dd94e7 |
children |
line wrap: on
line diff
--- a/src/DLD-FUNCTIONS/__qp__.cc Sun Oct 12 19:38:43 2008 +0200 +++ b/src/DLD-FUNCTIONS/__qp__.cc Fri Oct 31 07:52:55 2008 +0100 @@ -233,6 +233,9 @@ // FIXME -- still remain to handle the case of // non-full rank active set matrix. + // Computing the Y matrix (orthogonal to Z) + Y = Aact.pseudo_inverse (); + // Reduced Hessian Matrix Zt = Z.transpose (); Matrix rH = Zt * H * Z; @@ -327,7 +330,6 @@ // Computing the multipliers only for the inequality // constraints that are active. We do NOT compute // multipliers for the equality constraints. - Y = Aact.pseudo_inverse (); Matrix Yt = Y.transpose (); Yt = Yt.extract_n (n_eq, 0, n_act-n_eq, n); lambda_tmp = Yt * (g + H * p);