Go to the documentation of this file.
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
157 virtual int set_parameters(
double lambda=0.0,
double fourferm=0.0,
195 double &gap2,
double &gap3,
thermo <h);
201 double &gap2,
double &gap3,
thermo <h);
207 double &gap1,
double &gap2,
double &gap3,
208 thermo &qb,
double temper);
257 double up_default_mass;
258 double down_default_mass;
259 double strange_default_mass;
296 virtual const char *
type() {
return "eos_quark_njl"; }
308 double ms, m, mu, temper, limit;
326 #ifndef DOXYGEN_INTERNAL
349 double iqq(
double x,
const njtp &pa);
351 double ide(
double x,
const njtp &pa);
353 double ied(
double x,
const njtp &pa);
355 double ipr(
double x,
const njtp &pa);
364 #ifndef DOXYGEN_NO_O2NS
void njbag(quark &q)
Calculates the contribution to the bag constant from quark q.
int gapfunms(size_t nv, const ubvector &x, ubvector &y)
Calculates gap equations in y as a function of the constituent masses in x.
int B0fun(size_t nv, const ubvector &x, ubvector &y)
Used by calc_B0() to compute the bag constant.
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 at finite temperature.
struct o2scl::eos_quark_njl::njtp_s njtp
A structure for passing parameters to the integrands.
virtual int calc_p(quark &u, quark &d, quark &s, thermo <h)
Equation of state as a function of chemical potentials.
inte_qag_gsl def_it
The default integrator.
int gapfunqqT(size_t nv, const ubvector &x, ubvector &y)
Calculates gap equations in y as a function of the quark condensates in x.
double ied(double x, const njtp &pa)
The integrand for the energy density.
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.
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 set_quarks(quark &u, quark &d, quark &s)
Set the quark objects to use.
double G
The four-fermion coupling (in )
virtual const char * type()
Return string denoting type ("eos_quark_njl")
A structure for passing parameters to the integrands.
mroot_hybrids< mm_funct, boost::numeric::ublas::vector< double >, boost::numeric::ublas::matrix< double >, jac_funct > def_solver
The default solver.
double ide(double x, const njtp &pa)
The integrand for the density.
int gapfunmsT(size_t nv, const ubvector &x, ubvector &y)
Calculates gap equations in y as a function of the constituent masses in x.
double ipr(double x, const njtp &pa)
The integrand for the pressure.
double limit
Accuracy limit for Fermi integrals for finite temperature.
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.
virtual int set_solver(mroot< mm_funct, boost::numeric::ublas::vector< double >, jac_funct > &s)
Set solver to use in set_parameters()
inte * it
The integrator for finite temperature integrals.
double K
The 't Hooft six-fermion interaction coupling (in )
Quark matter equation of state base.
std::function< int(size_t, const boost::numeric::ublas::vector< double > &, boost::numeric::ublas::vector< double > &) > mm_funct
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.
virtual int set_parameters(double lambda=0.0, double fourferm=0.0, double sixferm=0.0)
Set the parameters and the bag constant B0.
std::function< int(size_t, boost::numeric::ublas::vector< double > &, size_t, boost::numeric::ublas::vector< double > &, boost::numeric::ublas::matrix< double > &) > jac_funct
quark * strange
The strange quark.
double cp_temp
The temperature for calc_temp_p()
double iqq(double x, const njtp &pa)
The integrand for the quark condensate.
mroot< mm_funct, boost::numeric::ublas::vector< double >, jac_funct > * solver
The solver to use for set_parameters()
virtual int set_inte(inte<> &i)
Set integration object.
quark * down
The down quark.
double B0
The bag constant (in )
double L
The momentum cutoff (in )
bool fromqq
Calculate from quark condensates if true (default true)
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).