Go to the documentation of this file.
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;
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);
88 th.
pr=-th.
ed+ne.n*ne.mu+pr.n*pr.mu;
95 #ifndef DOXYGEN_NO_O2NS
double deriv_x(double x, double y, std::string name) const
double deriv_y(double x, double y, std::string name) const
table3d t3d
The original EOS data.
double interp(double x, double y, std::string name) const
virtual int calc_e(fermion &ne, fermion &pr, thermo &th)
Equation of state as a function of density.
A hadronic EOS based on a function of the densities [abstract base].
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).