10   Mat         A, RHS, C, F, X, AX, spRHST;  in main()  local
156     PetscCall(MatCreate(PETSC_COMM_WORLD, &spRHST));  in main()
159       PetscCall(MatSetSizes(spRHST, nrhs, M, PETSC_DECIDE, PETSC_DECIDE));  in main()
161       PetscCall(MatSetSizes(spRHST, 0, 0, PETSC_DECIDE, PETSC_DECIDE));  in main()
163     PetscCall(MatSetType(spRHST, MATAIJ));  in main()
164     PetscCall(MatSetFromOptions(spRHST));  in main()
165     PetscCall(MatSetUp(spRHST));  in main()
168       for (i = 0; i < nrhs; i++) PetscCall(MatSetValues(spRHST, 1, &i, 1, &i, &v, INSERT_VALUES));  in main()
170     PetscCall(MatAssemblyBegin(spRHST, MAT_FINAL_ASSEMBLY));  in main()
171     PetscCall(MatAssemblyEnd(spRHST, MAT_FINAL_ASSEMBLY));  in main()
173     PetscCall(MatMatTransposeSolve(F, spRHST, X));  in main()
195       PetscCall(MatCreateTranspose(spRHST, &spRHS));  in main()
199       PetscCall(MatMumpsGetInverseTranspose(F, spRHST));  in main()
202         PetscCall(MatView(spRHST, PETSC_VIEWER_STDOUT_WORLD));  in main()
209     PetscCall(MatDestroy(&spRHST));  in main()