sdc_impl.hpp

Go to the documentation of this file.

#include "pfasst/controller/sdc.hpp"


namespace pfasst
{
  template<typename time>
  void SDC<time>::run()
  {
    auto sweeper = this->get_level(0);

    for (; this->get_time() < this->get_end_time(); this->advance_time()) {
      bool initial = this->get_step() == 0;
      for (this->set_iteration(0);
           this->get_iteration() < this->get_max_iterations();
           this->advance_iteration()) {
        bool predict = this->get_iteration() == 0;
        if (predict) {
          sweeper->predict(initial);
          sweeper->post_predict();
        } else {
          sweeper->sweep();
          sweeper->post_sweep();
        }
        if (sweeper->converged()) {
          break;
        }
      }
      sweeper->post_step();
      if (this->get_time() + this->get_step_size() < this->get_end_time()) {
        sweeper->advance();
      }
    }
  }
}  // ::pfasst