This version fixes the error where the network stalled before integrating, and the pointer/memory errors. The only remaining known bugs are (1) the partial equilibrium bug, and (2) it is likely that the QSS approximation isn't quite working correctly. I have also added to this version code that reads in a hydro profile file, and outputs that profile at the end to a gnu format file. The hydro profile is now stored in arrays, but we still have to implement the spline interpolator to interpolate the temperature and density from the hydro profile arrays. Also need to implement the spline interpolator for the partition function arrays to complete adding partition function corrections to the calculation.
This version fixes the error where the network stalled before integrating, and the pointer/memory errors. The only remaining known bugs are (1) the partial equilibrium bug, and (2) it is likely that the QSS approximation isn't quite working correctly. I have also added to this version code that reads in a hydro profile file, and outputs that profile at the end to a gnu format file. The hydro profile is now stored in arrays, but we still have to implement the spline interpolator to interpolate the temperature and density from the hydro profile arrays. Also need to implement the spline interpolator for the partition function arrays to complete adding partition function corrections to the calculation.