tboudreaux/GridFire
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

fix(atomicSpecies.h-->-species.h): added species.h includes for spesific species where needed

Browse Source
This commit is contained in:
Emily Boudreaux
2025-06-19 14:51:42 -04:00
parent 856ab51b4c
commit c3ff2420e9
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/network/private/netgraph.cpp
View File

@@ -296,7 +296,7 @@ namespace serif::network {
const double T9, const double rho) const { const double T9, const double rho) const {
const auto &constants = serif::constant::Constants::getInstance(); const auto &constants = serif::constant::Constants::getInstance();
constexpr auto u = constants.get("u"); // Atomic mass unit in g/mol const auto u = constants.get("u"); // Atomic mass unit in g/mol
const double k_reaction = reaction.calculate_rate(T9); // PERF: Consider precomputing all of these and putting them into an O(1) lookup table. const double k_reaction = reaction.calculate_rate(T9); // PERF: Consider precomputing all of these and putting them into an O(1) lookup table.
double reactant_product = 1.0; double reactant_product = 1.0;
@@ -327,7 +327,7 @@ namespace serif::network {
// Convert to number density // Convert to number density
double ni; double ni;
const double denominator = atomicMassAMU * u.value(); const double denominator = atomicMassAMU * u.value;
if (denominator > minThreshold) { if (denominator > minThreshold) {
ni = (Yi * rho) / (denominator); ni = (Yi * rho) / (denominator);
} else { } else {
@@ -341,7 +341,7 @@ namespace serif::network {
reactant_product /= static_cast<double>(std::tgamma(count + 1)); // Gamma function for factorial reactant_product /= static_cast<double>(std::tgamma(count + 1)); // Gamma function for factorial
} }
} }
constexpr double Na = constants.get("N_a").value(); // Avogadro's number in mol^-1 const double Na = constants.get("N_a").value; // Avogadro's number in mol^-1
const double molarCorrectionFactor = std::pow(Na, reaction.reactants().size() - 1); const double molarCorrectionFactor = std::pow(Na, reaction.reactants().size() - 1);
return (reactant_product * k_reaction) / molarCorrectionFactor; // reaction rate in per volume per time (particles/cm^3/s) return (reactant_product * k_reaction) / molarCorrectionFactor; // reaction rate in per volume per time (particles/cm^3/s)
} }