7   Mat       F, RHS, X, C1;  in TestMatrix()  local
81   PetscCall(MatCreate(PETSC_COMM_WORLD, &RHS));  in TestMatrix()
82   PetscCall(MatSetSizes(RHS, PETSC_DECIDE, PETSC_DECIDE, n, nrhs));  in TestMatrix()
83   PetscCall(MatSetType(RHS, MATSEQDENSE));  in TestMatrix()
84   PetscCall(MatSetUp(RHS));  in TestMatrix()
85   for (i = 0; i < nrhs; i++) PetscCall(MatSetValue(RHS, i, i, 1.0, INSERT_VALUES));  in TestMatrix()
86   PetscCall(MatAssemblyBegin(RHS, MAT_FINAL_ASSEMBLY));  in TestMatrix()
87   PetscCall(MatAssemblyEnd(RHS, MAT_FINAL_ASSEMBLY));  in TestMatrix()
88   PetscCall(MatDuplicate(RHS, MAT_DO_NOT_COPY_VALUES, &X));  in TestMatrix()
91     PetscCall(MatMatSolve(F, RHS, X));  in TestMatrix()
93     PetscCall(MatAXPY(C1, -1.0, RHS, SAME_NONZERO_PATTERN));  in TestMatrix()
98     PetscCall(MatMatSolveTranspose(F, RHS, X));  in TestMatrix()
100     PetscCall(MatAXPY(C1, -1.0, RHS, SAME_NONZERO_PATTERN));  in TestMatrix()
112   PetscCall(MatDestroy(&RHS));  in TestMatrix()