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: }