inherits: AMREX# Grid # grid contains scalars

inherits: AMREX # Grid # grid contains scalars

SCHEDULE TestProlongate_Test IN CCTK_INITIAL

SCHEDULE TestProlongate_Regrid AT CCTK_POSTINITIAL
{
  LANG: C
  WRITES: AMReX::regrid_error(everywhere)
} "Set up test grid"
SCHEDULE TestProlongate_Regrid AT CCTK_POSTSTEP

  WRITES: data

  WRITES: AMReX::regrid_error(everywhere)
} "Set up test grid"
SCHEDULE TestProlongate_Test AT CCTK_INITIAL
{
  LANG: C
  WRITES: data(interior)

SCHEDULE TestProlongate_Sync IN CCTK_INITIAL AFTER TestProlongate_Test

SCHEDULE TestProlongate_Sync AT CCTK_INITIAL AFTER TestProlongate_Test

#include <stddef.h>

#include <cstddef>
#include <cstdio>

}
void TestProlongate_Regrid(CCTK_ARGUMENTS)
{
  DECLARE_CCTK_ARGUMENTS;
  DECLARE_CCTK_PARAMETERS;
  if(cctk_iteration < 1)
    return;
  printf("Setting grid\n");
  Loop::loop_int<1, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    if(p.x > -0.5 && p.x < 0.5 &&
       p.y > -0.5 && p.y < 0.5 &&
       p.z > -0.5 && p.z < 0.5) {
      fprintf(stderr,"Setting error %g %g %g\n", p.x, p.y, p.z);
      regrid_error[p.idx] = 1e3;
    } else {
      regrid_error[p.idx] = 0.;
    }
  });
  fprintf(stderr, "done\n");

$ncells = 16

$ncells = 8
$blocking_factor = 2

Cactus::cctk_itlast = 0

Cactus::cctk_itlast = 4

#AMReX::blocking_factor_y = 1
#AMReX::blocking_factor_z = 1

AMReX::ghost_size = 3
AMReX::prolongation_order = 4
AMReX::verbose = yes

AMReX::blocking_factor_x = $blocking_factor
AMReX::blocking_factor_y = 1 # $blocking_factor
AMReX::blocking_factor_z = 1 # $blocking_factor

AMReX::regrid_every = 16

AMReX::regrid_every = 1

IO::out_every = 16   #TODO $ncells * 2 ** ($nlevels - 1) / 2

IO::out_every = 1  #TODO $ncells * 2 ** ($nlevels - 1) / 2