Numerical Enhancements for the Chiral Mean Field Model

A very successful model for the equation of state at large baryon densities that has been used to describe neutron stars, neutron star mergers, and heavy-ion collisions is the Chiral Mean Field model \cite{Dexheimer:2009hi}. Currently, the Chiral Mean Field model has a very slow execution time that lasts on the order of one month to construct 3D tables. Because heavy-ion collisions require 4D tables (baryon, strangeness, and electric charge are conserved) reduced running time is a high priority. In this work, we significantly reduce this time by employing program optimization techniques at the source, build, and compilation levels. Dynamic analysis was done by time profiling with Gprof and analyzing memory issues with Valgrind. For numerical integration, the QUADPACK library was linked. Finally, the execution was parallelized in constant temperature slices. Either for stellar mater or isospin symmetric matter, the execution time is reduced from months to hours.