Jacobi SVD fails to terminate on certain matrices containing NaN
Submitted by ren..@..mu.edu
Assigned to Nobody
Link to original bugzilla bug (#791)
Version: 3.2
Operating system: Linux
Description
Created attachment 450
program illustrating the problem
JacobiSVD fails to terminate on this 3x3 matrix of doubles:
0, 2.22045e-16, 0.5,
0, -0.5, 0,
nan, 0, 0;
and on this 4x4 matrix matrix of doubles:
1, 0, 0, 0,
0, 3, 1, 2e-308,
1, 0, 1, nan,
0, nan, nan, 0;
The attached program illustrates this behavior.
We expect termination because the documentation claims that "the computation is guaranteed to terminate in finite (and reasonable) time"
I can reproduce the problem in revision 5877:cd6ffdff46d2d3b9de945d81e20604c9648b83a1.
(This issue was found as part of the ASTAA robustness testing project.)
Attachment 450, "program illustrating the problem":
bug134a.cpp