[Octave-bug-tracker] [bug #54399] test eigs fail

[Dmitri A. Sergatskov]

Follow-up Comment #7, bug #54399 (project octave):

I did not notice that it was different tests that fail with openbls and with
reference blas. Attached two test scripts -- one fails with openblas, another
with ref blas. 
Essentially I used A = sin(a*magic(10)) variations to construct the matrix.
The values of "a" does not appear critical, it would fail every time I tried.


Openblas:


LD_PRELOAD=/usr/lib64/libopenblas.so ./run-octave -f -q
octave:1> eigs_fail_openblas
ans = ba8227df92ae
error: ASSERT errors for:  assert (Evector,Evector_f)

  Location  |  Observed  |  Expected  |  Reason
   (1,1)     0.040785+0i  0.040785+0i    Abs err 1.4572e-16 exceeds tol 0 by
1e-16
   (2,1)     -0.055363+0i -0.055363+0i   Abs err 2.7756e-17 exceeds tol 0 by
3e-17
   (3,1)     0.037331+0i  0.037331+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (4,1)      0.19551+0i   0.19551+0i    Abs err 1.3878e-16 exceeds tol 0 by
1e-16
   (5,1)      0.13734+0i   0.13734+0i    Abs err 8.3267e-17 exceeds tol 0 by
8e-17
   (6,1)     -0.040516+0i -0.040516+0i   Abs err 7.6328e-17 exceeds tol 0 by
8e-17
   (7,1)     0.045628+0i  0.045628+0i    Abs err 1.1102e-16 exceeds tol 0 by
1e-16
   (8,1)      -0.0456+0i   -0.0456+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (9,1)     -0.20177+0i  -0.20177+0i    Abs err 1.1102e-16 exceeds tol 0 by
1e-16
   (10,1)     -0.1282+0i   -0.1282+0i    Abs err 2.7756e-17 exceeds tol 0 by
3e-17
   (1,2)     -0.04529+0i  -0.04529+0i    Abs err 8.3267e-17 exceeds tol 0 by
8e-17
   (2,2)     -0.22421+0i  -0.22421+0i    Abs err 2.7756e-17 exceeds tol 0 by
3e-17
   (3,2)     0.081909+0i  0.081909+0i    Abs err 1.3878e-17 exceeds tol 0 by
1e-17
   (4,2)     -0.021274+0i -0.021274+0i   Abs err 1.3878e-16 exceeds tol 0 by
1e-16
   (5,2)      0.21574+0i   0.21574+0i    Abs err 2.7756e-17 exceeds tol 0 by
3e-17
   (8,2)     -0.074681+0i -0.074681+0i   Abs err 6.9389e-17 exceeds tol 0 by
7e-17
   (9,2)     0.030061+0i  0.030061+0i    Abs err 4.1633e-17 exceeds tol 0 by
4e-17
   (10,2)    -0.21714+0i  -0.21714+0i    Abs err 8.3267e-17 exceeds tol 0 by
8e-17
   (1,3)     0.27472-0.12463i 0.27472-0.12463i   Abs err 5.4672e-16 exceeds
tol 0 by 5e-16
   (2,3)     0.24391-0.08352i 0.24391-0.08352i   Abs err 4.003e-16 exceeds tol
0 by 4e-16
   (3,3)     -0.23805-0.022089i -0.23805-0.022089i   Abs err 1.3531e-16
exceeds tol 0 by 1e-16
   (4,3)     0.34398-0.24846i 0.34398-0.24846i   Abs err 1.1322e-15 exceeds
tol 0 by 1e-15
   (5,3)     -0.26944-0.14689i -0.26944-0.14689i   Abs err 4.4496e-16 exceeds
tol 0 by 4e-16
   (6,3)     -0.26338+0.11476i -0.26338+0.11476i   Abs err 2.0015e-16 exceeds
tol 0 by 2e-16
   (7,3)     -0.24263+0.071979i -0.24263+0.071979i   Abs err 5.6201e-16
exceeds tol 0 by 6e-16
   (8,3)     0.21314+0.046626i 0.21314+0.046626i   Abs err 6.9389e-17 exceeds
tol 0 by 7e-17
   (9,3)     -0.33649+0.2332i -0.33649+0.2332i   Abs err 8.6933e-16 exceeds
tol 0 by 9e-16
   (10,3)    0.26938+0.12755i 0.26938+0.12755i   Abs err 2.2204e-16 exceeds
tol 0 by 2e-16
   (1,4)     0.27472+0.12463i 0.27472+0.12463i   Abs err 5.4672e-16 exceeds
tol 0 by 5e-16
   (2,4)     0.24391+0.08352i 0.24391+0.08352i   Abs err 4.003e-16 exceeds tol
0 by 4e-16
   (3,4)     -0.23805+0.022089i -0.23805+0.022089i   Abs err 1.3531e-16
exceeds tol 0 by 1e-16
   (4,4)     0.34398+0.24846i 0.34398+0.24846i   Abs err 1.1322e-15 exceeds
tol 0 by 1e-15
   (5,4)     -0.26944+0.14689i -0.26944+0.14689i   Abs err 4.4496e-16 exceeds
tol 0 by 4e-16
   (6,4)     -0.26338-0.11476i -0.26338-0.11476i   Abs err 2.0015e-16 exceeds
tol 0 by 2e-16
   (7,4)     -0.24263-0.071979i -0.24263-0.071979i   Abs err 5.6201e-16
exceeds tol 0 by 6e-16
   (8,4)     0.21314-0.046626i 0.21314-0.046626i   Abs err 6.9389e-17 exceeds
tol 0 by 7e-17
   (9,4)     -0.33649-0.2332i -0.33649-0.2332i   Abs err 8.6933e-16 exceeds
tol 0 by 9e-16
   (10,4)    0.26938-0.12755i 0.26938-0.12755i   Abs err 2.2204e-16 exceeds
tol 0 by 2e-16


and reference blas:


 LD_PRELOAD=/usr/lib64/libblas.so ./run-octave -f -q
octave:1> eigs_fail_refblas
ans = ba8227df92ae
error: ASSERT errors for:  assert (Evector,Evector_f)

  Location  |  Observed  |  Expected  |  Reason
   (1,1)     -0.071017+0.074565i -0.071017+0.074566i   Abs err 3.3524e-07
exceeds tol 0 by 3e-07
   (2,1)     -0.17565-0.2761i -0.17565-0.2761i   Abs err 1.0653e-06 exceeds
tol 0 by 1e-06
   (3,1)     -0.0067953+0.006023i -0.0067953+0.006023i   Abs err 2.9562e-08
exceeds tol 0 by 3e-08
   (4,1)     -0.12045-0.13675i -0.12045-0.13675i   Abs err 5.9328e-07 exceeds
tol 0 by 6e-07
   (5,1)     -0.10914+0.12296i -0.10914+0.12296i   Abs err 5.3526e-07 exceeds
tol 0 by 5e-07
   (6,1)     0.045553+0.13595i 0.045553+0.13595i   Abs err 4.6679e-07 exceeds
tol 0 by 5e-07
   (7,1)     -0.14379-0.1448i -0.14379-0.1448i   Abs err 6.6435e-07 exceeds
tol 0 by 7e-07
   (8,1)     -0.035338+0.12118i -0.035337+0.12118i   Abs err 4.1094e-07
exceeds tol 0 by 4e-07
   (9,1)     -0.16312-0.09351i -0.16312-0.09351i   Abs err 6.1211e-07 exceeds
tol 0 by 6e-07
   (10,1)    -0.035689+0.15288i -0.035688+0.15288i   Abs err 5.111e-07 exceeds
tol 0 by 5e-07
   (1,2)     -0.015924+0.001396i 0.004222+0.015417i   Abs err 0.024545 exceeds
tol 0 by 0.02
   (2,2)     -0.002571+0.16368i 0.1615-0.026751i   Abs err 0.25137 exceeds tol
0 by 0.3
   (3,2)     0.11059+0.17509i 0.15248-0.14012i   Abs err 0.31798 exceeds tol 0
by 0.3
   (4,2)     -0.11432+0.079543i 0.09871+0.09825i   Abs err 0.21385 exceeds tol
0 by 0.2
   (5,2)     0.17136-0.05308i -0.082878-0.1591i   Abs err 0.27546 exceeds tol
0 by 0.3
   (6,2)     0.025526-0.088726i -0.091861-0.0092429i   Abs err 0.14177 exceeds
tol 0 by 0.1
   (7,2)     0.061377+0.11307i 0.10027-0.080615i   Abs err 0.19756 exceeds tol
0 by 0.2
   (8,2)     -0.057129-0.11445i -0.10239+0.076681i   Abs err 0.19642 exceeds
tol 0 by 0.2
   (9,2)     0.0041094+0.23767i 0.2331-0.046559i   Abs err 0.365 exceeds tol 0
by 0.4
   (10,2)    0.03398-0.20873i -0.21144+0.003907i   Abs err 0.32472 exceeds tol
0 by 0.3
   (1,3)     0.12455-0.0879i 0.097869-0.11688i   Abs err 0.039388 exceeds tol
0 by 0.04
   (2,3)     -0.043022+0.12201i -0.010325+0.12896i   Abs err 0.033428 exceeds
tol 0 by 0.03
   (3,3)     -0.019835+0.042526i -0.008277+0.046188i   Abs err 0.012124
exceeds tol 0 by 0.01
   (4,3)     0.010388+0.082815i 0.031259+0.077389i   Abs err 0.021565 exceeds
tol 0 by 0.02
   (5,3)     0.10949-0.048964i 0.093288-0.075382i   Abs err 0.030989 exceeds
tol 0 by 0.03
   (6,3)     -0.016475-0.15047i -0.054477-0.14122i   Abs err 0.03911 exceeds
tol 0 by 0.04
   (7,3)     0.0016431+0.22524i 0.059296+0.2173i   Abs err 0.058197 exceeds
tol 0 by 0.06
   (8,3)     0.015394+0.011785i 0.017899+0.0074471i   Abs err 0.0050091
exceeds tol 0 by 0.005
   (9,3)     0.028258+0.12216i 0.058612+0.11084i   Abs err 0.032396 exceeds
tol 0 by 0.03
   (10,3)    -0.040279-0.12763i -0.071634-0.11305i   Abs err 0.034579 exceeds
tol 0 by 0.03
   (1,4)     0.016914+0.18931i 0.016914+0.18931i   Abs err 2.163e-15 exceeds
tol 0 by 2e-15
   (2,4)     -0.072252+0.010435i -0.072252+0.010435i   Abs err 9.3353e-16
exceeds tol 0 by 9e-16
   (3,4)     -0.12413-0.15762i -0.12413-0.15762i   Abs err 2.2996e-15 exceeds
tol 0 by 2e-15
   (4,4)     -0.047428+0.013488i -0.047428+0.013488i   Abs err 6.2107e-16
exceeds tol 0 by 6e-16
   (5,4)     0.0093526+0.073147i 0.0093526+0.073147i   Abs err 8.0083e-16
exceeds tol 0 by 8e-16
   (6,4)     -0.051259+0.049822i -0.051259+0.049822i   Abs err 7.7595e-16
exceeds tol 0 by 8e-16
   (7,4)     -0.036039-0.06536i -0.036039-0.06536i   Abs err 7.5923e-16
exceeds tol 0 by 8e-16
   (8,4)     0.066386+0.064302i 0.066386+0.064302i   Abs err 1.1089e-15
exceeds tol 0 by 1e-15
   (9,4)     0.018931+0.0059828i 0.018931+0.0059828i   Abs err 2.541e-16
exceeds tol 0 by 3e-16
   (10,4)    0.0699+0.15279i 0.0699+0.15279i   Abs err 1.9502e-15 exceeds tol
0 by 2e-15


Dmitri.
-- 


(file #44644, file #44645)
    _______________________________________________________

Additional Item Attachment:

File name: eigs_fail_refblas.m            Size:0 KB
File name: eigs_fail_openblas.m           Size:0 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?54399>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/