All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Pages
Todo List
Global 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 error handler is called sometimes when calc_e() is used to compute pure neutron matter. This should be fixed.
Global 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.
Global 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
Global 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.
Global o2scl::eos_had_schematic::baryon_suscep (double n, double x)
This function is untested.
Global 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
  • Make sure that this class properly handles particles for which inc_rest_mass is true/false
  • 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().
  • Document load() file format.
  • Update reference list.
Global 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 formulae from Margueron02

Global o2scl::eos_had_skyrme::landau_neutron (double n0, double m, double &f0, double &g0, double &f1, double &g1)
This needs to be checked
Global 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.
Global 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?
Global 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_interp

It might be useful to exit more gracefully when non-finite values are obtained in interpolation, analogous to the err_nonconv mechanism elsewhere.

Create a sanity check where core_auxp is nonzero only if core_table is also nonzero. Alternatively, this complication is due to the fact that this class works in two ways: one where it reads a table (and adds a crust), and one where it reads in vectors (with no crust). Maybe these two modes of operation should be separated into two classes.

Read different vector types than just std::vector<>.

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_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).
Hosted at Get Object-oriented Scientific Computing
Lib at SourceForge.net. Fast, secure and Free Open Source software
downloads..