Method Run

examples/107_wasm/diffusion.cpp:142–167

}

void Solver::Run() {
  auto sem = m.GetSem();
  auto state = g_state;
  auto& s = *state;
  if (sem("init") && s.to_init_field) {
    Init(fcu, m);
  }
  if (sem("flux")) {
    for (auto f : m.Faces()) {
      ff_flux[f] = s.velocity.dot(m.GetSurface(f));
    }
  }
  if (sem()) {
    s.to_init_field = false;
    const Scal h = m.GetCellSize()[0];
    s.dt = std::min(
        0.125 * h * h / s.diffusion, 0.25 * h / s.velocity.abs().max());
    Step(fcu, s.dt, s.diffusion, ff_flux, m);
    Render(*g_canvas, fcu, m);
    m.Comm(&fcu);
    if (m.IsRoot()) {
      s.time += s.dt;
      ++s.step;
    }
  }
}

static void main_loop() {
  auto state = g_state;