9   Mat                A2, B2; // diag, offdiag of mat2  in main()  local
35   PetscCall(MatConvert(B, mat_type, MAT_INITIAL_MATRIX, &B2));  in main()
45   PetscCall(MatCreateMPIAIJWithSeqAIJ(PETSC_COMM_WORLD, M, N, A2, B2, garray_h, &mat2));  in main()
57   PetscCall(MatCreate(PETSC_COMM_SELF, &B2));  in main()
58   PetscCall(MatSetSizes(B2, m, N, m, N));  in main()
59   PetscCall(MatSetType(B2, mat_type));  in main()
60   PetscCall(MatSeqAIJSetPreallocation(B2, o_nz, NULL));  in main()
62 …for (int j = 0; j < bi[i + 1] - bi[i]; j++) PetscCall(MatSetValue(B2, i, garray[bj[bi[i] + j]], ba…  in main()
64   PetscCall(MatAssemblyBegin(B2, MAT_FINAL_ASSEMBLY));  in main()
65   PetscCall(MatAssemblyEnd(B2, MAT_FINAL_ASSEMBLY));  in main()
70   PetscCall(MatCreateMPIAIJWithSeqAIJ(PETSC_COMM_WORLD, M, N, A2, B2, NULL, &mat2));  in main()