[Octave-bug-tracker] [bug #62817] Regression in interpn, likely harmless

[Arun Giridhar]

Follow-up Comment #7, bug #62817 (project octave):

So! I ran this test a statistically large number of times and found it's not
100% repeatable!
Test:

for i in `seq 1 100`; do octave --eval "test interpn"; done >/tmp/tmp


Results: it fails about half the time, and within the failures about 1/3 is
for the imag part and 2/3 for the real part.

$ grep PASSES tmp | sort | uniq -c
     48 PASSES 18 out of 19 tests
     52 PASSES 19 out of 19 tests

$ grep ASSERT tmp | sort | uniq -c
     15 ASSERT errors for:  assert (imag (vi_complex),vi_imag)
     33 ASSERT errors for:  assert (real (vi_complex),vi_real)


The numerical errors are all of the order of eps or less:

$ cat tmp | grep exceeds | awk '{print $(NF-5)}' | sort -n | uniq -c
     34 1.1102e-16
      2 2.2204e-16
      1 2.7756e-17
     14 5.5511e-17


As you said, it's better solved by setting the seed (or not using random
numbers) than only by increasing the tolerance.


    _______________________________________________________

Reply to this item at:

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

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