28 #include "libmesh/equation_systems.h" 29 #include "libmesh/getpot.h" 30 #include "libmesh/exodusII_io.h" 31 #include "libmesh/mesh.h" 32 #include "libmesh/mesh_generation.h" 33 #include "libmesh/exact_solution.h" 34 #include "libmesh/string_to_enum.h" 35 #include "libmesh/enum_solver_package.h" 36 #include "libmesh/enum_norm_type.h" 39 #include "curl_curl_system.h" 40 #include "libmesh/diff_solver.h" 41 #include "libmesh/steady_solver.h" 42 #include "solution_function.h" 49 int main (
int argc,
char ** argv)
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);
130 std::vector<FunctionBase<Number> *> sols(1, &soln_func);
131 std::vector<FunctionBase<Gradient> *> grads(1, &soln_grad);
137 int extra_error_quadrature = infile(
"extra_error_quadrature", 2);
145 << exact_sol.
l2_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...
void attach_exact_values(const std::vector< FunctionBase< Number > *> &f)
Clone and attach arbitrary functors which compute the exact values of the EquationSystems' solutions ...
int main(int argc, char **argv)
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.
void print_info(std::ostream &os=libMesh::out) const
Prints information about the equation systems, by default to libMesh::out.
Real hcurl_error(std::string_view sys_name, std::string_view unknown_name)
void extra_quadrature_order(const int extraorder)
Increases or decreases the order of the quadrature rule used for numerical integration.
The LibMeshInit class, when constructed, initializes the dependent libraries (e.g.
The libMesh namespace provides an interface to certain functionality in the library.
void attach_exact_derivs(const std::vector< FunctionBase< Gradient > *> &g)
Clone and attach arbitrary functors which compute the exact gradients of the EquationSystems' solutio...
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 compute_error(std::string_view sys_name, std::string_view unknown_name)
Computes and stores the error in the solution value e = u-u_h, the gradient grad(e) = grad(u) - grad(...
void init(triangulateio &t)
Initializes the fields of t to nullptr/0 as necessary.
Real l2_error(std::string_view sys_name, std::string_view unknown_name)
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.
virtual void init()
Initialize all the systems.
virtual System & add_system(std::string_view system_type, std::string_view name)
Add the system of type system_type named name to the systems array.
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
Real error_norm(std::string_view sys_name, std::string_view unknown_name, const FEMNormType &norm)