[Octave-bug-tracker] [bug #65413] Missing function: pagesvd

[anonymous]

Follow-up Comment #9, bug #65413 (group octave):

OpenBLAS implementations are internally parallelized, but it doesn't help much
with small matrices. Turning off the internal OpenBLAS parallelization -
openblas_set_num_threads(1) - and using explicit OpenMP parallel loops means
each SVD is done on a single core.

For the use cases it is meant for, it's a big win.

>> A=rand(50,50,100000);
>> tic;[U S V]=pagesvd(A);toc
Elapsed time is 4.410267 seconds.

>> 
>> A=rand(50,50);
>> tic;for k=1:100000;[U S V]=svd(A);end;toc
Elapsed time is 64.020166 seconds.




    _______________________________________________________

Reply to this item at:

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

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