eschnett/cactusamrex - Change 2DKSL6DKZAIYQUJGDULORCKU5K4Z5Z3W4RIKQYDSLKMCNQNDZFBAC

AMR evolution should work now

Created by Erik Schnetter on July 24, 2019

2DKSL6DKZAIYQUJGDULORCKU5K4Z5Z3W4RIKQYDSLKMCNQNDZFBAC

Dependencies

In channels

main

master

Change contents

Insertion in AMReX/interface.ccl at line 4 [4.17976]
[4.18035]
[5.0]
```
INHERITS: IO
```
Insertion in AMReX/param.ccl at line 4 [4.17918]
[6.10]
[7.0]
```
BOOLEAN verbose "Verbose output"
{
} "no"
```
Replacement in AMReX/param.ccl at line 45 [4.17918]
B:BD[6.359] → [8.79:91]
```
  1:* :: ""
```
[6.359]
[9.0]
```
  1:* :: "must be a multiple of the blocking factor"
```
Replacement in AMReX/param.ccl at line 50 [4.17918]
B:BD[6.415] → [8.92:104]
```
  1:* :: ""
```
[6.415]
[9.7]
```
  1:* :: "must be a multiple of the blocking factor"
```

Insertion in AMReX/param.ccl at line 54 [4.17918]

[6.469]

{
  1:* :: "must be a multiple of the blocking factor"
} 128
CCTK_INT blocking_factor_x "Blocking factor"
{
  1:* :: ""
} 8
CCTK_INT blocking_factor_y "Blocking factor"

Replacement in AMReX/param.ccl at line 66 [4.17918]

∅:D[8.117] → [9.14:20]

B:BD[6.477] → [9.14:20]

} 128

[8.117]

[8.118]

} 8
CCTK_INT blocking_factor_z "Blocking factor"
{
  1:* :: "must be a power of 2"
} 8
CCTK_INT max_grid_size_x "Maximum grid size"
{
  1:* :: "must be a multiple of the blocking factor"
} 32
CCTK_INT max_grid_size_y "Maximum grid size"
{
  1:* :: "must be a multiple of the blocking factor"
} 32
CCTK_INT max_grid_size_z "Maximum grid size"
{
  1:* :: "must be a multiple of the blocking factor"
} 32
CCTK_INT max_tile_size_x "Maximum tile size"
{
  1:* :: ""
} 1024000
CCTK_INT max_tile_size_y "Maximum tile size"
{
  1:* :: ""
} 16
CCTK_INT max_tile_size_z "Maximum tile size"
{
  1:* :: ""
} 32

Insertion in AMReX/param.ccl at line 111 [4.17918]
[6.486]
[6.486]

Replacement in AMReX/param.ccl at line 119 [4.17918]

B:BD[2.7] → [2.7:65]

CCTK_INT out_every "How often to output" STEERABLE=always

[2.7]

[2.65]

CCTK_REAL dtfac "The standard timestep condition dt = dtfac*min(delta_space)"

Replacement in AMReX/param.ccl at line 121 [4.17918]
B:BD[2.67] → [2.67:124]
∅:D[2.124] → [6.544:548]
B:BD[6.544] → [6.544:548]
```
  0 :: "no output"
  1:* :: "every that many iterations"
} 1
```
[2.67]
[6.548]
```
  *:* :: ""
} 0.5 
```

Replacement in AMReX/param.ccl at line 126 [4.17918]

B:BD[6.551] → [6.551:629]

CCTK_REAL dtfac "The standard timestep condition dt = dtfac*min(delta_space)"

[6.551]

[6.629]

BOOLEAN out_tsv "Output in TSV (tab-separated values) format?" STEERABLE=always

Replacement in AMReX/param.ccl at line 128 [4.17918]
B:BD[6.631] → [6.631:650]
```
  *:* :: ""
} 0.5 
```
[6.631]
[6.650]
```
} no
```
Insertion in AMReX/param.ccl at line 140 [4.17918]
[7.186]
```
SHARES: IO
USES STRING out_dir
USES CCTK_INT out_every
```
Insertion in AMReX/src/driver.cxx at line 11 [4.15972]
[9.46]
[4.16090]
```
#include <AMReX_ParmParse.H>
```
Insertion in AMReX/src/driver.cxx at line 70 [4.15972]
[5.361]
[5.361]
```
  DECLARE_CCTK_PARAMETERS;
```

Replacement in AMReX/src/driver.cxx at line 76 [4.15972]

B:BD[5.463] → [5.463:505]

  CCTK_VINFO("ErrorEst level %d", level);

[5.463]

[5.505]

  if (verbose)
    CCTK_VINFO("ErrorEst level %d", level);

Replacement in AMReX/src/driver.cxx at line 104 [4.15972]

B:BD[5.1031] → [5.1031:1110]

  auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling({1024000, 16, 32});

[5.1031]

[5.1110]

  auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
      {max_tile_size_x, max_tile_size_y, max_tile_size_z});

Insertion in AMReX/src/driver.cxx at line 115 [4.15972]

[5.1411]

    // Note: This excludes ghosts, it's not the proper Cactus lsh

Replacement in AMReX/src/driver.cxx at line 119 [4.15972]

B:BD[5.1502] → [5.1502:1774]

    const Array4<CCTK_REAL> &ctagarr = groupdata.mfab.at(tl)->array(mfi);
    const CCTK_REAL *restrict ctags = ctagarr.ptr(imin.x, imin.y, imin.z, vi);
    const Array4<char> &atagarr = tags.array(mfi);
    char *restrict atags = atagarr.ptr(imin.x, imin.y, imin.z, vi);

[5.1502]

[5.1774]

    const Array4<CCTK_REAL> &vars = groupdata.mfab.at(tl)->array(mfi);
    const CCTK_REAL *restrict ptr = vars.ptr(imin.x, imin.y, imin.z, vi);
    const Array4<char> &tagarr = tags.array(mfi);

Insertion in AMReX/src/driver.cxx at line 126 [4.15972]
[5.1865]
[5.1865]

Replacement in AMReX/src/driver.cxx at line 132 [4.15972]

B:BD[5.2051] → [5.2051:2124]


          atags[idx] = ctags[idx] == 0.0 ? TagBox::CLEAR : TagBox::SET;

[5.2051]

[5.2124]

          tagarr(imin.x + i, imin.y + j, imin.z + k) =
              ptr[idx] == 0.0 ? TagBox::CLEAR : TagBox::SET;

Replacement in AMReX/src/driver.cxx at line 143 [4.15972]

B:BD[8.1479] → [8.1479:1534]

  CCTK_VINFO("MakeNewLevelFromScratch level %d", lev);

[8.1479]

[8.1534]

  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("MakeNewLevelFromScratch level %d", lev);

Replacement in AMReX/src/driver.cxx at line 151 [4.15972]

B:BD[8.1739] → [8.1739:1793]

  CCTK_VINFO("MakeNewLevelFromCoarse level %d", lev);

[8.1739]

[10.1244]

  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("MakeNewLevelFromCoarse level %d", lev);

Replacement in AMReX/src/driver.cxx at line 158 [4.15972]

B:BD[8.1931] → [8.1931:1974]

  CCTK_VINFO("RemakeLevel level %d", lev);

[8.1931]

[8.1974]

  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("RemakeLevel level %d", lev);

Replacement in AMReX/src/driver.cxx at line 164 [4.15972]

B:BD[8.2019] → [8.2019:2061]

  CCTK_VINFO("ClearLevel level %d", lev);

[8.2019]

[8.2061]

  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("ClearLevel level %d", lev);

Replacement in AMReX/src/driver.cxx at line 175 [4.15972]
B:BD[4.16301] → [11.333:358]
```
  CCTK_VINFO("Startup");
```
[4.16301]
[11.358]
```
  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("Startup");
```
Replacement in AMReX/src/driver.cxx at line 213 [4.15972]
∅:D[12.297] → [11.988:1013]
B:BD[11.988] → [11.988:1013]
```
  CCTK_VINFO("SetupGH");
```
[12.297]
[11.1013]
```
  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("SetupGH");
```
Replacement in AMReX/src/driver.cxx at line 233 [4.15972]
∅:D[6.681] → [11.1152:1176]
∅:D[12.333] → [11.1152:1176]
B:BD[11.1152] → [11.1152:1176]
```
  CCTK_VINFO("InitGH");
```
[6.681]
[11.1176]
```
  if (verbose)
    CCTK_VINFO("InitGH");
```

Insertion in AMReX/src/driver.cxx at line 252 [4.15972]

[4.16918]

[8.2065]

  // Set blocking factors via parameter table since AmrMesh needs to
  // know them when its constructor is running, but there are no
  // constructor arguments for them
  ParmParse pp;
  pp.add("amr.blocking_factor_x", blocking_factor_x);
  pp.add("amr.blocking_factor_y", blocking_factor_y);
  pp.add("amr.blocking_factor_z", blocking_factor_z);
  pp.add("amr.max_grid_size_x", max_grid_size_x);
  pp.add("amr.max_grid_size_y", max_grid_size_y);
  pp.add("amr.max_grid_size_z", max_grid_size_z);

Deletion in AMReX/src/driver.cxx at line 265 [4.15972]
∅:D[8.2182] → [10.1580:1622]
B:BD[10.1580] → [10.1580:1622]
B:BD[10.1622] → [8.2183:2306]
∅:D[10.1745] → [9.680:712]
∅:D[8.2306] → [9.680:712]
B:BD[9.680] → [9.680:712]
B:BD[9.712] → [8.2307:2356]
```
#warning "TODO: increase blocking factor"
  const int blocking_factor = 8;
  // const int blocking_factor = 1;
  ghext->amrcore->SetBlockingFactor(blocking_factor);
  const int max_grid_size = 32;
  ghext->amrcore->SetMaxGridSize(max_grid_size);
```

Replacement in AMReX/src/driver.cxx at line 266 [4.15972]

B:BD[4.17048] → [8.2357:2410]

∅:D[8.2410] → [10.1746:1846]

B:BD[4.17048] → [10.1746:1846]

B:BD[10.1846] → [8.2411:2460]

  int maxnumlevels = ghext->amrcore->maxLevel() + 1;
  for (int level = 0; level < maxnumlevels; ++level) {
    CCTK_VINFO("Geometry level %d:", level);
    cout << ghext->amrcore->Geom(level) << "\n";

[4.17048]

[10.1895]

  if (verbose) {
    int maxnumlevels = ghext->amrcore->maxLevel() + 1;
    for (int level = 0; level < maxnumlevels; ++level) {
      CCTK_VINFO("Geometry level %d:", level);
      cout << ghext->amrcore->Geom(level) << "\n";
    }

Replacement in AMReX/src/driver.cxx at line 289 [4.15972]

B:BD[5.2623] → [5.2623:2674]

  CCTK_VINFO("CreateRefinedGrid level %d", level);

[5.2623]

[4.17361]

  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("CreateRefinedGrid level %d", level);

Deletion in AMReX/src/driver.cxx at line 301 [4.15972]

B:BD[8.2814] → [8.2814:2918]

  cout << "CRG numlevels=" << numlevels << "\n";
  cout << "CRG maxnumlevels=" << maxnumlevels << "\n";

Replacement in AMReX/src/driver.cxx at line 310 [4.15972]

B:BD[5.3257] → [5.3257:3301]

  CCTK_VINFO("SetupLevel level %d", level);

[5.3257]

[7.827]

  if (verbose)
    CCTK_VINFO("SetupLevel level %d", level);

Replacement in AMReX/src/driver.cxx at line 346 [4.15972]

∅:D[12.483] → [13.83:158]

B:BD[11.1300] → [13.83:158]

  CCTK_VINFO("ScheduleTraverseGH [%d] %s", cctkGH->cctk_iteration, where);

[12.483]

[11.1336]

  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("ScheduleTraverseGH [%d] %s", cctkGH->cctk_iteration, where);

Replacement in AMReX/src/driver.cxx at line 358 [4.15972]
B:BD[4.17509] → [11.1680:1706]
```
  CCTK_VINFO("Shutdown");
```
[4.17509]
[11.1706]
```
  DECLARE_CCTK_PARAMETERS;
  if (verbose)
    CCTK_VINFO("Shutdown");
```
Insertion in AMReX/src/io.cxx at line 10 [14.107]
[14.243]
[14.243]
```
#include <AMReX_Orientation.H>
```
Insertion in AMReX/src/io.cxx at line 14 [14.107]
[14.294]
[14.294]
```
#include <fstream>
```

Insertion in AMReX/src/io.cxx at line 26 [14.107]

[14.489]


namespace {
// Convert a (direction, face) pair to an AMReX Orientation
Orientation orient(int d, int f) {
  return Orientation(d, Orientation::Side(f));
}
} // namespace
template <typename T>
void write_arrays(ostream &os, const cGH *restrict cctkGH, int level,
                  int component, const vector<const T *> ptrs,
                  const int *restrict ash, const int *restrict lbnd,
                  const int *restrict lsh, const int *restrict bbox,
                  const int *restrict nghostzones) {
  DECLARE_CCTK_ARGUMENTS;
  const int levfac = 1 << level;
  CCTK_REAL x0[dim], dx[dim];
  for (int d = 0; d < dim; ++d) {
    dx[d] = cctk_delta_space[d] / levfac;
    x0[d] = cctk_origin_space[d] + 0.5 * dx[d];
  }
  int imin[dim], imax[dim];
  for (int d = 0; d < dim; ++d) {
    imin[d] = cctk_bbox[2 * d] ? 0 : cctk_nghostzones[d];
    imax[d] = cctk_lsh[d] - (cctk_bbox[2 * d + 1] ? 0 : cctk_nghostzones[d]);
  }
  constexpr int di = 1;
  const int dj = di * ash[0];
  const int dk = dj * ash[1];
  for (int k = 0; k < lsh[2]; ++k) {
    for (int j = 0; j < lsh[1]; ++j) {
      for (int i = 0; i < lsh[0]; ++i) {
        const int idx = di * i + dj * j + dk * k;
        const int isghost = !(i >= imin[0] && i < imax[0] && j >= imin[1] &&
                              j < imax[1] && k >= imin[2] && k < imax[2]);
        T x = x0[0] + (lbnd[0] + i) * dx[0];
        T y = x0[1] + (lbnd[1] + j) * dx[1];
        T z = x0[2] + (lbnd[2] + k) * dx[2];
        os << cctk_iteration << "\t" << cctk_time << "\t" << level << "\t"
           << component << "\t" << isghost << "\t" << (lbnd[0] + i) << "\t"
           << (lbnd[1] + j) << "\t" << (lbnd[2] + k) << "\t" << x << "\t" << y
           << "\t" << z;
        for (const auto &ptr : ptrs)
          os << "\t" << ptr[idx];
        os << "\n";
      }
    }
  }
}
void WriteASCII(const cGH *restrict cctkGH, const string &filename, int gi,
                const vector<string> &varnames) {
  ostringstream buf;
  buf << filename << ".tsv";
  ofstream file(buf.str());

Insertion in AMReX/src/io.cxx at line 86 [14.107]

[14.490]

  // Output header
  file << "#"
       << "\t"
       << "time"
       << "\t"
       << "level"
       << "\t"
       << "component"
       << "\t"
       << "isghost"
       << "\t"
       << "i"
       << "\t"
       << "j"
       << "\t"
       << "k"
       << "\t"
       << "x"
       << "\t"
       << "y"
       << "\t"
       << "z";
  for (const auto &varname : varnames)
    file << "\t" << varname;
  file << "\n";
  for (const auto &leveldata : ghext->leveldata) {
    const auto &groupdata = leveldata.groupdata.at(gi);
    const int tl = 0;
    const MultiFab &mfab = *groupdata.mfab.at(tl);
    for (MFIter mfi(mfab); mfi.isValid(); ++mfi) {
      const Box &fbx = mfi.fabbox();       // allocated array
      const Box &bx = mfi.tilebox();       // current region (without ghosts)
      const Box &gbx = mfi.growntilebox(); // current region (with ghosts)
      const Dim3 imin = lbound(gbx);
      // Local shape
      int lsh[dim];
      for (int d = 0; d < dim; ++d)
        lsh[d] = gbx[orient(d, 1)] - gbx[orient(d, 0)] + 1;
      // Allocated shape
      int ash[dim];
      for (int d = 0; d < dim; ++d)
        ash[d] = fbx[orient(d, 1)] - fbx[orient(d, 0)] + 1;
      // Local extent
      int lbnd[dim];
      for (int d = 0; d < dim; ++d)
        lbnd[d] = gbx[orient(d, 0)];
      // Boundaries
      int bbox[2 * dim];
      for (int d = 0; d < dim; ++d)
        for (int f = 0; f < 2; ++f)
          bbox[2 * d + f] = bx[orient(d, f)] == gbx[orient(d, f)];
      // Number of ghost zones
      int nghostzones[dim];
      for (int d = 0; d < dim; ++d)
        nghostzones[d] = mfab.fb_nghost[d];
      const Array4<CCTK_REAL> &vars = groupdata.mfab.at(tl)->array(mfi);
      vector<const CCTK_REAL *> ptrs(groupdata.numvars);
      for (int vi = 0; vi < groupdata.numvars; ++vi)
        ptrs.at(vi) = vars.ptr(imin.x, imin.y, imin.z, vi);
      write_arrays(file, cctkGH, leveldata.level, mfi.index(), ptrs, ash, lbnd,
                   lsh, bbox, nghostzones);
    }
  }
  file.close();
}

Insertion in AMReX/src/io.cxx at line 164 [14.107]

[14.587]

[2.125]

  const bool is_root = CCTK_MyProc(nullptr) == 0;
  if (is_root)
    cout << "OutputGH: iteration " << cctk_iteration << ", time " << cctk_time
         << ", run time " << CCTK_RunTime() << " s\n";

Replacement in AMReX/src/io.cxx at line 184 [14.107]
B:BD[15.510] → [15.510:548]
```
      buf << "wavetoy/" << groupname;
```
[15.510]
[15.548]
```
      buf << out_dir << "/" << groupname;
```

Insertion in AMReX/src/io.cxx at line 211 [14.107]

[15.1672]

[16.1090]


      if (out_tsv) {
        ostringstream procbuf;
        procbuf << out_dir << "/" << groupname;
        procbuf << ".it" << setw(6) << setfill('0') << cctk_iteration;
        procbuf << ".p" << setw(4) << setfill('0') << CCTK_MyProc(nullptr);
        string procfilename = procbuf.str();

Insertion in AMReX/src/io.cxx at line 219 [14.107]

[16.1091]

[15.1673]

        WriteASCII(cctkGH, procfilename, gi, varnames);
      }

Replacement in AMReX/src/io.cxx at line 235 [14.107]

B:BD[17.692] → [17.692:800]

      CCTK_VINFO("maxabs(%s)=%g", CCTK_VarName(groupdata.firstvarindex + vi),
                 red.maxabs);

[17.692]

[18.0]

      if (is_root)
        cout << "  "
             << unique_ptr<char>(CCTK_FullName(groupdata.firstvarindex + vi))
                    .get()
             << ": maxabs=" << red.maxabs << "\n";

Insertion in AMReX/src/reduction.cxx at line 3 [17.870]
[17.918]
[17.918]
```
#include <cctk_Parameters.h>
```

Replacement in AMReX/src/reduction.cxx at line 20 [17.870]

B:BD[17.1208] → [10.3862:3933]

reduction<T> reduce_array(const Geometry &geom, const T *restrict var,

[17.1208]

[10.3933]

reduction<T> reduce_array(const Geometry &geom, const T *restrict ptr,

Replacement in AMReX/src/reduction.cxx at line 35 [17.870]

B:BD[17.1617] → [10.4133:4194]

        red = reduction<T>(red, reduction<T>(dV, var[idx]));

[17.1617]

[17.1674]

        red = reduction<T>(red, reduction<T>(dV, ptr[idx]));

Insertion in AMReX/src/reduction.cxx at line 43 [17.870]
[17.1763]
[17.1763]
```
  DECLARE_CCTK_PARAMETERS;
```

Replacement in AMReX/src/reduction.cxx at line 52 [17.870]

B:BD[17.2007] → [17.2007:2100]

      auto mfitinfo =
          MFItInfo().SetDynamic(true).EnableTiling({1024000, 16, 32});

[17.2007]

[17.2100]

      auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
          {max_tile_size_x, max_tile_size_y, max_tile_size_z});

Insertion in AMReX/src/reduction.cxx at line 62 [17.870]

[17.2408]

        // Note: This excludes ghosts, it's not the proper Cactus lsh

Replacement in AMReX/src/reduction.cxx at line 67 [17.870]

B:BD[17.2582] → [17.2582:2660]

        const CCTK_REAL *restrict var = vars.ptr(imin.x, imin.y, imin.z, vi);

[17.2582]

[17.2660]

        const CCTK_REAL *restrict ptr = vars.ptr(imin.x, imin.y, imin.z, vi);

Replacement in AMReX/src/reduction.cxx at line 69 [17.870]

B:BD[17.2661] → [10.4269:4319]

        red += reduce_array(geom, var, ash, lsh);

[17.2661]

[17.2705]

#warning "TODO: Skip refined regions"
        red += reduce_array(geom, ptr, ash, lsh);

Replacement in AMReX/src/schedule.cxx at line 102 [12.531]
B:BD[6.1302] → [6.1302:1334]
```
  for (int d = 0; d < dim; ++d)
```
[6.1302]
[10.5140]
```
  for (int d = 0; d < dim; ++d) {
```
Deletion in AMReX/src/schedule.cxx at line 104 [12.531]
∅:D[10.5196] → [12.1888:1920]
∅:D[6.1376] → [12.1888:1920]
∅:D[19.920] → [12.1888:1920]
B:BD[12.1888] → [12.1888:1920]
```
  for (int d = 0; d < dim; ++d)
```
Insertion in AMReX/src/schedule.cxx at line 105 [12.531]
[6.1418]
[6.1683]
```
  }
```

Insertion in AMReX/src/schedule.cxx at line 108 [12.531]

[6.1714]

[10.5197]

  // CCTK_REAL mindx = 1.0 / 0.0;
  // const int numlevels = ghext->amrcore->finestLevel() + 1;
  // for (int level = 0; level < numlevels; ++level) {
  //   const Geometry &geom = ghext->amrcore->Geom(level);
  //   const CCTK_REAL *restrict dx = geom.CellSize();
  //   for (int d = 0; d < dim; ++d)
  //     mindx = fmin(mindx, dx[d]);
  // }

Replacement in AMReX/src/schedule.cxx at line 117 [12.531]

B:BD[10.5228] → [8.4923:4982]

∅:D[8.4982] → [18.421:473]

B:BD[10.5287] → [18.421:473]

B:BD[18.473] → [8.4983:5039]

∅:D[8.5039] → [10.5394:5518]

B:BD[10.5394] → [10.5394:5518]

  const int numlevels = ghext->amrcore->finestLevel() + 1;
  for (int level = 0; level < numlevels; ++level) {
    const Geometry &geom = ghext->amrcore->Geom(level);
    const CCTK_REAL *restrict dx = geom.CellSize();
    for (int d = 0; d < dim; ++d)
      mindx = fmin(mindx, dx[d]);
  }

[10.5228]

[12.1959]

  for (int d = 0; d < dim; ++d)
    mindx = fmin(mindx, dx[d]);
  mindx = mindx / (1 << (max_num_levels - 1));

Replacement in AMReX/src/schedule.cxx at line 140 [12.531]
B:BD[12.2155] → [6.1877:1907]
```
  // TODO: Get this from mfab
```
[12.2155]
[12.2155]
```
  // TODO: Get this from mfab (mfab.fb_ghosts)
```

Replacement in AMReX/src/schedule.cxx at line 163 [12.531]

B:BD[12.2680] → [10.5781:5846]

  // Offset between this level's and the coarsest level's origin

[12.2680]

[10.5846]

  // Offset between this level's and the coarsest level's origin as multiple of
  // the grid spacing

Replacement in AMReX/src/schedule.cxx at line 192 [12.531]

B:BD[12.3361] → [20.517:584]

    cctkGH->cctk_lsh[d] = bx[orient(d, 1)] - bx[orient(d, 0)] + 1;

[12.3361]

[12.3424]

    cctkGH->cctk_lsh[d] = gbx[orient(d, 1)] - gbx[orient(d, 0)] + 1;

Replacement in AMReX/src/schedule.cxx at line 200 [12.531]

∅:D[9.1816] → [12.3594:3639]

B:BD[12.3594] → [12.3594:3639]

B:BD[12.3671] → [12.3671:3716]

    cctkGH->cctk_lbnd[d] = bx[orient(d, 0)];
    cctkGH->cctk_ubnd[d] = bx[orient(d, 1)];

[9.1816]

[9.1817]

    cctkGH->cctk_lbnd[d] = gbx[orient(d, 0)];
    cctkGH->cctk_ubnd[d] = gbx[orient(d, 1)];

Replacement in AMReX/src/schedule.cxx at line 207 [12.531]

B:BD[12.3797] → [20.655:731]

      cctkGH->cctk_bbox[2 * d + f] = bx[orient(d, f)] != gbx[orient(d, f)];

[12.3797]

[20.731]

      cctkGH->cctk_bbox[2 * d + f] = bx[orient(d, f)] == gbx[orient(d, f)];

Replacement in AMReX/src/schedule.cxx at line 210 [12.531]
B:BD[20.760] → [20.760:792]
```
  const Dim3 imin = lbound(bx);
```
[20.760]
[21.2429]
```
  const Dim3 imin = lbound(gbx);
```
Replacement in AMReX/src/schedule.cxx at line 299 [12.531]
B:BD[19.1147] → [19.1147:1181]
```
    for (int d = 0; d < dim; ++d)
```
[19.1147]
[19.1181]
```
    for (int d = 0; d < dim; ++d) {
```
Replacement in AMReX/src/schedule.cxx at line 301 [12.531]
B:BD[19.1224] → [19.1224:1258]
B:BD[19.1258] → [10.6057:6115]
∅:D[10.6115] → [19.1316:1350]
B:BD[19.1316] → [19.1316:1350]
```
    for (int d = 0; d < dim; ++d)
      x0[d] = cctkGH->cctk_origin_space[d] - 0.0 * dx[d];
    for (int d = 0; d < dim; ++d)
```
[19.1224]
[10.6116]
```
      x0[d] = cctkGH->cctk_origin_space[d];
```
Insertion in AMReX/src/schedule.cxx at line 303 [12.531]
[10.6154]
[19.1391]
```
    }
```

Replacement in AMReX/src/schedule.cxx at line 309 [12.531]

B:BD[19.1642] → [19.1642:1698]

    CCTK_VINFO("  dx=[%g,%g,%g]", dx[0], dx[1], dx[2]);

[19.1642]

[19.1698]

    CCTK_VINFO("  base dx=[%g,%g,%g]", dx[0], dx[1], dx[2]);

Insertion in AMReX/src/schedule.cxx at line 322 [12.531]
[5.4623]
[7.2545]
```
    CCTK_VINFO("Recovering from checkpoint...");
```
Insertion in AMReX/src/schedule.cxx at line 342 [12.531]
[5.5085]
[5.5085]
```
    CCTK_VINFO("Setting up initial conditions...");
```
Insertion in AMReX/src/schedule.cxx at line 346 [12.531]
[8.5164]
[5.5154]
```
      CCTK_VINFO("Initializing level %d...", current_level);
```

Deletion in AMReX/src/schedule.cxx at line 358 [12.531]

B:BD[5.5624] → [8.5301:5500]

      cout << "I old_numlevels=" << old_numlevels << "\n";
      cout << "I new_numlevels=" << new_numlevels << "\n";
      cout << "I max_numlevels=" << int(ghext->amrcore->maxLevel() + 1) << "\n";

Replacement in AMReX/src/schedule.cxx at line 359 [12.531]

B:BD[5.5853] → [8.5501:5543]

      if (new_numlevels <= old_numlevels)

[5.5853]

[5.5895]

      const bool did_create_new_level = new_numlevels > old_numlevels;
      if (!did_create_new_level)

Insertion in AMReX/src/schedule.cxx at line 365 [12.531]

[7.3122]

  CCTK_VINFO("Initialized %d levels", int(ghext->leveldata.size()));

Deletion in AMReX/src/schedule.cxx at line 432 [12.531]
B:BD[22.2105] → [6.2134:2170]
```
  // TODO: Get ghost_size from mfab
```
Insertion in AMReX/src/schedule.cxx at line 451 [12.531]
[7.5445]
[7.5445]
```
  CCTK_VINFO("Starting evolution...");
```
Insertion in AMReX/src/schedule.cxx at line 477 [12.531]
[7.6001]
[7.6001]
```
  CCTK_VINFO("Shutting down...");
```
Insertion in AMReX/src/schedule.cxx at line 489 [12.531]
[12.4810]
[12.4810]
```
  DECLARE_CCTK_PARAMETERS;
```

Replacement in AMReX/src/schedule.cxx at line 497 [12.531]

B:BD[12.4926] → [12.4926:5065]

  CCTK_VINFO("CallFunction [%d] %s: %s::%s", cctkGH->cctk_iteration,
             attribute->where, attribute->thorn, attribute->routine);

[12.4926]

[12.5065]

  if (verbose)
    CCTK_VINFO("CallFunction [%d] %s: %s::%s", cctkGH->cctk_iteration,
               attribute->where, attribute->thorn, attribute->routine);

Replacement in AMReX/src/schedule.cxx at line 513 [12.531]

∅:D[21.2939] → [23.3616:3713]

B:BD[23.3616] → [23.3616:3713]

        auto mfitinfo =
            MFItInfo().SetDynamic(true).EnableTiling({1024000, 16, 32});

[21.2939]

[23.3713]

        auto mfitinfo = MFItInfo().SetDynamic(true).EnableTiling(
            {max_tile_size_x, max_tile_size_y, max_tile_size_z});

Insertion in AMReX/src/schedule.cxx at line 567 [12.531]
[16.2313]
[16.2313]
```
#warning "TODO: during evolution, sync only one time level"
```

Replacement in AMReX/src/schedule.cxx at line 579 [12.531]

B:BD[18.752] → [18.752:873]

        // Refined level: Prolongate from next coarser level, and then
        // copy from adjacent boxes on same level

[18.752]

[18.873]

        // Refined level: Prolongate boundaries from next coarser
        // level, and copy from adjacent boxes on same level

Replacement in AMReX/src/schedule.cxx at line 589 [12.531]
B:BD[18.1187] → [18.1187:1216]
```
        CellBilinear interp;
```
[18.1187]
[18.1216]
```
        // CellBilinear interp;
```

Replacement in AMReX/src/schedule.cxx at line 594 [12.531]

B:BD[18.1461] → [18.1461:1618]

B:BD[18.1618] → [8.5613:5755]

∅:D[8.5755] → [18.1760:1792]

B:BD[18.1760] → [18.1760:1792]

        for (int tl = 0; tl < sync_tl; ++tl)
          InterpFromCoarseLevel(
              *groupdata.mfab.at(tl), 0.0, *coarsegroupdata.mfab.at(tl), 0, 0,
              groupdata.numvars, ghext->amrcore->Geom(level - 1),
              ghext->amrcore->Geom(level), cphysbc, 0, fphysbc, 0, reffact,
              &interp, bcs, 0);

[18.1461]

[18.1792]

        for (int tl = 0; tl < sync_tl; ++tl) {
          // InterpFromCoarseLevel(
          //     *groupdata.mfab.at(tl), 0.0, *coarsegroupdata.mfab.at(tl), 0,
          //     0, groupdata.numvars, ghext->amrcore->Geom(level - 1),
          //     ghext->amrcore->Geom(level), cphysbc, 0, fphysbc, 0, reffact,
          //     &interp, bcs, 0);
#warning "TODO: make copy of fine level"
          FillPatchTwoLevels(
              *groupdata.mfab.at(tl), 0.0, {&*coarsegroupdata.mfab.at(tl)},
              {0.0}, {&*groupdata.mfab.at(tl)}, {0.0}, 0, 0, groupdata.numvars,
              ghext->amrcore->Geom(level - 1), ghext->amrcore->Geom(level),
              cphysbc, 0, fphysbc, 0, reffact, &cell_bilinear_interp, bcs, 0);
        }

Insertion in AMReX/src/schedule.cxx at line 622 [12.531]
[3.832]
[3.832]
```
#warning "TODO: during evolution, restrict only one time level"
```
Replacement in AMReX/src/schedule.cxx at line 626 [12.531]
B:BD[3.956] → [3.956:1003]
```
    for (int tl = 0; tl < restrict_tl; ++tl) {
```
[3.956]
[3.1003]
```
    for (int tl = 0; tl < restrict_tl; ++tl)
```
Deletion in AMReX/src/schedule.cxx at line 631 [12.531]
B:BD[3.1253] → [3.1253:1259]
```
    }
```

Insertion in WaveToyAMReX/README at line 10 [4.14802]

[4.15080]


2. First-order equations:
Scalar wave equation:
  phi,tt = phi,xx
Introduce new variable:
  psi = phi,t
New set of equations:
  phi,t = psi
  psi,t = phi,tt = phi,xx
Discretization:
  phi_0 = phi_1 + dt psi_1
  psi_0 = psi_1 + dt phi_1,xx
3. Fancy equations:
Scalar wave equation:
  phi,tt = phi,xx
Direct second order discretization:
  phi_0 = 2 phi_1 - phi_2 + dt^2 phi_1,xx
Introduce new variable:
  psi = phi,t
New set of equations:
  phi,t = psi
  psi,t = phi,tt = phi,xx
New, equivalent discretization:
  phi_0 - phi_1 = dt psi_1
  phi_0 - 2 phi_1 + phi_2 = dt^2 phi_1,xx
  dt psi_1 - dt psi_2 = dt^2 phi_1,xx
  psi_1 = psi_2 + dt^2 phi_1,xx
  phi_0 = phi_1 + dt psi_1
  psi_0 = psi_1 + dt phi_0,xx
        = psi_1 + dt phi_1,xx + dt^2 psi_1,xx

Insertion in WaveToyAMReX/interface.ccl at line 12 [4.10196]
[7.6916]
[7.6916]
```
  # psi
```
Insertion in WaveToyAMReX/interface.ccl at line 15 [4.10196]
[24.1617]
[7.6956]
```
CCTK_REAL energy TYPE=gf
{
  eps
} "Energy density for wave equation"
```
Replacement in WaveToyAMReX/interface.ccl at line 22 [4.10196]
B:BD[20.2030] → [20.2030:2036]
```
  err
```
[20.2006]
[7.6979]
```
  phierr
  # psierr
```

Insertion in WaveToyAMReX/schedule.ccl at line 19 [4.9594]

[4.9911]

[24.1995]


SCHEDULE WaveToyAMReX_Energy AT analysis
{
  LANG: C
  SYNC: energy
} "Calculate energy density for the wave equation"

Insertion in WaveToyAMReX/src/wavetoy.cxx at line 6 [4.37]
[4.250]
[4.271]
```
#include <iostream>
```

Insertion in WaveToyAMReX/src/wavetoy.cxx at line 20 [4.37]

[4.559]


// Time derivative
template <typename T> auto timederiv(T f(T t, T x, T y, T z), T dt) {
  return [=](T t, T x, T y, T z) {
    // return (f(t + dt, x, y, z) - f(t, x, y, z)) / dt;
    return (f(t, x, y, z) - f(t - dt, x, y, z)) / dt;
  };
}

Replacement in WaveToyAMReX/src/wavetoy.cxx at line 31 [4.37]

B:BD[4.632] → [20.2101:2129]

∅:D[20.2129] → [4.660:724]

B:BD[4.660] → [4.660:724]

  T Lx = 2, Ly = 2, Lz = 2;
  T kx = 2 * M_PI / Lx, ky = 2 * M_PI / Ly, kz = 2 * M_PI / Lz;

[4.632]

[4.724]

  T kx = 2 * M_PI / 2, ky = 2 * M_PI / 2, kz = 2 * M_PI / 2;
  // T kx = 2 * M_PI / 2, ky = 0, kz = 0;

Replacement in WaveToyAMReX/src/wavetoy.cxx at line 51 [4.37]

B:BD[4.4712] → [16.2521:2607]

  for (int k = 0; k < cctk_lsh[2]; ++k) {
    for (int j = 0; j < cctk_lsh[1]; ++j) {

[4.4712]

[4.4979]

  int imin[dim], imax[dim];
  for (int d = 0; d < dim; ++d) {
    imin[d] = cctk_bbox[2 * d] ? 0 : cctk_nghostzones[d];
    imax[d] = cctk_lsh[d] - (cctk_bbox[2 * d + 1] ? 0 : cctk_nghostzones[d]);
  }
  for (int k = imin[2]; k < imax[2]; ++k) {
    for (int j = imin[1]; j < imax[1]; ++j) {

Replacement in WaveToyAMReX/src/wavetoy.cxx at line 60 [4.37]
B:BD[4.4996] → [20.2913:2959]
```
      for (int i = 0; i < cctk_lsh[0]; ++i) {
```
[4.4996]
[20.2959]
```
      for (int i = imin[0]; i < imax[0]; ++i) {
```
Insertion in WaveToyAMReX/src/wavetoy.cxx at line 67 [4.37]
[20.3279]
[24.3756]
```
        // psi[idx] = timederiv(standing, dt)(t, x, y, z);
```
Insertion in WaveToyAMReX/src/wavetoy.cxx at line 92 [4.37]
[5.7295]
[3.1267]
```
        // CCTK_REAL r = fabs(x);
```

Insertion in WaveToyAMReX/src/wavetoy.cxx at line 107 [4.37]

[10.6575]

  }
  int imin[dim], imax[dim];
  for (int d = 0; d < dim; ++d) {
    imin[d] = cctk_bbox[2 * d] ? 0 : cctk_nghostzones[d];
    imax[d] = cctk_lsh[d] - (cctk_bbox[2 * d + 1] ? 0 : cctk_nghostzones[d]);

Replacement in WaveToyAMReX/src/wavetoy.cxx at line 119 [4.37]

B:BD[4.6652] → [16.2619:2705]

  for (int k = 0; k < cctk_lsh[2]; ++k) {
    for (int j = 0; j < cctk_lsh[1]; ++j) {

[4.6652]

[4.6740]

  for (int k = imin[2]; k < imax[2]; ++k) {
    for (int j = imin[1]; j < imax[1]; ++j) {

Replacement in WaveToyAMReX/src/wavetoy.cxx at line 122 [4.37]
B:BD[4.6757] → [20.3527:3573]
```
      for (int i = 0; i < cctk_lsh[0]; ++i) {
```
[4.6757]
[20.3573]
```
      for (int i = imin[0]; i < imax[0]; ++i) {
```

Replacement in WaveToyAMReX/src/wavetoy.cxx at line 130 [4.37]

B:BD[24.5321] → [24.5321:5373]

        phi[idx] = -phi_p_p[idx] + 2 * phi_p[idx] +

[24.5321]

[24.5373]

        // CCTK_REAL ddx_psi =
        //     (psi_p[idx - di] - 2 * psi_p[idx] + psi_p[idx + di]) /
        //     sqr(dx[0]);
        // CCTK_REAL ddy_psi =
        //     (psi_p[idx - dj] - 2 * psi_p[idx] + psi_p[idx + dj]) /
        //     sqr(dx[1]);
        // CCTK_REAL ddz_psi =
        //     (psi_p[idx - dk] - 2 * psi_p[idx] + psi_p[idx + dk]) /
        //     sqr(dx[2]);
        // phi[idx] = phi_p[idx] + dt * psi_p[idx];
        // psi[idx] = psi_p[idx] + dt * (ddx_phi + ddy_phi + ddz_phi) +
        //            sqr(dt) * (ddx_psi + ddy_psi + ddz_psi);
        // phi[idx] = phi_p[idx] + dt * psi_p[idx];
        // psi[idx] = psi_p[idx] + dt * (ddx_phi + ddy_phi + ddz_phi);
        phi[idx] = 2 * phi_p[idx] - phi_p_p[idx] +

Insertion in WaveToyAMReX/src/wavetoy.cxx at line 146 [4.37]

[24.5433]

      }
    }
  }
}
extern "C" void WaveToyAMReX_Energy(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS;
  DECLARE_CCTK_PARAMETERS;
  const CCTK_REAL dt = CCTK_DELTA_TIME;
  CCTK_REAL dx[dim];
  for (int d = 0; d < dim; ++d)
    dx[d] = CCTK_DELTA_SPACE(d);
  int imin[dim], imax[dim];
  for (int d = 0; d < dim; ++d) {
    imin[d] = cctk_bbox[2 * d] ? 0 : cctk_nghostzones[d];
    imax[d] = cctk_lsh[d] - (cctk_bbox[2 * d + 1] ? 0 : cctk_nghostzones[d]);
  }
  constexpr int di = 1;
  const int dj = di * cctk_ash[0];
  const int dk = dj * cctk_ash[1];
  for (int k = imin[2]; k < imax[2]; ++k) {
    for (int j = imin[1]; j < imax[1]; ++j) {
#pragma omp simd
      for (int i = imin[0]; i < imax[0]; ++i) {
        const int idx = CCTK_GFINDEX3D(cctkGH, i, j, k);
        CCTK_REAL dt_phi = (phi[idx] - phi_p[idx]) / dt;
        CCTK_REAL dx_phi = (phi[idx + di] - phi[idx - di]) / (2 * dx[0]);
        CCTK_REAL dy_phi = (phi[idx + dj] - phi[idx - dj]) / (2 * dx[1]);
        CCTK_REAL dz_phi = (phi[idx + dk] - phi[idx - dk]) / (2 * dx[2]);
        eps[idx] = (sqr(dt_phi) + sqr(dx_phi) + sqr(dy_phi) + sqr(dz_phi)) / 2;

Insertion in WaveToyAMReX/src/wavetoy.cxx at line 190 [4.37]
[20.3664]
[10.6580]
```
  // const CCTK_REAL dt = CCTK_DELTA_TIME;
```

Replacement in WaveToyAMReX/src/wavetoy.cxx at line 205 [4.37]

B:BD[20.4140] → [20.4140:4192]

        err[idx] = phi[idx] - standing(t, x, y, z);

[20.4140]

[24.6606]

        phierr[idx] = phi[idx] - standing(t, x, y, z);
        // psierr[idx] = psi[idx] - timederiv(standing, dt)(t, x, y, z);

File addition: amrex-1.par (-rw-r--r--)

[25.18]

ActiveThorns = "
    AMReX
    IOUtil
    WaveToyAMReX
"
$maxlev = 1
$ncells = 128
Cactus::cctk_show_schedule = no
Cactus::cctk_itlast = $ncells * 2 ** ($maxlev - 1) * 2
AMReX::ncells_x = $ncells
AMReX::ncells_y = 1
AMReX::ncells_z = 1
AMReX::blocking_factor_x = 1
AMReX::blocking_factor_y = 1
AMReX::blocking_factor_z = 1
AMReX::max_grid_size_x = 1048576
AMReX::max_grid_size_y = 1048576
AMReX::max_grid_size_z = 1048576
AMReX::max_num_levels = $maxlev
AMReX::dtfac = 0.5 / 4
IO::out_dir = "wavetoy-1"
IO::out_every = $ncells * 2 ** ($maxlev - 1) / 2

File addition: amrex-2.par (-rw-r--r--)

[25.18]

ActiveThorns = "
    AMReX
    IOUtil
    WaveToyAMReX
"
$maxlev = 2
$ncells = 128
Cactus::cctk_show_schedule = no
Cactus::cctk_itlast = $ncells * 2 ** ($maxlev - 1) * 2
AMReX::ncells_x = $ncells
AMReX::ncells_y = 1
AMReX::ncells_z = 1
AMReX::blocking_factor_x = 1
AMReX::blocking_factor_y = 1
AMReX::blocking_factor_z = 1
AMReX::max_grid_size_x = 1048576
AMReX::max_grid_size_y = 1048576
AMReX::max_grid_size_z = 1048576
AMReX::max_num_levels = $maxlev
AMReX::dtfac = 0.5 / 2
IO::out_dir = "wavetoy-2"
IO::out_every = $ncells * 2 ** ($maxlev - 1) / 2

File addition: amrex-3.par (-rw-r--r--)

[25.18]

ActiveThorns = "
    AMReX
    IOUtil
    WaveToyAMReX
"
$maxlev = 3
$ncells = 128
Cactus::cctk_show_schedule = no
Cactus::cctk_itlast = $ncells * 2 ** ($maxlev - 1) * 2
AMReX::ncells_x = $ncells
AMReX::ncells_y = 1
AMReX::ncells_z = 1
AMReX::blocking_factor_x = 1
AMReX::blocking_factor_y = 1
AMReX::blocking_factor_z = 1
AMReX::max_grid_size_x = 1048576
AMReX::max_grid_size_y = 1048576
AMReX::max_grid_size_z = 1048576
AMReX::max_num_levels = $maxlev
AMReX::dtfac = 0.5
IO::out_dir = "wavetoy-3"
IO::out_every = $ncells * 2 ** ($maxlev - 1) / 2

Insertion in azure-pipelines/amrex.par at line 3 [26.12]
[26.40]
[26.40]
```
    IOUtil
```
Insertion in azure-pipelines/amrex.par at line 6 [26.12]
[26.59]
[26.59]
```
 
$nlevels = 3
$ncells = 128
```

Replacement in azure-pipelines/amrex.par at line 12 [26.12]

B:BD[10.6717] → [8.5803:5869]

AMReX::ncells_x = 128
AMReX::ncells_y = 128
AMReX::ncells_z = 128

[10.6717]

[10.6783]

Cactus::cctk_itlast = $ncells * 2 ** ($nlevels - 1) * 2
AMReX::ncells_x = $ncells
AMReX::ncells_y = $ncells
AMReX::ncells_z = $ncells

Replacement in azure-pipelines/amrex.par at line 18 [26.12]
B:BD[10.6784] → [8.5870:5896]
```
AMReX::max_num_levels = 3
```
[10.6784]
[8.5896]
```
AMReX::max_num_levels = $nlevels
```
Replacement in azure-pipelines/amrex.par at line 20 [26.12]
B:BD[8.5897] → [8.5897:5923]
```
Cactus::cctk_itlast = 128
```
[8.5897]
[15.1850]
```
AMReX::dtfac = 0.5
```
Replacement in azure-pipelines/amrex.par at line 22 [26.12]
B:BD[15.1851] → [27.27:50]
```
AMReX::out_every = 128
```
[15.1851]
```
IO::out_dir = "wavetoy"
IO::out_every = $ncells * 2 ** ($nlevels - 1) / 2
```

AMR evolution should work now

Dependencies

In channels

Change contents

Insertion in AMReX/interface.ccl at line 4 [4.17976]

Insertion in AMReX/param.ccl at line 4 [4.17918]

Replacement in AMReX/param.ccl at line 45 [4.17918]

Replacement in AMReX/param.ccl at line 50 [4.17918]

Insertion in AMReX/param.ccl at line 54 [4.17918]

Replacement in AMReX/param.ccl at line 66 [4.17918]

Insertion in AMReX/param.ccl at line 111 [4.17918]

Replacement in AMReX/param.ccl at line 119 [4.17918]

Replacement in AMReX/param.ccl at line 121 [4.17918]

Replacement in AMReX/param.ccl at line 126 [4.17918]

Replacement in AMReX/param.ccl at line 128 [4.17918]

Insertion in AMReX/param.ccl at line 140 [4.17918]

Insertion in AMReX/src/driver.cxx at line 11 [4.15972]

Insertion in AMReX/src/driver.cxx at line 70 [4.15972]

Replacement in AMReX/src/driver.cxx at line 76 [4.15972]

Replacement in AMReX/src/driver.cxx at line 104 [4.15972]

Insertion in AMReX/src/driver.cxx at line 115 [4.15972]

Replacement in AMReX/src/driver.cxx at line 119 [4.15972]

Insertion in AMReX/src/driver.cxx at line 126 [4.15972]

Replacement in AMReX/src/driver.cxx at line 132 [4.15972]

Replacement in AMReX/src/driver.cxx at line 143 [4.15972]

Replacement in AMReX/src/driver.cxx at line 151 [4.15972]

Replacement in AMReX/src/driver.cxx at line 158 [4.15972]

Replacement in AMReX/src/driver.cxx at line 164 [4.15972]

Replacement in AMReX/src/driver.cxx at line 175 [4.15972]

Replacement in AMReX/src/driver.cxx at line 213 [4.15972]

Replacement in AMReX/src/driver.cxx at line 233 [4.15972]

Insertion in AMReX/src/driver.cxx at line 252 [4.15972]

Deletion in AMReX/src/driver.cxx at line 265 [4.15972]

Replacement in AMReX/src/driver.cxx at line 266 [4.15972]

Replacement in AMReX/src/driver.cxx at line 289 [4.15972]

Deletion in AMReX/src/driver.cxx at line 301 [4.15972]

Replacement in AMReX/src/driver.cxx at line 310 [4.15972]

Replacement in AMReX/src/driver.cxx at line 346 [4.15972]

Replacement in AMReX/src/driver.cxx at line 358 [4.15972]

Insertion in AMReX/src/io.cxx at line 10 [14.107]

Insertion in AMReX/src/io.cxx at line 14 [14.107]

Insertion in AMReX/src/io.cxx at line 26 [14.107]

Insertion in AMReX/src/io.cxx at line 86 [14.107]

Insertion in AMReX/src/io.cxx at line 164 [14.107]

Replacement in AMReX/src/io.cxx at line 184 [14.107]

Insertion in AMReX/src/io.cxx at line 211 [14.107]

Insertion in AMReX/src/io.cxx at line 219 [14.107]

Replacement in AMReX/src/io.cxx at line 235 [14.107]

Insertion in AMReX/src/reduction.cxx at line 3 [17.870]

Replacement in AMReX/src/reduction.cxx at line 20 [17.870]

Replacement in AMReX/src/reduction.cxx at line 35 [17.870]

Insertion in AMReX/src/reduction.cxx at line 43 [17.870]

Replacement in AMReX/src/reduction.cxx at line 52 [17.870]

Insertion in AMReX/src/reduction.cxx at line 62 [17.870]

Replacement in AMReX/src/reduction.cxx at line 67 [17.870]

Replacement in AMReX/src/reduction.cxx at line 69 [17.870]

Replacement in AMReX/src/schedule.cxx at line 102 [12.531]

Deletion in AMReX/src/schedule.cxx at line 104 [12.531]

Insertion in AMReX/src/schedule.cxx at line 105 [12.531]

Insertion in AMReX/src/schedule.cxx at line 108 [12.531]

Replacement in AMReX/src/schedule.cxx at line 117 [12.531]

Replacement in AMReX/src/schedule.cxx at line 140 [12.531]

Replacement in AMReX/src/schedule.cxx at line 163 [12.531]

Replacement in AMReX/src/schedule.cxx at line 192 [12.531]

Replacement in AMReX/src/schedule.cxx at line 200 [12.531]

Replacement in AMReX/src/schedule.cxx at line 207 [12.531]

Replacement in AMReX/src/schedule.cxx at line 210 [12.531]

Replacement in AMReX/src/schedule.cxx at line 299 [12.531]

Replacement in AMReX/src/schedule.cxx at line 301 [12.531]

Insertion in AMReX/src/schedule.cxx at line 303 [12.531]

Replacement in AMReX/src/schedule.cxx at line 309 [12.531]

Insertion in AMReX/src/schedule.cxx at line 322 [12.531]

Insertion in AMReX/src/schedule.cxx at line 342 [12.531]

Insertion in AMReX/src/schedule.cxx at line 346 [12.531]

Deletion in AMReX/src/schedule.cxx at line 358 [12.531]

Replacement in AMReX/src/schedule.cxx at line 359 [12.531]

Insertion in AMReX/src/schedule.cxx at line 365 [12.531]

Deletion in AMReX/src/schedule.cxx at line 432 [12.531]

Insertion in AMReX/src/schedule.cxx at line 451 [12.531]

Insertion in AMReX/src/schedule.cxx at line 477 [12.531]