qp() in Octave 3.0.0 returns result egregiously violating input constraints
Joshua Redstone
redstone at gmail.com
Fri Apr 4 20:24:09 CDT 2008
Hi John,
I've attached a data file created with (on OSX):
save -binary -z "./example.dat" x0 H Q A B LB UB A_LB A_IN A_UB
The following should reproduce the problem for you:
load "example.dat"
[X, OBJ, INFO, LAMBDA] = qp(x0, H, Q, A, B, LB, UB, A_LB, A_IN, A_UB)
With my qp.m that I modified so the last arguement is maxiter, I get the
follow first few lines of output:
octave-3.0.0:5> [X, OBJ, INFO, LAMBDA] = qp(x0, H, Q, A, B, LB, UB, A_LB,
A_IN, A_UB, 187)
X =
0.17954
0.00000
0.01742
0.00304
0.16220
[... output truncated.... ]
octave-3.0.0:6> [X, OBJ, INFO, LAMBDA] = qp(x0, H, Q, A, B, LB, UB, A_LB,
A_IN, A_UB, 188)
X =
8.7626e+13
5.1887e+00
1.0476e+13
-9.8103e+13
2.0382e+15
-6.1036e+00
2.3327e+00
[.... output truncated ....]
I'm not sure how to scale down the problem or generate it with a script as
other datasets of this type don't exhibit this behavior. The numbers
originate from some empirical data sets.
Let me know if there anything else I can do to help.
Thanks,
Josh
On Fri, Apr 4, 2008 at 1:48 PM, Joshua Redstone <redstone at gmail.com> wrote:
> OK - I"ll get it in that form. I think I may have lost a race with the
> moderator and so you may see an old email come thru
> with an attachment in a less usefull format.
> Josh
>
>
> On Fri, Apr 4, 2008 at 1:32 PM, John W. Eaton <jwe at bevo.che.wisc.edu>
> wrote:
>
> > On 4-Apr-2008, Joshua Redstone wrote:
> >
> > | Hi Gabriele,
> > | Thanks for your resonse. The problematic example involves 100
> > variables.
> > | What format would be most convenient for me to send it in?
> > | I've set 'format long g' and printed out all the args to qp() that I
> > use,
> > | naming them as they are named in qp.m, and attached the result
> > | to this email (gzipped).
> >
> > How about in a format that can be read by Octave, along with the
> > command that you use to invoke the qp function? Make it so that we
> > can do
> >
> > load somefile.dat
> > qp (?)
> >
> > where you tell us what ? should be so that we can reproduce the
> > problem.
> >
> > Where do the numbers come from? Is it possible to generate them with
> > a script?
> >
> > Make it easy for us to reproduce the problem and we are more
> > likely to help debug it...
> >
> > jwe
> >
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://www.cae.wisc.edu/pipermail/bug-octave/attachments/20080404/7bdc3f5c/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: example.dat
Type: application/octet-stream
Size: 82018 bytes
Desc: not available
Url : https://www.cae.wisc.edu/pipermail/bug-octave/attachments/20080404/7bdc3f5c/attachment-0001.obj
More information about the Bug-octave
mailing list