# Recommendation for C/C++ library which offers Schur complement functions?

I need to find C/C++ libraries which offer function for computing Schur complement. I know about MUMPS and Pastix, but I need more of them to compare them in my research.

Do you have any experience with other libraries?

• What size are the matrices you're working with? What type of compute node(s)/computers are you running it on? If you're having to deal with distributed memory systems, the answers will be very different from a suggestion for your laptop. Jul 17 '15 at 2:47
• @NoseKnowsAll Very large matrices, it's meant to run on some nodes of our university supercomputer. But the sequential library is enough for now, its results will be send out separately. Jul 17 '15 at 8:23
• Probably one of the default libraries for you to look into for distributed computing would be ScaLapack: netlib.org/scalapack Jul 17 '15 at 15:29

I assume one could also directly use BLAS to compute a Schur complement $S = A - B*D^{-1}*C$ using routines for matrix products and inversion.