cl_EF_HEX8TB.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_EF_HEX8TB_HPP
#define BELFEM_CL_EF_HEX8TB_HPP
 
#include "nedelec/cl_EF_EdgeFunction.hpp"
#include "cl_FEM_Element.hpp"
 
namespace belfem
{
    namespace fem
    {
        class Element;

        class EF_HEX8TB : public EdgeFunction
        {
            Cell< Vector< real > > mF ;
            Cell< Matrix< real > > mFxi ;
 
            // edge directions
            real mS[ 4 ];
 
            // help vectors
            Vector< real > mU ; // midpoint of face {0,3,4,7} at xi = -1, then tangent t
            Vector< real > mV ; // midpoint of face {1,2,6,5} at xi = +1, then normal n
            Vector< real > mW ; // scratch, then binormal b

            Vector< real > mP ;

            Vector< real > mQ ;
 
            real mThickness  = BELFEM_QUIET_NAN ;
            real mWidth      = BELFEM_QUIET_NAN ;
            real mLength     = BELFEM_QUIET_NAN ;
            real mCrossSection = BELFEM_QUIET_NAN ;
 
            real mVolume  = BELFEM_QUIET_NAN ;
 
//------------------------------------------------------------------------------
        public:
//------------------------------------------------------------------------------

            EF_HEX8TB();

            ~EF_HEX8TB() override = default ;
 
//------------------------------------------------------------------------------

            void
            link( Element *aElement ) override;
 
//------------------------------------------------------------------------------
 
            void
            precompute( const Matrix< real > & aXi ) override ;
 
//------------------------------------------------------------------------------
 
            const Matrix <real> &
            E( const uint aIndex ) override ;
 
//------------------------------------------------------------------------------
 
            const Matrix <real> &
            C( const uint aIndex ) override ;
 
//------------------------------------------------------------------------------

            const Matrix< real > &
            G( const uint aIndex ) override ;
 
//------------------------------------------------------------------------------

            void
            update_nabla( const uint aIndex ) override ;
 
//------------------------------------------------------------------------------
        };
 
//------------------------------------------------------------------------------
 
        inline void
        EF_HEX8TB::update_nabla( const uint /* aIndex */ )
        {
            // nothing to do: mJ, mInvJ, mDetJ are element constants,
            // computed once in link()
        }
 
//------------------------------------------------------------------------------
    }
}
 
#endif //BELFEM_CL_EF_HEX8TB_HPP