Recursion.h

Go to the documentation of this file.

#pragma once

 
#include <concepts>
 
#include <EmrysFilter/Core/Math.h>
#include <EmrysFilter/Distributions/Params.h>
#include <EmrysFilter/Spaces/Space.h>
 
namespace Emrys
{

template <typename R>
concept Recursion = RealScalar<typename R::Scalar> && Space<typename R::SpaceType> && requires(R const & recursion) {
    typename R::StateType;
    typename R::ModelType;
    typename R::MethodType;
    { R::Reentrant } -> std::convertible_to<bool>;
    { recursion.model() } -> std::same_as<typename R::ModelType const &>;
    { recursion.method() } -> std::same_as<typename R::MethodType const &>;
};

template <typename R, typename Component_>
concept RecursionOver = Recursion<R> && std::same_as<typename R::StateType, Params<Component_>>;
 
namespace internal
{

struct UpdateAccess
{
    template <typename Update>
    [[nodiscard]] static auto &
    logLikelihood(Update & update)
    {
        return update.m_logLikelihood;
    }

    template <typename Update>
    [[nodiscard]] static auto &
    posteriorEnds(Update & update)
    {
        return update.m_posteriorEnds;
    }
};
 
}  // namespace internal
 
}  // namespace Emrys