> 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
