octave-bug-tracker
[Top][All Lists]
Advanced

[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/




reply via email to

[Prev in Thread] Current Thread [Next in Thread]