Global cfl_njl_eos::calc_eq_temp_p (quark &u, quark &d, quark &s, double &qq1, double &qq2, double &qq3, double &gap1, double &gap2, double &gap3, double mu3, double mu8, double &n3, double &n8, thermo &qb, double temper)
It surprises me that n3 is not -res[11]. Is there a sign error in the color densities?
Global cfl_njl_eos::gapped_eigenvalues (double m1, double m2, double lmom, double mu1, double mu2, double tdelta, double lam[4], double dldmu1[4], double dldmu2[4], double dldm1[4], double dldm2[4], double dldg[4])
In the code, the equal mass case seems to be commented out. Why?
Global ddc_eos::calc_eq_e (fermion &neu, fermion &p, double sig, double ome, double rho, double &f1, double &f2, double &f3, thermo &th)
Is the thermodynamic identity is satisfied even when the field equations are not solved? Check this.
Class rmf_eos
Check the formulas in the "Background" section There are two calc_e() functions that solve. One is specially designed to work without a good initial guess. Possibly the other calc_e() function should be similarly designed? Make sure that this class properly handles particles for which inc_rest_mass is true/false The error handler is called sometimes when calc_e() is used to compute pure neutron matter. This should be fixed.
Global rmf_eos::fix_saturation (double guess_cs=4.0, double guess_cw=3.0, double guess_b=0.001, double guess_c=-0.001)
Fix this for zm_mode=true Ensure solver is more robust
Global rmf_eos::fkprime_fields (double sig, double ome, double nb, double &k, double &kprime)
Does this work? Fix fkprime_fields() if it does not.
Global rmf_eos::n_charge
Should use hadronic_eos::proton_frac instead?
Global rmf_eos::check_naturalness (rmf_eos &re)
I may have ignored some signs in the above, which are unimportant for this application, but it would be good to fix them for posterity.
Global schematic_eos::set_a_from_mstar (double u_msom, double mnuc)
This was computed in schematic_sym.nb, which might be added to the documentation?
Class skyrme_eos
Make sure that this class properly handles particles for which inc_rest_mass is true/false What about the spin-orbit units? Need to write a function that calculates saturation density? Remove use of mnuc in calparfun()? The compressibility could probably use some simplification Make sure the finite-temperature part is properly tested The testing code doesn't work if err_mode is 2, probably because of problems in load().
Global skyrme_eos::calpar (double gt0=-10.0, double gt3=70.0, double galpha=0.2, double gt1=2.0, double gt2=-1.0)
Does this work for both 'a' and 'b' non-zero?
Global skyrme_eos::calpar (double gt0=-10.0, double gt3=70.0, double galpha=0.2, double gt1=2.0, double gt2=-1.0)
Compare to similar formulae from Margueron02
Global skyrme_eos::landau_neutron (double n0, double m, double &f0, double &g0, double &f1, double &g1)
This needs to be checked
Global skyrme_eos::landau_nuclear (double n0, double m, double &f0, double &g0, double &f0p, double &g0p, double &f1, double &g1, double &f1p, double &g1p)
This needs to be checked.
Class tov_interp_eos
Warn that the pressure in the low-density eos is not strictly increasing! (see at P=4.3e-10)
Class tov_solve
baryon mass doesn't work for fixed() (This may be fixed. We should make sure it's tested.) Combine maxoutsize and kmax? Document column naming issues Document surface gravity and redshift
Documentation generated with Doxygen and provided under the
GNU Free Documentation License. See License Information for details.
Project hosting provided by
,
O2scl Sourceforge Project Page