#include <mesh_builder.h>

Public Member Functions
	MeshBuilder ()
	This Object handles building a libMesh::UnstructuredMesh subclass. More...

	~MeshBuilder ()

SharedPtr< libMesh::UnstructuredMesh >	build (const GetPot &input, const libMesh::Parallel::Communicator &comm LIBMESH_CAN_DEFAULT_TO_COMMWORLD)
	Builds the libMesh::Mesh according to input options. More...

void	do_mesh_refinement_from_input (const GetPot &input, const libMesh::Parallel::Communicator &comm, libMesh::UnstructuredMesh &mesh) const
	Refine the mesh based on user input parameters. More...

Private Member Functions
void	generate_mesh (const std::string &mesh_build_type, const GetPot &input, libMesh::UnstructuredMesh *mesh)

template<typename T >
void	deprecated_option (const GetPot &input, const std::string &old_option, const std::string &new_option, const T &default_value, T &option_value) const
	Helper function for displaying deprecated warnings. More...

Detailed Description

Definition at line 41 of file mesh_builder.h.

Constructor & Destructor Documentation

GRINS::MeshBuilder::MeshBuilder ( )

inline

This Object handles building a libMesh::UnstructuredMesh subclass.

Based on runtime input, either a generic 1, 2, or 3-dimensional mesh is built; or is read from input from a specified file.

Definition at line 48 of file mesh_builder.h.

48 {};

GRINS::MeshBuilder::~MeshBuilder ( )

inline

Definition at line 49 of file mesh_builder.h.

49 {};

Member Function Documentation

SharedPtr< libMesh::UnstructuredMesh > GRINS::MeshBuilder::build	(	const GetPot &	input,
		const libMesh::Parallel::Communicator &comm	LIBMESH_CAN_DEFAULT_TO_COMMWORLD
	)

Builds the libMesh::Mesh according to input options.

Todo:: Can remove last 4 checks once mesh-options/mesh_option support is removed.

Todo:: Need to a check a GMSH meshes

Definition at line 46 of file mesh_builder.C.

Referenced by GRINSTesting::MeshBuilderTest::build_mesh(), GRINSTesting::RayfireTest::build_mesh(), GRINSTesting::RayfireTestAMR::build_mesh(), main(), and GRINSTesting::SystemHelper::setup_multiphysics_system().

   {
     // First check if the user has both old and new versions of mesh input
     if( input.have_section("mesh-options/") &&
         input.have_section("Mesh/") )
       {
         libmesh_error_msg("Error: Detected illegal simulataneous use of [mesh-options] and [Mesh] in input!");
       }
 
     // User needs to tell us if we are generating or reading a mesh
     // We infer this by checking and seeing if the use has a Mesh/Read
     // or a Mesh/Generation section
     if( !input.have_variable("mesh-options/mesh_option") &&
         !input.have_section("Mesh/Read/") &&
         !input.have_section("Mesh/Generation/") )
       {
         libmesh_error_msg("ERROR: Must specify either Mesh/Read or Mesh/Generation in input.");
       }
 
     // But you can't have it both ways
     if( input.have_section("Mesh/Read/") && input.have_section("Mesh/Generation/") )
       {
         libmesh_error_msg("ERROR: Can only specify one of Mesh/Read and Mesh/Generation");
       }
 
     // Are we generating the mesh or are we reading one in from a file?
     std::string mesh_build_type = "NULL";
     if( input.have_section("Mesh/Read/") )
       mesh_build_type = "read";
 
     else if( input.have_section("Mesh/Generation/") )
       mesh_build_type = "generate";
 
     this->deprecated_option<std::string>( input, "mesh-options/mesh_option", "Mesh/Read or Mesh/Generation", "DIE!", mesh_build_type);
 
     // Make sure the user gave a valid option
     if( mesh_build_type != std::string("generate") &&
         mesh_build_type != std::string("read") &&
         mesh_build_type != std::string("read_mesh_from_file") &&
         mesh_build_type != std::string("create_1D_mesh") &&
         mesh_build_type != std::string("create_2D_mesh") &&
         mesh_build_type != std::string("create_3D_mesh") )
       {
         std::string error = "ERROR: Invalid value of "+mesh_build_type+" for Mesh/type.\n";
         error += "       Valid values are: generate\n";
         error += "                         read\n";
         libmesh_error_msg(error);
       }
 
     // Create UnstructuredMesh object (defaults to dimension 1).
     libMesh::UnstructuredMesh* mesh;
 
     // Were we specifically asked to use a ParallelMesh or SerialMesh?
     {
       std::string mesh_class = input("Mesh/class", "default");
 
       this->deprecated_option<std::string>( input, "mesh-options/mesh_class", "Mesh/class", "default", mesh_class);
 
       if (mesh_class == "parallel")
         mesh = new libMesh::ParallelMesh(comm);
       else if (mesh_class == "serial")
         mesh = new libMesh::SerialMesh(comm);
       else if (mesh_class == "default")
         mesh = new libMesh::Mesh(comm);
       else
         {
           std::string error = "ERROR: Invalid class "+mesh_class+" input for Mesh/class.\n";
           error += "       Valid choices are: serial, parallel.\n";
           libmesh_error_msg(error);
         }
     }
 
     // Read mesh from file
     if(mesh_build_type =="read_mesh_from_file" /* This is deprecated */ ||
        mesh_build_type == "read" )
       {
         // Make sure the user set the filename to read
         if( !input.have_variable("mesh-options/mesh_filename") /* This is deprecated */ &&
             !input.have_variable("Mesh/Read/filename") )
           {
             libmesh_error_msg("ERROR: Must specify Mesh/Read/filename for reading mesh.");
           }
 
         std::string mesh_filename = input("Mesh/Read/filename", "DIE!");
 
         this->deprecated_option<std::string>( input, "mesh-options/mesh_filename", "Mesh/Read/filename", "DIE!", mesh_filename);
 
         // According to Roy Stogner, the only read format
         // that won't properly reset the dimension is gmsh.
         mesh->read(mesh_filename);
 
         // If we have a first order mesh file but we need second order
         // elements we should fix that.
         bool all_second_order = input("Mesh/all_second_order", false);
         if (all_second_order)
           mesh->all_second_order();
       }
 
     // Generate the mesh using built-in libMesh functions
     else if(mesh_build_type=="create_1D_mesh" /* This is deprecated */ ||
             mesh_build_type=="create_2D_mesh" /* This is deprecated */ ||
             mesh_build_type=="create_3D_mesh" /* This is deprecated */ ||
             mesh_build_type=="generate")
       {
         this->generate_mesh(mesh_build_type,input,mesh);
       }
 
     else
       {
         // Shouldn't have gotten here
         libmesh_error();
       }
 
     /* Only do the mesh refinement here if we don't have a restart file.
        Otherwise, we need to wait until we've read in the restart file.
        That is done in Simulation::check_for_restart */
     if( !input.have_variable("restart-options/restart_file") )
       {
         this->do_mesh_refinement_from_input( input, comm, *mesh );
       }
 
     return SharedPtr<libMesh::UnstructuredMesh>(mesh);
   }

template<typename T >

void GRINS::MeshBuilder::deprecated_option	(	const GetPot &	input,
		const std::string &	old_option,
		const std::string &	new_option,
		const T &	default_value,
		T &	option_value
	)		const

inlineprivate

Helper function for displaying deprecated warnings.

Definition at line 81 of file mesh_builder.h.

References grins_warning.

Referenced by do_mesh_refinement_from_input(), and generate_mesh().

   {
     if( input.have_variable(old_option) )
       {
         std::string warning = "WARNING: "+old_option+" is DEPRECATED.\n";
         warning += "         Please update to use "+new_option+".\n";
         grins_warning(warning);
 
         option_value = input(old_option, default_value);
       }
 
     return;
   }

void GRINS::MeshBuilder::do_mesh_refinement_from_input	(	const GetPot &	input,
		const libMesh::Parallel::Communicator &	comm,
		libMesh::UnstructuredMesh &	mesh
	)		const

Refine the mesh based on user input parameters.

There are several parameters that allow for the user to specify mesh refinements. They are factored out here because where we need to apply the refinements depends on the existence of a restart file.

Definition at line 355 of file mesh_builder.C.

References deprecated_option().

Referenced by GRINS::Simulation::init_restart().

   {
     std::string redistribution_function_string =
       input("Mesh/Redistribution/function", std::string("0"));
     this->deprecated_option<std::string>( input, "mesh-options/redistribute", "Mesh/Redistribution/function", "0", redistribution_function_string );
 
     if (redistribution_function_string != "0")
       {
         libMesh::ParsedFunction<libMesh::Real>
           redistribution_function(redistribution_function_string);
 
         libMesh::MeshTools::Modification::redistribute
           (mesh, redistribution_function);
 
         // Redistribution can create distortions *within* second-order
         // elements, which can then be magnified by refinement.  Let's
         // undistort everything by converting to first order and back
         // if necessary.
 
         // FIXME - this only works for meshes with uniform geometry
         // order equal to FIRST or (full-order) SECOND.
 
         const libMesh::Elem *elem = *mesh.elements_begin();
 
         if (elem->default_order() != libMesh::FIRST)
           {
             mesh.all_first_order();
             mesh.all_second_order();
           }
       }
 
     int uniformly_refine = input("Mesh/Refinement/uniformly_refine", 0);
     this->deprecated_option( input, "mesh-options/uniformly_refine", "Mesh/Refinement/uniformly_refine", 0, uniformly_refine );
 
     if( uniformly_refine > 0 )
       {
         libMesh::MeshRefinement(mesh).uniformly_refine(uniformly_refine);
       }
 
     std::string h_refinement_function_string =
       input("Mesh/Refinement/locally_h_refine", std::string("0"));
     this->deprecated_option<std::string>( input, "mesh-options/locally_h_refine", "Mesh/Refinement/locally_h_refine", "0", h_refinement_function_string );
 
     if (h_refinement_function_string != "0")
       {
         libMesh::ParsedFunction<libMesh::Real>
           h_refinement_function(h_refinement_function_string);
 
         libMesh::MeshRefinement mesh_refinement(mesh);
 
         libMesh::dof_id_type found_refinements = 0;
         do {
           found_refinements = 0;
           unsigned int max_level_refining = 0;
 
           libMesh::MeshBase::element_iterator elem_it =
             mesh.active_elements_begin();
           libMesh::MeshBase::element_iterator elem_end =
             mesh.active_elements_end();
           for (; elem_it != elem_end; ++elem_it)
             {
               libMesh::Elem *elem = *elem_it;
 
               const libMesh::Real refinement_val =
                 h_refinement_function(elem->centroid());
 
               const unsigned int n_refinements = refinement_val > 0 ?
                 refinement_val : 0;
 
               if (elem->level() - uniformly_refine < n_refinements)
                 {
                   elem->set_refinement_flag(libMesh::Elem::REFINE);
                   found_refinements++;
                   max_level_refining = std::max(max_level_refining,
                                                 elem->level());
                 }
             }
 
           comm.max(found_refinements);
           comm.max(max_level_refining);
 
           if (found_refinements)
             {
               std::cout << "Found up to " << found_refinements <<
                 " elements to refine on each processor," << std::endl;
               std::cout << "with max level " << max_level_refining << std::endl;
               mesh_refinement.refine_and_coarsen_elements();
 
               if( input.have_variable("restart-options/restart_file") )
                 {
                   std::cout << "Warning: it is known that locally_h_refine is broken when restarting." << std::endl
                             << "         and multiple refinement passes are done. We are forcibly" << std::endl
                             << "         limiting the refinement to one pass until this issue is resolved." << std::endl;
                   break;
                 }
             }
 
         } while(found_refinements);
 
       }
 
     return;
   }

void GRINS::MeshBuilder::generate_mesh	(	const std::string &	mesh_build_type,
		const GetPot &	input,
		libMesh::UnstructuredMesh *	mesh
	)

private

Definition at line 173 of file mesh_builder.C.

References deprecated_option().

   {
     unsigned int dimension = input("Mesh/Generation/dimension",0);
 
     if( !input.have_variable("mesh-options/mesh_option") /* Deprecated */ &&
         !input.have_variable("Mesh/Generation/dimension") )
       {
         libmesh_error_msg("ERROR: Must specify Mesh/Generation/dimension for generating mesh.");
       }
 
     /* Remove these once suport for mesh-options/mesh_option is removed */
     if( mesh_build_type == "create_1D_mesh" /* This is deprecated */ )
       dimension = 1;
 
     if( mesh_build_type=="create_2D_mesh" /* This is deprecated */ )
       dimension = 2;
 
     if( mesh_build_type=="create_3D_mesh" /* This is deprecated */ )
       dimension = 3;
 
     // Set the mesh dimension
     mesh->set_mesh_dimension(dimension);
 
     /* Now look for spatial extent of the grid that the user wants to generate. */
 
     libMesh::Real x_min = input("Mesh/Generation/x_min", 0.0);
     this->deprecated_option( input, "mesh-options/domain_x1_min", "Mesh/Generation/x_min", 0.0, x_min );
 
     libMesh::Real x_max = input("Mesh/Generation/x_max", 1.0);
     this->deprecated_option( input, "mesh-options/domain_x1_max", "Mesh/Generation/x_max", 1.0, x_max );
 
 
     /* We only check the y_{min,max} input if dimension is > 1 so that GetPot
        UFO detection will give us an error if we have this in the input file
        and are only using a 1D grid. */
     libMesh::Real y_min = 0.0;
     libMesh::Real y_max = 0.0;
 
     if( dimension > 1 )
       {
         y_min = input("Mesh/Generation/y_min", 0.0);
         this->deprecated_option( input, "mesh-options/domain_x2_min", "Mesh/Generation/y_min", 0.0, y_min );
 
         y_max = input("Mesh/Generation/y_max", 1.0);
         this->deprecated_option( input, "mesh-options/domain_x2_max", "Mesh/Generation/y_max", 1.0, y_max );
       }
 
     /* We only check the z_{min,max} input if dimension is > 2 so that GetPot
        UFO detection will give us an error if we have this in the input file
        and are only using a 1D or 2D grid. */
     libMesh::Real z_min = 0.0;
     libMesh::Real z_max = 0.0;
 
     if( dimension > 2 )
       {
         z_min = input("Mesh/Generation/z_min", 0.0);
         this->deprecated_option( input, "mesh-options/domain_x3_min", "Mesh/Generation/z_min", 0.0, z_min );
 
         z_max = input("Mesh/Generation/z_max", 1.0);
         this->deprecated_option( input, "mesh-options/domain_x3_max", "Mesh/Generation/z_max", 1.0, z_max );
       }
 
     /* Now check for the number of elements in each direction */
 
     // Make sure user gave us info about how many elements to use
     if( !input.have_variable("mesh-options/mesh_nx1") /* Deprecated */ &&
         !input.have_variable("Mesh/Generation/n_elems_x") )
       {
         libmesh_error_msg("ERROR: Must supply Mesh/Generation/n_elems_x for mesh generation.");
       }
 
     unsigned int n_elems_x = input("Mesh/Generation/n_elems_x", 0);
     this->deprecated_option<unsigned int>( input, "mesh-options/mesh_nx1", "Mesh/Generation/n_elems_x", 0, n_elems_x );
 
     /* We only check n_elems_y input if dimension is > 1 so that GetPot
        UFO detection will give us an error if we have this in the input file
        and are only using a 1D grid. */
     unsigned int n_elems_y = 0;
     if( dimension > 1 )
       {
         if( !input.have_variable("mesh-options/mesh_nx2") /* Deprecated */ &&
             !input.have_variable("Mesh/Generation/n_elems_y") )
           {
             libmesh_error_msg("ERROR: Must supply Mesh/Generation/n_elems_y for mesh generation.");
           }
 
         n_elems_y = input("Mesh/Generation/n_elems_y", 0);
         this->deprecated_option<unsigned int>( input, "mesh-options/mesh_nx2", "Mesh/Generation/n_elems_y", 0, n_elems_y );
       }
 
     /* We only check n_elems_z input if dimension is > 2 so that GetPot
        UFO detection will give us an error if we have this in the input file
        and are only using a 1D or 2D grid. */
     unsigned int n_elems_z = 0;
     if( dimension > 2 )
       {
         if( !input.have_variable("mesh-options/mesh_nx3") /* Deprecated */ &&
             !input.have_variable("Mesh/Generation/n_elems_z") )
           {
             libmesh_error_msg("ERROR: Must supply Mesh/Generation/n_elems_z for mesh generation.");
           }
 
         n_elems_z = input("Mesh/Generation/n_elems_z", 0);
         this->deprecated_option<unsigned int>( input, "mesh-options/mesh_nx3", "Mesh/Generation/n_elems_z", 0, n_elems_z );
       }
 
     /* Now grab the element_type the user wants for the mesh. */
 
     std::string element_type = input("Mesh/Generation/element_type", "default");
     this->deprecated_option<std::string>( input, "mesh-options/element_type", "Mesh/Generation/element_type", "default", element_type );
 
     /* Now generate the mesh. */
     if( dimension == 1 )
       {
         if(element_type=="default")
           {
             element_type = "EDGE3";
           }
 
         GRINSEnums::ElemType element_enum_type =
           libMesh::Utility::string_to_enum<GRINSEnums::ElemType>(element_type);
 
         libMesh::MeshTools::Generation::build_line(*mesh,
                                                    n_elems_x,
                                                    x_min,
                                                    x_max,
                                                    element_enum_type);
       }
 
     else if( dimension == 2 )
       {
         if(element_type=="default")
           {
             element_type = "TRI6";
           }
 
         GRINSEnums::ElemType element_enum_type =
           libMesh::Utility::string_to_enum<GRINSEnums::ElemType>(element_type);
 
         libMesh::MeshTools::Generation::build_square(*mesh,
                                                      n_elems_x,
                                                      n_elems_y,
                                                      x_min,
                                                      x_max,
                                                      y_min,
                                                      y_max,
                                                      element_enum_type);
       }
 
     else if( dimension == 3 )
       {
         if(element_type=="default")
           {
             element_type = "TET10";
           }
 
         GRINSEnums::ElemType element_enum_type =
           libMesh::Utility::string_to_enum<GRINSEnums::ElemType>(element_type);
 
         libMesh::MeshTools::Generation::build_cube(*mesh,
                                                    n_elems_x,
                                                    n_elems_y,
                                                    n_elems_z,
                                                    x_min,
                                                    x_max,
                                                    y_min,
                                                    y_max,
                                                    z_min,
                                                    z_max,
                                                    element_enum_type);
       }
 
     else
       {
         // This shouldn't have happened
         libmesh_error();
       }
 
     return;
   }

The documentation for this class was generated from the following files:

src/solver/include/grins/mesh_builder.h
src/solver/src/mesh_builder.C

Public Member Functions

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation