MAGMA Forum

In sgemm_fermi.cu, there is: #define magmablas_sgemm_fermi magmablas_sgemm extern "C" void magmablas_sgemm_fermi( char TRANSA, char TRANSB, int m , int n , int k , float alpha, const float *A, int lda, const float *B, int ldb, float beta, float *C, int ldc ) { // the actual sgemm code } Similarly, t...

John, I would say distributed use of MAGMA is of interest to us, but we don't currently have any code that does this, so you may have to wait a while for that functionality. Hooking MAGMA into Scalapack is one potential route. Another route could be using DAGue, which currently extends PLASMA to dis...

Alice, What I meant in my original reply is to use getrf and getrs -- or equivalently, gesv -- to solve Ax = b. getrf computes the LU factors, then getrs solves LUx = b using forward and back triangular solves. The inverse A^{-1} is never explicitly computed. Yes, getri takes the LU factors generate...

The routines are magma_*getrf() for LU, magma_*potrf() for Cholesky, magma_*qrtrf() for QR. These automatically do out-of-core computation, so the regular testing_*trf sample code should work, given a large matrix. An out-of-core triangular solve is not currently available in MAGMA, so you have to u...

Hi Alice, Usually that is caused by accessing the wrong memory, like a segfault in CPUs. Looks like I can replicate your problem, so it's a bug in MAGMA. I'll look into it. In the meantime, may I ask why are you computing the inverse? If you are trying to compute b = A^{-1} x, it is generally both f...

Hi John, Both of these goals -- running multiple problems simultaneously on a single GPU and running MAGMA distributed across several GPUs on several nodes -- are of interest to the MAGMA team. We are moving in that direction, with MAGMA 1.1 adding support for streams in the MAGMA BLAS and support f...

Those variables (blockIdx, etc.) are CUDA specific, as is the kernel launch <<< >>> notation. In order to compile CUDA specific code, you need to be use the CUDA nvcc compiler, instead of gcc. CUDA source code files should have an extension .cu instead of .c or .cpp, to differentiate them. The testi...

The omp_get_num_threads is from missing openmp, I believe from a quick web search. Unfortunately, the Intel MKL libraries vary between different versions. Using gcc and MKL, I currently use: LIB = -lmkl_gf_lp64 -lmkl_gnu_thread -lmkl_core -lpthread -lcublas -lm -fopenmp and LIBDIR = -L${MKLROOT}/lib...

Notice the Fortran compile command is gibberish: gfortran -O3 -DADD_ -x f95-cpp-input -Dmagma_devptr_t="integer(kind=../control/sizeptr.c: In function ‘main’: ../control/sizeptr.c:6:3: warning: format ‘%lu’ expects type ‘long unsigned int’, but argument 2 has type ‘unsigned int’ 4)" ... It should lo...

Most of the MAGMA routines are self-documented, which you can read in the files section of the doxygen documentation at http://icl.cs.utk.edu/magma/docs/ or in the source code itself. We don't yet have an updated manual. As for eigenvalues, there are syevd and heevd routines for the symmetric/Hermit...