Skip to content

Commit

Permalink
[0D] Fix type-casts and formatting
Browse files Browse the repository at this point in the history
  • Loading branch information
ischoegl committed Jan 15, 2022
1 parent abd4159 commit 0f23536
Showing 1 changed file with 13 additions and 12 deletions.
25 changes: 13 additions & 12 deletions src/zeroD/ReactorNet.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -269,16 +269,17 @@ void ReactorNet::eval(doublereal t, doublereal* y,
"eval(double t, double* ydot)");
m_have_deprecated_eval[n] = true;
} catch (NotImplementedError& err) {
(void)err; // suppress compiler warning
m_reactors[n]->eval(t, m_LHS.data() + m_start[n], m_RHS.data() + m_start[n]);
int yEnd = 0;
if(n == m_reactors.size()-1){
size_t yEnd = 0;
if (n == m_reactors.size() - 1) {
yEnd = m_RHS.size();
} else {
yEnd = m_start[n+1];
yEnd = m_start[n + 1];
}
for (size_t i = m_start[n]; i < yEnd; i++) {
ydot[i] = m_RHS[i] / m_LHS[i];
}
for (int i = m_start[n]; i < yEnd; i++) {
ydot[i] = m_RHS[i]/m_LHS[i];
}
}
m_reactors[n]->resetSensitivity(p);
}
Expand All @@ -290,15 +291,15 @@ void ReactorNet::eval(doublereal t, doublereal* y,
m_reactors[n]->evalEqs(t, y + m_start[n], ydot + m_start[n], p);
} else {
m_reactors[n]->eval(t, m_LHS.data() + m_start[n], m_RHS.data() + m_start[n]);
int yEnd = 0;
if (n == m_reactors.size()-1) {
size_t yEnd = 0;
if (n == m_reactors.size() - 1) {
yEnd = m_RHS.size();
} else {
yEnd = m_start[n+1];
yEnd = m_start[n + 1];
}
for (size_t i = m_start[n]; i < yEnd; i++) {
ydot[i] = m_RHS[i] / m_LHS[i];
}
for (int i = m_start[n]; i < yEnd; i++) {
ydot[i] = m_RHS[i]/m_LHS[i];
}
}
m_reactors[n]->resetSensitivity(p);
}
Expand Down

0 comments on commit 0f23536

Please sign in to comment.