FaceRule.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "Element.h"
 
namespace MeshLib
{
 
class FaceRule
{
public:
    static const unsigned dimension = 2u;

    static const Element* getFace(const Element* e, unsigned i) { return e->getEdge(i); }

    static const unsigned n_faces = 0;

    static bool testElementNodeOrder(Element const& e);

    static Eigen::Vector3d getFirstSurfaceVector(Element const& e);

    static Eigen::Vector3d getSecondSurfaceVector(Element const& e);

    static Eigen::Vector3d getSurfaceNormal(Element const& e);
 
protected:
    template <typename ElementRule>
    static unsigned identifyFace(Node const* const* element_nodes,
                                 Node const* nodes[ElementRule::dimension])
    {
        for (unsigned i = 0; i < ElementRule::n_edges; i++)
        {
            unsigned flag(0);
            constexpr std::size_t n = sizeof(ElementRule::edge_nodes[0]) /
                                      sizeof(ElementRule::edge_nodes[0][0]);
            for (unsigned j = 0; j < n; j++)
            {
                for (unsigned k = 0; k < ElementRule::dimension; k++)
                {
                    if (element_nodes[ElementRule::edge_nodes[i][j]] ==
                        nodes[k])
                    {
                        flag++;
                    }
                }
            }
            if (flag == ElementRule::dimension)
            {
                return i;
            }
        }
        return std::numeric_limits<unsigned>::max();
    }
 
}; /* class */
 
}  // namespace MeshLib