[petsc-users] small dense solves

[Andrew T Barker]

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