BUG: reactionSensitivity - corrected integral calculations

This commit is contained in:
Andrew Heather
2017-06-28 14:07:57 +01:00
parent 367095f666
commit bb2041ef69

View File

@ -110,12 +110,12 @@ calculateSpeciesRR
if (RR[0] > 0.0) if (RR[0] > 0.0)
{ {
production_[speciei][reactioni] = RR[0]; production_[speciei][reactioni] = RR[0];
productionInt_[speciei][reactioni] =+ dt*RR[0]; productionInt_[speciei][reactioni] += dt*RR[0];
} }
else if (RR[0] < 0.0) else if (RR[0] < 0.0)
{ {
consumption_[speciei][reactioni] = RR[0]; consumption_[speciei][reactioni] = RR[0];
consumptionInt_[speciei][reactioni] =+ dt*RR[0]; consumptionInt_[speciei][reactioni] += dt*RR[0];
} }
else else
{ {