CellRule.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 "MeshLib/Node.h"
 
namespace MeshLib
{
 
class Element;
 
class CellRule
{
public:
    static const unsigned dimension = 3u;

    static bool testElementNodeOrder(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_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();
    }
}; /* class */
 
}  // namespace MeshLib