As one part of a larger parallel code, I want to do many small, dense, local linear solves. Small means 6 by 6 or 9 by 9. Building a KSP / Mat for this seems like overkill.
I could probably link to LAPACK and use dgesv_(), but I'm looking for a way to use "Petsc's LAPACK" if that makes sense. A quick look at $PETSC_DIR/include/petscblaslapack* didn't help me. What's the best way to do this? Thanks, Andrew
