HexRule20.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 "HexRule.h"
 
namespace MeshLib
{

class HexRule20 : public HexRule
{
public:
    static const unsigned n_all_nodes = 20u;

    static const CellType cell_type = CellType::HEX20;

    constexpr static const unsigned face_nodes[6][8] = {
        {0, 3, 2, 1, 11, 10, 9, 8},    // Face 0
        {0, 1, 5, 4, 8, 17, 12, 16},   // Face 1
        {1, 2, 6, 5, 9, 18, 13, 17},   // Face 2
        {2, 3, 7, 6, 10, 19, 14, 18},  // Face 3
        {3, 0, 4, 7, 11, 16, 15, 19},  // Face 4
        {4, 5, 6, 7, 12, 13, 14, 15}   // Face 5
    };

    constexpr static const unsigned edge_nodes[12][3] = {
        {0, 1, 8},   // Edge 0
        {1, 2, 9},   // Edge 1
        {2, 3, 10},  // Edge 2
        {0, 3, 11},  // Edge 3
        {4, 5, 12},  // Edge 4
        {5, 6, 13},  // Edge 5
        {6, 7, 14},  // Edge 6
        {4, 7, 15},  // Edge 7
        {0, 4, 16},  // Edge 8
        {1, 5, 17},  // Edge 9
        {2, 6, 18},  // Edge 10
        {3, 7, 19}   // Edge 11
    };

    using EdgeReturn = MeshLib::QuadraticEdgeReturn;

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

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