[Octave-bug-tracker] [bug #60552] BIST for lscov.m tolerance failure with OpenBLAS (v0.3.13–v0.3.15)

[Mike Miller]

URL:
  <https://savannah.gnu.org/bugs/?60552>

                 Summary: BIST for lscov.m tolerance failure with OpenBLAS
(v0.3.13–v0.3.15)
                 Project: GNU Octave
            Submitted by: mtmiller
            Submitted on: Fri 07 May 2021 12:17:28 PM PDT
                Category: Test Suite
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Inaccurate Result
                  Status: None
             Assigned to: None
         Originator Name: 
        Originator Email: 
             Open/Closed: Open
                 Release: dev
         Discussion Lock: Any
        Operating System: GNU/Linux

    _______________________________________________________

Details:

In Octave version 6.2.0 or version 7.0.0 development branch, with OpenBLAS
0.3.13 (on Debian native system):


>> test lscov.m
***** test
 ## Adapted from example in Matlab documentation
 x1 = [.2 .5 .6 .8 1.0 1.1]';
 x2 = [.1 .3 .4 .9 1.1 1.4]';
 X = [ones(size(x1)) x1 x2];
 y = [.17 .26 .28 .23 .27 .34]';
 [b, se_b, mse, S] = lscov(X, y);
 assert (b, [0.1203 0.3284 -0.1312]', 1E-4);
 assert (se_b, [0.0643 0.2267 0.1488]', 1E-4);
 assert (mse, 0.0015, 1E-4);
 assert (S, [0.0041 -0.0130 0.0075; -0.0130 0.0514 -0.0328; 0.0075 -0.0328
0.0221], 1E-4);
 w = [1 1 1 1 1 .1]';
 [bw, sew_b, msew] = lscov (X, y, w);
 assert (bw, [0.1046 0.4614 -0.2621]', 1E-4);
 assert (sew_b, [0.0309 0.1152 0.0814]', 1E-4);
 assert (msew, 3.4741e-004, -1E-4);
 V = .2*ones (length (x1)) + .8*diag (ones (size (x1)));
 [bg, sew_b, mseg] = lscov (X, y, V);
 assert (bg, [0.1203 0.3284 -0.1312]', 1E-4);
 assert (sew_b, [0.0672 0.2267 0.1488]', 1E-4);
 assert (mseg, 0.0019, 1E-4);
 y2 = [y 2*y];
 [b2, se_b2, mse2, S2] = lscov (X, y2);
 assert (b2, [b 2*b], 2*eps);
 assert (se_b2, [se_b 2*se_b], 2*eps);
 assert (mse2, [mse 4*mse], eps);
 assert (S2(:, :, 1), S, eps);
 assert (S2(:, :, 2), 4*S, eps);
!!!!! test failed
ASSERT errors for:  assert (S2 (:, :, 2),4 * S,eps)

  Location  |  Observed  |  Expected  |  Reason
   (2,2)       0.20555      0.20555      Abs err 2.7756e-16 exceeds tol
2.2204e-16 by 6e-17


In Octave version 6.2.0 or version 7.0.0 development branch, with OpenBLAS
0.3.15, running in Flatpak container, same test, slightly different result:


!!!!! test failed
ASSERT errors for:  assert (S2 (:, :, 2),4 * S,eps)

  Location  |  Observed  |  Expected  |  Reason
   (2,2)       0.20555      0.20555      Abs err 3.0531e-16 exceeds tol
2.2204e-16 by 8e-17





    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?60552>

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