193     Mat pP, hP;  in main()  local
204     PetscCall(MatPtAP(C, B, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &hP));  in main()
205     PetscCall(MatPtAP(C, B, MAT_REUSE_MATRIX, PETSC_DETERMINE, &hP));  in main()
206     PetscCall(MatPtAPMultEqual(C, B, hP, 10, &flg));  in main()
211     PetscCall(MatAXPY(hP, -1., pP, DIFFERENT_NONZERO_PATTERN));  in main()
212     PetscCall(MatHasOperation(hP, MATOP_NORM, &flg));  in main()
214       PetscCall(MatConvert(hP, MATAIJ, MAT_INPLACE_MATRIX, &hP));  in main()
216     PetscCall(MatNorm(hP, NORM_INFINITY, &err));  in main()
217 …PetscCheck(err / norm <= PETSC_SMALL, PetscObjectComm((PetscObject)hP), PETSC_ERR_PLIB, "Error Mat…  in main()
220     PetscCall(MatDestroy(&hP));  in main()
224     PetscCall(MatPtAP(A, B, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &hP));  in main()
225     PetscCall(MatPtAP(A, B, MAT_REUSE_MATRIX, PETSC_DETERMINE, &hP));  in main()
226     PetscCall(MatPtAPMultEqual(A, B, hP, 10, &flg));  in main()
228     PetscCall(MatDestroy(&hP));  in main()