26   Mat           A, F, RHS, X, C1;  in main()  local
64   PetscCall(ComputeRHSMatrix(m, nrhs, &RHS));  in main()
65   PetscCall(MatDuplicate(RHS, MAT_DO_NOT_COPY_VALUES, &X));  in main()
116   PetscCall(MatMatSolve(F, RHS, X));  in main()
118   PetscCall(MatAXPY(C1, -1.0, RHS, SAME_NONZERO_PATTERN));  in main()
128   PetscCall(MatDestroy(&RHS));  in main()
153   Mat          RHS;  in ComputeRHSMatrix()  local
158   PetscCall(MatCreate(PETSC_COMM_WORLD, &RHS));  in ComputeRHSMatrix()
159   PetscCall(MatSetSizes(RHS, m, PETSC_DECIDE, PETSC_DECIDE, nrhs));  in ComputeRHSMatrix()
160   PetscCall(MatSetType(RHS, MATSEQDENSE));  in ComputeRHSMatrix()
161   PetscCall(MatSetUp(RHS));  in ComputeRHSMatrix()
165   PetscCall(MatDenseGetArray(RHS, &array));  in ComputeRHSMatrix()
175   PetscCall(MatDenseRestoreArray(RHS, &array));  in ComputeRHSMatrix()
176   PetscCall(MatAssemblyBegin(RHS, MAT_FINAL_ASSEMBLY));  in ComputeRHSMatrix()
177   PetscCall(MatAssemblyEnd(RHS, MAT_FINAL_ASSEMBLY));  in ComputeRHSMatrix()
178   *C = RHS;  in ComputeRHSMatrix()