chemicalpotential.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2020-
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
 
#ifndef SRC_INCLUDE_SMASH_CHEMICALPOTENTIAL_H_
#define SRC_INCLUDE_SMASH_CHEMICALPOTENTIAL_H_
 
#include <gsl/gsl_multiroots.h>
#include <gsl/gsl_vector.h>
 
#include <smash/integrate.h>
 
namespace smash {
 
class ChemicalPotentialSolver {
 public:
  static double density_one_species(double degeneracy, double mass,
                                    double temperature,
                                    double effective_chemical_potential,
                                    double statistics, Integrator* integrator);
 
  struct ParametersForChemPotSolver {
    double degeneracy;
 
    double mass;
 
    double number_density;
 
    double temperature;
 
    double statistics;
 
    Integrator* integrator;
  };
 
  static double root_equation_effective_chemical_potential(
      double degeneracy, double mass, double number_density, double temperature,
      double effective_chemical_potential, double statistics,
      Integrator* integrator);
 
  static int root_equation_effective_chemical_potential_for_GSL(
      const gsl_vector* roots_array, void* parameters, gsl_vector* function);
 
  static void print_state_effective_chemical_potential(
      unsigned int iter, gsl_multiroot_fsolver* solver);
 
  static int find_effective_chemical_potential(
      double degeneracy, double mass, double number_density, double temperature,
      double statistics, double mu_initial_guess, double solution_precision,
      Integrator* integrator, double* effective_chemical_potential);
 
  double effective_chemical_potential(double degeneracy, double mass,
                                      double number_density, double temperature,
                                      double statistics,
                                      double solution_precision);
 
 private:
  Integrator integrator_ = Integrator(1E6);
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_CHEMICALPOTENTIAL_H_