cl_ODE.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_ODE_HPP
#define BELFEM_CL_ODE_HPP
 
#include "typedefs.hpp"
#include "cl_Vector.hpp"
#include "en_ODE_Status.hpp"
namespace belfem
{
    namespace ode
    {
//------------------------------------------------------------------------------

        class ODE
        {
            const uint mDimension;
 
//------------------------------------------------------------------------------
        public:
//------------------------------------------------------------------------------
 
            ODE( const uint aDimension ) :
                mDimension( aDimension )
            {};
 
//------------------------------------------------------------------------------
 
            virtual ~ODE() = default;
 
//------------------------------------------------------------------------------
 
            virtual void
            compute(
                    const real           & aT,
                    const Vector< real > & aY,
                          Vector< real > & adYdT ) = 0;
 
//------------------------------------------------------------------------------

            inline uint
            dimension() const
            {
                return mDimension;
            }
 
 
//------------------------------------------------------------------------------
 
            virtual real
            check_events(
                    const real           & aT,
                    const Vector< real > & aY )
            {
                return 0.0 ;
            }
 
        };
    }
}
#endif //BELFEM_CL_ODE_HPP