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

[Dmitri A. Sergatskov]

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

                 Summary: test eigs fail
                 Project: GNU Octave
            Submitted by: dasergatskov
            Submitted on: Sat 28 Jul 2018 02:38:40 AM UTC
                Category: Test Suite
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: None
                  Status: None
             Assigned to: None
         Originator Name: 
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: dev
        Operating System: Any

    _______________________________________________________

Details:

It appears that after the last check-in:

hg head
changeset:   25697:91f0416c2ba7
bookmark:    @
tag:         tip
parent:      25696:3b6691ff0f59
parent:      25695:cfc32b131755
user:        John W. Eaton <address@hidden>
date:        Fri Jul 27 18:15:17 2018 -0400
summary:     maint: Merge stable to default.

changeset:   25695:cfc32b131755
branch:      stable
user:        John W. Eaton <address@hidden>
date:        Fri Jul 27 13:21:50 2018 -0400
summary:     avoid questionable cast

test eigs fil with:



octave:1> test eigs
***** testif HAVE_ARPACK
 A = rand (10);
 B = rand (10);
 B = B * B';
 opts.v0 = (1:10)';
 [Evector, Evalues] = eigs (A, B, 4, "SM", opts);
 Afun = @(x) A \ x;
 [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.039202+0i  0.039202+0i    Abs err 8.3267e-17 exceeds tol 0 by
8e-17
   (2,1)     -0.47728+0i  -0.47728+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (3,1)     -0.26688+0i  -0.26688+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (4,1)      0.24262+0i   0.24262+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (5,1)     -0.30222+0i  -0.30222+0i    Abs err 1.6653e-16 exceeds tol 0 by
2e-16
   (6,1)      0.3733+0i    0.3733+0i     Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (7,1)      0.20615+0i   0.20615+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (8,1)      0.24472+0i   0.24472+0i    Abs err 2.7756e-16 exceeds tol 0 by
3e-16
   (9,1)     -0.032564+0i -0.032564+0i   Abs err 8.3267e-17 exceeds tol 0 by
8e-17
   (10,1)    0.067874+0i  0.067874+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (1,2)     -0.14068+0i  -0.14068+0i    Abs err 6.9389e-16 exceeds tol 0 by
7e-16
   (2,2)     -0.47604+0i  -0.47604+0i    Abs err 8.3267e-16 exceeds tol 0 by
8e-16
   (3,2)      0.12838+0i   0.12838+0i    Abs err 1.1102e-16 exceeds tol 0 by
1e-16
   (4,2)     -0.10735+0i  -0.10735+0i    Abs err 5.5511e-17 exceeds tol 0 by
6e-17
   (5,2)     -0.0013885+0i -0.0013885+0i   Abs err 1.1178e-15 exceeds tol 0 by
1e-15
   (6,2)     0.098606+0i  0.098606+0i    Abs err 4.1633e-16 exceeds tol 0 by
4e-16
   (7,2)      0.20902+0i   0.20902+0i    Abs err 5.8287e-16 exceeds tol 0 by
6e-16
   (8,2)      0.23101+0i   0.23101+0i    Abs err 1.3323e-15 exceeds tol 0 by
1e-15
   (9,2)     0.052558+0i  0.052558+0i    Abs err 9.4369e-16 exceeds tol 0 by
9e-16
   (10,2)    -0.37364+0i  -0.37364+0i    Abs err 7.7716e-16 exceeds tol 0 by
8e-16
   (1,3)     -0.31312+0i  -0.31312+0i    Abs err 4.4409e-16 exceeds tol 0 by
4e-16
   (2,3)     -0.55234+0i  -0.55234+0i    Abs err 2.2204e-16 exceeds tol 0 by
2e-16
   (3,3)     0.082344+0i  0.082344+0i    Abs err 4.5797e-16 exceeds tol 0 by
5e-16
   (4,3)     -0.080656+0i -0.080656+0i   Abs err 1.3878e-16 exceeds tol 0 by
1e-16
   (5,3)      0.3594+0i    0.3594+0i     Abs err 3.8858e-16 exceeds tol 0 by
4e-16
   (6,3)      0.22037+0i   0.22037+0i    Abs err 6.1062e-16 exceeds tol 0 by
6e-16
   (7,3)      0.35982+0i   0.35982+0i    Abs err 3.8858e-16 exceeds tol 0 by
4e-16
   (8,3)      0.27437+0i   0.27437+0i    Abs err 4.996e-16 exceeds tol 0 by
5e-16
   (9,3)      0.39752+0i   0.39752+0i    Abs err 1.0547e-15 exceeds tol 0 by
1e-15
   (10,3)    -0.50035+0i  -0.50035+0i    Abs err 1.1102e-16 exceeds tol 0 by
1e-16
   (1,4)     -0.084833+0.40857i 0.40857+0.084833i   Abs err 0.59012 exceeds
tol 0 by 0.6
   (2,4)     0.14708-0.48518i -0.48518-0.14708i   Abs err 0.71698 exceeds tol
0 by 0.7
   (3,4)     -0.29952-0.10353i -0.10353+0.29952i   Abs err 0.44817 exceeds tol
0 by 0.4
   (4,4)     -0.028404+0.22032i 0.22032+0.028404i   Abs err 0.31416 exceeds
tol 0 by 0.3
   (5,4)     -0.059085-0.18242i -0.18242+0.059085i   Abs err 0.27118 exceeds
tol 0 by 0.3
   (6,4)     -0.13868+0.15719i 0.15719+0.13868i   Abs err 0.29645 exceeds tol
0 by 0.3
   (7,4)     -0.32197+0.42888i 0.42888+0.32197i   Abs err 0.75843 exceeds tol
0 by 0.8
   (8,4)     0.11036-0.7379i -0.7379-0.11036i   Abs err 1.0552 exceeds tol 0
by 1
   (9,4)     0.062964+0.40427i 0.40427-0.062964i   Abs err 0.57862 exceeds tol
0 by 0.6
   (10,4)    0.60332-0.30705i -0.30705-0.60332i   Abs err 0.95737 exceeds tol
0 by 1
octave:2> 



I think assert needs some tolerance.

Dmitri.
-- 





    _______________________________________________________

Reply to this item at:

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

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