potentials.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2014-2019
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */

#ifndef SRC_INCLUDE_POTENTIALS_H_
#define SRC_INCLUDE_POTENTIALS_H_

#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;

  std::pair<double, int> force_scale(const ParticleType &data) const;

  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_; }

 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_POTENTIALS_H_