25   Mat          A, A_save, B, AT, ATT, BT, BTT, P, R, C, C1;  in main()  local
187     PetscCall(MatCreate(PETSC_COMM_WORLD, &P));  in main()
188     PetscCall(MatSetSizes(P, PETSC_DECIDE, PETSC_DECIDE, M, PN));  in main()
189     PetscCall(MatSetType(P, mattype));  in main()
190     PetscCall(MatSeqAIJSetPreallocation(P, nzp, NULL));  in main()
191     PetscCall(MatMPIAIJSetPreallocation(P, nzp, NULL, nzp, NULL));  in main()
192     PetscCall(MatGetOwnershipRange(P, &rstart, &rend));  in main()
199       PetscCall(MatSetValues(P, 1, &i, nzp, idxn, a, ADD_VALUES));  in main()
201     PetscCall(MatAssemblyBegin(P, MAT_FINAL_ASSEMBLY));  in main()
202     PetscCall(MatAssemblyEnd(P, MAT_FINAL_ASSEMBLY));  in main()
205     PetscCall(MatTranspose(P, MAT_INITIAL_MATRIX, &R));  in main()
209       PetscCall(MatTranspose(P, MAT_REUSE_MATRIX, &R));  in main()
215     PetscCall(MatTransposeMatMult(P, B, MAT_INITIAL_MATRIX, fill, &C));  in main()
219     PetscCall(MatTransposeMatMult(P, B, MAT_REUSE_MATRIX, fill, &C));  in main()
252       PetscCall(MatMatMult(B, P, MAT_INITIAL_MATRIX, fill, &C1));  in main()
258     PetscCall(MatDestroy(&P));  in main()
274     PetscCall(MatCreate(PETSC_COMM_WORLD, &P));  in main()
275     PetscCall(MatSetSizes(P, PETSC_DECIDE, PETSC_DECIDE, N, PN));  in main()
276     PetscCall(MatSetType(P, mattype));  in main()
277     PetscCall(MatSeqAIJSetPreallocation(P, nzp, NULL));  in main()
278     PetscCall(MatMPIAIJSetPreallocation(P, nzp, NULL, nzp, NULL));  in main()
280     PetscCall(MatGetOwnershipRange(P, &rstart, &rend));  in main()
286       PetscCall(MatSetValues(P, 1, &i, nzp, idxn, a, ADD_VALUES));  in main()
288     PetscCall(MatAssemblyBegin(P, MAT_FINAL_ASSEMBLY));  in main()
289     PetscCall(MatAssemblyEnd(P, MAT_FINAL_ASSEMBLY));  in main()
292     PetscCall(MatGetSize(P, &pM, &pN));  in main()
293     PetscCall(MatGetLocalSize(P, &pm, &pn));  in main()
294     PetscCall(MatPtAP(A, P, MAT_INITIAL_MATRIX, fill, &C));  in main()
301       PetscCall(MatPtAP(A, P, MAT_REUSE_MATRIX, fill, &C));  in main()
308       PetscCall(MatConvert(P, MATSEQDENSE, MAT_INITIAL_MATRIX, &Pdense));  in main()
311       PetscCall(MatConvert(P, MATMPIDENSE, MAT_INITIAL_MATRIX, &Pdense));  in main()
354       PetscCall(MatTranspose(P, MAT_INITIAL_MATRIX, &R));  in main()
379       PetscCall(MatTranspose(P, MAT_INITIAL_MATRIX, &R));  in main()
380       PetscCall(MatMatMatMult(R, A, P, MAT_INITIAL_MATRIX, 2.0, &RAP));  in main()
381       PetscCall(MatMatMatMult(R, A, P, MAT_REUSE_MATRIX, 2.0, &RAP));  in main()
390     PetscCall(MatGetLocalSize(P, &m, &n));  in main()
402       PetscCall(MatMult(P, x, v1));  in main()
405       PetscCall(MatMultTranspose(P, v2, v3)); /* v3 = Pt*A*P*x */  in main()
417     PetscCall(MatDestroy(&P));  in main()