We need to solve linear systems with both dense matrices and tri-diagonal matrices on a cluster and I have been reading the SLUG for a week. My understanding is that to solve a linear system with Ax=b using ScaLAPACK, I need to store the entire matrix A in a node, convert A to the 2-D block cyclic form, and pass the submatrices to PDGESV_.
Writing a script to convert a matrix to its 2-D block-cyclic form is not difficult, but could be error-prone. My question is, is there a C/FORTRAN code that handle this task? Could such a function be made as a part of ScaLapack auxilliary functions?
Thank you for any comment/pointer.