eschnett/cactusamrex - Change E3MBKFT4GEFDAGZQQW4OROY5F6FWC46G6MRH54GDYTGO7O5YSRIAC

Make HydroToy work with multiple levels

Created by Erik Schnetter on October 10, 2019

E3MBKFT4GEFDAGZQQW4OROY5F6FWC46G6MRH54GDYTGO7O5YSRIAC

Dependencies

In channels

master

main

Change contents

Replacement in AMReX/src/driver.cxx at line 417 [6.15972]

B:BD[5.2340] → [5.2340:2424]

          valid.at(tl).valid_int = false;
          valid.at(tl).valid_bnd = false;

[5.2340]

[5.2424]

          // do nothing

Replacement in AMReX/src/driver.cxx at line 448 [6.15972]
B:BD[2.479] → [2.479:516]
```
            valid.at(tl).valid_int =
```
[2.479]
[2.516]
```
            valid.at(vi).valid_int =
```

Replacement in AMReX/src/driver.cxx at line 453 [6.15972]

B:BD[2.762] → [2.762:806]

            valid.at(tl).valid_bnd = false;

[2.762]

[2.806]

            valid.at(vi).valid_bnd = false;

Replacement in AMReX/src/schedule.cxx at line 792 [8.531]

B:BD[7.7955] → [7.7955:8031]

  for (int level = int(ghext->leveldata.size()) - 2; level >= 0; --level) {

[7.7955]

[9.247]

  for (int level = int(ghext->leveldata.size()) - 2; level >= 0; --level)

Replacement in AMReX/src/schedule.cxx at line 794 [8.531]

B:BD[9.268] → [7.8032:8133]

    current_level = level;
    CCTK_Traverse(cctkGH, "CCTK_POSTRESTRICT");
  }
  current_level = -1;

[9.268]

  CCTK_Traverse(cctkGH, "CCTK_POSTRESTRICT");

Insertion in AMReX/src/schedule.cxx at line 875 [8.531]
[5.5608]
[5.5608]
```
  assert(current_level == -1);
```
Insertion in AMReX/src/schedule.cxx at line 899 [8.531]
[10.2105]
[11.2656]
```
  assert(current_level == -1);
```
Insertion in AMReX/src/schedule.cxx at line 958 [8.531]
[12.874]
[12.874]
```
      assert(current_level == -1);
```

Replacement in AMReX/src/schedule.cxx at line 975 [8.531]

B:BD[7.8634] → [7.8634:8712]

    for (int level = int(ghext->leveldata.size()) - 2; level >= 0; --level) {

[7.8634]

[9.360]

    for (int level = int(ghext->leveldata.size()) - 2; level >= 0; --level)

Replacement in AMReX/src/schedule.cxx at line 977 [8.531]

B:BD[9.383] → [7.8713:8822]

      current_level = level;
      CCTK_Traverse(cctkGH, "CCTK_POSTRESTRICT");
    }
    current_level = -1;

[9.383]

[13.5723]

    CCTK_Traverse(cctkGH, "CCTK_POSTRESTRICT");

Replacement in AMReX/src/schedule.cxx at line 1167 [8.531]

B:BD[14.392] → [11.3098:3152]

  for (auto &restrict leveldata : ghext->leveldata) {

[14.268]

[15.4086]

  const int min_level = current_level == -1 ? 0 : current_level;
  const int max_level =
      current_level == -1 ? ghext->leveldata.size() : current_level + 1;
  for (int level = min_level; level < max_level; ++level) {
    auto &restrict leveldata = ghext->leveldata.at(level);

Replacement in AMReX/src/schedule.cxx at line 1287 [8.531]

B:BD[16.161] → [5.8748:8792]

      // Only restrict valid grid functions

[16.161]

[5.8792]

      // Only restrict valid grid functions. Restriction only uses the interior.

Deletion in AMReX/src/schedule.cxx at line 1291 [8.531]

B:BD[3.2992] → [3.2992:3063]

                       !finegroupdata.valid.at(tl).at(vi).valid_bnd &&

Deletion in AMReX/src/schedule.cxx at line 1304 [8.531]
B:BD[3.3185] → [3.3185:3249]
```
                 finegroupdata.valid.at(tl).at(vi).valid_bnd &&
```

Insertion in AMReX/src/schedule.cxx at line 1305 [8.531]

[5.9558]

[3.3250]

          poison_invalid(fineleveldata, finegroupdata, vi, tl);

Insertion in AMReX/src/schedule.cxx at line 1307 [8.531]
[3.3311]
[3.3311]
```
          poison_invalid(leveldata, groupdata, vi, tl);
```

Replacement in AMReX/src/schedule.cxx at line 1332 [8.531]

B:BD[5.10729] → [3.3817:3937]

              finegroupdata.valid.at(tl).at(vi).valid_int &&
              finegroupdata.valid.at(tl).at(vi).valid_bnd;

[5.10729]

[5.10851]

              finegroupdata.valid.at(tl).at(vi).valid_int;

Replacement in HydroToyAMReX/par/hydro-shocktube.par at line 20 [4.22]
B:BD[4.310] → [4.310:331]
```
AMReX::verbose = yes
```
[4.310]
[4.331]
```
AMReX::verbose = no
```
Replacement in HydroToyAMReX/par/hydro-sphericalshock.par at line 10 [4.1031]
B:BD[4.1141] → [4.1141:1154]
```
$nlevels = 1
```
[4.1141]
[4.1154]
```
$nlevels = 3
```
Replacement in HydroToyAMReX/par/hydro-sphericalshock.par at line 36 [4.1031]
B:BD[4.1593] → [4.1593:1629]
```
AMReX::regrid_error_threshold = 0.1
```
[4.1593]
[4.1629]
```
AMReX::regrid_error_threshold = 0.01
```
Insertion in HydroToyAMReX/schedule.ccl at line 6 [18.6006]
[17.2929]
[18.6054]
Insertion in HydroToyAMReX/schedule.ccl at line 9 [18.6006]
[17.2931]
[17.2931]
```
# Initial conditions
```

Replacement in HydroToyAMReX/schedule.ccl at line 18 [18.6006]

B:BD[17.2967] → [17.2967:3140]

# SCHEDULE HydroToyAMReX_EstimateError AT postinitial
# {
#   LANG: C
#   WRITES: AMReX::regrid_error(interior)
# } "Estimate local error for regridding initial conditions"

[17.2967]

[17.3140]

SCHEDULE HydroToyAMReX_Pressure AT initial AFTER HydroToyAMReX_Boundaries
{
  LANG: C
  READS: conserved(everywhere)
  WRITES: primitive(everywhere)
} "Calculate pressure"
SCHEDULE HydroToyAMReX_Fluxes AT initial AFTER HydroToyAMReX_Pressure
{
  LANG: C
  READS: conserved(everywhere)
  READS: primitive(everywhere)
  WRITES: flux_x(interior) flux_y(interior) flux_z(interior)
  SYNC: flux_x flux_y flux_z
} "Calculate the hydro fluxes"
SCHEDULE HydroToyAMReX_EstimateError AT postinitial AFTER HydroToyAMReX_Boundaries
{
  LANG: C
  READS: conserved(everywhere)
  WRITES: AMReX::regrid_error(interior)
} "Estimate local error for regridding initial conditions"

Insertion in HydroToyAMReX/schedule.ccl at line 43 [18.6006]
[17.3141]
[17.3141]
```
# Regridding
```

Insertion in HydroToyAMReX/schedule.ccl at line 45 [18.6006]

[17.3142]

SCHEDULE HydroToyAMReX_Boundaries AT postregrid
{
  LANG: C
  SYNC: conserved
} "Apply boundary conditions after regridding"

Replacement in HydroToyAMReX/schedule.ccl at line 51 [18.6006]

B:BD[17.3143] → [17.3143:3187]

SCHEDULE HydroToyAMReX_Pressure AT poststep

[17.3143]

[17.3187]

SCHEDULE HydroToyAMReX_Pressure AT postregrid AFTER HydroToyAMReX_Boundaries

Replacement in HydroToyAMReX/schedule.ccl at line 58 [18.6006]

B:BD[18.6149] → [17.3286:3357]

SCHEDULE HydroToyAMReX_Fluxes AT poststep AFTER HydroToyAMReX_Pressure

[18.6149]

[17.3357]

SCHEDULE HydroToyAMReX_Fluxes AT postregrid AFTER HydroToyAMReX_Pressure

Replacement in HydroToyAMReX/schedule.ccl at line 67 [18.6006]
B:BD[17.3553] → [4.2194:2236]
```
SCHEDULE HydroToyAMReX_Output AT poststep
```
[17.3553]
[4.2236]
```
# Dependent quantities
SCHEDULE HydroToyAMReX_Pressure AT poststep
```
Replacement in HydroToyAMReX/schedule.ccl at line 75 [18.6006]
B:BD[4.2279] → [4.2279:2300]
```
} "Output grid data"
```
[4.2279]
[4.2300]
```
  WRITES: primitive(everywhere)
} "Calculate pressure"
```

Replacement in HydroToyAMReX/schedule.ccl at line 78 [18.6006]

∅:D[4.2301] → [17.3553:3604]

B:BD[17.3553] → [17.3553:3604]

# SCHEDULE HydroToyAMReX_EstimateError AT poststep

[4.2301]

[17.3604]

SCHEDULE HydroToyAMReX_Fluxes AT poststep AFTER HydroToyAMReX_Pressure
{
  LANG: C
  READS: conserved(everywhere)
  READS: primitive(everywhere)
  WRITES: flux_x(interior) flux_y(interior) flux_z(interior)
  SYNC: flux_x flux_y flux_z
} "Calculate the hydro fluxes"
# SCHEDULE HydroToyAMReX_Output AT poststep

Replacement in HydroToyAMReX/schedule.ccl at line 90 [18.6006]

B:BD[17.3620] → [17.3620:3721]

#   WRITES: AMReX::regrid_error(interior)
# } "Estimate local error for regridding during evolution"

[17.3620]

[17.3721]

#   READS: conserved(everywhere)
# } "Output grid data"
SCHEDULE HydroToyAMReX_EstimateError AT poststep
{
  LANG: C
  READS: conserved(everywhere)
  WRITES: AMReX::regrid_error(interior)
} "Estimate local error for regridding during evolution"

Insertion in HydroToyAMReX/schedule.ccl at line 102 [18.6006]
[17.3723]
[17.3723]
```
# Time stepping
```
Insertion in HydroToyAMReX/schedule.ccl at line 111 [18.6006]
[18.6230]
[4.2302]

Replacement in HydroToyAMReX/schedule.ccl at line 113 [18.6006]

B:BD[4.2303] → [4.2303:3264]

# SCHEDULE HydroToyAMReX_Evolve1 AT evol
# {
#   LANG: C
#   READS: conserved_p(everywhere)
#   READS: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)
#   WRITES: conserved(everywhere)
# } "Evolve the hydro equations, step 1"
# 
# SCHEDULE HydroToyAMReX_Pressure AT evol AFTER HydroToyAMReX_Evolve1
# {
#   LANG: C
#   READS: conserved(everywhere)
#   WRITES: primitive(everywhere)
# } "Calculate pressure"
# 
# SCHEDULE HydroToyAMReX_Fluxes AT evol AFTER HydroToyAMReX_Pressure
# {
#   LANG: C
#   READS: conserved(everywhere)
#   READS: primitive(everywhere)
#   WRITES: flux_x(interior) flux_y(interior) flux_z(interior)
#   SYNC: flux_x flux_y flux_z
# } "Calculate the hydro fluxes"
# 
# SCHEDULE HydroToyAMReX_Evolve AT evol AFTER HydroToyAMReX_Fluxes
# {
#   LANG: C
#   READS: conserved_p(everywhere)
#   READS: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)
#   WRITES: conserved(everywhere)
# } "Evolve the hydro equations, step 2"

[4.2303]


# Restriction
SCHEDULE HydroToyAMReX_Boundaries AT postrestrict
{
  LANG: C
  SYNC: conserved
} "Apply boundary conditions after restricting"

Deletion in HydroToyAMReX/src/hydrotoy.cxx at line 15 [18.6253]
B:BD[18.6479] → [18.6479:6561]
```
////////////////////////////////////////////////////////////////////////////////
```

Insertion in HydroToyAMReX/src/hydrotoy.cxx at line 85 [18.6253]

[17.4652]

}
extern "C" void HydroToyAMReX_Boundaries(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS;
  // do nothing

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 193 [18.6253]
B:BD[17.8519] → [17.8519:8558]
```
    auto mkflux{[&](auto &u, auto f) {
```
[17.8519]
[4.4430]
```
    auto calcflux{[&](auto &u, auto f) {
```

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 205 [18.6253]

B:BD[17.8828] → [17.8828:8904]

    fxrho_(p.I) = mkflux(rho_, [&](auto I) { return rho_(I) * velx_(I); });

[17.8828]

[17.8904]

    fxrho_(p.I) = calcflux(rho_, [&](auto I) { return rho_(I) * velx_(I); });

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 207 [18.6253]

B:BD[17.8905] → [17.8905:9162]

    fxmomx_(p.I) =
        mkflux(momx_, [&](auto I) { return momx_(I) * velx_(I) + press_(I); });
    fxmomy_(p.I) = mkflux(momy_, [&](auto I) { return momy_(I) * velx_(I); });
    fxmomz_(p.I) = mkflux(momz_, [&](auto I) { return momz_(I) * velx_(I); });

[17.8905]

[17.9162]

    fxmomx_(p.I) = calcflux(
        momx_, [&](auto I) { return momx_(I) * velx_(I) + press_(I); });
    fxmomy_(p.I) = calcflux(momy_, [&](auto I) { return momy_(I) * velx_(I); });
    fxmomz_(p.I) = calcflux(momz_, [&](auto I) { return momz_(I) * velx_(I); });

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 212 [18.6253]
B:BD[17.9163] → [17.9163:9190]
```
    fxetot_(p.I) = mkflux(
```
[17.9163]
[17.9190]
```
    fxetot_(p.I) = calcflux(
```
Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 217 [18.6253]
B:BD[17.9338] → [17.9338:9377]
```
    auto mkflux{[&](auto &u, auto f) {
```
[17.9338]
[4.4650]
```
    auto calcflux{[&](auto &u, auto f) {
```

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 229 [18.6253]

B:BD[17.9647] → [17.9647:9723]

    fyrho_(p.I) = mkflux(rho_, [&](auto I) { return rho_(I) * vely_(I); });

[17.9647]

[17.9723]

    fyrho_(p.I) = calcflux(rho_, [&](auto I) { return rho_(I) * vely_(I); });

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 231 [18.6253]

B:BD[17.9724] → [17.9724:9981]

    fymomx_(p.I) = mkflux(momx_, [&](auto I) { return momx_(I) * vely_(I); });
    fymomy_(p.I) =
        mkflux(momy_, [&](auto I) { return momy_(I) * vely_(I) + press_(I); });
    fymomz_(p.I) = mkflux(momz_, [&](auto I) { return momz_(I) * vely_(I); });

[17.9724]

[17.9981]

    fymomx_(p.I) = calcflux(momx_, [&](auto I) { return momx_(I) * vely_(I); });
    fymomy_(p.I) = calcflux(
        momy_, [&](auto I) { return momy_(I) * vely_(I) + press_(I); });
    fymomz_(p.I) = calcflux(momz_, [&](auto I) { return momz_(I) * vely_(I); });

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 236 [18.6253]
B:BD[17.9982] → [17.9982:10009]
```
    fyetot_(p.I) = mkflux(
```
[17.9982]
[17.10009]
```
    fyetot_(p.I) = calcflux(
```
Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 241 [18.6253]
B:BD[17.10157] → [17.10157:10196]
```
    auto mkflux{[&](auto &u, auto f) {
```
[17.10157]
[4.4870]
```
    auto calcflux{[&](auto &u, auto f) {
```

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 253 [18.6253]

B:BD[17.10466] → [17.10466:10542]

    fzrho_(p.I) = mkflux(rho_, [&](auto I) { return rho_(I) * velz_(I); });

[17.10466]

[17.10542]

    fzrho_(p.I) = calcflux(rho_, [&](auto I) { return rho_(I) * velz_(I); });

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 255 [18.6253]

B:BD[17.10543] → [17.10543:10800]

    fzmomx_(p.I) = mkflux(momx_, [&](auto I) { return momx_(I) * velz_(I); });
    fzmomy_(p.I) = mkflux(momy_, [&](auto I) { return momy_(I) * velz_(I); });
    fzmomz_(p.I) =
        mkflux(momz_, [&](auto I) { return momz_(I) * velz_(I) + press_(I); });

[17.10543]

[17.10800]

    fzmomx_(p.I) = calcflux(momx_, [&](auto I) { return momx_(I) * velz_(I); });
    fzmomy_(p.I) = calcflux(momy_, [&](auto I) { return momy_(I) * velz_(I); });
    fzmomz_(p.I) = calcflux(
        momz_, [&](auto I) { return momz_(I) * velz_(I) + press_(I); });

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 260 [18.6253]
B:BD[17.10801] → [17.10801:10828]
```
    fzetot_(p.I) = mkflux(
```
[17.10801]
[17.10828]
```
    fzetot_(p.I) = calcflux(
```
Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 265 [18.6253]
B:BD[4.5091] → [4.5091:5147]
```
extern "C" void HydroToyAMReX_Evolve1(CCTK_ARGUMENTS) {
```
[4.5091]
[4.5147]
```
extern "C" void HydroToyAMReX_Evolve(CCTK_ARGUMENTS) {
```

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 310 [18.6253]

B:BD[4.7259] → [4.7259:7490]

    auto mkupdate{[&](auto &fx, auto &fy, auto &fz) {
      return 0.5 * dt *
             ((fx(p.I + p.DI(0)) - fx(p.I)) / dx +
              (fy(p.I + p.DI(1)) - fy(p.I)) / dy +
              (fz(p.I + p.DI(2)) - fz(p.I)) / dz);

[4.7259]

[4.7490]

    auto calcupdate{[&](auto &fx, auto &fy, auto &fz) {
      return dt * ((fx(p.I + p.DI(0)) - fx(p.I)) / dx +
                   (fy(p.I + p.DI(1)) - fy(p.I)) / dy +
                   (fz(p.I + p.DI(2)) - fz(p.I)) / dz);

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 316 [18.6253]

B:BD[4.7499] → [4.7499:7563]

    rho_(p.I) = rho_p_(p.I) - mkupdate(fxrho_, fyrho_, fzrho_);

[4.7499]

[4.7563]

    rho_(p.I) = rho_p_(p.I) - calcupdate(fxrho_, fyrho_, fzrho_);

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 318 [18.6253]

B:BD[4.7564] → [4.7564:7771]

    momx_(p.I) = momx_p_(p.I) - mkupdate(fxmomx_, fymomx_, fzmomx_);
    momy_(p.I) = momy_p_(p.I) - mkupdate(fxmomy_, fymomy_, fzmomy_);
    momz_(p.I) = momz_p_(p.I) - mkupdate(fxmomz_, fymomz_, fzmomz_);

[4.7564]

[4.7771]

    momx_(p.I) = momx_p_(p.I) - calcupdate(fxmomx_, fymomx_, fzmomx_);
    momy_(p.I) = momy_p_(p.I) - calcupdate(fxmomy_, fymomy_, fzmomy_);
    momz_(p.I) = momz_p_(p.I) - calcupdate(fxmomz_, fymomz_, fzmomz_);

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 322 [18.6253]

B:BD[4.7772] → [4.7772:7841]

    etot_(p.I) = etot_p_(p.I) - mkupdate(fxetot_, fyetot_, fzetot_);

[4.7772]

[4.7841]

    etot_(p.I) = etot_p_(p.I) - calcupdate(fxetot_, fyetot_, fzetot_);

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 326 [18.6253]
∅:D[4.7850] → [18.7058:7113]
B:BD[18.7058] → [18.7058:7113]
```
extern "C" void HydroToyAMReX_Evolve(CCTK_ARGUMENTS) {
```
[4.7850]
[18.7113]
```
extern "C" void HydroToyAMReX_EstimateError(CCTK_ARGUMENTS) {
```

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 330 [18.6253]

B:BD[18.7167] → [18.7167:7207]

B:BD[18.7207] → [19.1816:1948]

  const CCTK_REAL dt = CCTK_DELTA_TIME;
  const CCTK_REAL dx = CCTK_DELTA_SPACE(0);
  const CCTK_REAL dy = CCTK_DELTA_SPACE(1);
  const CCTK_REAL dz = CCTK_DELTA_SPACE(2);

[18.7167]

[18.7339]

  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> rho_(cctkGH, rho);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> momx_(cctkGH, momx);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> momy_(cctkGH, momy);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> momz_(cctkGH, momz);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> etot_(cctkGH, etot);

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 336 [18.6253]

B:BD[18.7340] → [17.10904:11950]

  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> rho_p_(cctkGH, rho_p);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> momx_p_(cctkGH, momx_p);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> momy_p_(cctkGH, momy_p);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> momz_p_(cctkGH, momz_p);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> etot_p_(cctkGH, etot_p);
  const Loop::GF3D<const CCTK_REAL, 0, 1, 1> fxrho_(cctkGH, fxrho);
  const Loop::GF3D<const CCTK_REAL, 0, 1, 1> fxmomx_(cctkGH, fxmomx);
  const Loop::GF3D<const CCTK_REAL, 0, 1, 1> fxmomy_(cctkGH, fxmomy);
  const Loop::GF3D<const CCTK_REAL, 0, 1, 1> fxmomz_(cctkGH, fxmomz);
  const Loop::GF3D<const CCTK_REAL, 0, 1, 1> fxetot_(cctkGH, fxetot);
  const Loop::GF3D<const CCTK_REAL, 1, 0, 1> fyrho_(cctkGH, fyrho);
  const Loop::GF3D<const CCTK_REAL, 1, 0, 1> fymomx_(cctkGH, fymomx);
  const Loop::GF3D<const CCTK_REAL, 1, 0, 1> fymomy_(cctkGH, fymomy);
  const Loop::GF3D<const CCTK_REAL, 1, 0, 1> fymomz_(cctkGH, fymomz);
  const Loop::GF3D<const CCTK_REAL, 1, 0, 1> fyetot_(cctkGH, fyetot);

[18.7340]

[17.11950]

  const Loop::GF3D<CCTK_REAL, 1, 1, 1> regrid_error_(cctkGH, regrid_error);

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 338 [18.6253]

B:BD[17.11951] → [17.11951:13011]

  const Loop::GF3D<const CCTK_REAL, 1, 1, 0> fzrho_(cctkGH, fzrho);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 0> fzmomx_(cctkGH, fzmomx);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 0> fzmomy_(cctkGH, fzmomy);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 0> fzmomz_(cctkGH, fzmomz);
  const Loop::GF3D<const CCTK_REAL, 1, 1, 0> fzetot_(cctkGH, fzetot);
  const Loop::GF3D<CCTK_REAL, 1, 1, 1> rho_(cctkGH, rho);
  const Loop::GF3D<CCTK_REAL, 1, 1, 1> momx_(cctkGH, momx);
  const Loop::GF3D<CCTK_REAL, 1, 1, 1> momy_(cctkGH, momy);
  const Loop::GF3D<CCTK_REAL, 1, 1, 1> momz_(cctkGH, momz);
  const Loop::GF3D<CCTK_REAL, 1, 1, 1> etot_(cctkGH, etot);
  // Transport
  // dt rho + d_i (rho vel^i) = 0
  // dt mom_j + d_i (mom_j vel^i) = 0
  // dt etot + d_i (etot vel^i) = 0
  Loop::loop_all<1, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    auto mkupdate{[&](auto &fx, auto &fy, auto &fz) {
      return dt * ((fx(p.I + p.DI(0)) - fx(p.I)) / dx +
                   (fy(p.I + p.DI(1)) - fy(p.I)) / dy +
                   (fz(p.I + p.DI(2)) - fz(p.I)) / dz);

[17.11951]

[17.13011]

  Loop::loop_int<1, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    auto calcerr{[&](auto &var_) {
      CCTK_REAL err{0};
      for (int d = 0; d < dim; ++d) {
        auto varm = var_(p.I - p.DI(d));
        auto var0 = var_(p.I);
        auto varp = var_(p.I + p.DI(d));
        err = fmax(err, fabs(varm - 2 * var0 + varp));
      }
      return err;

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 350 [18.6253]

B:BD[18.7435] → [17.13020:13084]

∅:D[19.2168] → [18.7754:7755]

∅:D[20.17559] → [18.7754:7755]

∅:D[17.13084] → [18.7754:7755]

B:BD[18.7754] → [18.7754:7755]

B:BD[18.7755] → [17.13085:13292]

∅:D[17.13292] → [18.7925:7926]

B:BD[18.7925] → [18.7925:7926]

B:BD[18.7926] → [17.13293:13362]

    rho_(p.I) = rho_p_(p.I) - mkupdate(fxrho_, fyrho_, fzrho_);
    momx_(p.I) = momx_p_(p.I) - mkupdate(fxmomx_, fymomx_, fzmomx_);
    momy_(p.I) = momy_p_(p.I) - mkupdate(fxmomy_, fymomy_, fzmomy_);
    momz_(p.I) = momz_p_(p.I) - mkupdate(fxmomz_, fymomz_, fzmomz_);
    etot_(p.I) = etot_p_(p.I) - mkupdate(fxetot_, fyetot_, fzetot_);

[18.7435]

[18.8055]

    regrid_error_(p.I) =
        fmax(fmax(fmax(fmax(calcerr(rho_), calcerr(momx_)), calcerr(momy_)),
                  calcerr(momz_)),
             calcerr(etot_));

Deletion in HydroToyAMReX/src/hydrotoy.cxx at line 357 [18.6253]

B:BD[18.8064] → [17.13363:13583]

B:BD[17.13583] → [4.7851:7854]

extern "C" void HydroToyAMReX_EstimateError(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS;
  DECLARE_CCTK_PARAMETERS;
  Loop::loop_int<1, 1, 1>(
      cctkGH, [&](const Loop::PointDesc &p) { regrid_error[p.idx] = 0.0; });
}

Make HydroToy work with multiple levels

Dependencies

In channels

Change contents

Replacement in AMReX/src/driver.cxx at line 417 [6.15972]

Replacement in AMReX/src/driver.cxx at line 448 [6.15972]

Replacement in AMReX/src/driver.cxx at line 453 [6.15972]

Replacement in AMReX/src/schedule.cxx at line 792 [8.531]

Replacement in AMReX/src/schedule.cxx at line 794 [8.531]

Insertion in AMReX/src/schedule.cxx at line 875 [8.531]

Insertion in AMReX/src/schedule.cxx at line 899 [8.531]

Insertion in AMReX/src/schedule.cxx at line 958 [8.531]

Replacement in AMReX/src/schedule.cxx at line 975 [8.531]

Replacement in AMReX/src/schedule.cxx at line 977 [8.531]

Replacement in AMReX/src/schedule.cxx at line 1167 [8.531]

Replacement in AMReX/src/schedule.cxx at line 1287 [8.531]

Deletion in AMReX/src/schedule.cxx at line 1291 [8.531]

Deletion in AMReX/src/schedule.cxx at line 1304 [8.531]

Insertion in AMReX/src/schedule.cxx at line 1305 [8.531]

Insertion in AMReX/src/schedule.cxx at line 1307 [8.531]

Replacement in AMReX/src/schedule.cxx at line 1332 [8.531]

Replacement in HydroToyAMReX/par/hydro-shocktube.par at line 20 [4.22]

Replacement in HydroToyAMReX/par/hydro-sphericalshock.par at line 10 [4.1031]

Replacement in HydroToyAMReX/par/hydro-sphericalshock.par at line 36 [4.1031]

Insertion in HydroToyAMReX/schedule.ccl at line 6 [18.6006]

Insertion in HydroToyAMReX/schedule.ccl at line 9 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 18 [18.6006]

Insertion in HydroToyAMReX/schedule.ccl at line 43 [18.6006]

Insertion in HydroToyAMReX/schedule.ccl at line 45 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 51 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 58 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 67 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 75 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 78 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 90 [18.6006]

Insertion in HydroToyAMReX/schedule.ccl at line 102 [18.6006]

Insertion in HydroToyAMReX/schedule.ccl at line 111 [18.6006]

Replacement in HydroToyAMReX/schedule.ccl at line 113 [18.6006]

Deletion in HydroToyAMReX/src/hydrotoy.cxx at line 15 [18.6253]

Insertion in HydroToyAMReX/src/hydrotoy.cxx at line 85 [18.6253]

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 193 [18.6253]

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 205 [18.6253]

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 207 [18.6253]

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 212 [18.6253]

Replacement in HydroToyAMReX/src/hydrotoy.cxx at line 217 [18.6253]