10   Mat         mat;         /* matrix */  in main()  local
31   PetscCall(MatCreate(PETSC_COMM_SELF, &mat));  in main()
32   PetscCall(MatSetType(mat, MATSEQAIJ));  in main()
33   PetscCall(MatSetSizes(mat, n, n, n, n));  in main()
34   PetscCall(MatSetFromOptions(mat));  in main()
35   PetscCall(MatSeqAIJSetPreallocation(mat, 3, NULL));  in main()
36   PetscCall(MatSeqBAIJSetPreallocation(mat, 1, 3, NULL));  in main()
37   PetscCall(MatSeqSBAIJSetPreallocation(mat, 1, 3, NULL));  in main()
38   PetscCall(MatSeqSELLSetPreallocation(mat, 3, NULL));  in main()
46     PetscCall(MatSetValues(mat, 1, &i, 3, col, value, INSERT_VALUES));  in main()
51   PetscCall(MatSetValues(mat, 1, &i, 2, col, value, INSERT_VALUES));  in main()
57   PetscCall(MatSetValues(mat, 1, &i, 2, col, value, INSERT_VALUES));  in main()
58   PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY));  in main()
59   PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY));  in main()
62   PetscCall(MatMult(mat, ustar, b));  in main()
68   PetscCall(PCSetOperators(pc, mat, mat));  in main()
75   PetscCall(PCSetOperators(pc, mat, mat));  in main()
92   PetscCall(MatDestroy(&mat));  in main()