I really wish column-major matrix order was never invented. It very quickly stops making sense after two dimensions and now I have to deal with it when interfacing with the cuSolver functions.
The Nvidia Cuda cuBlas and cuSolver library functions all accept matrices in column major order. I've been very successful up to now in adapting them for row-major matrices by moving arguments passed into them around, but with the factorization functions I am out of luck and will have to do an explicit transpose somewhere probably.
I am not at all familiar with how LU and QR decomposition work, so I am not sure whether I just need to transpose the input or whether I will also need to do the output.