libMesh::ReferenceElem Namespace Reference

This namespace implements singleton reference elements for each fundamental element type supported by libMesh. More...

Functions
const Elem &	get (const ElemType type_in)

Detailed Description

This namespace implements singleton reference elements for each fundamental element type supported by libMesh.

Author

Benjamin S. Kirk

Date

2013 Namespace providing access to reference geometric element types.

Function Documentation

get()

const Elem & libMesh::ReferenceElem::get

(

const ElemType

type_in

)

Returns

A constant reference to the reference element of the user-requested type.

Definition at line 254 of file reference_elem.C.

References libMesh::Utility::enum_to_string(), libMesh::INVALID_ELEM, libMesh::QUAD4, libMesh::QUAD8, libMesh::QUAD9, libMesh::QUADSHELL4, libMesh::QUADSHELL8, libMesh::QUADSHELL9, libMesh::TRI3, and libMesh::TRISHELL3.

Referenced by NodalNeighborsTest::do_test(), main(), libMesh::Elem::reference_elem(), SideVertexAverageNormalTest::testEdge2(), SideVertexAverageNormalTest::testEdge3(), VolumeTest::testEdge4Invertible(), WhichNodeAmITest::testHexes(), SideVertexAverageNormalTest::testHexes(), WhichNodeAmITest::testPrisms(), SideVertexAverageNormalTest::testPrisms(), VolumeTest::testPyramid5TrueCentroid(), WhichNodeAmITest::testPyramids(), SideVertexAverageNormalTest::testPyramids(), VolumeTest::testQuad4TrueCentroid(), SideVertexAverageNormalTest::testQuads(), WhichNodeAmITest::testTets(), SideVertexAverageNormalTest::testTets(), VolumeTest::testTri3TrueCentroid(), SideVertexAverageNormalTest::testTris(), and MeshSmootherTest::testVariationalSmoother().

{
  ElemType base_type = type_in;

  // For shell elements, use non shell type as the base type
  if (type_in == TRISHELL3)
    base_type = TRI3;

  if (type_in == QUADSHELL4)
    base_type = QUAD4;

  if (type_in == QUADSHELL8)
    base_type = QUAD8;

  if (type_in == QUADSHELL9)
    base_type = QUAD9;

  init_ref_elem_table();

  // Throw an error if the user asked for an ElemType that we don't
  // have a reference element for.
  libmesh_error_msg_if(ref_elem_map[base_type] == nullptr || type_in == INVALID_ELEM,
                       "No reference elem data available for ElemType " << type_in
                       << " = " << Utility::enum_to_string(type_in) << ".");

  return *ref_elem_map[base_type];
}