oomph::MeshChecker Namespace Reference

Functions
template<class GEOM_ELEMENT_BASE , class ELEMENT >
void	assert_geometric_element (const unsigned &dim, const unsigned &nnode_1d=0)

Detailed Description

//////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////// Namespace with helper function to check element type in mesh constructors (say).

Function Documentation

assert_geometric_element()

template<class GEOM_ELEMENT_BASE , class ELEMENT >

void oomph::MeshChecker::assert_geometric_element	(	const unsigned &	dim,
		const unsigned &	nnode_1d = 0
	)

Helper function to assert that finite element of type ELEMENT can be cast to base class of type GEOM_ELEMENT_BASE – useful to avoid confusion if a mesh that was written for a specific element type (e.g. a QElement) is used with another one (e.g. a TElement. First argument specifies the required spatial dimension of the element (i.e. the number of local coordinates). The optional second argument specifies the required nnode_1d (i.e. the number of nodes along a 1D element edge). Can be omitted if the mesh can handle any number in which case this test is skipped.

    {
      // Only do tests in paranoia mode
#ifndef PARANOID
      return;
#endif
 
      // Instantiate element
      ELEMENT* el_pt = new ELEMENT;
 
      // Can we cast to required geometric element base type
      if (dynamic_cast<GEOM_ELEMENT_BASE*>(el_pt) == 0)
      {
        std::stringstream error_message;
        error_message << "You have specified an illegal element type! Element "
                         "is of type \n\n"
                      << typeid(el_pt).name()
                      << "\n\nand cannot be cast to type \n\n "
                      << typeid(GEOM_ELEMENT_BASE).name() << "\n\n";
        throw OomphLibError(error_message.str(),
                            OOMPH_CURRENT_FUNCTION,
                            OOMPH_EXCEPTION_LOCATION);
      }
 
      // Does the dimension match?
      if (dim != el_pt->dim())
      {
        std::stringstream error_message;
        error_message << "You have specified an illegal element type! Element "
                         "is of type \n\n"
                      << typeid(el_pt).name()
                      << "\n\nand has dimension = " << el_pt->dim()
                      << " but we need dim = " << dim << std::endl;
        throw OomphLibError(error_message.str(),
                            OOMPH_CURRENT_FUNCTION,
                            OOMPH_EXCEPTION_LOCATION);
      }
 
      // Does nnode_1d match?
      if (nnode_1d != 0)
      {
        if (nnode_1d != el_pt->nnode_1d())
        {
          std::stringstream error_message;
          error_message << "You have specified an illegal element type! "
                           "Element is of type \n\n"
                        << typeid(el_pt).name()
                        << "\n\nand has nnode_1d = " << el_pt->nnode_1d()
                        << " but we need nnode_1d = " << nnode_1d << std::endl;
          throw OomphLibError(error_message.str(),
                              OOMPH_CURRENT_FUNCTION,
                              OOMPH_EXCEPTION_LOCATION);
        }
      }
 
      // Clean up
      delete el_pt;
    }

References plotDoE::name, OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.