cl_FEM_DofMgr_SideSetData.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_FEM_DOFMGR_SIDESETDATA_HPP
#define BELFEM_CL_FEM_DOFMGR_SIDESETDATA_HPP
 
#include "typedefs.hpp"
#include "cl_Cell.hpp"
#include "cl_Map.hpp"
#include "cl_Vector.hpp"
#include "commtools.hpp"
#include "cl_IWG.hpp"
#include "cl_FEM_SideSet.hpp"
 
namespace belfem
{
    class Mesh;
 
    namespace fem
    {
 
        class Kernel;
 
        class DofManager;
 
        namespace dofmgr
        {
            class DofData ;
 
            class SideSetData
            {
                DofManager * mParent;

                Kernel * mKernel;

                Mesh * mMesh;
 
                // my rank
                const proc_t mCommRank;
 
                Cell< SideSet * > mSideSets;
                Map< id_t, SideSet * > mSideSetMap;
 
                // an empty sidset that is called if this set does not exist
                // on current proc
                SideSet * mEmptySideset = nullptr;
 
                // counts how many nodes on all sidesets experience convection
                index_t mNumberOfConvectionNodes ;
 
                Cell< Vector< index_t > > mConvectionNodes ;
 
//-----------------------------------------------------------------------------
            public:
//-----------------------------------------------------------------------------
 
                SideSetData( DofManager * aParent ) ;
 
//-----------------------------------------------------------------------------
 
                ~SideSetData();
 
//-----------------------------------------------------------------------------
 
                void
                create_sidesets();
 
//-----------------------------------------------------------------------------

                SideSet *
                sideset( const id_t aID ) ;
 
//-----------------------------------------------------------------------------

                Cell< SideSet * > &
                sidesets() ;
 
//-----------------------------------------------------------------------------

                 bool
                 sideset_exists( const id_t aID ) const;
 
//-----------------------------------------------------------------------------

                 void
                 collect_wetted_sidesets();
 
//-----------------------------------------------------------------------------
 
                 void
                 count_wetted_nodes();
 
//-----------------------------------------------------------------------------
 
                 index_t
                 number_of_wetted_nodes() const ;
 
//-----------------------------------------------------------------------------

                void
                create_alpha_fields();
 
 
//------------------------------------------------------------------------------

                void
                set_boundary_conditions() ;
 
//-----------------------------------------------------------------------------

                void
                reset();
 
//----------------------------------------------------------------------------
            };
 
//----------------------------------------------------------------------------
 
            inline SideSet *
            SideSetData::sideset( const id_t aID )
            {
                if ( mSideSetMap.key_exists( aID ) )
                {
                    return mSideSetMap( aID );
                }
                else
                {
                    mEmptySideset->set_mesh_id( aID );
                    return mEmptySideset ;
                }
            }
 
//----------------------------------------------------------------------------
 
            inline Cell< SideSet * > &
            SideSetData::sidesets()
            {
                return mSideSets ;
            }
 
//------------------------------------------------------------------------------
 
            inline bool
            SideSetData::sideset_exists( const id_t aID ) const
            {
                return mSideSetMap.key_exists( aID );
            }
 
//------------------------------------------------------------------------------
 
            inline index_t
            SideSetData::number_of_wetted_nodes() const
            {
                return mNumberOfConvectionNodes ;
            }
 
//------------------------------------------------------------------------------
 
        } /* end namespace dofmgr */
    } /* end namespace fem */
} /* end namespace belfem */
 
#endif //BELFEM_CL_FEM_DOFMGR_SIDESETDATA_HPP