It might be worth coding up direct differentiation, or differentiating the eff results, as these may succeed more generally.
This class will have difficulty with extremely degenerate or extremely non-degnerate systems. Fix this.
Create a more intelligent method for dealing with bad initial guesses for the chemical potential in calc_density().
Include explicit zero-temperature calculation, maybe by making this a child of fermion_zerot or by making a new fermion_deriv_zerot?
There is also a closed form for the derivatives of massless fermions with pairs at finite temperature in Constantiou et al. 2014 which could be implemented here.
Also calibrate massless fermions?
Convert into separate class?
The expressions which appear in in the integrand functions density_fun(), etc. could likely be improved, especially in the case where o2scl::part::inc_rest_mass is false
. There should not be a need to check if ret
is finite.
It appears this class doesn't compute the uncertainty in the chemical potential or density with calc_density(). This could be fixed.
I'd like to change the lower limit on the entropy integration, but the value in the code at the moment (stored in ll
) makes bm_part2.cpp worse.
The function pair_mu() should set the antiparticle integrators as done in fermion_deriv_rel.
Create a caching and more intelligent search system for the table. The table is sorted by A and then N, so we could probably just copy the search routine from mnmsk_mass, which is sorted by Z and then N (some code written for this, but it doesn't work yet).
Should m_neut and m_prot be set to the neutron and proton masses from the table by default?
Warn about malformed combinations like Carbon-5
Right now, n4
is interpreted incorrectly as Nitrogen-4, rather than the tetraneutron.
Interpret strings with the full name rather than just the abbreviation.
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).
Hosted at
.