entity-toolkit · haykh · Jan 27, 2026 · Jan 29, 2026 · Jan 30, 2026 · Jan 30, 2026
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -47,9 +47,6 @@ set(pgen
     ${default_pgen}
     CACHE STRING "Problem generator")
 
-set(gui
-    ${default_gui}
-    CACHE BOOL "Use GUI [nttiny]")
 set(output
     ${default_output}
     CACHE BOOL "Enable output")
@@ -70,7 +67,7 @@ if(${DEBUG} STREQUAL "OFF")
   set(CMAKE_BUILD_TYPE
       Release
       CACHE STRING "CMake build type")
-  set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -DNDEBUG")
+  set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -DNDEBUG -O3")
 else()
   set(CMAKE_BUILD_TYPE
       Debug
@@ -170,21 +167,29 @@ endif()
 
 link_libraries(${DEPENDENCIES})
 
+set(SRC_DIR ${CMAKE_CURRENT_SOURCE_DIR}/src/)
+add_subdirectory(${SRC_DIR}/global ${CMAKE_CURRENT_BINARY_DIR}/global)
+add_subdirectory(${SRC_DIR}/metrics ${CMAKE_CURRENT_BINARY_DIR}/metrics)
+add_subdirectory(${SRC_DIR}/kernels ${CMAKE_CURRENT_BINARY_DIR}/kernels)
+add_subdirectory(${SRC_DIR}/archetypes ${CMAKE_CURRENT_BINARY_DIR}/archetypes)
+add_subdirectory(${SRC_DIR}/framework ${CMAKE_CURRENT_BINARY_DIR}/framework)
+add_subdirectory(${SRC_DIR}/output ${CMAKE_CURRENT_BINARY_DIR}/output)
+
+# ------------------------------- Main source ------------------------------ #
+if(NOT ${pgen} STREQUAL ${default_pgen})
+  set_problem_generator(${pgen})
+  add_subdirectory(${SRC_DIR}/engines ${CMAKE_CURRENT_BINARY_DIR}/engines)
+  add_subdirectory(${SRC_DIR} ${CMAKE_CURRENT_BINARY_DIR}/src)
+endif()
+
 if(TESTS)
   # ---------------------------------- Tests --------------------------------- #
   include(${CMAKE_CURRENT_SOURCE_DIR}/cmake/tests.cmake)
-elseif(BENCHMARK)
+endif()
+
+if(BENCHMARK)
   # ------------------------------ Benchmark --------------------------------- #
   include(${CMAKE_CURRENT_SOURCE_DIR}/cmake/benchmark.cmake)
-else()
-  # ----------------------------------- GUI ---------------------------------- #
-  if(${gui})
-    find_or_fetch_dependency(nttiny FALSE QUIET)
-  endif()
-
-  # ------------------------------- Main source ------------------------------ #
-  set_problem_generator(${pgen})
-  add_subdirectory(${CMAKE_CURRENT_SOURCE_DIR}/src src)
 endif()
 
 include(${CMAKE_CURRENT_SOURCE_DIR}/cmake/report.cmake)
diff --git a/cmake/benchmark.cmake b/cmake/benchmark.cmake
@@ -2,17 +2,6 @@
 
 set(SRC_DIR ${CMAKE_CURRENT_SOURCE_DIR}/src)
 
-add_subdirectory(${SRC_DIR}/global ${CMAKE_CURRENT_BINARY_DIR}/global)
-add_subdirectory(${SRC_DIR}/metrics ${CMAKE_CURRENT_BINARY_DIR}/metrics)
-add_subdirectory(${SRC_DIR}/kernels ${CMAKE_CURRENT_BINARY_DIR}/kernels)
-add_subdirectory(${SRC_DIR}/archetypes ${CMAKE_CURRENT_BINARY_DIR}/archetypes)
-add_subdirectory(${SRC_DIR}/framework ${CMAKE_CURRENT_BINARY_DIR}/framework)
-
-if(${output})
-  add_subdirectory(${SRC_DIR}/output ${CMAKE_CURRENT_BINARY_DIR}/output)
-  add_subdirectory(${SRC_DIR}/checkpoint ${CMAKE_CURRENT_BINARY_DIR}/checkpoint)
-endif()
-
 set(exec benchmark.xc)
 set(src ${CMAKE_CURRENT_SOURCE_DIR}/benchmark/benchmark.cpp)
 

diff --git a/cmake/report.cmake b/cmake/report.cmake
@@ -8,7 +8,9 @@ if(${PGEN_FOUND})
     "${Blue}"
     PGEN_REPORT
     0)
-elseif(${TESTS})
+endif()
+
+if(${TESTS})
   set(TEST_NAMES "")
   foreach(test_dir IN LISTS TEST_DIRECTORIES)
     get_property(
@@ -18,14 +20,38 @@ elseif(${TESTS})
     list(APPEND TEST_NAMES ${LOCAL_TEST_NAMES})
   endforeach()
   printchoices(
-    "Test cases"
+    "Tests"
+    "TESTS"
+    "${ON_OFF_VALUES}"
+    "ON"
+    "OFF"
+    "${Green}"
+    TESTS_REPORT_1
+    46)
+  printchoices(
+    ""
     ""
     "${TEST_NAMES}"
     ""
     "${ColorReset}"
     ""
-    TESTS_REPORT
+    TESTS_REPORT_2
     0)
+  # remove only first line of TESTS_REPORT_2
+  string(REPLACE "\n" ";" TESTS_REPORT_2_LIST "${TESTS_REPORT_2}")
+  list(REMOVE_AT TESTS_REPORT_2_LIST 0)
+  string(REPLACE ";" "\n" TESTS_REPORT_2 "${TESTS_REPORT_2_LIST}")
+  set(TESTS_REPORT "${TESTS_REPORT_1}\n${TESTS_REPORT_2}")
+else()
+  printchoices(
+    "Tests"
+    "TESTS"
+    "${ON_OFF_VALUES}"
+    "OFF"
+    "OFF"
+    "${Green}"
+    TESTS_REPORT
+    46)
 endif()
 
 printchoices(
@@ -120,9 +146,8 @@ string(APPEND REPORT_TEXT ${DASHED_LINE_SYMBOL} "\n" "Configurations" "\n")
 
 if(${PGEN_FOUND})
   string(APPEND REPORT_TEXT "  " ${PGEN_REPORT} "\n")
-elseif(${TESTS})
-  string(APPEND REPORT_TEXT "  " ${TESTS_REPORT} "\n")
 endif()
+string(APPEND REPORT_TEXT "  " ${TESTS_REPORT} "\n")
 
 string(
   APPEND

diff --git a/cmake/tests.cmake b/cmake/tests.cmake
@@ -3,13 +3,6 @@ enable_testing()
 
 set(SRC_DIR ${CMAKE_CURRENT_SOURCE_DIR}/src)
 
-add_subdirectory(${SRC_DIR}/global ${CMAKE_CURRENT_BINARY_DIR}/global)
-add_subdirectory(${SRC_DIR}/metrics ${CMAKE_CURRENT_BINARY_DIR}/metrics)
-add_subdirectory(${SRC_DIR}/kernels ${CMAKE_CURRENT_BINARY_DIR}/kernels)
-add_subdirectory(${SRC_DIR}/archetypes ${CMAKE_CURRENT_BINARY_DIR}/archetypes)
-add_subdirectory(${SRC_DIR}/framework ${CMAKE_CURRENT_BINARY_DIR}/framework)
-add_subdirectory(${SRC_DIR}/output ${CMAKE_CURRENT_BINARY_DIR}/output)
-
 set(TEST_DIRECTORIES "")
 
 if(NOT ${mpi})

diff --git a/dev/nix/kokkos.nix b/dev/nix/kokkos.nix
@@ -7,26 +7,26 @@
 
 let
   name = "kokkos";
-  pversion = "5.0.0";
+  pversion = "5.0.1";
   compilerPkgs = {
     "HIP" = with pkgs.rocmPackages; [
-      llvm.llvm
+      clang
       rocm-core
       clr
       rocthrust
       rocprim
       rocminfo
       rocm-smi
-      pkgs.llvmPackages_19.clang-tools
+      pkgs.clang-tools
     ];
     "CUDA" = with pkgs.cudaPackages; [
-      llvmPackages_19.clang-tools
+      pkgs.clang-tools
       cudatoolkit
       cuda_cudart
       pkgs.gcc13
     ];
     "NONE" = [
-      pkgs.llvmPackages_19.clang-tools
+      pkgs.clang-tools
       pkgs.gcc13
     ];
   };
@@ -41,7 +41,7 @@ let
       "-D Kokkos_ENABLE_HIP=ON"
       "-D Kokkos_ARCH_${getArch { }}=ON"
       "-D GPU_TARGETS=${builtins.replaceStrings [ "amd_" ] [ "" ] (pkgs.lib.toLower (getArch { }))}"
-      "-D CMAKE_CXX_COMPILER=hipcc"
+      "-D CMAKE_CXX_COMPILER=clang++"
     ];
     "CUDA" = [
       "-D Kokkos_ENABLE_CUDA=ON"
@@ -57,7 +57,7 @@ pkgs.stdenv.mkDerivation rec {
   src = pkgs.fetchgit {
     url = "https://github.com/kokkos/kokkos/";
     rev = "${pversion}";
-    sha256 = "sha256-C4DarqnEcdF3+19TPqcM0A9bcQSkKTJkB8b7OkzC7T8=";
+    sha256 = "sha256-ChpwGBwE7sNovjdAM/iCeOqqwGufKxAh5vQ3qK6aFBU=";
   };
 
   nativeBuildInputs = with pkgs; [

diff --git a/dev/nix/shell.nix b/dev/nix/shell.nix
@@ -29,8 +29,8 @@ let
         CC = "gcc";
       };
       HIP = {
-        CXX = "hipcc";
-        CC = "hipcc";
+        CXX = "clang++";
+        CC = "clang";
       };
       CUDA = { };
     };

diff --git a/input.example.toml b/input.example.toml
@@ -230,6 +230,18 @@
       #   @required
       photon_species = ""
 
+      # @inferred:
+      # - nominal_probability
+      #     @brief: Nominal probability of the emission for a particle with `gamma * beta = 1`, charge-to-mass = `q0 / m0`
+      #     @type: float
+      #     @from: `.gamma_qed`, `.photon_weight`, `...drag.synchrotron.gamma_rad`, `scales.omegaB0`, `algorithms.timestep.dt`
+      #     @value: `0.1 * omegaB0 * dt * (gamma_qed / gamma_rad)^2 / photon_weight`
+      # - nominal_photon_energy
+      #     @brief: Nominal energy of the emitted photon for a particle with `gamma * beta = 1`, mass = `m0`
+      #     @type: float
+      #     @from: `.gamma_qed`
+      #     @value: `(1 / gamma_qed)^2`
+
     [radiation.emission.compton]
       # Gamma-factor of a particle emitting inverse Compton photons at energy `m0 c^2` in fiducial magnetic field `B0`
       #   @type: float [> 1.0]
@@ -248,6 +260,18 @@
       #   @required
       photon_species = ""
 
+      # @inferred:
+      # - nominal_probability
+      #     @brief: Nominal probability of the emission for a particle with `gamma * beta = 1`, charge-to-mass = `q0 / m0`
+      #     @type: float
+      #     @from: `.gamma_qed`, `.photon_weight`, `...drag.compton.gamma_rad`, `scales.omegaB0`, `algorithms.timestep.dt`
+      #     @value: `0.1 * omegaB0 * dt * (gamma_qed / gamma_rad)^2 / photon_weight`
+      # - nominal_photon_energy
+      #     @brief: Nominal energy of the emitted photon for a particle with `gamma * beta = 1`, mass = `m0`
+      #     @type: float
+      #     @from: `.gamma_qed`
+      #     @value: `(1 / gamma_qed)^2`
+
 [algorithms]
   # Number of current smoothing passes
   #   @type: ushort [>= 0]

diff --git a/pgens/accretion/pgen.hpp b/pgens/accretion/pgen.hpp
@@ -239,7 +239,7 @@ namespace user {
 
     inline void InitPrtls(Domain<S, M>& local_domain) {
       const auto energy_dist  = arch::Maxwellian<S, M>(local_domain.mesh.metric,
-                                                      local_domain.random_pool,
+                                                      local_domain.random_pool(),
                                                       temperature);
       const auto spatial_dist = PointDistribution<S, M>(xi_min,
                                                         xi_max,
@@ -260,7 +260,7 @@ namespace user {
 
     void CustomPostStep(std::size_t, long double time, Domain<S, M>& local_domain) {
       const auto energy_dist  = arch::Maxwellian<S, M>(local_domain.mesh.metric,
-                                                      local_domain.random_pool,
+                                                      local_domain.random_pool(),
                                                       temperature);
       const auto spatial_dist = PointDistribution<S, M>(xi_min,
                                                         xi_max,

diff --git a/pgens/reconnection/pgen.hpp b/pgens/reconnection/pgen.hpp
@@ -216,7 +216,7 @@ namespace user {
 
       // current layer
       auto       edist_cs = arch::Maxwellian<S, M>(local_domain.mesh.metric,
-                                             local_domain.random_pool,
+                                             local_domain.random_pool(),
                                              cs_temperature,
                                                    { ZERO, ZERO, cs_drift_u });
       const auto sdist_cs = CurrentLayer<S, M>(local_domain.mesh.metric,
@@ -243,7 +243,7 @@ namespace user {
       }
 
       const auto energy_dist = arch::Maxwellian<S, M>(domain.mesh.metric,
-                                                      domain.random_pool,
+                                                      domain.random_pool(),
                                                       bg_temperature);
 
       const auto dx = domain.mesh.metric.template sqrt_h_<1, 1>({});

diff --git a/pgens/turbulence/pgen.hpp b/pgens/turbulence/pgen.hpp
@@ -401,7 +401,7 @@ namespace user {
 
       // particle escape (resample velocities)
       const auto energy_dist = arch::Maxwellian<S, M>(domain.mesh.metric,
-                                                      domain.random_pool,
+                                                      domain.random_pool(),
                                                       temperature);
       for (const auto& sp : { 0, 1 }) {
         if (domain.species[sp].npld_r() > 1) {

diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt
@@ -24,15 +24,6 @@
 
 set(SRC_DIR ${CMAKE_CURRENT_SOURCE_DIR})
 
-# dependencies
-add_subdirectory(${SRC_DIR}/global ${CMAKE_CURRENT_BINARY_DIR}/global)
-add_subdirectory(${SRC_DIR}/metrics ${CMAKE_CURRENT_BINARY_DIR}/metrics)
-add_subdirectory(${SRC_DIR}/kernels ${CMAKE_CURRENT_BINARY_DIR}/kernels)
-add_subdirectory(${SRC_DIR}/archetypes ${CMAKE_CURRENT_BINARY_DIR}/archetypes)
-add_subdirectory(${SRC_DIR}/framework ${CMAKE_CURRENT_BINARY_DIR}/framework)
-add_subdirectory(${SRC_DIR}/engines ${CMAKE_CURRENT_BINARY_DIR}/engines)
-add_subdirectory(${SRC_DIR}/output ${CMAKE_CURRENT_BINARY_DIR}/output)
-
 set(ENTITY ${PROJECT_NAME}.xc)
 set(SOURCES ${SRC_DIR}/entity.cpp)
 

diff --git a/src/archetypes/particle_injector.h b/src/archetypes/particle_injector.h
@@ -257,7 +257,7 @@ namespace arch {
                              energy_dists.first,
                              energy_dists.second,
                              ONE / params.template get<real_t>("scales.V0"),
-                             domain.random_pool));
+                             domain.random_pool()));
       domain.species[species.first - 1].set_npart(
         domain.species[species.first - 1].npart() + nparticles);
       domain.species[species.second - 1].set_npart(
@@ -372,7 +372,7 @@ namespace arch {
         energy_dists.second,
         spatial_dist,
         ONE / params.template get<real_t>("scales.V0"),
-        domain.random_pool);
+        domain.random_pool());
       Kokkos::parallel_for("InjectNonUniformNumberDensity",
                            cell_range,
                            injector_kernel);

diff --git a/src/archetypes/tests/pgen.cpp b/src/archetypes/tests/pgen.cpp
@@ -67,7 +67,7 @@ struct CustomPgen : public arch::ProblemGenerator<S, M> {
     return init_flds;
   }
 
-  auto FixFieldsConst(const bc_in&, const ntt::em&) const
+  auto FixFieldsConst(simtime_t, const bc_in&, ntt::em) const
     -> std::pair<real_t, bool> {
     return { ZERO, false };
   }

diff --git a/src/archetypes/traits.h b/src/archetypes/traits.h
@@ -114,11 +114,12 @@ namespace arch {
       };
 
       template <class PG>
-      concept HasFixFieldsConst = requires(const PG&      pgen,
-                                           const bc_in&   bc,
-                                           const ntt::em& comp) {
+      concept HasFixFieldsConst = requires(const PG&    pgen,
+                                           simtime_t    time,
+                                           const bc_in& bc,
+                                           ntt::em      comp) {
         {
-          pgen.FixFieldsConst(bc, comp)
+          pgen.FixFieldsConst(time, bc, comp)
         } -> std::convertible_to<std::pair<real_t, bool>>;
       };
 

diff --git a/src/archetypes/utils.h b/src/archetypes/utils.h
@@ -52,11 +52,11 @@ namespace arch {
     const auto temperature_2 = temperatures.second / mass_2;
 
     const auto maxwellian_1 = arch::Maxwellian<S, M>(domain.mesh.metric,
-                                                     domain.random_pool,
+                                                     domain.random_pool(),
                                                      temperature_1,
                                                      drift_four_vels.first);
     const auto maxwellian_2 = arch::Maxwellian<S, M>(domain.mesh.metric,
-                                                     domain.random_pool,
+                                                     domain.random_pool(),
                                                      temperature_2,
                                                      drift_four_vels.second);
 

diff --git a/src/engines/engine.hpp b/src/engines/engine.hpp
@@ -120,6 +120,13 @@ namespace ntt {
     virtual void step_forward(timer::Timers&, Domain<S, M>&) = 0;
 
     void run();
+
+    auto engineParams() const -> prm::Parameters {
+      auto parameters = prm::Parameters {};
+      parameters.set("dt", static_cast<real_t>(dt));
+      parameters.set("time", static_cast<simtime_t>(time));
+      return parameters;
+    }
   };
 
   template <SimEngine::type S, class M>