eschnett/cactusamrex - Change PGEO55PNCMJMUNKGEFPXH4RPRM7VJ2XR6IOO3ZKLK4KDXB77XE7QC

CarpetX: Allow different numbers of ghost zones in different directions

Created by Erik Schnetter on August 5, 2020

PGEO55PNCMJMUNKGEFPXH4RPRM7VJ2XR6IOO3ZKLK4KDXB77XE7QC

Dependencies

In channels

main

Change contents

Insertion in CarpetX/param.ccl at line 183 [4.17918]

[3.530]

{
  -1 :: "use ghost_size_[xyz]"
  0:* :: ""
} -1
CCTK_INT ghost_size_x "Number of ghost zones"
{
  0:* :: ""
} 1
CCTK_INT ghost_size_y "Number of ghost zones"
{
  0:* :: ""
} 1
CCTK_INT ghost_size_z "Number of ghost zones"

Replacement in CarpetX/src/driver.cxx at line 202 [4.15972]

B:BD[5.3382] → [5.3382:3438]

    nghostzones = {ghost_size, ghost_size, ghost_size};

[5.3382]

[5.3438]

    nghostzones = {ghost_size >= 0 ? ghost_size : ghost_size_x,
                   ghost_size >= 0 ? ghost_size : ghost_size_y,
                   ghost_size >= 0 ? ghost_size : ghost_size_z};

Replacement in CarpetX/src/driver.cxx at line 706 [4.15972]

B:BD[2.6550] → [2.6550:6625]

  leveldata.fab = make_unique<amrex::FabArrayBase>(ba, dm, 1, ghost_size);

[2.6550]

[2.6625]

  const amrex::IntVect nghostzones = {
      ghost_size >= 0 ? ghost_size : ghost_size_x,
      ghost_size >= 0 ? ghost_size : ghost_size_y,
      ghost_size >= 0 ? ghost_size : ghost_size_z};
  leveldata.fab = make_unique<amrex::FabArrayBase>(ba, dm, 1, nghostzones);

Replacement in CarpetX/src/driver.cxx at line 946 [4.15972]

B:BD[5.4900] → [6.4626:4701]

  leveldata.fab = make_unique<amrex::FabArrayBase>(ba, dm, 1, ghost_size);

[5.4900]

[6.4701]

  const amrex::IntVect nghostzones = {
      ghost_size >= 0 ? ghost_size : ghost_size_x,
      ghost_size >= 0 ? ghost_size : ghost_size_y,
      ghost_size >= 0 ? ghost_size : ghost_size_z};
  leveldata.fab = make_unique<amrex::FabArrayBase>(ba, dm, 1, nghostzones);

Replacement in CarpetX/src/schedule.cxx at line 872 [7.531]

∅:D[3.1907] → [7.2155:2187]

∅:D[8.9502] → [7.2155:2187]

B:BD[7.2155] → [7.2155:2187]

B:BD[7.2187] → [3.1908:1954]

  for (int d = 0; d < dim; ++d)
    cctkGH->cctk_nghostzones[d] = ghost_size;

[8.9502]

[9.7220]

  cctkGH->cctk_nghostzones[0] = ghost_size >= 0 ? ghost_size : ghost_size_x;
  cctkGH->cctk_nghostzones[1] = ghost_size >= 0 ? ghost_size : ghost_size_y;
  cctkGH->cctk_nghostzones[2] = ghost_size >= 0 ? ghost_size : ghost_size_z;