| 467 | } |
| 468 | |
| 469 | void ReactorNet::evalDae(double t, double* y, double* ydot, double* p, double* residual) |
| 470 | { |
| 471 | m_time = t; |
| 472 | updateState(y); |
| 473 | for (size_t n = 0; n < m_reactors.size(); n++) { |
| 474 | m_reactors[n]->applySensitivity(p); |
| 475 | m_reactors[n]->evalDae(t, y, ydot, residual); |
| 476 | m_reactors[n]->resetSensitivity(p); |
| 477 | } |
| 478 | checkFinite("ydot", ydot, m_nv); |
| 479 | } |
| 480 | |
| 481 | void ReactorNet::getConstraints(double* constraints) |
| 482 | { |
nothing calls this directly
no test coverage detected