potentials.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2014-2020
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_POTENTIALS_H_
#define SRC_INCLUDE_SMASH_POTENTIALS_H_
 
#include <tuple>
#include <utility>
#include <vector>
 
#include "configuration.h"
#include "density.h"
#include "forwarddeclarations.h"
#include "particledata.h"
#include "threevector.h"
 
namespace smash {
 
class Potentials {
 public:
  Potentials(Configuration conf, const DensityParameters &parameters);
  virtual ~Potentials();
 
  double skyrme_pot(const double baryon_density) const;
 
  double symmetry_pot(const double baryon_isospin_density,
                      const double baryon_density) const;
 
  double symmetry_S(const double baryon_density) const;
 
  double potential(const ThreeVector &r, const ParticleList &plist,
                   const ParticleType &acts_on) const;
 
  static std::pair<double, int> force_scale(const ParticleType &data);
 
  std::pair<ThreeVector, ThreeVector> skyrme_force(
      const double density, const ThreeVector grad_rho, const ThreeVector dj_dt,
      const ThreeVector rot_j) const;
 
  std::pair<ThreeVector, ThreeVector> symmetry_force(
      const double rhoI3, const ThreeVector grad_rhoI3,
      const ThreeVector djI3_dt, const ThreeVector rot_jI3, const double rhoB,
      const ThreeVector grad_rhoB, const ThreeVector djB_dt,
      const ThreeVector rot_jB) const;
  virtual std::tuple<ThreeVector, ThreeVector, ThreeVector, ThreeVector>
  all_forces(const ThreeVector &r, const ParticleList &plist) const;
 
  virtual bool use_skyrme() const { return use_skyrme_; }
  virtual bool use_symmetry() const { return use_symmetry_; }
 
  double skyrme_a() const { return skyrme_a_; }
  double skyrme_b() const { return skyrme_b_; }
  double skyrme_tau() const { return skyrme_tau_; }
 
 private:
  const DensityParameters param_;
 
  bool use_skyrme_;
 
  bool use_symmetry_;
 
  double skyrme_a_;
 
  double skyrme_b_;
 
  double skyrme_tau_;
 
  double symmetry_S_Pot_;
 
  bool symmetry_is_rhoB_dependent_ = false;
  double symmetry_gamma_;
 
  double dVsym_drhoI3(const double rhoB, const double rhoI3) const;
 
  double dVsym_drhoB(const double rhoB, const double rhoI3) const;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_POTENTIALS_H_