Search found 32 matches

by hsahasra
Tue Apr 22, 2014 5:38 pm
Forum: User discussion
Topic: zgemm for matrices that don't fit
Replies: 2
Views: 2042

zgemm for matrices that don't fit

Hi, I need to multiply matrices of size 12000x12000. 3 of these matrices don't fit in a GPU so I can't use cuBLAS. I'm not able to use magmablas_zgemm either as it works on matrices in GPU memory. Is there a method in MAGMA or some other library which takes matrices in RAM larger than GPU memory and...
by hsahasra
Tue Apr 22, 2014 5:34 pm
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8629

Re: Unable to run zgetri_gpu

Thanks Mark. I'll start from this code and see what's wrong with mine. There must be some small bug somewhere. Will post when I find something.
by hsahasra
Wed Apr 16, 2014 9:25 pm
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8629

Re: Unable to run zgetri_gpu

Hi Mark,

The matrices which are inverting properly are of sizes 80x80, 320x320 and 1280x1280. The inversion fails when I try to do it for a matrix of size 720x720. Is this problem related to the number of rows/columns being divisible by 16/32?

Harshad
by hsahasra
Wed Apr 16, 2014 7:11 pm
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8629

Re: Unable to run zgetri_gpu

Sorry, I didn't see that you said invert and multiply x=A^{-1}*b. In my code I'm actually inverting the matrix and not solving a linear equation. I can also implicitly invert a matrix using zgesv, by setting the RHS to an identity matrix. So my question was, do you know of a way of creating an ident...
by hsahasra
Tue Apr 15, 2014 1:02 pm
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8629

Re: Unable to run zgetri_gpu

Hi Mark, Yes, I'm multiplying by sizeof(cuDoubleComplex). I tried to use zgesv, but I get the same error. How do you usually set the rhs for zgesv while inverting a matrix? I don't understand why the code works for 80x80 matrices and fails for larger matrices. Here's what I currently do: //Memory al...
by hsahasra
Tue Apr 15, 2014 1:18 am
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8629

Re: Unable to run zgetri_gpu

The only difference in my code is that I'm using cudaMalloc instead of magma_malloc. Should it matter?
by hsahasra
Tue Apr 15, 2014 12:59 am
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8629

Unable to run zgetri_gpu

Hi, I'm having some weird problems using zgetrf+zgetri while inverting a matrix. I'm doing exactly what is being done in the testing_zgetri_gpu example. I get the same error if I use magma_zgesv_gpu with an identity RHS. I'm able to run my code and get correct answers for small (80x80) matrices. But...
by hsahasra
Thu Apr 10, 2014 2:30 pm
Forum: User discussion
Topic: how to avoid copying to/from the GPU
Replies: 4
Views: 2231

Re: how to avoid copying to/from the GPU

Yes, transfer the starting data once then successively use magma_*_gpu calls. _gpu calls take GPU data pointers. You can also have a look at CUBLAS and CUSPARSE libraries.

http://docs.nvidia.com/cuda/cublas/index.html
http://docs.nvidia.com/cuda/cusparse/
by hsahasra
Wed Oct 09, 2013 6:41 pm
Forum: User discussion
Topic: Running MAGMA across several GPUs on several nodes
Replies: 2
Views: 2703

Re: Running MAGMA across several GPUs on several nodes

We need to call LU factorization, symmetric or nonsymmetric eigenvalues and mat-mat multiplication. I am working on adding an interface for MAGMA in PETSc library. I have completed work on LU factorization. Now I'm working on adding xgetrf_gpu+magma_getrs_gpu and xgemm. I can't find any tutorial for...
by hsahasra
Tue Aug 13, 2013 1:33 pm
Forum: User discussion
Topic: std::complex<double> to cuComplexDouble
Replies: 2
Views: 6772

Re: std::complex<double> to cuComplexDouble

Thanks, this works!