cl_EF_TRI6.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_TRI6_HPP
#define BELFEM_CL_EF_TRI6_HPP
 
#include "nedelec/cl_EF_EdgeFunction.hpp"
#include "cl_FEM_Element.hpp"
 
namespace belfem
{
    namespace fem
    {
        class Element;
 
        class EF_TRI6 : public EdgeFunction
        {
            // node coordinates in matrix
            Matrix< real > mNodeCoords ;
 
            // node coordinates in vectors
            Vector< real > mX;
            Vector< real > mY;
 
            real mS[ 3 ];
 
            real mNablaXi[ 2 ];
            real mNablaEta[ 2 ];

            Matrix< real > mNxi ;
            Matrix< real > mNeta ;
 
            // help coefficients for shape function
            Matrix< real > mG ;
            Matrix< real > mGxi ;
            Matrix< real > mGeta ;
 
            Matrix< real > mH ;
            Matrix< real > mHxi ;
            Matrix< real > mHeta ;
 
            real mExi[ 16 ] ;
            real mEeta[ 16 ] ;
 
            // additional work vector
            real mW[ 16 ];
 
            // contracted second derivatives of the quadratic geometry map,
            // S^(pq)_m ( IF-Voigt rows xx, yy, xy ), constant per element;
            // filled on the CURVED branch of link(), NaN-filled on the
            // straight branch ( the hess channel of G must never read it
            // there — the affine map has S = 0 and the channel is skipped )
            real mCurv[ 3 ][ 2 ];
 
            uint mLastJ = BELFEM_UINT_MAX ;
            uint mLastNabla = BELFEM_UINT_MAX ;
            void
            ( EF_TRI6::*mFunInterpolation )( const uint aIndex );
 
            void
            ( EF_TRI6::*mFunCurl )( const uint aIndex );
 
            void
            ( EF_TRI6::*mFunGrad )( const uint aIndex );
 
//------------------------------------------------------------------------------
        public:
//------------------------------------------------------------------------------

            EF_TRI6();
 
//------------------------------------------------------------------------------

            ~EF_TRI6() override = default;
 
//------------------------------------------------------------------------------
 
            void
            precompute( const Matrix< real > & aXi ) override;
 
//------------------------------------------------------------------------------

            void
            link( Element * aElement ) override;
 
//------------------------------------------------------------------------------

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

            const Matrix< real > &
            C( const uint aIndex = 0 ) override;
 
//------------------------------------------------------------------------------

            const Matrix< real > &
            G( const uint aIndex = 0 ) override ;
 
//------------------------------------------------------------------------------
        private:
//------------------------------------------------------------------------------
 
            void
            E_curved( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            C_curved( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            C_straight( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            compute_J( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            compute_nabla( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            compute_E( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            compute_E_xi( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            compute_E_eta( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            compute_C();
 
//------------------------------------------------------------------------------

            void
            G_term1( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            G_straight( const uint aIndex );
 
//------------------------------------------------------------------------------
 
            void
            G_curved( const uint aIndex );
 
//------------------------------------------------------------------------------
        };
 
 
//------------------------------------------------------------------------------
 
        inline const Matrix< real > &
        EF_TRI6::E( const uint aIndex )
        {
            ( this->*mFunInterpolation )( aIndex );
            return mE ;
        }
 
//------------------------------------------------------------------------------
 
        inline const Matrix< real > &
        EF_TRI6::C( const uint aIndex )
        {
            ( this->*mFunCurl )( aIndex );
            return mC ;
        }
 
//------------------------------------------------------------------------------
    } /* end namespace fem */
}  /* end namespace belfem */
 
#endif //BELFEM_CL_EF_TRI6_HPP