inverting independent, small dense matrices

[Shao-Ching Huang]

Hi

I need to invert a bunch of independent, dense 4x4 matrices on each
MPI process. (This is separate from the main computation, a parallel
Ax=b solve.) I am thinking of using LAPACK's DGESV to invert these 4x4
matrices.

In my code, is it necessary to call LAPACK directly? Or does PETSc
already have a convenient interface to LAPACK?

The size 4x4 may be increased to 8x8 or so later, so I am not trying
to use the closed-form inversion formula.

Thanks,

Shao-Ching