bugGNU Octave - Bugs: bug #37742, sqp test failure with ATLAS on...

 
 

bug #37742: sqp test failure with ATLAS on Debian amd64

Submitter:  Mike Miller <mtmiller>
Submitted:  Thu 15 Nov 2012 03:42:16 AM UTC
   
 
Category:  Libraries Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  Fixed Assigned to:  mtmiller
Originator Name:  Open/Closed:  * Closed
Release:  * dev Operating System:  * GNU/Linux
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Sat 08 Dec 2012 11:21:20 PM UTC, comment #1: 

Fixed by increasing the tolerance on the failing test. This same test tolerance had already been increased a few times over the years.

http://hg.savannah.gnu.org/hgweb/octave/rev/534a2c881f45

Mike Miller <mtmiller>
Group Member
Thu 15 Nov 2012 03:42:16 AM UTC, original submission:  

I'm not using sqp, but I have been getting a consistent failure when running make check, affects both stable and default branches for me. Error is slightly over the tolerance specified in the test case.

The gotcha is it only fails when libblas.so.3 is configured to point to ATLAS, and only on Debian 64-bit (testing or unstable). This is the Debian libatlas3-base package, not a self-compiled ATLAS.

This could conceivably be a bug with the atlas library in Debian, but the result of the test as shown below is so close and the tolerance of 5*sqrt(eps) looks arbitrary to me, I think this could be resolved by bumping up the tolerance here.


$ readlink -f /usr/lib/libblas.so.3
/usr/lib/atlas-base/atlas/libblas.so.3.0
$ ./run-octave -qf --eval="test sqp"
  ***** test

 x0 = [-1.8; 1.7; 1.9; -0.8; -0.8];

 [x, obj, info, iter, nf, lambda] = sqp (x0, @__phi, @__g, []);

 x_opt = [-1.717143501952599;
           1.595709610928535;
           1.827245880097156;
          -0.763643103133572;
          -0.763643068453300];

 obj_opt = 0.0539498477702739;

 assert (x, x_opt, 5*sqrt (eps));
 assert (obj, obj_opt, sqrt (eps));
!!!!! test failed
assert (x,x_opt,5 * sqrt (eps)) expected
  -1.71714
   1.59571
   1.82725
  -0.76364
  -0.76364
but got
  -1.71714
   1.59571
   1.82725
  -0.76364
  -0.76364
maximum absolute error 9.34422e-08 exceeds tolerance 7.45058e-08


Mike Miller <mtmiller>
Group Member

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach Files:
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by mtmiller (Submitted the item)
  •  

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    Follow 3 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2012-12-08 mtmiller StatusNone Fixed
        Assigned toNone mtmiller
        Open/ClosedOpen Closed

    Back to the top

    Powered by Savane 3.13-758e.
    Corresponding source code