26 #ifndef O2SCL_GOGNY_EOS_H 27 #define O2SCL_GOGNY_EOS_H 30 #include <o2scl/constants.h> 31 #include <o2scl/eos_had_base.h> 32 #include <o2scl/table3d.h> 34 #ifndef DOXYGEN_NO_O2NS 65 O2SCL_ERR(
"No data loaded in eos_had_gogny::calc_e().",
70 double asym=(ne.
n-pr.
n)/(ne.
n+pr.
n);
78 double E=t3d.
interp(rho,asym,
"enneut")/hc;
79 th.
ed=E*rho+ne.
n*ne.
m+pr.
n*pr.
m;
82 double dEdrho=t3d.
deriv_x(rho,asym,
"enneut")/hc;
83 double dEda=t3d.
deriv_y(rho,asym,
"enneut")/hc;
84 ne.
mu=ne.
m+E+ne.
n*(dEdrho+dEda*2.0*pr.
n/rho/rho);
85 pr.
mu=pr.
m+E+pr.
n*(dEdrho-dEda*2.0*ne.
n/rho/rho);
95 #ifndef DOXYGEN_NO_O2NS
double deriv_x(double x, double y, std::string name)
table3d t3d
The original EOS data.
A hadronic EOS based on a function of the densities [abstract base].
double deriv_y(double x, double y, std::string name)
virtual int calc_e(fermion &ne, fermion &pr, thermo &th)
Equation of state as a function of density.
double interp(double x, double y, std::string name)