2: #include <petscmat.h>
6: PetscErrorCode MatPtAP_Basic(Mat A,Mat P,MatReuse scall,PetscReal fill,Mat *C)
7: {
11: if (scall == MAT_INITIAL_MATRIX){
12: MatPtAPSymbolic(A,P,fill,C);
13: }
14: MatPtAPNumeric(A,P,*C);
15: return(0);
16: }