void FUNCTION GetLoopBoxAll(
  CCTK_POINTER_TO_CONST IN cctkGH,
  CCTK_INT IN size,
  CCTK_INT ARRAY OUT loop_min,
  CCTK_INT ARRAY OUT loop_max)
PROVIDES FUNCTION GetLoopBoxAll WITH CarpetX_GetLoopBoxAll LANGUAGE C
void FUNCTION GetLoopBoxInt(
  CCTK_POINTER_TO_CONST IN cctkGH,
  CCTK_INT IN size,
  CCTK_INT ARRAY OUT loop_min,
  CCTK_INT ARRAY OUT loop_max)
PROVIDES FUNCTION GetLoopBoxInt WITH CarpetX_GetLoopBoxInt LANGUAGE C

#include <cctk_Loop.h>

#include "loop.hxx"

extern "C" void CarpetX_GetLoopBoxAll(const void *restrict const cctkGH_,
                                      const CCTK_INT size,
                                      CCTK_INT *restrict const loop_min,
                                      CCTK_INT *restrict const loop_max) {
  const cGH *restrict const cctkGH = static_cast<const cGH *>(cctkGH_);
  assert(size == dim);
  const GridDescBase grid(cctkGH);
  vect<int, dim> imin, imax;
  // Assumd cell centred box
  grid.box_all<1, 1, 1>(grid.nghostzones, imin, imax);
  for (int d = 0; d < dim; ++d)
    loop_min[d] = imin[d];
  for (int d = 0; d < dim; ++d)
    loop_max[d] = imax[d];
}
extern "C" void CarpetX_GetLoopBoxInt(const void *restrict const cctkGH_,
                                      const CCTK_INT size,
                                      CCTK_INT *restrict const loop_min,
                                      CCTK_INT *restrict const loop_max) {
  const cGH *restrict const cctkGH = static_cast<const cGH *>(cctkGH_);
  assert(size == dim);
  const GridDescBase grid(cctkGH);
  vect<int, dim> imin, imax;
  // Assumd cell centred box
  grid.box_int<1, 1, 1>(grid.nghostzones, imin, imax);
  for (int d = 0; d < dim; ++d)
    loop_min[d] = imin[d];
  for (int d = 0; d < dim; ++d)
    loop_max[d] = imax[d];
}