| 285 | } |
| 286 | |
| 287 | Eigen::SparseMatrix<double> ReactorNet::steadyJacobian(double rdt) |
| 288 | { |
| 289 | if (!m_init) { |
| 290 | initialize(); |
| 291 | } else if (!m_integrator_init) { |
| 292 | reinitialize(); |
| 293 | } |
| 294 | vector<double> y0(neq()); |
| 295 | vector<double> y1(neq()); |
| 296 | getState(y0.data()); |
| 297 | SteadyReactorSolver solver(this, y0.data()); |
| 298 | solver.evalJacobian(y0.data()); |
| 299 | if (rdt) { |
| 300 | solver.linearSolver()->updateTransient(rdt, solver.transientMask().data()); |
| 301 | } |
| 302 | return std::dynamic_pointer_cast<EigenSparseJacobian>(solver.linearSolver())->jacobian(); |
| 303 | } |
| 304 | |
| 305 | void ReactorNet::getEstimate(double time, int k, double* yest) |
| 306 | { |
nothing calls this directly
no test coverage detected