I've installed the new realesed of MAGMA library 1.1. The device on which I'm working is a Tesla M2050 with compute capability 2.0.I'm facing problem using the function magma_dgetri_gpu for computing the inverse of a matrix, using LU factorization, that give me a strange error 'can not bind to texture', which I'm not able to resolve.I've noticed that also the testing program provided with MAGMA library failed with the same error.

N CPU GFlop/s GPU GFlop/s ||R||_F / ||A||_F========================================================can not bind to texture can not bind to texture can not bind to texture can not bind to texture can not bind to texture

Have you ever encountered this problem? Have you got any idea about this message error??I'm working with other functions of the MAGMA 1.1 without any problems.

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 faster and more accurate to solve Ax=b by using gesv to compute the LU factors and solve LUx = b.

N CPU GFlop/s GPU GFlop/s ||R||_F / ||A||_F========================================================can not bind to texture can not bind to texture .......can not bind to texture can not bind to texture 2048 23.90 80.40 1.177364e-01

When I run the whole test case it runs past that size and only crashes later:

Dear all, Thank you. John, I faced the same problem when I run the test case without optional argument and when I set the matrix dimension with -N parameter. In the first case it crashes later. Mark, I've taken a look to magma documentation about magma_dgetri_gpu, it computes the inverse of a matrix using LU factorization and then it solves a system of the form LUx=b, as you suggested. In facts this function required on entry the factor L and U of LU decomposition instead of the original matrix.

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.