Pynucastro

* Switch to using math.fsum() instead of sum() for better python
3.12 compatibility (720)

* Some modernization of the AMReX-Astro C++ code output (709, 714,
717)

* Documentation improvements, including references and contributing
instructions, and NSE table (669, 679, 695, 697, 706, 707)

* Clean-up up some plot interfaces to always take rho, T, comp in
that order (693)

* Performance improvements (680, 685, 692)

* A new constants module was added to provide fundamental constants
(688)

* A new NumpyNetwork was added for use with the rate reduction
algorithm (684)

* Neutrino cooling routines based on Itoh et al. 1996 were added
(628)

* A new method for generating an NSE table was added (676)

* A new NSENetwork class was added for dealing with NSE (675)

* Moved some screening utilities out of RateCollection (666)

* A PythonNetwork now write out a function to compute the energy
release (640)

* `RateCollection.evaluate_energy_generation()` can now return the
neutrino losses from weak rates separately. (665)

* `AmrexAstroCxxNetwork` now capitalizes the species names to be
consistent with other AMReX-Astro Microphysics networks (652)

* The pynucastro project configuration was updated and reorganized
to pyproject.toml (657, 660, 661)

* `Library.get_rate_by_name` no longer returns None if one of
the input rate names is not found (659)

* A C++ indexing error at the boundary of a weak rate table was
fixed. (654)

* The C++ AmrexAstroCxxNetwork for evaluating weak rate neutrino
losses was cleaned up to fix a bounds error (655)

* We now explicitly check in the C++ code if a weak rate table
if found and contains data (654)

* a RateCollection now enforces that there are no duplicate rates. A
`find_duplicate_links()` method has been added to `Library` to allow
one to construct a `Library` without duplicates and then create a
`RateCollection` or other network from that. (651)

* Some documentation updates (643, 650)

* The positron rates in the Langanke weak rate tables have been
merged with the corresponding electron rates and the n <-> p rate
was added (648)

* The `Composition.bin_as()` method now takes a list of nuclei
to exclude from being binned into (642)

* `RateCollection.evaluate_ydots()` now can take a filter to
only evaluate a subset of the rates. (641)

* add `eval_zbar()` to Composition (632)

* fix `get_rate_by_name` to work with "pp" reactions (632)

* created a method to reduce a Composition from one set
of nuclei to another based on the nuclei masses and
charge number (625)

* switch AmrexAstroCxxNetwork to do bilinear interpolation
in terms of log(T) and log(rhoY) (592, 611)

* tabular rates in python now do linear interpolation (602)

* added an example of creating a custom rate (615)

* Rate now calls `_set_q()` to set the Q value if it is not
passed in (617)

* added a TableInterpolator that works both for interactive
python and PythonNetwork (612, 610, 609)

* added a RateCollection method to find duplicate links (601)

* python networks with tabular rates now copy the table files (605)

* added a `get_nuclei_in_range` method to return a range of
nuclei (593)

* we now do a binary search in the C++ partition function
interpolation (581)

* added a simple C++ network (SimpleCxxNetwork) (591, 585)

* added the weak rates from Langanke 2001 (536)

* cleaned up partition functions in C++ (578, 573, 569, 565)

* converted the Suzuki tabular rates to be in terms of log (550)

* fixed a bounds issue in C++ table interpolation (566)

* eliminated a variable length array in the C++ table interpolation
(567)

* added rate indices to the C++ networks (563)

* added a network reduction algorithm (529, 528, 527, 526, 525,
523)

* added a molar fraction method to Composition (546)

* added examples of interfacing with Julia (539)

* added a code of conduct (504)

* added gamma heating to tabular weak rates (502)

* add functions to numerically evaluate the Jacobian (467)

* Fix C++ compilation with partition functions and approximate
rates and derived rates (407, 494)

* The Explorer class now takes all of the options as
RateCollection.plot() (496)

* Rates edges can be highlighted via a user-defined function in the
RateCollection plot (486)

* All of the plot routines in RateCollection now return a
matplotlib Figure object (490)

* add get_rate_by_nuclei and add_rate to Library() (484)

* RateCollection.remove_rates() can now act on a single rate and a
new method, add_rates(), will add a rate directly to the
collection.

* Fixed printing of TabularLibrary (479)

* Documentation updates (473)

* added get_rate_by_name to RateCollection (468)

* added tabular rates to PythonNetwork (451)

* fixed NSE state for network without protons (463)

fix a setup.py install issue for the C++ templates

-- Fix some setup.py metadata for PyPI (453)

-- Change how tabular rate tables are stored and read (450)

-- Support historial ReacLib format (449)