m_default.hpp

#ifndef __DEFAULT_MODEL_IMPLEMENTATION_HPP__
#define __DEFAULT_MODEL_IMPLEMENTATION_HPP__
 
#include <mc/alias.hpp>
#include <mc/traits.hpp>
#include <optional>
 
/***
DEFAULT MODEL IMPLEMENTATION
***/
 
struct DefaultModel
{
  enum class particle_var : int
  {
    a = 0,
  };
  static constexpr std::size_t n_var = 1;
  static constexpr std::size_t n_c = 1;
  static constexpr std::string_view name = "default_dynamic";
  using uniform_weight = std::true_type; // Using type alias
  using Self = DefaultModel;
  using FloatType = float;
  using SelfParticle = MC::ParticlesModel<Self::n_var, Self::FloatType>;
  using SelfContribs = MC::ParticlesModel<Self::n_c, Self::FloatType>;
  using Config = std::nullopt_t;
  static constexpr bool uniform_weigth = false;
 
  KOKKOS_INLINE_FUNCTION static void
  init([[maybe_unused]] const MC::pool_type& random_pool,
       [[maybe_unused]] std::size_t idx,
       [[maybe_unused]] const SelfParticle& arr)
  {
  }
 
  KOKKOS_INLINE_FUNCTION static double
  mass([[maybe_unused]] std::size_t idx,
       [[maybe_unused]] const SelfParticle& arr)
  {
    return 1.;
  }
 
  KOKKOS_INLINE_FUNCTION static MC::Status
  update([[maybe_unused]] const MC::pool_type& random_pool,
         [[maybe_unused]] FloatType d_t,
         [[maybe_unused]] std::size_t idx,
         [[maybe_unused]] const SelfParticle& arr,
         [[maybe_unused]] const SelfContribs& contribs_arr,
         [[maybe_unused]] const std::size_t position_index,
         [[maybe_unused]] const MC::LocalConcentration& c)
  {
 
    return MC::Status::Idle;
  }
 
  KOKKOS_INLINE_FUNCTION static void
  division([[maybe_unused]] const MC::pool_type& random_pool,
           [[maybe_unused]] std::size_t idx,
           [[maybe_unused]] std::size_t idx2,
           [[maybe_unused]] const SelfParticle& arr,
           [[maybe_unused]] const SelfParticle& buffer_arr)
  {
  }
};
 
CHECK_MODEL(DefaultModel)
 
/*
  Dynamic Model (with size not known at compile time)
*/
 
struct DynamicDefaultModel
{
  enum class particle_var : int
  {
    a = 0,
  };
 
  static constexpr std::size_t n_var = 1;
  static constexpr std::size_t n_c = 2;
  static constexpr std::string_view name = "default_dynamic";
  using uniform_weight = std::true_type; // Using type alias
  using Self = DynamicDefaultModel;
  using FloatType = float;
  using SelfParticle = MC::DynParticlesModel<FloatType>;
  using SelfContribs = MC::DynParticlesContribs<FloatType>;
  using Config = float;
  static constexpr bool uniform_weigth = false;
 
  static Config
  get_config(const std::size_t size)
  {
    (void)size;
    return 1.;
  };
 
  KOKKOS_INLINE_FUNCTION static void
  init([[maybe_unused]] const MC::pool_type& random_pool,
       [[maybe_unused]] std::size_t idx,
       [[maybe_unused]] const SelfParticle& arr,
       [[maybe_unused]] const Config& config)
  {
  }
 
  KOKKOS_INLINE_FUNCTION static double
  mass([[maybe_unused]] std::size_t idx,
       [[maybe_unused]] const SelfParticle& arr)
  {
    return 1.;
  }
 
  KOKKOS_INLINE_FUNCTION static MC::Status
  update([[maybe_unused]] const MC::pool_type& random_pool,
         [[maybe_unused]] FloatType d_t,
         [[maybe_unused]] std::size_t idx,
         [[maybe_unused]] const SelfParticle& arr,
         [[maybe_unused]] const SelfContribs& contribs_arr,
         [[maybe_unused]] const std::size_t position_index,
         [[maybe_unused]] const MC::LocalConcentration& c)
  {
 
    return MC::Status::Idle;
  }
 
  MC::ContribIndexBounds static get_bounds()
  {
    return { 0, 0 };
  }
 
  KOKKOS_INLINE_FUNCTION static void
  division([[maybe_unused]] const MC::pool_type& random_pool,
           [[maybe_unused]] std::size_t idx,
           [[maybe_unused]] std::size_t idx2,
           [[maybe_unused]] const SelfParticle& arr,
           [[maybe_unused]] const SelfParticle& buffer_arr)
  {
  }
};
 
CHECK_MODEL(DynamicDefaultModel)
 
#endif