Search found 32 matches

by hsahasra
Fri Jun 27, 2014 1:14 pm
Forum: User discussion
Topic: creating libmagma.a in 'lib'.
Replies: 5
Views: 3140

Re: creating libmagma.a in 'lib'.

You need to add -fPIC in all the compiler flags in make.inc file. For example: OPTS = -fPIC -O3 -DADD_ -Wall -openmp -DMAGMA_WITH_MKL -DMAGMA_SETAFFINITY F77OPTS = -fPIC -O3 -DADD_ -warn all FOPTS = -fPIC -O3 -DADD_ -warn all NVOPTS = -O3 -DADD_ -Xcompiler "-fno-strict-aliasing -fPIC"
by hsahasra
Tue Jun 24, 2014 10:02 am
Forum: User discussion
Topic: Can I call magma functions inside CUDA GPU kernel
Replies: 1
Views: 2290

Re: Can I call magma functions inside CUDA GPU kernel

You can call kernels from within kernels by using dynamic parallelism: http://developer.download.nvidia.com/as ... n_CUDA.pdf

Should work for MAGMA kernels too
by hsahasra
Fri Jun 20, 2014 10:56 am
Forum: User discussion
Topic: Possible bug in ztranspose-v2.cu
Replies: 3
Views: 2811

Re: Possible bug in ztranspose-v2.cu

Just doing this probably fixes the error, not sure if it still gives the correct result. if (ibx+inx+ZSIZE_1SHARED < m) { if (t2 < n) { sA[iny+0][inx] = A[0*lda]; if (t2+ 8 < n) { sA[iny+8][inx] = A[8*lda]; if (t2 + 16 < n) { sA[iny+16][inx] = A[16*lda]; if (t2 + 24 < n) { sA[iny+24][inx] = A[24*lda...
by hsahasra
Fri Jun 20, 2014 10:03 am
Forum: User discussion
Topic: Possible bug in ztranspose-v2.cu
Replies: 3
Views: 2811

Possible bug in ztranspose-v2.cu

Hi, I'm having a problem with zgetri_gpu or zgesv_gpu when they perform LU factorization at zgetrf_gpu.cpp:167. I checked that all the pointers I pass to zgesv/zgetri are fine. What is more weird is that there is no such error when only the tests for zgesv/zgetri are run. This only happens when othe...
by hsahasra
Sun May 04, 2014 10:07 pm
Forum: User discussion
Topic: MAGMA with CUDA MPS
Replies: 1
Views: 2480

MAGMA with CUDA MPS

Hi, I'm inverting matrices using MAGMA (zgetri_gpu). To use the GPU and host properly, I want to run multiple processes accessing a single GPU on the BlueWaters supercomputer. I'm using the CUDA MPS for sharing the GPU. When a single process uses a GPU, the program runs properly. When multiple proce...
by hsahasra
Sun May 04, 2014 11:31 am
Forum: User discussion
Topic: Why does magma_dgemm silently fail?
Replies: 3
Views: 3007

Re: Why does magma_dgemm silently fail?

You can correctly check for errors when the stream is synchronized. For example: cudaStream_t stream; . . cublasStatus_t cublas_error = cublasDgemm(stream, ....); <check cublas_error> cudaError_t cuda_error = cudaStreamSynchronize(stream); <check cuda_error> for MAGMA, just set the blas kernel strea...
by hsahasra
Sun May 04, 2014 11:25 am
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8630

Re: Unable to run zgetri_gpu

Hi Mark, I think I'm closer to the bug in my code and possibly a bug in MAGMA. I'm using magma_zgetrf_gpu + magma_zgetri_gpu in several iterations on different matrices. What I found was that if I copy those different matrices to a fixed memory on GPU and then invert the matrix, the inversion runs p...
by hsahasra
Wed Apr 23, 2014 2:38 pm
Forum: User discussion
Topic: Parallel MAGMA function call
Replies: 2
Views: 2278

Re: Parallel MAGMA function call

Can you specify more details about the next magma_sgemm call? Why are you using gemm for matrix vector multiplication?

What do you mean by extracting vectors? You could just use the starting pointer eg. v1 = v[0], v2 = v[10], .... and length = 10 to specify each vector.
by hsahasra
Tue Apr 22, 2014 7:59 pm
Forum: User discussion
Topic: zgemm for matrices that don't fit
Replies: 2
Views: 2042

Re: zgemm for matrices that don't fit

Thanks!
by hsahasra
Tue Apr 22, 2014 6:08 pm
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8630

Re: Unable to run zgetri_gpu

test_zgetri.c is working fine, which means there is some bug in my code.