[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #54399] test eigs fail
From: |
Dmitri A. Sergatskov |
Subject: |
[Octave-bug-tracker] [bug #54399] test eigs fail |
Date: |
Sat, 28 Jul 2018 11:28:40 -0400 (EDT) |
User-agent: |
Mozilla/5.0 (X11; Fedora; Linux x86_64; rv:61.0) Gecko/20100101 Firefox/61.0 |
Follow-up Comment #4, bug #54399 (project octave):
Hmm... This is only with openblas. With reference blas it still fails, but no
swapping:
LD_PRELOAD=/usr/lib64/libblas.so ./run-octave -f -q
octave:1> test eigs
***** testif HAVE_ARPACK
A = rand (10) + 1i * rand (10);
B = rand (10) + 1i * rand (10);
B = B * B';
opts.v0 = (1:10)';
[Evector, Evalues] = eigs (A, B, 4, "SM", opts);
Afun = @(x) A \ x;
opts.isreal = false;
[Evector_f, Evalues_f] = eigs (Afun, 10, B, 4, "SM", opts);
assert (Evector, Evector_f);
assert (Evalues, Evalues_f);
!!!!! test failed
ASSERT errors for: assert (Evector,Evector_f)
Location | Observed | Expected | Reason
(1,1) 0.44932+0.36053i 0.062793+0.57265i Abs err 0.44091 exceeds tol
0 by 0.4
(2,1) -0.2014+0.061636i -0.186-0.098823i Abs err 0.1612 exceeds tol 0
by 0.2
(3,1) 0.11407+0.016089i 0.069286+0.092035i Abs err 0.088168 exceeds
tol 0 by 0.09
(4,1) 0.0063383+0.01503i -0.0061453+0.015109i Abs err 0.012484
exceeds tol 0 by 0.01
(5,1) -0.21947-0.20332i -0.011424-0.29896i Abs err 0.22898 exceeds
tol 0 by 0.2
(6,1) 0.078061-0.03096i 0.077091+0.033304i Abs err 0.064272 exceeds
tol 0 by 0.06
(7,1) 0.21379-0.029145i 0.17178+0.13056i Abs err 0.16514 exceeds tol
0 by 0.2
(8,1) 0.026375-0.17993i 0.14588-0.10859i Abs err 0.13918 exceeds tol
0 by 0.1
(9,1) -0.26341-0.039268i -0.15849-0.21402i Abs err 0.20383 exceeds
tol 0 by 0.2
(10,1) 0.022875-0.31838i 0.2413-0.20896i Abs err 0.2443 exceeds tol 0
by 0.2
(1,2) -0.097689+0.26252i -0.082802+0.26758i Abs err 0.015726 exceeds
tol 0 by 0.02
(2,2) -0.28431+0.16468i -0.27462+0.18038i Abs err 0.018447 exceeds
tol 0 by 0.02
(3,2) -0.059822+0.18621i -0.049277+0.18928i Abs err 0.010981 exceeds
tol 0 by 0.01
(4,2) 0.19295-0.019498i 0.19155-0.030296i Abs err 0.010888 exceeds
tol 0 by 0.01
(5,2) 0.29759+0.010432i 0.29771-0.0062855i Abs err 0.016718 exceeds
tol 0 by 0.02
(6,2) -0.28995+0.021831i -0.28827+0.038069i Abs err 0.016325 exceeds
tol 0 by 0.02
(7,2) -0.091281+0.2275i -0.078369+0.23227i Abs err 0.013763 exceeds
tol 0 by 0.01
(8,2) -0.043801-0.23869i -0.057128-0.23586i Abs err 0.013625 exceeds
tol 0 by 0.01
(9,2) 0.16326-0.35502i 0.14308-0.36362i Abs err 0.021939 exceeds tol
0 by 0.02
(10,2) 0.31081-0.18609i 0.29988-0.20324i Abs err 0.020339 exceeds tol
0 by 0.02
(1,3) -0.18758-0.014342i -0.071058-0.17419i Abs err 0.1978 exceeds
tol 0 by 0.2
(2,3) 0.026441+0.1237i -0.098819+0.078972i Abs err 0.13301 exceeds
tol 0 by 0.1
(3,3) 0.0013691+0.078902i -0.06996+0.036511i Abs err 0.082975 exceeds
tol 0 by 0.08
(4,3) -0.044882+0.037667i -0.053762-0.023299i Abs err 0.061609
exceeds tol 0 by 0.06
(5,3) 0.015823-0.01741i 0.022648+0.0063664i Abs err 0.024736 exceeds
tol 0 by 0.02
(6,3) -0.17024-0.12416i 0.03492-0.20779i Abs err 0.22155 exceeds tol
0 by 0.2
(7,3) -0.18329+0.01315i -0.09373-0.15806i Abs err 0.19321 exceeds tol
0 by 0.2
(8,3) 0.042603-0.059632i 0.072389+0.011437i Abs err 0.077058 exceeds
tol 0 by 0.08
(9,3) 0.10692-0.077796i 0.1174+0.060844i Abs err 0.13904 exceeds tol
0 by 0.1
(10,3) 0.053564+0.089147i -0.055781+0.087777i Abs err 0.10935 exceeds
tol 0 by 0.1
(1,4) 0.10283-0.28385i 0.10283-0.28385i Abs err 1.5927e-15 exceeds
tol 0 by 2e-15
(2,4) 0.048707+0.26316i 0.048707+0.26316i Abs err 1.4289e-15 exceeds
tol 0 by 1e-15
(3,4) 0.13125+0.10726i 0.13125+0.10726i Abs err 9.4532e-16 exceeds
tol 0 by 9e-16
(4,4) 0.20488+0.03423i 0.20488+0.03423i Abs err 1.1271e-15 exceeds
tol 0 by 1e-15
(5,4) -0.10285+0.11154i -0.10285+0.11154i Abs err 8.2803e-16 exceeds
tol 0 by 8e-16
(6,4) -0.26098-0.37498i -0.26098-0.37498i Abs err 2.2651e-15 exceeds
tol 0 by 2e-15
(7,4) -0.35516+0.012622i -0.35516+0.012622i Abs err 1.8629e-15
exceeds tol 0 by 2e-15
(8,4) -0.029845-0.049972i -0.029845-0.049972i Abs err 3.1101e-16
exceeds tol 0 by 3e-16
(9,4) 0.45968+0.13373i 0.45968+0.13373i Abs err 2.5853e-15 exceeds
tol 0 by 3e-15
(10,4) -0.10302+0.075966i -0.10302+0.075966i Abs err 6.5506e-16
exceeds tol 0 by 7e-16
Dmitri.
--
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?54399>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
- [Octave-bug-tracker] [bug #54399] test eigs fail, Dmitri A. Sergatskov, 2018/07/27
- [Octave-bug-tracker] [bug #54399] test eigs fail, Rik, 2018/07/28
- [Octave-bug-tracker] [bug #54399] test eigs fail, Dmitri A. Sergatskov, 2018/07/28
- [Octave-bug-tracker] [bug #54399] test eigs fail, Marco Caliari, 2018/07/28
- [Octave-bug-tracker] [bug #54399] test eigs fail,
Dmitri A. Sergatskov <=
- [Octave-bug-tracker] [bug #54399] test eigs fail, Dmitri A. Sergatskov, 2018/07/28
- [Octave-bug-tracker] [bug #54399] test eigs fail, Marco Caliari, 2018/07/28
- [Octave-bug-tracker] [bug #54399] test eigs fail, Dmitri A. Sergatskov, 2018/07/28
- [Octave-bug-tracker] [bug #54399] test eigs fail, Marco Caliari, 2018/07/29
- [Octave-bug-tracker] [bug #54399] test eigs fail, Dmitri A. Sergatskov, 2018/07/29
- [Octave-bug-tracker] [bug #54399] test eigs fail, Marco Caliari, 2018/07/30
- [Octave-bug-tracker] [bug #54399] test eigs fail, Dmitri A. Sergatskov, 2018/07/30
- [Octave-bug-tracker] [bug #54399] test eigs fail, Marco Caliari, 2018/07/30
- [Octave-bug-tracker] [bug #54399] test eigs fail, Rik, 2018/07/30
- [Octave-bug-tracker] [bug #54399] test eigs fail, Marco Caliari, 2018/07/30