Class eos_leptons_multip (o2scl)

O2scl : Class List

class eos_leptons_multip : public o2scl::eos_leptons

Internal particle objects [protected]

fermion_ld eld

Electron in long double precision.

fermion_cdf25 ecdf25

Electron in 25-digit precision.

fermion_ld muld

Muon in long double precision.

fermion_cdf25 mucdf25

Muon in 25-digit precision.

fermion_ld tauld

Tau in long double precision.

fermion_cdf25 taucdf25

Tau in 25-digit precision.

Unit conversion objects to set the lepton masses

convert_units<long double> cu_ld

Long double precision unit conversion object.

convert_units<cpp_dec_float_25> cu_cdf25

25-digit precision unit conversion object

fermion_rel_ld frel_ld

Relativistic fermion thermodynamics in long double precision.

fermion_rel_cdf25 frel_cdf25

Relativistic fermion thermodynamics in 25 digit precision.

virtual int electron_density(double T)

Compute electron thermodynamics from the electron density.

Because this function uses fermion_rel::pair_density(), the current electron chemical potential (from e) is used as an initial guess.

The temperature should be in units of \( 1/\mathrm{fm} \) .

Note

This internal function is designed to be used when the flag fermion_rel::err_nonconv for the frel object is set to false . If this function fails, it returns a non-zero value rather than calling the error handler. This allows it to be used inside pair_density_eq_fun(). Thus this function is protected instead of public.

virtual int pair_density_eq_fun(size_t nv, const ubvector &x, ubvector &y, double T, double nq)

Function to solve for pair_density_eq()

The temperature should be in units of \( 1/\mathrm{fm} \) .

The value x is used for the electron density or chemical potential, depending on the value of pde_from_density .

The value

\[ \frac{\left(n_e + n_{\mu} - n_q\right)}{|n_q|} \]

is stored in y.

Note

This internal function presumes that include_muons is true (otherwise there’s nothing to solve) and that the fermion_rel::err_nonconv flag in the frel object is set to false .

eos_leptons_multip()

inline virtual ~eos_leptons_multip()