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.