[petsc-users] question about small matrices

[Povolotskyi, Mykhailo]

Dear Petsc developers,

in my application I have to solve millions of linear and non-linear 
systems with small matrices (2x2, 3x3,..., 10x10).

I consider them as dense, and use SNES with KSP method PREONLY, and LU 
preconditioner.

I found that when KSPSolve is called only 25% of time is spend in 
lapack, the rest is PETSc overhead.

I know how to call lapack directly to solve a linear system.

Question: is it possible to call lapack directly in the SNES solver to 
avoid the KSPSolve overhead?

Thank you,

Michael.