I am new to using cusolver with dense matrices and would like to know if there is information on how to use this library to implement a particular version of LU Decomposition. I would like to be able to pass a matrix and vector to solve using LU Decomposition whereby the vector is actually made up of multiple columns rather than a single column. Is this possible and if so can anyone provide a link/example code to this?
solution of a dense linear system using LU factorization is supported using the routines cusolverDnXgetrf for factorization followed by cusolverDnXgetrs for the solve. Multiple RHS are supported. A sample code is here. Yes, I’m aware that the sample code doesn’t demonstrate multiple RHS. That should be a fairly trivial adaptation.