cl_EF_LINE3.hpp

Go to the documentation of this file.

//
// Created by christian on 8/11/22.
//
 
#ifndef BELFEM_CL_EF_LINE3_HPP
#define BELFEM_CL_EF_LINE3_HPP
 
#include "nedelec/cl_EF_EdgeFunction.hpp"
#include "cl_FEM_Element.hpp"
#include "fn_polyval.hpp"
#include "fn_dot.hpp"
#include "fn_norm.hpp"
 
namespace belfem
{
    namespace fem
    {
        class Element;

        class EF_LINE3 : public EdgeFunction
        {
            // coordinates of edge
            Vector< real > mX ;
            Vector< real > mY ;
 
            // coordinates of edge dofs
            Vector< real > mXi;
 
            // scaling parameters
            Cell< Vector< real > > mNxi ;
            Cell< Matrix< real > > mE ;
 
            // coefficients for edge function on dof 0
            Vector< real > mC0 ;
 
            // coefficients for edge function on dof 1
            Vector< real > mC1 ;
 
            // shape function of volume element
            InterpolationFunction * mVolumeLagrange = nullptr ;
            Cell< Matrix< real > >  mVolumeNxi ;
 
//------------------------------------------------------------------------------
        public:
//------------------------------------------------------------------------------

            EF_LINE3();
 
//------------------------------------------------------------------------------

            ~EF_LINE3() override;
 
//------------------------------------------------------------------------------

            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:
//------------------------------------------------------------------------------
 
            real
            compute_det_J( const uint aIndex );
 
//------------------------------------------------------------------------------
        };
 
//------------------------------------------------------------------------------
 
        inline const Matrix< real > &
        EF_LINE3::E( const uint aIndex )
        {
            mAbsDetJ = this->compute_det_J( aIndex );
            return mE( aIndex );
        }
 
//--------------------------------------------------------------------------
    }
}
#endif //BELFEM_CL_EF_LINE3_HPP