SVD on Magma

[jacquesdutoit]

Hi

I've had a look at the SVD code, and if I understand it correctly, the routine takes a matrix resident on the host and returns a result on the host. This seems to suggest that the GPU is used to accelerate certain parts of the algorithm, but that a lot of the work is done on the host.

Are there any plans to make a version of SVD which takes a matrix resident on the GPU and returns a result resident on the GPU?

Regards
Jacques

[fletchjp]

Stan

Thank you for that. I have been silent for a long time pursuing other activities. I am now coming back to this again.

Best wishes

John

[fletchjp]

Just to add to the old thread some tests with Magma 1.2.1 and OpenBLAS.

Code: Select all

~/magma-1.2.1/testing$ ./testing_dgesvd
device 0: GeForce GTX 460, 1400.0 MHz clock, 2047.2 MB memory, capability 2.1

Usage:
  testing_dgesvd -M 1024 -N 1024

  N     CPU Time(s)    GPU Time(s)     ||R||_F / ||A||_F
==========================================================
 1024      12.19          12.26         2.184792e-15
 2048     151.06         149.38         2.239185e-15
 3072     529.26         523.26         2.523813e-15
 4032     1098.63         1082.79         2.940447e-15
 5184     2648.70         2615.76         3.537992e-15

Any gain here seems to be the better BLAS. it still doesn't multithread.