22 #include "libmesh/libmesh.h" 23 #include "libmesh/mesh.h" 24 #include "libmesh/equation_systems.h" 25 #include "libmesh/enum_norm_type.h" 31 for (
unsigned int k = 0; k != sys.
n_vars(); ++k)
34 ", vector " << vecname <<
52 int main(
int argc,
char ** argv)
60 <<
" mesh solution" << std::endl;
72 libMesh::out <<
"Loaded solution " << argv[2] << std::endl;
77 for (
unsigned int i = 0; i != es.n_systems(); ++i)
79 System & sys = es.get_system(i);
82 for (
unsigned int j = 0; j != sys.
n_vectors(); ++j)
This is the EquationSystems class.
virtual void read(const std::string &name, void *mesh_data=nullptr, bool skip_renumber_nodes_and_elements=false, bool skip_find_neighbors=false)=0
Interfaces for reading/writing a mesh to/from a file.
std::streamsize precision() const
Get the associated write precision.
The LibMeshInit class, when constructed, initializes the dependent libraries (e.g.
The libMesh namespace provides an interface to certain functionality in the library.
int main(int argc, char **argv)
unsigned int n_vectors() const
Manages consistently variables, degrees of freedom, and coefficient vectors.
void print_info(std::ostream &os=libMesh::out, const unsigned int verbosity=0, const bool global=true) const
Prints relevant information about the mesh.
std::unique_ptr< NumericVector< Number > > solution
Data structure to hold solution values.
Real calculate_norm(const NumericVector< Number > &v, unsigned int var, FEMNormType norm_type, std::set< unsigned int > *skip_dimensions=nullptr) const
void init(triangulateio &t)
Initializes the fields of t to nullptr/0 as necessary.
const std::string & variable_name(const unsigned int i) const
void output_norms(const System &sys, const NumericVector< Number > &vec, const std::string &vecname)
const std::string & vector_name(const unsigned int vec_num) const
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const std::string & name() const
unsigned int n_vars() const
The Mesh class is a thin wrapper, around the ReplicatedMesh class by default.
const NumericVector< Number > & get_vector(std::string_view vec_name) const