[petsc-users] small dense solves

[Satish Balay]

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
>