I was running magma_zgetrf_gpu (from magma_zgesv_gpu) and found a memory error. I am calling the routine magma_zgetrf_gpu with n=m=1024.
In line 107, nb is set to 32. In line 108, s is set to 32. Then, in line 219, we have
- Code: Select all
lapackf77_zgetrf( &rows, &nb0, work, &lddwork, ipiv+s*nb, &iinfo);
The magma_int_t array "ipiv" is supposed to be length n. The above code segment will reference ipiv out of bounds.
P.S. This is in magma_1.1. So, if it is fixed in magma_1.2, I apologize!