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()); } } }