void EventDriven::initOSIRhs()
{
// === initialization for OneStepIntegrators ===
OSI::TYPES osiType = (*_allOSI->begin())->getType();
for (OSIIterator itosi = _allOSI->begin(); itosi != _allOSI->end(); ++itosi)
{
//Check whether OSIs used are of the same type
if ((*itosi)->getType() != osiType)
RuntimeException::selfThrow("OSIs used must be of the same type");
// perform the initialization
DynamicalSystemsGraph::VIterator dsi, dsend;
SP::DynamicalSystemsGraph osiDSGraph = (*itosi)->dynamicalSystemsGraph();
for (std11::tie(dsi, dsend) = osiDSGraph->vertices(); dsi != dsend; ++dsi)
{
if (!(*itosi)->checkOSI(dsi)) continue;
SP::DynamicalSystem ds = osiDSGraph->bundle(*dsi);
// Initialize right-hand side
ds->initRhs(startingTime());
}
}
}
