inverting independent, small dense matrices

[Shao-Ching Huang]

>      For PETSc's own internal use we have written custom inverters for 
> small matrices (these are used by
> the factorization routines for BAIJ matrices).  They are called 
> Kernel_A_gets_inverse_A_* where * is the block size.
> The prototypes for the functions are in src/inline/ilu.h and the actual 
> source code in files in src/mat/impls/baij/seq/dgef*.c
> 

Beautiful! I will try it out.

Thank you.

Shao-Ching