PETSc just uses native lapack. There is nothing special here. petscblaslapack have c-prototypes for the native blas/lapack - so there is some error checking in petsc usage of blas/lapack.
Satish On Mon, 21 Jun 2010, Andrew T Barker wrote: > 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 >
