23   Mat               A, Ad, F, Fd, X, Xd, B;  in main()  local
62   PetscCall(MatCreateDense(PETSC_COMM_WORLD, m, PETSC_DECIDE, M, K, NULL, &Xd));  in main()
84     PetscCall(MatMatSolve(Fd, B, Xd));  in main()
86     PetscCall(MatViewFromOptions(Xd, NULL, "-Xd"));  in main()
87     PetscCall(MatAXPY(Xd, -1.0, X, SAME_NONZERO_PATTERN));  in main()
88     PetscCall(MatNorm(Xd, NORM_INFINITY, &norm));  in main()
89     PetscCall(MatViewFromOptions(Xd, NULL, "-MatMatSolve"));  in main()
93       PetscCall(MatMatSolveTranspose(Fd, B, Xd));  in main()
94       PetscCall(MatAXPY(Xd, -1.0, X, SAME_NONZERO_PATTERN));  in main()
95       PetscCall(MatNorm(Xd, NORM_INFINITY, &norm));  in main()
96       PetscCall(MatViewFromOptions(Xd, NULL, "-MatMatSolveTranspose"));  in main()
101     PetscCall(MatDenseGetColumnVecWrite(Xd, 0, &xd));  in main()
106     PetscCall(MatViewFromOptions(Xd, NULL, "-MatSolve"));  in main()
113       PetscCall(MatViewFromOptions(Xd, NULL, "-MatSolveTranspose"));  in main()
116     PetscCall(MatDenseRestoreColumnVecWrite(Xd, 0, &xd));  in main()
122   PetscCall(MatDestroy(&Xd));  in main()