2: #include <petsc-private/matimpl.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: }