ALGLIB forum
It is currently Tue May 28, 2024 12:06 pm

All times are UTC

Forum rules

1. This forum can be used for discussion of both ALGLIB-related and general numerical analysis questions
2. This forum is English-only - postings in other languages will be removed.

Post new topic Reply to topic  [ 1 post ] 
Author Message
 Post subject: False positive in test_c, minqp, Linearly Constrained - LC.
PostPosted: Tue Oct 06, 2020 8:46 pm 

Joined: Mon Nov 20, 2017 11:14 pm
Posts: 21
In testminqpunit_generallcqptest(), under "SMALL-SCALE TESTS: many tests for small N's", with solvertype = 1, under "Convex optimization problem with known answer", with n = 8 and akind = 1:

gtrial = b + a × (x1 - xorigin) + rep.lagbc + rep.laglc × rawc.

For i ∈ [0,n): the test is for |gtrial[i]| ≤ 0.001, but is only flagging if |gtrial[i]| > 0.001. However, gtrial[i] may be NAN, because rep.laglc occasionally has components equal to NAN. So it is possible for the test to fail, and yet fail to be flagged, resulting in an "ok" for the test - a false positive.

It may be the case that your algorithm allows rep.laglc to have NAN values and treats them as infinitesimals, but this should be examined more carefully. The false positive occurs about 5% of the time that the minqp test in test_c is run.

More generally, most of the tests in test_c flag on reverse logic for real-number comparisons, and may therefore yield false positives if NAN's occur. In our local version, we eliminated the possibility by making all the tests positive logic on real number comparisons. So far, this is the only false positive that's been uncovered, as a result.

Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

All times are UTC

Who is online

Users browsing this forum: No registered users and 1 guest

You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group