MeshLib::CellRule Class Reference

Detailed Description

Definition at line 20 of file CellRule.h.

#include <CellRule.h>

Inheritance diagram for MeshLib::CellRule:

Static Public Member Functions
static bool	testElementNodeOrder (Element const &e)

Static Public Attributes
static const unsigned	dimension = 3u
	Constant: Dimension of this mesh element.

Static Protected Member Functions
template<typename ElementRule >
static unsigned	identifyFace (Node const *const *element_nodes, Node const *nodes[ElementRule::dimension])
	Returns the ID of a face given an array of nodes.

Member Function Documentation

identifyFace()

template<typename ElementRule >

static unsigned MeshLib::CellRule::identifyFace ( Node const *const * element_nodes, Node const * nodes[ElementRule::dimension] )

inlinestaticprotected

Returns the ID of a face given an array of nodes.

Definition at line 39 of file CellRule.h.

    {
        for (unsigned i = 0; i < ElementRule::n_faces; i++)
        {
            unsigned flag(0);
            constexpr std::size_t n = sizeof(ElementRule::face_nodes[0]) /
                                      sizeof(ElementRule::face_nodes[0][0]);
            for (unsigned j = 0; j < n; j++)
            {
                for (unsigned k = 0; k < ElementRule::dimension; k++)
                {
                    if (ElementRule::face_nodes[i][j] != 99 &&
                        element_nodes[ElementRule::face_nodes[i][j]] ==
                            nodes[k])
                    {
                        flag++;
                    }
                }
            }
            if (flag == ElementRule::dimension)
            {
                return i;
            }
        }
        return std::numeric_limits<unsigned>::max();
    }

testElementNodeOrder()

bool MeshLib::CellRule::testElementNodeOrder ( Element const & e )

static

Checks if the node order of an element is correct by testing surface normals. For 3D elements true is returned if the normals of all faces points away from the centre of the element. Note: This method might give wrong results if something else is wrong with the element (non-planar faces, non-convex geometry, possibly zero volume) which causes the calculated center of gravity to lie outside of the actual element

Definition at line 18 of file CellRule.cpp.

{
    auto const cc = getCenterOfGravity(e).asEigenVector3d();
    const unsigned nFaces(e.getNumberOfFaces());
    for (unsigned j = 0; j < nFaces; ++j)
    {
        MeshLib::Element const* const face(e.getFace(j));
        // Node 1 is checked below because that way all nodes are used for the
        // test at some point, while for node 0 at least one node in every
        // element type would be used for checking twice and one wouldn't be
        // checked at all. (based on the definition of the _face_nodes variable)
        auto const& x = face->getNode(1)->asEigenVector3d();
        Eigen::Vector3d const cx = x - cc;
        const double s = FaceRule::getSurfaceNormal(*face).dot(cx);
        delete face;
        if (s >= 0)
        {
            return false;
        }
    }
    return true;
}

References MathLib::Point3d::asEigenVector3d(), MeshLib::getCenterOfGravity(), MeshLib::Element::getFace(), MeshLib::Element::getNode(), MeshLib::Element::getNumberOfFaces(), and MeshLib::FaceRule::getSurfaceNormal().

Member Data Documentation

dimension

const unsigned MeshLib::CellRule::dimension = 3u

static

Constant: Dimension of this mesh element.

Definition at line 24 of file CellRule.h.

---

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

• MeshLib/Elements/CellRule.h
• MeshLib/Elements/CellRule.cpp