magma_dsposv_gpu does not work when NRHS >1

[egyptwangja]

Hi,

magma_dsposv_gpu does not seem to work when NRHS > 1. In the example, testing_magma_dsposv_gpu.cpp, when I set NRHS = 2 on line 147 and line 278 (two places where NRHS is set), I get garbage values for the output X (res_ in this example). I tried manually assigning numbers in matrix A and b (As and Bs in this example) but I still get the same results.
Am I missing something?

Thanks!

[rnath]

Can you tell me about the matrix size and LDA ( Leading dimensions )
Thanks

[egyptwangja]

Hi,

I used the example in testing_magma_dsposv_gpu.cpp which runs through sizes of N = {1024,2048,3072,4032,5184,6016,7040,8064,9088,10112} where N=LDB = LDX = LDA = N. I also tried just for the first size 1024, and it didn't work either. I tried for a very small, 2X2 matrix, so that I could manually set the elements of A, B and verify the results in MATLAB, but it was the same issue.

Thanks!

[rnath]

Thank you for the information.
I am trying to isolate the problem. I will get back to you soon.

Thank you.