BOOLEAN use_subcycling_wip "Use subcycling in time"
{
} no  

#include <cassert>

template <typename T> constexpr T ipow(T x, int n) {
  assert(n >= 0);
  T r = 1;
  while (n > 0) {
    if (n % 2)
      r *= x;
    n /= 2;
    x *= x;
  }
  return r;
}

  const int timereffact = use_subcycling_wip ? 2 : 1;
  if (level == 0) {
    // We are creating the coarsest level
    leveldata.iteration = 0;
    leveldata.delta_iteration = leveldata.coarse_delta_iteration;
    leveldata.time = 0.0;
    leveldata.delta_time = 1.0;
  } else {
    // We are creating a new refined level
    auto &coarseleveldata = ghext->leveldata.at(level - 1);
    leveldata.iteration = coarseleveldata.iteration;
    assert(coarseleveldata.delta_iteration % timereffact == 0);
    leveldata.delta_iteration = coarseleveldata.delta_iteration / timereffact;
    leveldata.time = coarseleveldata.time;
    leveldata.delta_time = coarseleveldata.delta_time / timereffact;
  }

  // We assume that this level is at the same time as the next coarser level
  const int timereffact = use_subcycling_wip ? 2 : 1;
  // leveldata.iteration = coarseleveldata.iteration;
  // assert(coarseleveldata.delta_iteration % timereffact == 0);
  // leveldata.delta_iteration = coarseleveldata.delta_iteration / timereffact;
  // leveldata.time = coarseleveldata.time;
  // leveldata.delta_time = coarseleveldata.delta_time / timereffact;
  assert(leveldata.iteration == coarseleveldata.iteration);
  assert(leveldata.delta_iteration ==
         coarseleveldata.delta_iteration / timereffact);
  assert(leveldata.time == coarseleveldata.time);
  assert(leveldata.delta_time == coarseleveldata.delta_time / timereffact);

#include <cstdint>

    // For subcycling in time, there really should be one copy of each
    // integrated grid scalar per level. We don't do that yet; insted,
    // we assume that grid scalars only hold "analysis" data.

    // Iteration and time at which this level is valid
    constexpr static int64_t coarse_delta_iteration = 1LL << 32;
    int64_t iteration, delta_iteration;
    CCTK_REAL time, delta_time;