43 double solution(const mfem::Vector &pt) {  in solution()
52 double rhs(const mfem::Vector &pt) {  in rhs()
83   mfem::OptionsParser args(argc, argv);  in main()
104   mfem::Mesh *mesh = new mfem::Mesh(mesh_file, 1, 1);  in main()
117     mesh->SetCurvature(1, false, -1, mfem::Ordering::byNODES);  in main()
120     mesh->SetCurvature(order, false, -1, mfem::Ordering::byNODES);  in main()
126   mfem::FiniteElementCollection *fec     = new mfem::H1_FECollection(order, dim);  in main()
127   mfem::FiniteElementSpace      *fespace = new mfem::FiniteElementSpace(mesh, fec);  in main()
132   mfem::FunctionCoefficient sol_coeff(solution);  in main()
133   mfem::Array<int>          ess_tdof_list;  in main()
134   mfem::GridFunction        sol(fespace);  in main()
136     mfem::Array<int> ess_bdr(mesh->bdr_attributes.Max());  in main()
143   mfem::LinearForm          b(fespace);  in main()
144   mfem::FunctionCoefficient rhs_coeff(rhs);  in main()
145   b.AddDomainIntegrator(new mfem::DomainLFIntegrator(rhs_coeff));  in main()
151   mfem::Operator *D;  in main()
152   mfem::Vector    X, B;  in main()
156   mfem::CGSolver cg;  in main()
182     mfem::socketstream sol_sock(vishost, visport);  in main()