[Octave-bug-tracker] [bug #31742] sqp: incorrect order of removal of Inf bounds

[Rik]

Follow-up Comment #9, bug #31742 (project octave):

I credited you for the patch in the more important ChangeLog file kept within
Mercurial.  I reverted the ",1" indexing change I made in this changeset
(http://hg.savannah.gnu.org/hgweb/octave/rev/c767bb1afa03).

I think you need to keep the tmp_bnds variable to handle the case of
different numbers of bounds.


lb = [1; 2; 3];
ub = [10; Inf; 20];

There is no upper constraint on the second variable, which leads to
__sqp_lb__ being of size 3x1 and __sqp_ub__ beof of size 2x1 and the
comparison "any (__sqp_lb__ > __sqp_ub__)" would fail.



This could be handled by calculating the index and not applying it until
after the comparison above.  However, I think the code would be less clear. 
It is easier to understand calculating the index and removing the inapplicable
constraints as one block of code.

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?31742>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/