Method simulate

src/Solver/Simulator.cpp:79–175

void seissol::Simulator::simulate(seissol::SeisSol& seissolInstance) {
  SCOREP_USER_REGION( "simulate", SCOREP_USER_REGION_TYPE_FUNCTION )

  auto* faultOutputManager = seissolInstance.timeManager().getFaultOutputManager();
  faultOutputManager->writePickpointOutput(0.0, 0.0);

  Stopwatch simulationStopwatch;
  simulationStopwatch.start();

  Stopwatch computeStopwatch;
  Stopwatch ioStopwatch;

  ioStopwatch.start();

  // Set start time (required for checkpointing)
  seissolInstance.timeManager().setInitialTimes(m_currentTime);

  double timeTolerance = seissolInstance.timeManager().getTimeTolerance();

  // Write initial wave field snapshot
  if (m_currentTime == 0.0) {
    Modules::callHook<ModuleHook::SimulationStart>();
  }

  // intialize wave field and checkpoint time
  Modules::setSimulationStartTime(m_currentTime);

  // derive next synchronization time
  double upcomingTime = m_finalTime;
  // NOTE: This will not call the module specific implementation of the synchronization hook
  // since the current time is the simulation start time. We only use this function here to
  // get correct upcoming time. To be on the safe side, we use zero time tolerance.
  upcomingTime = std::min( upcomingTime, Modules::callSyncHook(m_currentTime, 0.0) );

  double lastSplit = 0;

  ioStopwatch.pause();

  Stopwatch::print("Time spent for initial IO:", ioStopwatch.split(), seissol::MPI::mpi.comm());

  while( m_finalTime > m_currentTime + timeTolerance ) {
    if (upcomingTime < m_currentTime + timeTolerance) {
      logError() << "Simulator did not advance in time from" << m_currentTime << "to" << upcomingTime;
    }
    if (m_abort) {
        logInfo(seissol::MPI::mpi.rank()) << "Aborting simulation.";
        break; 
    }

    // update the DOFs
    computeStopwatch.start();
    seissolInstance.timeManager().advanceInTime( upcomingTime );
    computeStopwatch.pause();

    ioStopwatch.start();

    // update current time
    m_currentTime = upcomingTime;