fix(network): added missing negative sign in energy calculation

This commit is contained in:
2025-03-21 14:44:35 -04:00
parent 3c657ca85a
commit 328b2cfe8d
2 changed files with 2 additions and 2 deletions

View File

@@ -421,7 +421,7 @@ namespace nnApprox8{
for (int i=0; i<Net::niso; i++) {
enuc += Net::mion[i]*dydt[i];
}
dydt[Net::iener] = enuc*avo*clight*clight;
dydt[Net::iener] = -enuc*avo*clight*clight;
}
nuclearNetwork::NetOut Approx8Network::evaluate(const nuclearNetwork::NetIn &netIn) {

View File

@@ -46,5 +46,5 @@ TEST_F(approx8Test, evaluate) {
EXPECT_DOUBLE_EQ(netOut.composition[nnApprox8::Net::ih1], 0.50166260916650918);
EXPECT_DOUBLE_EQ(netOut.composition[nnApprox8::Net::ihe4],0.48172270591286032);
EXPECT_DOUBLE_EQ(netOut.energy, -1.6433049870528356e+18);
EXPECT_DOUBLE_EQ(netOut.energy, 1.6433049870528356e+18);
}