%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                       HISTORY of PRIMAT VERSIONS
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0.1.0 Release version (February 2018). PRIMAT is described in arXiv:1801.08023, or Physics Reports, 04, (2018) 005.


0.1.1 Correction of incorrect implementation of analytic reactions (for piecewise functions)


0.1.2 (May 2019)
- Reimplementation of the way analytic functions must be given. Simpler and more compact and more transparent.
- Removes the possibility of setting the maximum number of equations. Now the equations used are all of the external file, plus the analytic ones.
- Possibility of writing a reaction with identical species as a+a or 2a in the external file.
- Possibility of adding global factors to some selected rates in order to assess how a reevaluation of a rate could affect final abundances.


0.2.0 (Nov 2020)
- New dpg (d+p <-> He3 + g) rate from LUNA experiment, Mossa et al. 2020, and included in the reaction file BBNRatesAC2021.dat. 
- Reaction Be7(d,p) is now in the tabulated file and not analytically given.
- Decoupling of neutrinos with oscillations (precomputed with NEVO, see [Froustey et al. 2020]).
- Constants reevaluated with PDG2020 data (notably G_Newton).
- QED corrections to some nuclear rates (see Pitrou&Pospelov 2020).
- A new option to use only reduced nuclear network of 12 nuclear reactions (much faster)
- QED at order (e^3) included in plasma thermodynamics (reduces Neff by 0.001, see e.g. [Bennett et al. 2019])
- be7np DAACV04 replaced by de Souza et al. 2020
- tdn & hdp DAACV04 replaced by de Souza et al. 2019a,b
- be7da from Rijal et al. 2019 added in table & removed from notebook

0.2.1 
-New dpg (d + p <-> He3 + g) rate (Moscoso et al. 2021 ) in new reaction file BBNRatesAC2022.dat. 

0.2.2:
-small modification to make it compatible with Mathematica 13.
-neutron lifetime updated with PDG 2023.

0.2.3 (Sept 2024)
- Possibility of adding spectral distortions in analytic form (of the mu-type and Y-type). This is beyond the standard model :
 Mu type is not really necessary since a Fermi-Dirac with chemical potential can be chosen when asking for DegenerateNeutrinos. 
 Y-type distortion is in principle not possible because neutrino scatter before decoupling and reach near equilibrium.

******** major changes **********

0.3.0 (Oct 2024)
- PRIMAT now looks more like a package even though it is not yet a Package. 
When loading PRIMAT-Main.m from an external notebook, this now only loads definitions (and takes 0.3 seconds), but no computation is performed. The user can then modify the various options about physics or numerics, and then call for the function which either perform one run (RunNumericalIntegrals) or a MonteCarlo estimation of errors (RunPRIMATMonteCarlo), and then read the results. Examples are provided in the Examples folder. It has been tested under a variety of combination of options, but please reports situations in which this still generates bugs. 
- It now takes less than a second for the small network of reactions (12 nuclear reactions) and 16 seconds for the full network (423 reactions), provided the weak rates are precomputed and stored in the CSV folder (which takes place the first time).
- There is a PythonInterface folder which explains how to call PRIMAT from python and how to read its results and modify its options. See readme file in that folder.
- The most crucial functions for performing numerical integration is compiled and stored on disk. Hence after the first run, the loss in compiling this function is avoided by reading the compiled functions which are stored in the lib folder. If this creates a problem, an option allows to turn off this feature
- A $Parthenope option has been added to replace the rates d+d <-> n+He3 and d+d <-> t+p by the ones used in Parthenope3.0 (rates given in Gariazzo et al. 2021, but inferred from Parthenope fortran code)
- Names of folder have changed (MonteCarlo -> Results, Interpolations -> CSV) and the files with nuclear rates are now in the Rates folder.
- Rate for d+a <-> Li6 + g replaced by Trezzi et al. 2017. New reaction file is BBNRatesAC2024.dat.