dynamicfluidfinder.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2022-2025
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_DYNAMICFLUIDFINDER_H_
#define SRC_INCLUDE_SMASH_DYNAMICFLUIDFINDER_H_
 
#include <limits>
#include <map>
#include <vector>
 
#include "actionfinderfactory.h"
#include "icparameters.h"
#include "input_keys.h"
 
namespace smash {
 
class DynamicFluidizationFinder : public ActionFinderInterface {
 public:
  DynamicFluidizationFinder(
      const RectangularLattice<EnergyMomentumTensor> &lattice,
      const std::map<int32_t, double> &background,
      const InitialConditionParameters &ic_params)
      : energy_density_lattice_{lattice},
        background_{background},
        energy_density_threshold_{ic_params.energy_density_threshold.value()},
        min_time_{ic_params.min_time.value()},
        max_time_{ic_params.max_time.value()},
        formation_time_fraction_{ic_params.formation_time_fraction.value()},
        smearing_kernel_at_0_{ic_params.smearing_kernel_at_0.value()},
        fluid_cells_{ic_params.num_fluid_cells.value()},
        fluidizable_processes_{ic_params.fluidizable_processes.value()},
        delay_initial_elastic_{ic_params.delay_initial_elastic.value()} {};
 
  ActionList find_actions_in_cell(
      const ParticleList &search_list, double dt, double gcell_vol,
      const std::vector<FourVector> &beam_momentum) const override;
 
  ActionList find_actions_with_neighbors(
      const ParticleList &, const ParticleList &, const double,
      const std::vector<FourVector> &) const override {
    return {};
  }
 
  ActionList find_actions_with_surrounding_particles(
      const ParticleList &, const Particles &, double,
      const std::vector<FourVector> &) const override {
    return {};
  }
 
  ActionList find_final_actions(const Particles &search_list) const override;
 
  bool above_threshold(const ParticleData &pdata) const;
 
  bool is_process_fluidizable(const HistoryData &history) const;
 
 private:
  const RectangularLattice<EnergyMomentumTensor> &energy_density_lattice_;
  const std::map<int32_t, double> &background_;
  const double energy_density_threshold_ = smash_NaN<double>;
  const double min_time_ = smash_NaN<double>;
  const double max_time_ = smash_NaN<double>;
  const double formation_time_fraction_ = smash_NaN<double>;
  const double smearing_kernel_at_0_ = smash_NaN<double>;
  const int fluid_cells_ = smash_NaN<int>;
  const FluidizableProcessesBitSet fluidizable_processes_;
  const bool delay_initial_elastic_ = true;
 
  mutable int particles_in_core_ = 0;
  mutable double energy_in_core_ = 0.;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_DYNAMICFLUIDFINDER_H_