7   Mat         pA, P, aijP;  in main()  local
27   PetscCall(MatCreate(PETSC_COMM_SELF, &pA));  in main()
28   PetscCall(MatSetSizes(pA, 3, 3, 3, 3));  in main()
29   PetscCall(MatSetType(pA, MATSEQAIJ));  in main()
30   PetscCall(MatSetUp(pA));  in main()
31   PetscCall(MatSetOption(pA, MAT_IGNORE_ZERO_ENTRIES, PETSC_TRUE));  in main()
32   PetscCall(MatSetValues(pA, 3, pij, 3, pij, pa, ADD_VALUES));  in main()
33   PetscCall(MatAssemblyBegin(pA, MAT_FINAL_ASSEMBLY));  in main()
34   PetscCall(MatAssemblyEnd(pA, MAT_FINAL_ASSEMBLY));  in main()
35   PetscCall(MatCreateMAIJ(pA, 3, &P));  in main()
36   PetscCall(MatDestroy(&pA));  in main()