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