events.hpp

#ifndef __MC_EVENTS_HPP__
#define __MC_EVENTS_HPP__
 
#include "biocma_cst_config.hpp"
#include <Kokkos_Core.hpp>
#include <common/execinfo.hpp>
#include <cstddef>
#include <span>
 
namespace MC
{
  enum class EventType : char
  {
    NewParticle = 0, 
    Exit,            
    Move,            
    Death,           
    Overflow,
    ChangeWeight, 
    __COUNT__
  };
 
  constexpr size_t number_event_type = static_cast<std::size_t>(
      EventType::__COUNT__); //< Number of different events

  template <EventType event>
  KOKKOS_INLINE_FUNCTION consteval size_t event_index()
  {
    return static_cast<size_t>(event);
  }

  struct EventContainer
  {
 
    // Use SharedHostPinnedSpace as this struct is meant to be small enough to
    // be shared between Host and Device According to SharedHostPinnedSpace
    // documentation, the size of this data can fit into one cache line so
    // transfer is not a botteneck
    Kokkos::View<std::size_t[number_event_type], Kokkos::SharedSpace> // NOLINT
        _events; // NOLINT(cppcoreguidelines-avoid-c-arrays,modernize-avoid-c-arrays)

    EventContainer() : _events("events")
    {
 
      Kokkos::deep_copy(_events, 0); // Ensure all event to 0 occurence
    }

    [[nodiscard]] std::span<std::size_t> get_span() const
    {
      // As we use SharedHostPinnedSpace, we can deal with classic std
      // containers to use host manipulation
      return {_events.data(), number_event_type};
    }

    KOKKOS_INLINE_FUNCTION void clear() const
    {
      Kokkos::deep_copy(_events, 0);
    }

    template <EventType event> [[nodiscard]] constexpr std::size_t get() const
    {
      static_assert(event != EventType::__COUNT__,
                    "Count is not a valid event");
      return _events(event_index<event>());
    }

    template <EventType event>
    KOKKOS_FORCEINLINE_FUNCTION constexpr void incr() const
    {
      static_assert(event != EventType::__COUNT__,
                    "Count is not a valid event");
      Kokkos::atomic_add(&_events[event_index<event>()], 1);
    }
 
    template <EventType event>
    KOKKOS_FORCEINLINE_FUNCTION constexpr void add(size_t val) const
    {
      static_assert(event != EventType::__COUNT__,
                    "Count is not a valid event");
      Kokkos::atomic_add(&_events[event_index<event>()], val);
    }
 
    template <EventType event>
    KOKKOS_FORCEINLINE_FUNCTION constexpr void wrap_incr() const
    {
      if constexpr (AutoGenerated::FlagCompileTime::enable_event_counter)
      {
        incr<event>();
      }
    }

    static EventContainer reduce(std::span<std::size_t> _data);
 
    template <class Archive> void save(Archive& ar) const
    {
      std::array<std::size_t, number_event_type> array{};
 
      auto rd = std::span<std::size_t>(_events.data(), number_event_type);
 
      std::copy(rd.begin(), rd.end(), array.begin());
      assert(rd[0] == _events[0] && rd[0] == array[0]);
 
      ar(array);
    }
 
    template <class Archive> void load(Archive& ar)
    {
 
      std::array<std::size_t, number_event_type> array{};
      ar(array);
 
      auto rd = std::span<std::size_t>(_events.data(), number_event_type);
 
      std::copy(array.begin(), array.end(), rd.begin());
      assert(rd[0] == _events[0]);
    }
  };
 
}; // namespace MC
 
#endif //__MC_EVENTS_HPP__