26 #ifndef O2SCL_NAMBUJL_EOS_H 27 #define O2SCL_NAMBUJL_EOS_H 30 #include <o2scl/constants.h> 31 #include <o2scl/eos_quark.h> 32 #include <o2scl/quark.h> 33 #include <o2scl/mm_funct.h> 34 #include <o2scl/mroot_hybrids.h> 35 #include <o2scl/inte.h> 36 #include <o2scl/inte_qag_gsl.h> 38 #ifndef DOXYGEN_NO_O2NS 152 virtual int set_parameters(
double lambda=0.0,
double fourferm=0.0,
190 double &gap2,
double &gap3,
thermo <h);
196 double &gap2,
double &gap3,
thermo <h);
202 double &gap1,
double &gap2,
double &gap3,
203 thermo &qb,
double temper);
212 int gapfunms(
size_t nv,
const ubvector &x, ubvector &y);
221 int gapfunqq(
size_t nv,
const ubvector &x, ubvector &y);
230 int gapfunmsT(
size_t nv,
const ubvector &x, ubvector &y);
239 int gapfunqqT(
size_t nv,
const ubvector &x, ubvector &y);
249 double up_default_mass;
250 double down_default_mass;
251 double strange_default_mass;
286 virtual const char *
type() {
return "eos_quark_njl"; }
298 double ms, m, mu, temper, limit;
316 #ifndef DOXYGEN_INTERNAL 326 int B0fun(
size_t nv,
const ubvector &x, ubvector &y);
339 double iqq(
double x,
const njtp &pa);
341 double ide(
double x,
const njtp &pa);
343 double ied(
double x,
const njtp &pa);
345 double ipr(
double x,
const njtp &pa);
354 #ifndef DOXYGEN_NO_O2NS virtual const char * type()
Return string denoting type ("eos_quark_njl")
inte * it
The integrator for finite temperature integrals.
int gapfunms(size_t nv, const ubvector &x, ubvector &y)
Calculates gap equations in y as a function of the constituent masses in x.
virtual int set_inte(inte<> &i)
Set integration object.
void njbag(quark &q)
Calculates the contribution to the bag constant from quark q.
std::function< int(size_t, const boost::numeric::ublas::vector< double > &, boost::numeric::ublas::vector< double > &) > mm_funct11
quark * strange
The strange quark.
int set_quarks(quark &u, quark &d, quark &s)
Set the quark objects to use.
double limit
Accuracy limit for Fermi integrals for finite temperature.
quark * down
The down quark.
struct o2scl::eos_quark_njl::njtp_s njtp
A structure for passing parameters to the integrands.
Nambu Jona-Lasinio EOS at zero temperature.
double G
The four-fermion coupling.
Quark matter equation of state base.
double K
The 't Hooft six-fermion interaction coupling.
double cp_temp
The temperature for calc_temp_p()
inte_qag_gsl def_it
The default integrator.
std::function< int(size_t, boost::numeric::ublas::vector< double > &, size_t, boost::numeric::ublas::vector< double > &, boost::numeric::ublas::matrix< double > &) > jac_funct11
double B0
The bag constant.
double ide(double x, const njtp &pa)
The integrand for the density.
virtual int calc_temp_p(quark &u, quark &d, quark &s, double T, thermo &th)
Equation of state as a function of chemical potentials at finite temperature.
double L
The momentum cutoff.
virtual int calc_eq_e(quark &u, quark &d, quark &s, double &gap1, double &gap2, double &gap3, thermo <h)
Equation of state and gap equations as a function of the densities.
double ipr(double x, const njtp &pa)
The integrand for the pressure.
virtual int calc_p(quark &u, quark &d, quark &s, thermo <h)
Equation of state as a function of chemical potentials.
virtual int set_solver(mroot< mm_funct11, boost::numeric::ublas::vector< double >, jac_funct11 > &s)
Set solver to use in set_parameters()
double iqq(double x, const njtp &pa)
The integrand for the quark condensate.
int gapfunqqT(size_t nv, const ubvector &x, ubvector &y)
Calculates gap equations in y as a function of the quark condensates in x.
A structure for passing parameters to the integrands.
virtual int set_parameters(double lambda=0.0, double fourferm=0.0, double sixferm=0.0)
Set the parameters and the bag constant B0.
bool fromqq
Calculate from quark condensates if true (default true)
int B0fun(size_t nv, const ubvector &x, ubvector &y)
Used by calc_B0() to compute the bag constant.
mroot< mm_funct11, boost::numeric::ublas::vector< double >, jac_funct11 > * solver
The solver to use for set_parameters()
double ied(double x, const njtp &pa)
The integrand for the energy density.
int gapfunqq(size_t nv, const ubvector &x, ubvector &y)
Calculates gap equations in y as a function of the quark condensates in x.
int gapfunmsT(size_t nv, const ubvector &x, ubvector &y)
Calculates gap equations in y as a function of the constituent masses in x.
mroot_hybrids< mm_funct11, boost::numeric::ublas::vector< double >, boost::numeric::ublas::matrix< double >, jac_funct11 > def_solver
The default solver.
int calc_eq_temp_p(quark &tu, quark &td, quark &ts, double &gap1, double &gap2, double &gap3, thermo &qb, double temper)
Equation of state and gap equations as a function of chemical potentials.
virtual int calc_eq_p(quark &u, quark &d, quark &s, double &gap1, double &gap2, double &gap3, thermo <h)
Equation of state and gap equations as a function of chemical potential.