inverting independent, small dense matrices

Shao-Ching Huang schuang at
Fri Feb 13 15:49:50 CST 2009

>      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.


More information about the petsc-users mailing list