25   Mat          A, A_save, B, C, P, C1, R;  in main()  local
92     PetscCall(MatProductCreate(A, B, NULL, &C));  in main()
93     PetscCall(MatSetOptionsPrefix(C, "AB_"));  in main()
94     PetscCall(MatProductSetType(C, MATPRODUCT_AB));  in main()
95     PetscCall(MatProductSetAlgorithm(C, MATPRODUCTALGORITHMDEFAULT));  in main()
96     PetscCall(MatProductSetFill(C, PETSC_DETERMINE));  in main()
97     PetscCall(MatProductSetFromOptions(C));  in main()
99     PetscCall(MatHasOperation(C, MATOP_PRODUCTSYMBOLIC, &flg));  in main()
100     PetscCall(MatProductSymbolic(C));  in main()
101     PetscCall(MatProductNumeric(C));  in main()
102     PetscCall(MatMatMultEqual(A, B, C, mcheck, &flg));  in main()
108     PetscCall(MatProductNumeric(C));  in main()
110     PetscCall(MatMatMultEqual(A, B, C, mcheck, &flg));  in main()
112     PetscCall(MatDestroy(&C));  in main()
115     PetscCall(MatMatMult(A, B, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &C));  in main()
122       PetscCall(MatMatMult(A, B, MAT_REUSE_MATRIX, PETSC_DETERMINE, &C));  in main()
124     PetscCall(MatMatMultEqual(A, B, C, mcheck, &flg));  in main()
129     PetscCall(MatProductClear(C));  in main()
130     PetscCall(MatDestroy(&C));  in main()
136       PetscCall(MatMatMult(A, B, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &C));  in main()
137       PetscCall(MatDestroy(&C));  in main()
173     PetscCall(MatProductCreate(P, B, NULL, &C));  in main()
174     PetscCall(MatSetOptionsPrefix(C, "AtB_"));  in main()
175     PetscCall(MatProductSetType(C, MATPRODUCT_AtB));  in main()
176     PetscCall(MatProductSetAlgorithm(C, MATPRODUCTALGORITHMDEFAULT));  in main()
177     PetscCall(MatProductSetFill(C, PETSC_DETERMINE));  in main()
178     PetscCall(MatProductSetFromOptions(C));  in main()
179     PetscCall(MatHasOperation(C, MATOP_PRODUCTSYMBOLIC, &flg));  in main()
181       PetscCall(MatProductSymbolic(C));                        /* equivalent to MatSetUp() */  in main()
182 …PetscCall(MatSetOption(C, MAT_USE_INODES, PETSC_FALSE)); /* illustrate how to call MatSetOption() …  in main()
183       PetscCall(MatProductNumeric(C));  in main()
184       PetscCall(MatProductNumeric(C)); /* test reuse symbolic C */  in main()
186       PetscCall(MatTransposeMatMultEqual(P, B, C, mcheck, &flg));  in main()
188       PetscCall(MatDestroy(&C));  in main()
191       PetscCall(MatTransposeMatMult(P, B, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &C));  in main()
192       PetscCall(MatTransposeMatMult(P, B, MAT_REUSE_MATRIX, PETSC_DETERMINE, &C));  in main()
193       PetscCall(MatGetInfo(C, MAT_GLOBAL_SUM, &info));  in main()
194       PetscCall(MatProductClear(C));  in main()
198       PetscCall(MatNormDifference(C, C1, &norm));  in main()
203       PetscCall(MatDuplicate(C, MAT_COPY_VALUES, &C1));  in main()
208     PetscCall(MatDestroy(&C));  in main()
217       PetscCall(MatMatTransposeMult(B, R, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &C));  in main()
218       PetscCall(MatSetOptionsPrefix(C, "ABt_")); /* enable '-ABt_' for matrix C */  in main()
219       PetscCall(MatGetInfo(C, MAT_GLOBAL_SUM, &info));  in main()
222       PetscCall(MatMatTransposeMult(B, R, MAT_REUSE_MATRIX, PETSC_DETERMINE, &C));  in main()
226       PetscCall(MatNormDifference(C, C1, &norm));  in main()
229       PetscCall(MatDestroy(&C));  in main()
239     PetscCall(MatProductCreate(A, P, NULL, &C));  in main()
240     PetscCall(MatSetOptionsPrefix(C, "PtAP_"));  in main()
241     PetscCall(MatProductSetType(C, MATPRODUCT_PtAP));  in main()
242     PetscCall(MatProductSetAlgorithm(C, MATPRODUCTALGORITHMDEFAULT));  in main()
243     PetscCall(MatProductSetFill(C, PETSC_DETERMINE));  in main()
244     PetscCall(MatProductSetFromOptions(C));  in main()
245     PetscCall(MatProductSymbolic(C));  in main()
246     PetscCall(MatProductNumeric(C));  in main()
247     PetscCall(MatPtAPMultEqual(A, P, C, mcheck, &flg));  in main()
249     PetscCall(MatProductNumeric(C)); /* reuse symbolic C */  in main()
251     PetscCall(MatPtAPMultEqual(A, P, C, mcheck, &flg));  in main()
253     PetscCall(MatDestroy(&C));  in main()
256     PetscCall(MatPtAP(A, P, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &C));  in main()
257     PetscCall(MatPtAPMultEqual(A, P, C, mcheck, &flg));  in main()
265       PetscCall(MatPtAP(A, P, MAT_REUSE_MATRIX, PETSC_DETERMINE, &C));  in main()
267     PetscCall(MatPtAPMultEqual(A, P, C, mcheck, &flg));  in main()
294         PetscCall(MatNormDifference(C, RARt, &norm));  in main()
303     PetscCall(MatDestroy(&C));  in main()