[Octave-bug-tracker] [bug #65178] [octave forge] (statistics) pdist2 BIST failures

[Dmitri A. Sergatskov]

Follow-up Comment #3, bug#65178 (group octave):

I have 3 failures on aarch64 (raspberry pi 4) Ubuntu 23.10:

 ***** test
 d = [1.2247, 2.4495, 4.4159; 2.1213, 4.2426, 6.1237];
 i = [2, 2, 2; 1, 4, 4];
 [D, I] = pdist2 (y, x, "mahalanobis", "smallest", 2);
 assert ({D, I}, {d, i}, 1e-4);
!!!!! test failed
ASSERT errors for:  assert (cond {i},expected {i},tol)

  Location  |  Observed  |  Expected  |  Reason
   (2,2)          1            4         Abs err 3 exceeds tol 0.0001 by 3
shared variables 
    x =

       1   1   1
       2   2   2
       3   3   3

    y =

       0   0   0
       1   2   3
       0   2   4
       4   7   1

    xx =

       1   2   3
       4   5   6
       7   8   9
       3   2   1

***** test
 d = [3.2404, 4.8990, 6.8191; 2.7386, 4.2426, 6.3640];
 i = [3, 3, 3; 4, 1, 1];
 [D, I] = pdist2 (y, x, "mahalanobis", "largest", 2);
 assert ({D, I}, {d, i}, 1e-4);
!!!!! test failed
ASSERT errors for:  assert (cond {i},expected {i},tol)

  Location  |  Observed  |  Expected  |  Reason
   (2,2)          4            1         Abs err 3 exceeds tol 0.0001 by 3
shared variables 
    x =

       1   1   1
       2   2   2
       3   3   3

    y =

       0   0   0
       1   2   3
       0   2   4
       4   7   1

    xx =

       1   2   3
       4   5   6
       7   8   9
       3   2   1

***** test
 yy = [1 2 3;5 6 7;9 5 1];
 d = [0, 8.4853, 18.0416; 2.4495, 10.0995, 19.4808; ...
      2.4495, 10.6771, 19.7104; 2.4495, 10.6771, 20.4573];
 i = [2, 2, 2; 1, 4, 4; 4, 1, 1; 3, 3, 3];
 [D, I] = pdist2 (y, yy, "mahalanobis", "smallest", 4);
 assert ({D, I}, {d, i}, 1e-4);
!!!!! test failed
ASSERT errors for:  assert (cond {i},expected {i},tol)

  Location  |  Observed  |  Expected  |  Reason
   (3,1)          3            4         Abs err 1 exceeds tol 0.0001 by 1
   (4,1)          4            3         Abs err 1 exceeds tol 0.0001 by 1
shared variables 
    x =

       1   1   1
       2   2   2
       3   3   3

    y =

       0   0   0
       1   2   3
       0   2   4
       4   7   1

    xx =

       1   2   3
       4   5   6
       7   8   9
       3   2   1



The test passes on the same Ubuntu 23.10 but on x86_64 (Intel(R) Core(TM)2 Duo
CPU     T9900 ).

So my guess is that some of the data is marginal and the results depends on FP
unit in use (core-II uses SSE and my other systems have AVX).

Dmitri.
-- 



    _______________________________________________________

Reply to this item at:

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

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