dynamicfluidfinder.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2022-2024
 *      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()},
        fluid_cells_{ic_params.num_fluid_cells.value()},
        fluidizable_processes_{ic_params.fluidizable_processes.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 &, bool) const override {
    return {};
  }
 
  bool above_threshold(const ParticleData &pdata) const;
 
  bool is_process_fluidizable(const ProcessType &type) const;
 
 private:
  const RectangularLattice<EnergyMomentumTensor> &energy_density_lattice_;
  const std::map<int32_t, double> &background_;
  const double energy_density_threshold_ = NAN;
  const double min_time_ = NAN;
  const double max_time_ = NAN;
  const double formation_time_fraction_ = NAN;
  const int fluid_cells_ = std::numeric_limits<int>::quiet_NaN();
  const FluidizableProcessesBitSet fluidizable_processes_;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_DYNAMICFLUIDFINDER_H_