update changes

CHANGES.md

@@ -3,6 +3,19 @@
   * A null screening routine was added to disable screening for any network at
     compile time. (#992)
 
+  * An option to disable the clipping of species in the VODE integration
+    was added (`integrator.do_species_clip`) (#989)
+
+  * A unit test for C++ partition functions was added (#980)
+
+  * An EOS for primordial chemistry was added (#981)
+
+  * A C++ version of the Powell's hybrid non linear system solver was
+    added (#976)
+
+  * screening in the approximate rates in pynucastro nets was fixed
+    (#978)
+
 # 22.09
 
   * An NSE solver was added (#963)
@@ -12,14 +25,15 @@
 
 # 22.08
 
-  * The subch network was renamed `subch_approx` and the subch2 network was renamed
-    `subch_full` (#947)
+  * The subch network was renamed `subch_approx` and the subch2
+    network was renamed `subch_full` (#947)
 
 # 22.07
 
-  * Two new screening formulations have been added for reaction rates, based
-    on Chugunov, DeWitt, and Yakovlev 2007 and Chugunov and DeWitt 2009.  These
-    can be used with any network by setting SCREEN_METHOD at compile time.(#887)
+  * Two new screening formulations have been added for reaction rates,
+    based on Chugunov, DeWitt, and Yakovlev 2007 and Chugunov and
+    DeWitt 2009.  These can be used with any network by setting
+    SCREEN_METHOD at compile time.(#887)
 
 # 22.06
 
@@ -33,8 +47,8 @@
 
 # 22.04
 
-  * aprox networks now use a templated C++ righthand side formulation that
-    builds the ODE system at compile time. (#802)
+  * aprox networks now use a templated C++ righthand side formulation
+    that builds the ODE system at compile time. (#802)
 
   * pynucastro networks were regenerated to take advantage of recent
     optimizations (#901)
@@ -81,7 +95,7 @@
   * With Intel compilers, logical runtime parameters in Fortran
     were not being correctly cast to int (#789)
 
-  * Simplified-SDC now works with Fortran nets (#786) 
+  * Simplified-SDC now works with Fortran nets (#786)
 
 # 21.09
 
@@ -101,8 +115,8 @@
   * the numerical Jacobian now uses a more robust one-sided difference
     algorithm (#660, #728)
 
-  * for simplified SDC, we now only integrate (rho X, rho e), and no longer
-    integrate (rho E) (#710, #712, #717)
+  * for simplified SDC, we now only integrate (rho X, rho e), and no
+    longer integrate (rho E) (#710, #712, #717)
 
   * for the NSE bailout, we can now relax the conditions needed to enter
     NSE after a failed burn (#702)
@@ -148,9 +162,10 @@
 
 # 21.03
 
-   * We now integrate internal energy (e) directly instead of integrating
-     temperature (T) for the thermodynamic evolution. T is obtained from e
-     with an EOS call when needed to evaluate the rates. (#496)
+   * We now integrate internal energy (e) directly instead of
+     integrating temperature (T) for the thermodynamic evolution. T is
+     obtained from e with an EOS call when needed to evaluate the
+     rates. (#496)
 
    * simplified-SDC can be used with the NSE table in aprox19 now
      (#423, #497)
@@ -171,8 +186,8 @@
 
    * Microphysics now requires C++17 (gcc >= 7, CUDA >= 11). (#485)
 
-   * The BS integrator was removed.  This was Fortran only, doesn't support
-     SDC integration, and not well used. (#488)
+   * The BS integrator was removed.  This was Fortran only, doesn't
+     support SDC integration, and not well used. (#488)
 
 # 20.12
 
@@ -207,8 +222,8 @@
    * A new single-zone EOS test routine was created as unit_test/eos_cell
      (#382)
 
-   * The gamma_law eos (not gamma_law_general) now fills the sound speed,
-     entropy, and derivatives for more inputs (#374)
+   * The gamma_law eos (not gamma_law_general) now fills the sound
+     speed, entropy, and derivatives for more inputs (#374)
 
    * The rprox network now has screening (#377)
 
@@ -226,8 +241,8 @@
    * The VBDF integrator was removed (#348)
 
    * VODE can now reject an internal timestep that has any abundance
-     change by more than a factor of 2, or an abundance < 0 or > 1,
-     as well as timesteps where the temperature ends up negative. (#350)
+     change by more than a factor of 2, or an abundance < 0 or > 1, as
+     well as timesteps where the temperature ends up negative. (#350)
 
 # 20.07