| 257 | } |
| 258 | |
| 259 | double ReactorNet::step() |
| 260 | { |
| 261 | if (!m_init) { |
| 262 | initialize(); |
| 263 | } else if (!m_integrator_init) { |
| 264 | reinitialize(); |
| 265 | } |
| 266 | m_time = m_integ->step(m_time + 1.0); |
| 267 | updateState(m_integ->solution()); |
| 268 | return m_time; |
| 269 | } |
| 270 | |
| 271 | void ReactorNet::solveSteady(int loglevel) |
| 272 | { |