adding more quantities to the eos test, including derivatives

This commit is contained in:
Aaron Dotter
2025-03-06 14:36:07 -05:00
parent 0455ef6b92
commit da7f5e5022
2 changed files with 40 additions and 18 deletions

View File

@@ -13,6 +13,7 @@
#include <iomanip>
#include <sstream>
#include <string>
#include <format>
/**
* @brief 2D array template alias.
@@ -193,18 +194,24 @@ namespace helmholtz
friend std::ostream& operator<<(std::ostream& os, const helmholtz::EOS& eos) {
os << "EOS Data:\n" << std::setw(20) << std::left;
os << " Electron Fraction: " << eos.ye << "\n";
os << " Electron Chemical Potential: " << eos.etaele << "\n";
os << " Electron Number Density: " << eos.xnefer << "\n";
os << " Total Pressure: " << eos.ptot << "\n";
os << " Gas Pressure: " << eos.pgas << "\n";
os << " Radiation Pressure: " << eos.prad << "\n";
os << " Total Energy: " << eos.etot << "\n";
os << " Gas Energy: " << eos.egas << "\n";
os << " Radiation Energy: " << eos.erad << "\n";
os << " Total Entropy: " << eos.stot << "\n";
os << " Gas Entropy: " << eos.sgas << "\n";
os << " Radiation Entropy: " << eos.srad;
os << " Electron Fraction: " << std::format("{0:24.16e}",eos.ye) << "\n";
os << " Electron Chemical Potential: " << std::format("{0:24.16e}",eos.etaele) << "\n";
os << " Electron Number Density: " << std::format("{0:24.16e}",eos.xnefer) << "\n";
os << " Total Pressure: " << std::format("{0:24.16e}",eos.ptot) << "\n";
os << " dPres/dRho: " << std::format("{0:24.16e}",eos.dpresdd) << "\n";
os << " dPres/dT: " << std::format("{0:24.16e}",eos.dpresdt) << "\n";
os << " Gas Pressure: " << std::format("{0:24.16e}",eos.pgas) << "\n";
os << " Radiation Pressure: " << std::format("{0:24.16e}",eos.prad) << "\n";
os << " Total Energy: " << std::format("{0:24.16e}",eos.etot) << "\n";
os << " dEner/dRho: " << std::format("{0:24.16e}",eos.denerdd) << "\n";
os << " dEner/dT: " << std::format("{0:24.16e}",eos.denerdt) << "\n";
os << " Gas Energy: " << std::format("{0:24.16e}",eos.egas) << "\n";
os << " Radiation Energy: " << std::format("{0:24.16e}",eos.erad) << "\n";
os << " Total Entropy: " << std::format("{0:24.16e}",eos.stot) << "\n";
os << " dEntr/dRho: " << std::format("{0:24.16e}",eos.dentrdd) << "\n";
os << " dEntr/dT: " << std::format("{0:24.16e}",eos.dentrdt) << "\n";
os << " Gas Entropy: " << std::format("{0:24.16e}",eos.sgas) << "\n";
os << " Radiation Entropy: " << std::format("{0:24.16e}",eos.srad);
return os;
}
};
@@ -356,4 +363,4 @@ namespace helmholtz
}
#endif // HELM_H
#endif // HELM_H