bugGNU Octave - Bugs: bug #31731, Test failed for sqp

 
 

bug #31731: Test failed for sqp

Submitter:  Marco Caliari <caliari>
Submitted:  Mon 22 Nov 2010 03:49:09 PM UTC
   
 
Category:  None Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  Works For Me Assigned to:  None
Originator Name:  Open/Closed:  * Closed
Release:  * 3.3.53 Operating System:  * GNU/Linux
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Thu 02 Dec 2010 06:34:45 PM UTC, comment #3: 

Since you can reach the correct result by adjusting the tolerance, and this seems a satisfactory outcome to you, I'm closing this report.

Rik <rik5>
Group administrator
Tue 23 Nov 2010 08:31:39 AM UTC, comment #2: 

I'm on a linux system, i686.
With the value x_opt found in the test, I see


phi(x_opt) = 5.39498477702741e-02
max (g(x_opt)) = 3.55271367880050e-15


while with the result x I get (with tol = eps)


phi(x) = 5.39498477702744e-02
max (g(x_opt)) = 5.32907051820075e-15


(so, both x_opt and x are double precision accurate), but


norm (x_opt-x, inf) = 1.80132234683938e-08


The result I get is satisfactory for me, but I cannot understand the difference.

Marco

Marco Caliari <caliari>
Group Member
Tue 23 Nov 2010 03:09:19 AM UTC, comment #1: 

I can't reproduce this issue.  I ran the tests for sqp 500 times and they never failed.  This is on a Linux system, Kubuntu 9.04, x86_64.

I did change the documentation to reflect the actual tolerance value.  See this changeset: http://hg.savannah.gnu.org/hgweb/octave/rev/646063a22a35

Rik <rik5>
Group administrator
Mon 22 Nov 2010 03:49:09 PM UTC, original submission:  

Dear all,

I noticed that the first test fails (for me) in sqp. I can fix it if I use 7*sqrt(eps) instead of 5*sqrt(eps) in the assert. It is not dramatic, but I find a little bit strange this difference in the result (around single precision) w.r.t. the author of the test.

Moreover, in the documentation, it is written that tol = eps is used by default (while tol = sqrt (eps) is used).

Marco

Marco Caliari <caliari>
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 rik5 (Posted a comment)
  • -email is unavailable- added by caliari (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 2 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2010-12-02 rik5 Open/ClosedOpen Closed
    2010-11-23 rik5 StatusNone Works For Me

    Back to the top

    Powered by Savane 3.13-02a9.
    Corresponding source code