cl_JcFunction_UserDefined.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_JCFUNCTION_USERDEFINED_HPP
#define BELFEM_CL_JCFUNCTION_USERDEFINED_HPP
 
#include "cl_JcFunction.hpp"
#include "cl_Material.hpp"
 
namespace belfem
{
    namespace material
    {
//------------------------------------------------------------------------------

        class JcFunctionUserDefined : public JcFunction
        {
            Material * mMaterial ;  
 
            MatFunc2  * mUserJcFunction2 = nullptr ;  
            MatFunc3  * mUserJcFunction3 = nullptr ;  
 
        public:

            JcFunctionUserDefined( Material * Mat ) : mMaterial( Mat ) {};
 
            JcFunctionUserDefined( const JcFunctionUserDefined & ) = delete;
            JcFunctionUserDefined & operator=( const JcFunctionUserDefined & ) = delete;
 
            ~JcFunctionUserDefined() override = default;

            real
            eval( const real normB, const real angle ) const override ;

            real
            eval( const real normB, const real angle, const real T ) const override ;

            void
            set_jc_function( MatFunc2 * Function )
            {
                mUserJcFunction2 = Function ;
                this->set_dependency( JcParameter::normB );
                this->set_dependency( JcParameter::angleNxB );
            }

            void
            set_jc_function( MatFunc3 * Function )
            {
                mUserJcFunction3 = Function ;
                this->set_dependency( JcParameter::normB );
                this->set_dependency( JcParameter::angleNxB );
                this->set_dependency( JcParameter::T );
            }
 
        private:

            real
            jc_user_2( const real normB, const real angle ) const
            {
                return mUserJcFunction2( mMaterial, normB, angle );
            }

            real
            jc_user_3( const real normB, const real angle, const real T ) const
            {
                return mUserJcFunction3( mMaterial, normB, angle, T );
            }
        };
 
        inline real
        JcFunctionUserDefined::eval( const real normB, const real angle ) const
        {
            BELFEM_ERROR( mUserJcFunction2 != nullptr,
                "2-parameter user function is not set" );
 
            return jc_user_2( normB, angle );
        }
 
        inline real
        JcFunctionUserDefined::eval( const real normB, const real angle, const real T ) const
        {
            BELFEM_ERROR( mUserJcFunction3 != nullptr,
                "3-parameter user function is not set" );
 
            return jc_user_3( normB, angle, T );
        }
 
//------------------------------------------------------------------------------
    }
}
#endif //BELFEM_CL_JCFUNCTION_USERDEFINED_HPP