43 double solution(const mfem::Vector &pt) {  in solution()
61   mfem::OptionsParser args(argc, argv);  in main()
82   mfem::Mesh *mesh = new mfem::Mesh(mesh_file, 1, 1);  in main()
95     mesh->SetCurvature(1, false, -1, mfem::Ordering::byNODES);  in main()
98     mesh->SetCurvature(order, false, -1, mfem::Ordering::byNODES);  in main()
104   mfem::FiniteElementCollection *fec     = new mfem::H1_FECollection(order, dim);  in main()
105   mfem::FiniteElementSpace      *fespace = new mfem::FiniteElementSpace(mesh, fec);  in main()
111   mfem::LinearForm          b(fespace);  in main()
112   mfem::FunctionCoefficient sol_coeff(solution);  in main()
113   b.AddDomainIntegrator(new mfem::DomainLFIntegrator(sol_coeff));  in main()
120   mfem::CGSolver cg;  in main()
130   mfem::GridFunction sol(fespace);  in main()
148     mfem::socketstream sol_sock(vishost, visport);  in main()