Search found 9 matches

by mrader1248
Wed Dec 06, 2017 4:42 am
Forum: User discussion
Topic: QDWH-SVD
Replies: 2
Views: 1494

Re: QDWH-SVD

Thank you for your reply! This paper is not included in the MAGMA publications site? http://icl.cs.utk.edu/magma/pubs/index.html Is there also a preprint available somewhere? (Unfortunately, my university has no subscription for this journal.) The highlights section of your publication says: Achieve...
by mrader1248
Tue Dec 05, 2017 6:58 am
Forum: User discussion
Topic: QDWH-SVD
Replies: 2
Views: 1494

QDWH-SVD

Hello, I stumbled over the paper "A HIGH PERFORMANCE QDWH-SVD SOLVER USING HARDWARE ACCELERATORS" by D. Sukkari et al. [1], where it is shown that SVDs can be performed much faster than it is done within MAGMA by first performing a polar decomposition (via QR decompositions) and then a hermitian EVD...
by mrader1248
Mon Oct 16, 2017 4:07 am
Forum: User discussion
Topic: MAGMA + pycuda + my own CUDA kernels
Replies: 11
Views: 5076

Re: MAGMA + pycuda + my own CUDA kernels

Just as a reminder: I want to obtain both Q and R. When I use magma_zgeqrf2_gpu, I have direct access to R, but there is no matching function to restore Q: magma_zungqr and magma_zungqr2 both require A to be in host memory and magma_zungqr_gpu requires the dT array which I don't get from magma_zgeqr...
by mrader1248
Fri Oct 06, 2017 6:15 am
Forum: User discussion
Topic: MAGMA + pycuda + my own CUDA kernels
Replies: 11
Views: 5076

Re: MAGMA + pycuda + my own CUDA kernels

After some testing, I finally understood, what you meant in the beginning with "You probably want to use magma_zgeqrf2_gpu". When I increase the dimensions of my matrices, R is no longer stored in the memory containing A initally. I want to be able to restore R (and also Q). In http://icl.cs.utk.edu...
by mrader1248
Wed Oct 04, 2017 5:12 am
Forum: User discussion
Topic: MAGMA + pycuda + my own CUDA kernels
Replies: 11
Views: 5076

Re: MAGMA + pycuda + my own CUDA kernels

Thank you very much for your help! Finally everything works. :)
by mrader1248
Mon Oct 02, 2017 8:27 am
Forum: User discussion
Topic: MAGMA + pycuda + my own CUDA kernels
Replies: 11
Views: 5076

Re: MAGMA + pycuda + my own CUDA kernels

I guess, I finally managed to get it working (for details see my example below). I took the final idea from this post: https://stackoverflow.com/a/10415966/8683723 Apparently one also has to call cudaFree(0) right after magma_setdevice. import ctypes import numpy as np import sys import pycuda.autoi...
by mrader1248
Mon Oct 02, 2017 4:58 am
Forum: User discussion
Topic: MAGMA + pycuda + my own CUDA kernels
Replies: 11
Views: 5076

Re: MAGMA + pycuda + my own CUDA kernels

The answers seems to be given here: https://stackoverflow.com/a/20553411/8683723 This is how I am interpreting it: 1. import pycuda.autoinit 2. Call magma_setdevice and don't call magma_init. 3. Now you are free to do whatever you like. But now I'm wondering if MAGMA will work without calling magma_...
by mrader1248
Fri Sep 29, 2017 4:37 am
Forum: User discussion
Topic: MAGMA + pycuda + my own CUDA kernels
Replies: 11
Views: 5076

Re: MAGMA + pycuda + my own CUDA kernels

Thank you very much for your comments! Regarding the int32/int64 thing: I actually took the long long type from the scikit-cuda binding. However, I also tested int32 and int64 explicitly with the same result. https://github.com/lebedov/scikit-cuda/blob/master/skcuda/magma.py Thank you for pointing o...
by mrader1248
Thu Sep 28, 2017 5:09 am
Forum: User discussion
Topic: MAGMA + pycuda + my own CUDA kernels
Replies: 11
Views: 5076

MAGMA + pycuda + my own CUDA kernels

Hey all, currently I have a problem when using MAGMA, pycuda and my own CUDA kernels together and I am not sure whether this is an issue with MAGMA, pycuda or my knowledge of CUDA. In the code attached below, I am using CUDA via pycuda and also load MAGMA to perform a QR decomposition. In the end I ...