26 #ifndef O2SCL_NSE_EOS_H
27 #define O2SCL_NSE_EOS_H
29 #include <o2scl/classical.h>
30 #include <o2scl/constants.h>
31 #include <o2scl/nucdist.h>
32 #include <o2scl/mm_funct.h>
33 #include <o2scl/mroot_hybrids.h>
35 #ifndef DOXYGEN_NO_O2NS
72 #ifndef DOXYGEN_INTERNAL
79 std::vector<nucleus> *ndp;
113 void calc_mu(
double mun,
double mup,
double T,
114 double &nb,
double &Ye,
thermo &th, std::vector<nucleus> &nd);
127 int calc_density(
double nb,
double Ye,
double T,
double &mun,
128 double &mup,
thermo &th, std::vector<nucleus> &nd);
142 #ifndef DOXYGEN_NO_O2NS
void calc_mu(double mun, double mup, double T, double &nb, double &Ye, thermo &th, std::vector< nucleus > &nd)
Calculate the equation of state as a function of the chemical potentials.
bool err_nonconv
If true, call the error handler if calc_density() does not converge (default true) ...
mroot< mm_funct11, ubvector, jac_funct11 > * root
Solver.
mroot_hybrids< mm_funct11, ubvector, ubmatrix, jac_funct11 > def_root
Default solver.
void set_mroot(mroot< mm_funct11, ubvector, jac_funct11 > &rp)
Set the solver for use in computing the chemical potentials.
int solve_fun(size_t nv, const ubvector &x, ubvector &y, solve_parms &pa)
Function to solve for baryon and charge conservation.
classical cla
Compute particle properties assuming classical thermodynamics.
int calc_density(double nb, double Ye, double T, double &mun, double &mup, thermo &th, std::vector< nucleus > &nd)
Calculate the equation of state as a function of the densities.
Equation of state for nuclei in statistical equilibrium.