DGESV doesn't use magma_dgetrs?

Open discussion for MAGMA library (Matrix Algebra on GPU and Multicore Architectures)

DGESV doesn't use magma_dgetrs?

Postby jah87 » Wed May 30, 2012 10:30 am

I was looking through the source for the magma_dgesv routine and noticed this:

Code: Select all
magma_dgetrf( n, n, A, lda, ipiv, info );
   94     if ( *info == 0 ) {
   95         lapackf77_dgetrs( MagmaNoTransStr, &n, &nrhs, A, &lda, ipiv, B, &ldb, info );
   96     }

There is not a magma_dgetrs routine being implemented?
Posts: 21
Joined: Tue May 01, 2012 1:54 pm

Re: DGESV doesn't use magma_dgetrs?

Postby mgates3 » Wed May 30, 2012 3:14 pm

There is a GPU interface, magma_dgetrs_gpu, which is used in magma_dgesv_gpu. For the CPU interface, magma_dgesv, if the matrix is not already on the GPU, it can be slower to transfer it, perform getrs, and transfer the results back, since getrs does not parallelize well. We may change magma_dgesv internally to use magma_dgetrf_gpu and magma_dgetrs_gpu, but need to deal with cases when the matrix does not fit on the GPU (non-resident) or using multiple GPUs.

Posts: 782
Joined: Fri Jan 06, 2012 2:13 pm

Return to User discussion

Who is online

Users browsing this forum: No registered users and 4 guests