cl_Vertex.hpp

Go to the documentation of this file.

/*
 * BELFEM -- The Berkeley Lab Finite Element Framework
 * Copyright (c) 2026, The Regents of the University of California,
 * through Lawrence Berkeley National Laboratory (subject to receipt of any required
 * approvals from the U.S. Dept. of Energy).  All rights reserved.
 *
 * Developers: Christian Messe, Gregory Giard
 *
 * See the top-level LICENSE file for the complete license and disclaimer.
 */
 
#ifndef BELFEM_CL_VERTEX_HPP
#define BELFEM_CL_VERTEX_HPP
 
#include "typedefs.hpp"
#include "assert.hpp"
#include "cl_Mesh_Basis.hpp"
#include "Mesh_Enums.hpp"
 
namespace belfem
{
    namespace mesh
    {
        // forward declarations for nodes and edges
        class Node ;
        class Element;
        class Facet;
        class Edge ;
 
        class Vertex : public Basis
        {
        protected:
            Node ** mNodes = nullptr ;

            Edge ** mEdges = nullptr ;

            Face ** mFaces = nullptr ;

            Facet ** mFacets = nullptr ;

            Element ** mElements = nullptr ;

            uint16_t mNodeCounter = 0 ;
            uint16_t mNodeCapacity = 0 ;

            uint16_t mEdgeCounter  = 0 ;
            uint16_t mEdgeCapacity = 0 ;

            uint16_t mFaceCounter = 0 ;
            uint16_t mFaceCapacity = 0 ;

            uint16_t mFacetCounter = 0 ;
            uint16_t mFacetCapacity = 0 ;

            uint16_t mElementCounter  = 0 ;
            uint16_t mElementCapacity = 0 ;
 
//------------------------------------------------------------------------------
        public:
//------------------------------------------------------------------------------
 
            Vertex();
 
            ~Vertex() override;
 
//------------------------------------------------------------------------------

            EntityType
            entity_type() const override ;
 
//------------------------------------------------------------------------------
//  nodes
//------------------------------------------------------------------------------
 
            void
            reset_node_container();
 
            void
            increment_node_counter();
 
            void
            allocate_node_container( const uint aCounter=0 );
 
            void
            add_node( Node * aNode );
 
            void
            insert_node( Node * aNode, const uint aIndex );
 
            uint
            number_of_nodes() const override ;
 
            Node *
            node( const uint aIndex ) override;
 
            const Node *
            node( const uint aIndex ) const override ;
 
//------------------------------------------------------------------------------
//  edges
//------------------------------------------------------------------------------
 
            void
            reset_edge_container();
 
            void
            increment_edge_counter();
 
            void
            allocate_edge_container( const uint aCounter=0 );
 
            void
            add_edge( Edge * aEdge );
 
            void
            insert_edge( Edge * aEdge, const uint aIndex );
 
            uint
            number_of_edges() const override ;
 
            Edge *
            edge( const uint aIndex ) override;
 
            const Edge *
            edge( const uint aIndex ) const override ;
 
//------------------------------------------------------------------------------
//  faces
//------------------------------------------------------------------------------
 
            void
            reset_face_container();
 
            void
            increment_face_counter();
 
            void
            allocate_face_container( const uint aCounter=0 );
 
            void
            add_face( Face * aFace );
 
            uint
            number_of_faces() const override ;
 
            Face *
            face( const uint aIndex ) override;
 
            const Face *
            face( const uint aIndex ) const override ;
 
//------------------------------------------------------------------------------
//  facets
//------------------------------------------------------------------------------
 
            void
            reset_facet_container();
 
            void
            increment_facet_counter();
 
            void
            allocate_facet_container( const uint aNumFacets=0 );
 
            void
            add_facet( Facet * aFacet );
 
            virtual uint
            number_of_facets() const ;
 
            uint16_t
            facet_capacity() const ;
 
            virtual Facet *
            facet( const uint aIndex );
 
            virtual const Facet *
            facet( const uint aIndex ) const ;
 
//------------------------------------------------------------------------------
//  elements
//------------------------------------------------------------------------------
 
            void
            reset_element_container() override;
 
            void
            increment_element_counter();
 
            void
            allocate_element_container( const uint aCounter=0 );
 
            void
            add_element( Element * aElement );
 
            uint
            number_of_elements() const override ;
 
            Element *
            element( const uint aIndex ) override;
 
            virtual const Element *
            element( const uint aIndex ) const ;
 
//------------------------------------------------------------------------------
 
            virtual void
            flag_nodes( const uint8_t aIndex = 0 );
 
            virtual void
            unflag_nodes( const uint8_t aIndex = 0 );
 
//------------------------------------------------------------------------------
        protected:
//------------------------------------------------------------------------------

            void
            delete_containers();
 
//------------------------------------------------------------------------------

            size_t
            array_memory() const ;
 
//------------------------------------------------------------------------------
        };
 
 
//------------------------------------------------------------------------------
 
        inline void
        Vertex::increment_node_counter()
        {
            BELFEM_ASSERT(  mNodeCounter < 65535, "Node counter overflow." );
            ++mNodeCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline void
        Vertex::increment_edge_counter()
        {
            BELFEM_ASSERT(  mEdgeCounter < 65535, "Edge counter overflow." );
            ++mEdgeCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline void
        Vertex::increment_face_counter()
        {
            BELFEM_ASSERT(  mFaceCounter < 65535, "Face counter overflow." );
            ++mFaceCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline void
        Vertex::increment_facet_counter()
        {
            BELFEM_ASSERT(  mFacetCounter < 65535, "Facet counter overflow." );
            ++mFacetCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline void
        Vertex::increment_element_counter()
        {
            BELFEM_ASSERT(  mElementCounter < 65535, "Element counter overflow." );
            ++mElementCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline uint
        Vertex::number_of_nodes() const
        {
            return mNodeCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline size_t
        Vertex::array_memory() const
        {
            size_t tMem = this->number_of_vertices() * sizeof( graph::Vertex * );
 
            tMem += this->number_of_dofs() * sizeof( graph::Vertex * );
 
            // sources and weights are parallel arrays with the same slot count
            tMem += this->number_of_sources() * ( sizeof( Basis * ) + sizeof( real ) );
 
            tMem += mNodeCapacity    * sizeof( Node * );
            tMem += mEdgeCapacity    * sizeof( Edge * );
            tMem += mFaceCapacity    * sizeof( Face * );
            tMem += mFacetCapacity   * sizeof( Facet * );
            tMem += mElementCapacity * sizeof( Element * );
 
            return tMem ;
        }
 
//------------------------------------------------------------------------------
 
        inline uint
        Vertex::number_of_edges() const
        {
            return mEdgeCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline uint
        Vertex::number_of_facets() const
        {
            return mFacetCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline uint16_t
        Vertex::facet_capacity() const
        {
            return mFacetCapacity ;
        }
 
//------------------------------------------------------------------------------
 
        inline uint
        Vertex::number_of_elements() const
        {
            return mElementCounter ;
        }
 
//------------------------------------------------------------------------------
 
        inline Node *
        Vertex::node( const uint aIndex )
        {
            BELFEM_ASSERT( aIndex < ( uint ) mNodeCapacity,
                          "Node index %u is out of bounds for element %u, must be < %u.",
                          ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mNodeCapacity );
 
            return mNodes[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline const Node *
        Vertex::node( const uint aIndex ) const
        {
            BELFEM_ASSERT( aIndex < ( uint ) mNodeCapacity,
                          "Node index %u is out of bounds for element %u, must be < %u.",
                                  ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mNodeCapacity );
 
            return mNodes[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline Edge *
        Vertex::edge( const uint aIndex )
        {
            BELFEM_ASSERT( aIndex < ( uint ) mEdgeCapacity,
                          "Edge index %u is out of bounds for element %u, must be < %u.",
                                  ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mEdgeCapacity );
 
            return mEdges[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline const Edge *
        Vertex::edge( const uint aIndex ) const
        {
            BELFEM_ASSERT( aIndex < ( uint ) mEdgeCapacity,
                          "Edge index %u is out of bounds for element %u, must be < %u.",
                                  ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mEdgeCapacity );
 
            return mEdges[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline uint
        Vertex::number_of_faces() const
        {
            return mFaceCounter ;
        }
 
        inline Face *
        Vertex::face( const uint aIndex )
        {
            BELFEM_ASSERT( aIndex < ( uint ) mFaceCapacity,
                           "Face index %u is out of bounds for face %u, must be < %u.",
                           ( unsigned int ) aIndex,
                           ( unsigned int ) this->id(),
                           ( unsigned int ) mFaceCapacity );
 
            return mFaces[ aIndex ];
        }
 
        inline const Face *
        Vertex::face( const uint aIndex ) const
        {
            BELFEM_ASSERT( aIndex < ( uint ) mFaceCapacity,
                           "Face index %u is out of bounds for face %u, must be < %u.",
                           ( unsigned int ) aIndex,
                           ( unsigned int ) this->id(),
                           ( unsigned int ) mFaceCapacity );
 
            return mFaces[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline Facet *
        Vertex::facet( const uint aIndex )
        {
            BELFEM_ASSERT( aIndex < ( uint ) mFacetCapacity,
                          "Facet index %u is out of bounds for element %u, must be < %u.",
                                  ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mFacetCapacity );
 
            return mFacets[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline const Facet *
        Vertex::facet( const uint aIndex ) const
        {
            BELFEM_ASSERT( aIndex < ( uint ) mFacetCapacity,
                          "Facet index %u is out of bounds for element %u, must be < %u.",
                                  ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mFacetCapacity );
 
            return mFacets[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline Element *
        Vertex::element( const uint aIndex )
        {
            BELFEM_ASSERT( aIndex < ( uint ) mElementCapacity,
                          "Element index %u is out of bounds for element %u, must be < %u.",
                                  ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mElementCapacity );
 
            return mElements[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline const Element *
        Vertex::element( const uint aIndex ) const
        {
            BELFEM_ASSERT( aIndex < ( uint ) mElementCapacity,
                          "Element index %u is out of bounds for element %u, must be < %u.",
                                  ( unsigned int ) aIndex,
                          ( unsigned int ) this->id(),
                          ( unsigned int ) mElementCapacity );
 
            return mElements[ aIndex ];
        }
 
//------------------------------------------------------------------------------
 
        inline EntityType
        Vertex::entity_type() const
        {
            BELFEM_ERROR( false, "invalid call to abstract class");
            return EntityType::UNDEFINED ;
        }
 
//------------------------------------------------------------------------------
    }
}
#endif //BELFEM_CL_VERTEX_HPP