PyramidRule5.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 "EdgeReturn.h"
#include "Element.h"
#include "PyramidRule.h"
 
namespace MeshLib
{

class PyramidRule5 : public PyramidRule
{
public:
    static const unsigned n_all_nodes = 5u;

    static const CellType cell_type = CellType::PYRAMID5;

    constexpr static const unsigned face_nodes[5][4] = {
        {0, 1, 4, 99},  // Face 0
        {1, 2, 4, 99},  // Face 1
        {2, 3, 4, 99},  // Face 2
        {3, 0, 4, 99},  // Face 3
        {0, 3, 2, 1}    // Face 4
    };

    constexpr static const unsigned edge_nodes[8][2] = {
        {0, 1},  // Edge 0
        {1, 2},  // Edge 1
        {2, 3},  // Edge 2
        {0, 3},  // Edge 3
        {0, 4},  // Edge 4
        {1, 4},  // Edge 5
        {2, 4},  // Edge 6
        {3, 4}   // Edge 7
    };

    static const unsigned n_face_nodes[5];

    using EdgeReturn = MeshLib::LinearEdgeReturn;

    static const Element* getFace(const Element* e, unsigned i);

    static unsigned identifyFace(Node const* const* element_nodes,
                                 Node const* nodes[3])
    {
        return CellRule::identifyFace<PyramidRule5>(element_nodes, nodes);
    }
}; /* class */
 
}  // namespace MeshLib