[Octave-bug-tracker] [bug #58926] Octave gives wrong results with intel-mkl when diagonalizing large matrices

[Dmitri A. Sergatskov]

Follow-up Comment #30, bug #58926 (project octave):

Speaking of benchmarks, on Ryzen i got best results whith openMP
interface and limiting number of threds to the number of CPU cores
(16 in my case):

 OMP_NUM_THREADS=16 LD_PRELOAD=./libopenblas.so octave -q -f
octave:1> ii=1:500;
octave:2> c = sin(ii' + ii.^2);
octave:3> whos c
Variables in the current scope:

   Attr Name        Size                     Bytes  Class
   ==== ====        ====                     =====  ===== 
        c         500x500                  2000000  double

Total is 250000 elements using 2000000 bytes

octave:4> tic; g = eig(c); toc
Elapsed time is 0.118513 seconds.
octave:5> max(real(g))
ans =  16.915

Dmitri.
-- 


    _______________________________________________________

Reply to this item at:

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

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