56 "--enable-petsc, --enable-trilinos, or --enable-eigen");
59 GetPot infile(
"vector_fe_ex3.in");
62 const unsigned int grid_size = infile(
"grid_size", 2);
65 libmesh_example_requires(3 <= LIBMESH_DIM,
"2D/3D support");
74 std::string elem_str =
78 libmesh_error_msg_if(elem_str !=
"TRI6" && elem_str !=
"TRI7" && elem_str !=
"QUAD8" && elem_str !=
"QUAD9",
79 "You selected: " << elem_str <<
80 " but this example must be run with TRI6, TRI7, QUAD8, or QUAD9.");
87 Utility::string_to_enum<ElemType>(elem_str));
101 system.
time_solver = std::make_unique<SteadySolver>(system);
104 equation_systems.init();
108 solver.
quiet = infile(
"solver_quiet",
true);
120 equation_systems.print_info();
130 std::vector<FunctionBase<Number> *> sols(1, &soln_func);
131 std::vector<FunctionBase<Gradient> *> grads(1, &soln_grad);
133 exact_sol.attach_exact_values(sols);
134 exact_sol.attach_exact_derivs(grads);
137 int extra_error_quadrature = infile(
"extra_error_quadrature", 2);
138 exact_sol.extra_quadrature_order(extra_error_quadrature);
141 exact_sol.compute_error(
"CurlCurl",
"u");
145 << exact_sol.l2_error(
"CurlCurl",
"u")
151 << exact_sol.hcurl_error(
"CurlCurl",
"u")
154 #ifdef LIBMESH_HAVE_EXODUS_API 159 #endif // #ifdef LIBMESH_HAVE_EXODUS_API This class handles the computation of the L2 and/or H1 error for the Systems in the EquationSystems o...
This is the EquationSystems class.
virtual void write_equation_systems(const std::string &, const EquationSystems &, const std::set< std::string > *system_names=nullptr)
This method implements writing a mesh with data to a specified file where the data is taken from the ...
bool quiet
The DiffSolver should not print anything to libMesh::out unless quiet is set to false; default is tru...
The ExodusII_IO class implements reading meshes in the ExodusII file format from Sandia National Labs...
unsigned int max_nonlinear_iterations
The DiffSolver should exit in failure if max_nonlinear_iterations is exceeded and continue_after_max_...
Real absolute_residual_tolerance
The DiffSolver should exit after the residual is reduced to either less than absolute_residual_tolera...
std::unique_ptr< TimeSolver > time_solver
A pointer to the solver object we're going to use.
The LibMeshInit class, when constructed, initializes the dependent libraries (e.g.
unsigned int variable_number(std::string_view var) const
SolverPackage default_solver_package()
This is a generic class that defines a solver to handle ImplicitSystem classes, including NonlinearIm...
T command_line_value(const std::string &, T)
unsigned int max_linear_iterations
Each linear solver step should exit after max_linear_iterations is exceeded.
void print_info(std::ostream &os=libMesh::out, const unsigned int verbosity=0, const bool global=true) const
Prints relevant information about the mesh.
void init(triangulateio &t)
Initializes the fields of t to nullptr/0 as necessary.
FEMSystem, TimeSolver and NewtonSolver will handle most tasks, but we must specify element residuals...
bool verbose
The DiffSolver may print a lot more to libMesh::out if verbose is set to true; default is false...
virtual void solve() override
Invokes the solver associated with the system.
The Mesh class is a thin wrapper, around the ReplicatedMesh class by default.
double initial_linear_tolerance
Any required linear solves will at first be done with this tolerance; the DiffSolver may tighten the ...
Real relative_residual_tolerance
Real relative_step_tolerance