cl_IF_InterpolationFunctionTemplate.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_IF_INTERPOLATIONFUNCTIONTEMPLATE_HPP
#define BELFEM_CL_IF_INTERPOLATIONFUNCTIONTEMPLATE_HPP
 
#include "cl_IF_InterpolationFunction.hpp"
#include "assert.hpp"
 
namespace belfem
{
    namespace fem
    {
//------------------------------------------------------------------------------
        template< GeometryType G, InterpolationType T, uint D, uint B >
        class InterpolationFunctionTemplate : public InterpolationFunction
        {
//------------------------------------------------------------------------------
        public:
//------------------------------------------------------------------------------

            InterpolationFunctionTemplate() = default;
 
//------------------------------------------------------------------------------

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

            void
            param_coords( Matrix< real > & aXiHat ) const override
            {
                BELFEM_ERROR( false,
                "illegal call to template function InterpolationFunctionTemplate::param_coords()" );
            }
 
//------------------------------------------------------------------------------

            void
            N( const Vector< real > & aXi,
                     Matrix< real > & aN  ) const override
            {
                BELFEM_ERROR( false,
                "illegal call to template function InterpolationFunctionTemplate::N()" );
            }
 
//------------------------------------------------------------------------------

            void
            dNdXi( const Vector< real > & aXi,
                       Matrix< real > & adNdXi  ) const override
            {
                BELFEM_ERROR( false,
                "illegal call to template function InterpolationFunctionTemplate::adNdXi()" );
            }
 
//------------------------------------------------------------------------------

            void
            d2NdXi2( const Vector< real > & aXi,
                           Matrix< real > & ad2NdXi2  ) const override
            {
                BELFEM_ERROR( false,
                 "illegal call to template function InterpolationFunctionTemplate::ad2NdXi2()" );
            }
 
//------------------------------------------------------------------------------

            uint
            number_of_bases() const override
            {
                return B;
            }
 
//------------------------------------------------------------------------------

            uint
            number_of_dimensions() const override
            {
                return D;
            }
 
//------------------------------------------------------------------------------

            InterpolationOrder
            interpolation_order() const override
            {
                BELFEM_ERROR( false,
                             "illegal call to template function InterpolationFunctionTemplate::interpolation_order()" );
 
                return InterpolationOrder::UNDEFINED;
            }
 
//------------------------------------------------------------------------------

            InterpolationType
            interpolation_type() const override
            {
                return T;
            }
 
//------------------------------------------------------------------------------

            GeometryType
            geometry_type() const override
            {
                return G;
            }
 
//------------------------------------------------------------------------------

            ElementType
            element_type() const override
            {
                BELFEM_ERROR( false,
                             "illegal call to template function InterpolationFunctionTemplate::element_type()" );
 
                return ElementType::UNDEFINED;
            }
 
//------------------------------------------------------------------------------
 
        };
 
//------------------------------------------------------------------------------
    }
}
#endif //BELFEM_CL_IF_INTERPOLATIONFUNCTIONTEMPLATE_HPP