When I solve SLE after matrix factorization (attached) with function sgetrs_gpu and nrhs>1, I get NaN in solution vector. If nrhs=1, answer is true. dgetrf.cpp is code for solution attached test matrix.
What am I doing wrong?
mgates3 wrote:There are a few things to improve in the code.
* You allocate A twice, once with malloc and once with magma_malloc_cpu, leading to a memory leak.
* Check whether info != 0 after all LAPACK and MAGMA functions that take info.
* Verify that A, b, ipiv, A_ptr, and b_ptr are not NULL after allocating them.
* There is no need to copy A_ptr back to A after getrf, unless you need it for some other purpose.
Users browsing this forum: No registered users and 5 guests