qp() in Octave 3.0.0 returns result egregiously violating input constraints
Ben Abbott
bpabbott at mac.com
Sun Apr 6 12:12:06 CDT 2008
On Apr 5, 2008, at 5:58 PM, Thomas Treichl wrote:
> Ben Abbott schrieb:
>> I've built octave using Apple's lapack/blas (i.e. vecLib), so the
>> problem isn't directly related to the lapack/blas, correct?
>
> True. Then what about the difference between your PPC and our IAs?
> Just an idea...
ok, I pulled sources to my intel and ppc based Macs (each running OSX
10.5.2) and built both branches on each.
The "release-3-0-x" fails on the intel box.
> load example.dat
> [X, OBJ, INFO, LAMBDA] = qp (X0, H, Q, A, B, LB, UB, A_LB, A_IN,
A_UB);
> result = 0.5*X'*H*X + X'*Q
result = 5.4445e+58
The qp example appears to work on release-3-0-x running on the ppc.
result = 2.7676e-06
I get identical answers running the default branch on both my intel
and ppc based macs.
Ben
More information about the Bug-octave
mailing list