1: #define PETSCMAT_DLL
3: #include petscmat.h
7: PetscErrorCode MatPtAP_Basic(Mat A,Mat P,MatReuse scall,PetscReal fill,Mat *C)
8: {
12: if (scall == MAT_INITIAL_MATRIX){
13: MatPtAPSymbolic(A,P,fill,C);
14: }
15: MatPtAPNumeric(A,P,*C);
16: return(0);
17: }