[Octave-bug-tracker] [bug #45387] Matrix product treats NaN times zero as zero

[Marco Caliari]

Follow-up Comment #6, bug #45387 (project octave):

It is reference BLAS. If you look at line 281 here

http://www.netlib.org/lapack/explore-html/dc/da8/dgemv_8f_source.html

you see that if the entry of the vector is zero, the multiplication is
skipped. I tried the attached fortran file and I confirm the correct behavior
for openblas and ATLAS and the wrong for reference.

I think that something should be written, otherwise a Linux user may have the
impression that 1) octave cannot properly handle NaN and 2) (much more
important IMHO) octave is much slower than Matlab.

Marco

(file #34369)
    _______________________________________________________

Additional Item Attachment:

File name: main_dgemv.f                   Size:0 KB


    _______________________________________________________

Reply to this item at:

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

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/