15   Mat                 A, mC, C;  in main()  local
59   PetscCall(MatProductCreate(A, P, NULL, &mC));  in main()
60   PetscCall(MatProductSetType(mC, MATPRODUCT_PtAP));  in main()
61   PetscCall(MatProductSetAlgorithm(mC, MATPRODUCTALGORITHMDEFAULT));  in main()
62   PetscCall(MatProductSetFill(mC, PETSC_DEFAULT));  in main()
63   PetscCall(MatProductSetFromOptions(mC));  in main()
64   PetscCall(MatProductSymbolic(mC));  in main()
65   PetscCall(MatProductNumeric(mC));  in main()
66   PetscCall(MatProductNumeric(mC));  in main()
67   PetscCall(MatProductGetAlgorithm(mC, &alg));  in main()
71   PetscCall(MatEqual(C, mC, &flg));  in main()
73   PetscCall(MatDestroy(&mC));  in main()
76   PetscCall(MatPtAP(A, P, MAT_INITIAL_MATRIX, 1., &mC));  in main()
77   PetscCall(MatPtAP(A, P, MAT_REUSE_MATRIX, 1., &mC));  in main()
80   PetscCall(MatEqual(C, mC, &flg));  in main()
82   PetscCall(MatDestroy(&mC));  in main()