Todo List
Member o2scl::eos_had_ddc::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 o2scl::eos_had_rmf
  • The functions fcomp_fields(), fkprime_fields(), and fesym_fields() are not quite correct if the neutron and proton masses are different. For this reason, they are currently unused by saturation().
  • The fix_saturation() and calc_cr() functions use mnuc, and should be modified to allow different neutron and proton masses.
  • Check the formulas in the "Background" section
  • Make sure that this class properly handles particles for which inc_rest_mass is true/false
  • The calc_e() function fails to converge at lower densities. See the testing code which has trouble with NL3 and RAPR
Member o2scl::eos_had_rmf::check_naturalness (eos_had_rmf &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.
Member o2scl::eos_had_rmf::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
Member o2scl::eos_had_rmf::fkprime_fields (double sig, double ome, double nb, double &k, double &kprime)
This function, o2scl::eos_had_rmf::fkprime_fields() is currently untested.
Member o2scl::eos_had_schematic::baryon_suscep (double n, double x)
This function is untested.
Member o2scl::eos_had_schematic::set_a_from_mstar (double u_msom, double mnuc)
This was computed in schematic_sym.nb, which might be added to the documentation?
Class o2scl::eos_had_skyrme
  • Convert W0 to b4 and b4p everywhere
  • Remove use of mnuc in calparfun()?
  • Document o2scl_hdf::skyrme_load() file format.
  • Update reference list.
Member o2scl::eos_had_skyrme::alt_params_set (double Crr00, double Crr10, double Crr0D, double Crr1D, double Crt0, double Crt1, double CrDr0, double CrDr1, double CrnJ0, double CrnJ1, double alpha2)
These expressions are not exactly the same as those in Bender03, so I need to find out why and make this more clear.
Member o2scl::eos_had_skyrme::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?

Compare to similar formulas in Margueron02

Member o2scl::eos_had_skyrme::landau_neutron (double n0, double m, double &f0, double &g0, double &f1, double &g1)
This needs to be checked
Member o2scl::eos_had_skyrme::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 o2scl::eos_nse_full
I don't think inc_lept_phot=false works because then all WS cells have infinite size because of no electrons. For the moment, this variable is protected to discourage the user from changing it.
Class o2scl::eos_nstar_rot_interp
Document what the nstar_rot code requirements are for the low-density part of the EOS.
Member o2scl::eos_quark_cfl::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?
Member o2scl::eos_quark_cfl::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?
Class o2scl::eos_sn_ls
There are still a few points for which the electron/photon EOS seems to be off, but this may be the result of small inaccuracies from finite-differencing the LS table.
Class o2scl::eos_sn_stos
Add the T=0 and Ye=0 data to this class. Separate tables for these cases have been released, but I don't think this class can read them yet.
Class o2scl::eos_tov_buchdahl
Fix base EOS functions
Class o2scl::nstar_cold
Maybe calc_eos() doesn't really need a thermo object since it can use the default EOS thermo object and we can get rid of the class variables h, hb, and l
Class o2scl::nstar_rot

Better documentation is needed everywhere.

Test the resize() function

Class o2scl::nucleus_rmf

Better documentation

Convert energies() to use EOS and possibly replace sigma_rhs() and related functions by the associated field equation method of eos_had_rmf.

Document hw=3.923+23.265/cbrt(atot);

I believe currently the center of mass correction for the binding energy per nucleon is not done and has to be added in after the fact

Class o2scl::nucmass_ldrop
12/4/14: This doesn't gracefully handle negative values of n0 and n1 as then the neutron and proton densities become negative. This needs to be addressed. For now, there is a fix at line 246 in nucmass_ldrop.cpp .
Class o2scl::nucmass_ldrop_skin

This is based on LPRL, but it's a little different in Lattimer and Swesty. I should document what the difference is.

The testing could be updated.

Documentation generated with Doxygen. Provided under the GNU Free Documentation License (see License Information).