eschnett/cactusamrex - Change W6W3EC56H4PPREAK7DMX3DAHYKDAHSASGDP63DSU55TRCBBGJXBAC

Hydro: Finish implementation

Created by Erik Schnetter on November 12, 2019

W6W3EC56H4PPREAK7DMX3DAHYKDAHSASGDP63DSU55TRCBBGJXBAC

Dependencies

[2] Y7FSCHC3E3KIGVSB2KQRIPKBD2HL6JPUNBBSVZQAJKPXJKKRVNQQC

In channels

master

main

Change contents

Replacement in Hydro/interface.ccl at line 20 [2.407]

B:BD[2.694] → [2.694:756]

CCTK_REAL conserved TYPE=gf TIMELEVELS=2 TAGS='index={1 1 1}'

[2.694]

[2.756]

CCTK_REAL primitive TYPE=gf TAGS='index={1 1 1} checkpoint="no" restrict="no"'

Deletion in Hydro/interface.ccl at line 23 [2.407]

B:BD[2.764] → [2.764:1123]

  momx momy momz
  etot
} "Conserved variables"
CCTK_REAL conserved_rhs TYPE=gf TAGS='index={1 1 1} fluxes="HydroToyCarpetX::flux_x HydroToyCarpetX::flux_y HydroToyCarpetX::flux_z"'
{
  dtrho
  dtmomx dtmomy dtmomz
  dtetot
} "Time derivatives of conserved variables"
CCTK_REAL primitive TYPE=gf TAGS='index={1 1 1} checkpoint="no" restrict="no"'
{
  press

Insertion in Hydro/interface.ccl at line 25 [2.407]
[2.1147]
[2.1147]
```
  press
```
Replacement in Hydro/interface.ccl at line 30 [2.407]
B:BD[2.1270] → [2.1270:1278]
```
  fxrho
```
[2.1270]
[2.1278]
```
  fxdens
```
Replacement in Hydro/interface.ccl at line 37 [2.407]
B:BD[2.1435] → [2.1435:1443]
```
  fyrho
```
[2.1435]
[2.1443]
```
  fydens
```
Replacement in Hydro/interface.ccl at line 44 [2.407]
B:BD[2.1600] → [2.1600:1608]
```
  fzrho
```
[2.1600]
[2.1608]
```
  fzdens
```

Insertion in Hydro/interface.ccl at line 48 [2.407]

[2.1666]


CCTK_REAL conserved_rhs TYPE=gf TAGS='index={1 1 1} checkpoint="no" fluxes="Hydro::flux_x Hydro::flux_y Hydro::flux_z"'
{
  dtdens
  dtmomx dtmomy dtmomz
  dtetot
} "Time derivatives of conserved variables"
CCTK_REAL conserved TYPE=gf TAGS='index={1 1 1} rhs="Hydro::conserved_rhs"'
{
  dens
  momx momy momz
  etot
} "Conserved variables"

File addition: par (drwxr-xr-x)
[2.8]

File addition: hydro-sphericalshock.par (-rw-r--r--)

[0.468]

ActiveThorns = "
    CarpetX
    Formaline
    IOUtil
    Hydro
    ODESolvers
    SystemTopology
    TimerReport
"
 
$nlevels = 3
$ncells = 32
Cactus::cctk_show_schedule = no
# Cactus::terminate = "iteration"
# Cactus::cctk_itlast = 0
Cactus::terminate = "time"
Cactus::cctk_final_time = 1.0
CarpetX::verbose = no
CarpetX::xmin = -1.0
CarpetX::ymin = -1.0
CarpetX::zmin = -1.0
CarpetX::xmax = +1.0
CarpetX::ymax = +1.0
CarpetX::zmax = +1.0
CarpetX::ncells_x = $ncells
CarpetX::ncells_y = $ncells
CarpetX::ncells_z = $ncells
CarpetX::max_num_levels = $nlevels
CarpetX::regrid_every = 16
CarpetX::regrid_error_threshold = 0.01
CarpetX::prolongation_type = "conservative"
CarpetX::do_reflux = no
CarpetX::do_restrict = no
CarpetX::restrict_during_sync = yes
CarpetX::dtfac = 0.25
Hydro::setup = "spherical shock"
Hydro::output_every = 1
IO::out_dir = $parfile
IO::out_every = $ncells * 2 ** ($nlevels - 1) / 32
CarpetX::out_tsv = no
TimerReport::out_every = $ncells * 2 ** ($nlevels - 1) / 32
TimerReport::out_filename = "TimerReport"
TimerReport::output_all_timers_together = yes
TimerReport::output_all_timers_readable = yes
TimerReport::n_top_timers = 50

Deletion in Hydro/schedule.ccl at line 3 [2.2254]
B:BD[2.2295] → [2.2295:2340]
```
STORAGE: conserved[2]
STORAGE: conserved_rhs
```
Insertion in Hydro/schedule.ccl at line 5 [2.2254]
[2.2389]
[2.2389]
```
STORAGE: conserved_rhs
STORAGE: conserved
```

Replacement in Hydro/schedule.ccl at line 15 [2.2254]

B:BD[2.2463] → [2.2463:2530]

  WRITES: conserved(interior)
  INVALIDATES: primitive(everywhere)

[2.2463]

[2.2530]

  WRITES: primitive(interior)
  INVALIDATES: conserved(everywhere)

Replacement in Hydro/schedule.ccl at line 18 [2.2254]
B:BD[2.2602] → [2.2602:2620]
```
  SYNC: conserved
```
[2.2602]
[2.2620]
```
  INVALIDATES: conserved_rhs(everywhere)
  SYNC: primitive
```

Replacement in Hydro/schedule.ccl at line 22 [2.2254]

B:BD[2.2657] → [2.2657:2701]

SCHEDULE Hydro_EstimateError AT postinitial

[2.2657]

[2.2701]

SCHEDULE Hydro_Prim2Con AT postinitial
{
  LANG: C
  READS: primitive(everywhere)
  WRITES: conserved(everywhere)
  INVALIDATES: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)
  INVALIDATES: conserved_rhs(everywhere)
} "Calculate conserved variables"
SCHEDULE Hydro_EstimateError AT postinitial AFTER Hydro_Prim2Con

Insertion in Hydro/schedule.ccl at line 46 [2.2254]
[2.2932]
[2.2932]
```
  INVALIDATES: primitive(everywhere)
```

Replacement in Hydro/schedule.ccl at line 65 [2.2254]

B:BD[2.3369] → [2.3369:3396]

    } "Calculate pressure"

[2.3369]

[2.3396]

      INVALIDATES: conserved_rhs(everywhere)
    } "Calculate primitive variables"

Insertion in Hydro/schedule.ccl at line 74 [2.2254]
[2.3620]
[2.3620]
```
      INVALIDATES: conserved_rhs(everywhere)
```

Insertion in Hydro/schedule.ccl at line 76 [2.2254]

[2.3655]


    SCHEDULE Hydro_RHS IN Hydro_OutputGroup AFTER Hydro_Fluxes
    {
      LANG: C
      READS: flux_x(interior) flux_y(interior) flux_z(interior)
      WRITES: conserved_rhs(interior)
      SYNC: conserved_rhs           # sync and restrict
    } "Calculate RHS of the hydro equations"

Replacement in Hydro/schedule.ccl at line 97 [2.2254]
B:BD[2.3861] → [2.3861:3893]
```
SCHEDULE Hydro_Con2prim AT evol
```
[2.3861]
[2.3893]
```
SCHEDULE Hydro_Con2prim IN ODESolvers_RHS
```
Insertion in Hydro/schedule.ccl at line 103 [2.2254]
[2.4040]
[2.4040]
```
  INVALIDATES: conserved_rhs(everywhere)
```

Replacement in Hydro/schedule.ccl at line 106 [2.2254]

B:BD[2.4064] → [2.4064:4115]

SCHEDULE Hydro_Fluxes AT evol AFTER Hydro_Con2prim

[2.4064]

[2.4115]

SCHEDULE Hydro_Fluxes IN ODESolvers_RHS AFTER Hydro_Con2prim

Insertion in Hydro/schedule.ccl at line 112 [2.2254]
[2.4250]
[2.4250]
```
  INVALIDATES: conserved_rhs(everywhere)
```

Replacement in Hydro/schedule.ccl at line 116 [2.2254]

B:BD[2.4325] → [2.4325:4371]

SCHEDULE Hydro_RHS AT evol AFTER Hydro_Fluxes

[2.4325]

[2.4371]

SCHEDULE Hydro_RHS IN ODESolvers_RHS AFTER Hydro_Fluxes

Replacement in Hydro/src/con2prim.cxx at line 36 [2.4919]

B:BD[2.5859] → [2.5859:5911]

      for (int i = imin[0]; i < imax[0]; i += VS) {

[2.5859]

[2.5911]

      for (int i = imin[0]; i < imax[0]; i += vsize) {

Replacement in Hydro/src/con2prim.cxx at line 39 [2.4919]

B:BD[2.5957] → [2.5957:6003]

        CCTK_REALVEC rho1 = vloadu(rho[ind]);

[2.5957]

[2.6003]

        CCTK_REALVEC dens1 = vloadu(dens[ind]);

Replacement in Hydro/src/con2prim.cxx at line 45 [2.4919]

B:BD[2.6196] → [2.6196:6250]

        CCTK_REALVEC rho_inv = CCTK_REAL(1.0) / rho1;

[2.6196]

[2.6250]

        CCTK_REALVEC rho1 = dens1;
        CCTK_REALVEC rho1_inv = CCTK_REAL(1.0) / rho1;

Replacement in Hydro/src/con2prim.cxx at line 49 [2.4919]

B:BD[2.6251] → [2.6251:6424]

        CCTK_REALVEC ekin = CCTK_REAL(0.5) * rho_inv *
                            sqrt(pow2(momx1) + pow2(momy1) + pow2(momz1));
        CCTK_REALVEC eint1 = etot1 - ekin;

[2.6251]

[2.6424]

        CCTK_REALVEC ekin1 = CCTK_REAL(0.5) * rho1_inv *
                             sqrt(pow2(momx1) + pow2(momy1) + pow2(momz1));
        CCTK_REALVEC eint1 = etot1 - ekin1;

Replacement in Hydro/src/con2prim.cxx at line 57 [2.4919]

B:BD[2.6575] → [2.6575:6713]

        CCTK_REALVEC velx1 = rho_inv * momx1;
        CCTK_REALVEC vely1 = rho_inv * momy1;
        CCTK_REALVEC velz1 = rho_inv * momz1;

[2.6575]

[2.6713]

        CCTK_REALVEC velx1 = rho1_inv * momx1;
        CCTK_REALVEC vely1 = rho1_inv * momy1;
        CCTK_REALVEC velz1 = rho1_inv * momz1;

Replacement in Hydro/src/con2prim.cxx at line 61 [2.4919]

B:BD[2.6714] → [2.6714:6778]

        press1.storeu_partial(press[ind], i, imin[0], imax[0]);

[2.6714]

[2.6778]

        rho1.storeu_partial(rho[ind], i, imin[0], imax[0]);

Insertion in Hydro/src/con2prim.cxx at line 66 [2.4919]

[2.7026]

        press1.storeu_partial(press[ind], i, imin[0], imax[0]);

Replacement in Hydro/src/defs.hxx at line 14 [2.7080]

B:BD[2.7235] → [2.7235:7287]

template <typename T> T pow2(T x) { return x * x; }

[2.7235]

[2.7287]

template <typename T> inline CCTK_ATTRIBUTE_ALWAYS_INLINE T pow2(T x) {
  return x * x;
}

Replacement in Hydro/src/defs.hxx at line 24 [2.7080]

B:BD[2.7413] → [2.7413:7482]

template <typename T> inline T fmax5(T x0, T x1, T x2, T x3, T x4) {

[2.7413]

[2.7482]

template <typename T>
inline CCTK_ATTRIBUTE_ALWAYS_INLINE T fmax5(T x0, T x1, T x2, T x3, T x4) {

Replacement in Hydro/src/defs.hxx at line 34 [2.7080]
B:BD[2.7647] → [2.7647:7695]
```
constexpr int VS = vecprops<CCTK_REAL>::size();
```
[2.7647]
[2.7695]
```
constexpr int vsize = vecprops<CCTK_REAL>::size();
```

Replacement in Hydro/src/defs.hxx at line 42 [2.7080]

B:BD[2.7878] → [2.7878:7939]

  array<CCTK_REAL, VS> iota_;
  for (int i = 0; i < VS; ++i)

[2.7878]

[2.7939]

  array<CCTK_REAL, vsize> iota_;
  for (int i = 0; i < vsize; ++i)

Replacement in Hydro/src/estimateerror.cxx at line 39 [2.8162]

B:BD[2.9327] → [2.9327:9379]

      for (int i = imin[0]; i < imax[0]; i += VS) {

[2.9327]

[2.9379]

      for (int i = imin[0]; i < imax[0]; i += vsize) {

Replacement in Hydro/src/estimateerror.cxx at line 56 [2.8162]

B:BD[2.10104] → [2.10104:10181]

            fmax5(calcerr(rho), calcerr(momx), calcerr(momy), calcerr(momz),

[2.10104]

[2.10181]

            fmax5(calcerr(dens), calcerr(momx), calcerr(momy), calcerr(momz),

Replacement in Hydro/src/fluxes.cxx at line 63 [2.10349]

B:BD[2.12284] → [2.12284:12345]

                          CCTK_REAL *restrict const fdirrho,

[2.12284]

[2.12345]

                          CCTK_REAL *restrict const fdirdens,

Insertion in Hydro/src/fluxes.cxx at line 78 [2.10349]

[2.13083]

    const ptrdiff_t offdir = nghostzones[0] * didir + nghostzones[1] * djdir +
                             nghostzones[2] * dkdir;

Replacement in Hydro/src/fluxes.cxx at line 83 [2.10349]

B:BD[2.13190] → [2.13190:13250]

        for (int i = imindir[0]; i < imaxdir[0]; i += VS) {

[2.13190]

[2.13250]

        for (int i = imindir[0]; i < imaxdir[0]; i += vsize) {

Replacement in Hydro/src/fluxes.cxx at line 85 [2.10349]

B:BD[2.13297] → [2.13297:13353]

          ptrdiff_t inddir = i + djdir * j + dkdir * k;

[2.13297]

[2.13353]

          ptrdiff_t inddir = i + djdir * j + dkdir * k - offdir;

Replacement in Hydro/src/fluxes.cxx at line 88 [2.10349]

B:BD[2.13406] → [2.13406:13481]

          array<CCTK_REALVEC, 2> rho2, momx2, momy2, momz2, etot2, press2,

[2.13406]

[2.13481]

          array<CCTK_REALVEC, 2> dens2, momx2, momy2, momz2, etot2, press2,

Replacement in Hydro/src/fluxes.cxx at line 91 [2.10349]

B:BD[2.13544] → [2.13544:13601]

            rho2[n] = vloadu(rho[ind + (n - 1) * ddir]);

[2.13544]

[2.13601]

            dens2[n] = vloadu(dens[ind + (n - 1) * ddir]);

Insertion in Hydro/src/fluxes.cxx at line 100 [2.10349]
[2.13974]
[2.13974]
```
#warning "TODO: This calculates the fluxes twice; change this"
```

Replacement in Hydro/src/fluxes.cxx at line 102 [2.10349]

B:BD[2.14012] → [2.14012:14088]

          array<CCTK_REALVEC, 2> fdirrho2, fdirmomx2, fdirmomy2, fdirmomz2,

[2.14012]

[2.14088]

          array<CCTK_REALVEC, 2> fdirdens2, fdirmomx2, fdirmomy2, fdirmomz2,

Replacement in Hydro/src/fluxes.cxx at line 105 [2.10349]

B:BD[2.14153] → [2.14153:14201]

            fdirrho2[n] = rho2[n] * veldir2[n];

[2.14153]

[2.14201]

            fdirdens2[n] = dens2[n] * veldir2[n];

Replacement in Hydro/src/fluxes.cxx at line 113 [2.10349]

B:BD[2.14550] → [2.14550:14613]

          CCTK_REALVEC fdirrho1 = dAdir * llf(rho2, fdirrho2);

[2.14550]

[2.14613]

          CCTK_REALVEC fdirdens1 = dAdir * llf(dens2, fdirdens2);

Replacement in Hydro/src/fluxes.cxx at line 120 [2.10349]

B:BD[2.14904] → [2.14904:14983]

          fdirrho1.storeu_partial(fdirrho[inddir], i, imindir[0], imaxdir[0]);

[2.14904]

[2.14983]

          fdirdens1.storeu_partial(fdirdens[inddir], i, imindir[0], imaxdir[0]);

Replacement in Hydro/src/fluxes.cxx at line 130 [2.10349]

B:BD[2.15338] → [2.15338:15557]

  calcflux(0, di, dAx, nx, velx, fxrho, fxmomx, fxmomy, fxmomz, fxetot);
  calcflux(1, dj, dAy, ny, vely, fyrho, fymomx, fymomy, fymomz, fyetot);
  calcflux(2, dk, dAz, nz, velz, fzrho, fzmomx, fzmomy, fzmomz, fzetot);

[2.15338]

[2.15557]

  calcflux(0, di, dAx, nx, velx, fxdens, fxmomx, fxmomy, fxmomz, fxetot);
  calcflux(1, dj, dAy, ny, vely, fydens, fymomx, fymomy, fymomz, fyetot);
  calcflux(2, dk, dAz, nz, velz, fzdens, fzmomx, fzmomy, fzmomz, fzetot);

Insertion in Hydro/src/initialize.cxx at line 26 [2.15599]

[2.16231]

  const Loop::vect<int, dim> lbnd{{cctk_lbnd[0], cctk_lbnd[1], cctk_lbnd[2]}};

Replacement in Hydro/src/initialize.cxx at line 48 [2.15599]

B:BD[2.17049] → [2.17049:17103]

        for (int i = imin[0]; i < imax[0]; i += VS) {

[2.17049]

[2.17103]

        for (int i = imin[0]; i < imax[0]; i += vsize) {

Replacement in Hydro/src/initialize.cxx at line 51 [2.15599]

B:BD[2.17196] → [2.17196:17384]

          CCTK_REALVEC momx1 = CCTK_REAL(0.0);
          CCTK_REALVEC momy1 = CCTK_REAL(0.0);
          CCTK_REALVEC momz1 = CCTK_REAL(0.0);
          CCTK_REALVEC etot1 = CCTK_REAL(1.0);

[2.17196]

[2.17384]

          CCTK_REALVEC velx1 = CCTK_REAL(0.0);
          CCTK_REALVEC vely1 = CCTK_REAL(0.0);
          CCTK_REALVEC velz1 = CCTK_REAL(0.0);
          CCTK_REALVEC eint1 = CCTK_REAL(1.0);

Replacement in Hydro/src/initialize.cxx at line 56 [2.15599]

B:BD[2.17446] → [2.17446:17702]

          momx1.storeu_partial(momx[ind], i, imin[0], imax[0]);
          momy1.storeu_partial(momy[ind], i, imin[0], imax[0]);
          momz1.storeu_partial(momz[ind], i, imin[0], imax[0]);
          etot1.storeu_partial(etot[ind], i, imin[0], imax[0]);

[2.17446]

[2.17702]

          velx1.storeu_partial(velx[ind], i, imin[0], imax[0]);
          vely1.storeu_partial(vely[ind], i, imin[0], imax[0]);
          velz1.storeu_partial(velz[ind], i, imin[0], imax[0]);
          eint1.storeu_partial(eint[ind], i, imin[0], imax[0]);

Replacement in Hydro/src/initialize.cxx at line 68 [2.15599]

B:BD[2.17871] → [2.17871:17925]

        for (int i = imin[0]; i < imax[0]; i += VS) {

[2.17871]

[2.17925]

        for (int i = imin[0]; i < imax[0]; i += vsize) {

Replacement in Hydro/src/initialize.cxx at line 70 [2.15599]

B:BD[2.17972] → [2.17972:18048]

          CCTK_REALVEC x1 = x0[0] + (i + CCTK_REAL(0.5) + viota()) * dx[0];

[2.17972]

[2.18048]

          CCTK_REALVEC x1 = x0[0] + (lbnd[0] + i) * dx[0] + viota() * dx[0];

Replacement in Hydro/src/initialize.cxx at line 72 [2.15599]

B:BD[2.18094] → [2.18094:18165]

          CCTK_REALVEC momx1 = CCTK_REAL(0.0) + CCTK_REAL(amplitude) *

[2.18094]

[2.18165]

          CCTK_REALVEC velx1 = CCTK_REAL(0.0) + CCTK_REAL(amplitude) *

Replacement in Hydro/src/initialize.cxx at line 74 [2.15599]

B:BD[2.18245] → [2.18245:18420]

          CCTK_REALVEC momy1 = CCTK_REAL(0.0);
          CCTK_REALVEC momz1 = CCTK_REAL(0.0);
          CCTK_REALVEC etot1 = CCTK_REAL(1.0); // should add kinetic energy here

[2.18245]

[2.18420]

          CCTK_REALVEC vely1 = CCTK_REAL(0.0);
          CCTK_REALVEC velz1 = CCTK_REAL(0.0);
          CCTK_REALVEC eint1 = CCTK_REAL(1.0);

Replacement in Hydro/src/initialize.cxx at line 78 [2.15599]

B:BD[2.18482] → [2.18482:18738]

          momx1.storeu_partial(momx[ind], i, imin[0], imax[0]);
          momy1.storeu_partial(momy[ind], i, imin[0], imax[0]);
          momz1.storeu_partial(momz[ind], i, imin[0], imax[0]);
          etot1.storeu_partial(etot[ind], i, imin[0], imax[0]);

[2.18482]

[2.18738]

          velx1.storeu_partial(velx[ind], i, imin[0], imax[0]);
          vely1.storeu_partial(vely[ind], i, imin[0], imax[0]);
          velz1.storeu_partial(velz[ind], i, imin[0], imax[0]);
          eint1.storeu_partial(eint[ind], i, imin[0], imax[0]);

Replacement in Hydro/src/initialize.cxx at line 90 [2.15599]

B:BD[2.18907] → [2.18907:18961]

        for (int i = imin[0]; i < imax[0]; i += VS) {

[2.18907]

[2.18961]

        for (int i = imin[0]; i < imax[0]; i += vsize) {

Replacement in Hydro/src/initialize.cxx at line 92 [2.15599]

B:BD[2.19008] → [2.19008:19084]

          CCTK_REALVEC x1 = x0[0] + (i + CCTK_REAL(0.5) + viota()) * dx[0];

[2.19008]

[2.19084]

          CCTK_REALVEC x1 = x0[0] + (lbnd[0] + i) * dx[0] + viota() * dx[0];

Replacement in Hydro/src/initialize.cxx at line 95 [2.15599]

B:BD[2.19155] → [2.19155:19347]

          CCTK_REALVEC momx1l = CCTK_REAL(0.0);
          CCTK_REALVEC momy1l = CCTK_REAL(0.0);
          CCTK_REALVEC momz1l = CCTK_REAL(0.0);
          CCTK_REALVEC etot1l = CCTK_REAL(2.0);

[2.19155]

[2.19347]

          CCTK_REALVEC velx1l = CCTK_REAL(0.0);
          CCTK_REALVEC vely1l = CCTK_REAL(0.0);
          CCTK_REALVEC velz1l = CCTK_REAL(0.0);
          CCTK_REALVEC eint1l = CCTK_REAL(2.0);

Replacement in Hydro/src/initialize.cxx at line 101 [2.15599]

B:BD[2.19419] → [2.19419:19611]

          CCTK_REALVEC momx1r = CCTK_REAL(0.0);
          CCTK_REALVEC momy1r = CCTK_REAL(0.0);
          CCTK_REALVEC momz1r = CCTK_REAL(0.0);
          CCTK_REALVEC etot1r = CCTK_REAL(1.0);

[2.19419]

[2.19611]

          CCTK_REALVEC velx1r = CCTK_REAL(0.0);
          CCTK_REALVEC vely1r = CCTK_REAL(0.0);
          CCTK_REALVEC velz1r = CCTK_REAL(0.0);
          CCTK_REALVEC eint1r = CCTK_REAL(1.0);

Replacement in Hydro/src/initialize.cxx at line 107 [2.15599]

B:BD[2.19705] → [2.19705:20013]

          CCTK_REALVEC momx1 = ifthen(x1 <= CCTK_REAL(0.0), momx1l, momx1r);
          CCTK_REALVEC momy1 = ifthen(x1 <= CCTK_REAL(0.0), momy1l, momy1r);
          CCTK_REALVEC momz1 = ifthen(x1 <= CCTK_REAL(0.0), momz1l, momz1r);
          CCTK_REALVEC etot1 = ifthen(x1 <= CCTK_REAL(0.0), etot1l, etot1r);

[2.19705]

[2.20013]

          CCTK_REALVEC velx1 = ifthen(x1 <= CCTK_REAL(0.0), velx1l, velx1r);
          CCTK_REALVEC vely1 = ifthen(x1 <= CCTK_REAL(0.0), vely1l, vely1r);
          CCTK_REALVEC velz1 = ifthen(x1 <= CCTK_REAL(0.0), velz1l, velz1r);
          CCTK_REALVEC eint1 = ifthen(x1 <= CCTK_REAL(0.0), eint1l, eint1r);

Replacement in Hydro/src/initialize.cxx at line 112 [2.15599]

B:BD[2.20075] → [2.20075:20331]

          momx1.storeu_partial(momx[ind], i, imin[0], imax[0]);
          momy1.storeu_partial(momy[ind], i, imin[0], imax[0]);
          momz1.storeu_partial(momz[ind], i, imin[0], imax[0]);
          etot1.storeu_partial(etot[ind], i, imin[0], imax[0]);

[2.20075]

[2.20331]

          velx1.storeu_partial(velx[ind], i, imin[0], imax[0]);
          vely1.storeu_partial(vely[ind], i, imin[0], imax[0]);
          velz1.storeu_partial(velz[ind], i, imin[0], imax[0]);
          eint1.storeu_partial(eint[ind], i, imin[0], imax[0]);

Replacement in Hydro/src/initialize.cxx at line 124 [2.15599]

B:BD[2.20505] → [2.20505:20559]

        for (int i = imin[0]; i < imax[0]; i += VS) {

[2.20505]

[2.20559]

        for (int i = imin[0]; i < imax[0]; i += vsize) {

Replacement in Hydro/src/initialize.cxx at line 126 [2.15599]

B:BD[2.20606] → [2.20606:20874]

          CCTK_REALVEC x1 = x0[0] + (i + CCTK_REAL(0.5) + viota()) * dx[0];
          CCTK_REALVEC y1 = x0[1] + (j + CCTK_REAL(0.5)) * dx[1];
          CCTK_REALVEC z1 = x0[2] + (k + CCTK_REAL(0.5)) * dx[2];
          CCTK_REALVEC r2 = pow2(x1) + pow2(y1) + pow2(z1);

[2.20606]

[2.20874]

          CCTK_REALVEC x1 = x0[0] + (lbnd[0] + i) * dx[0] + viota() * dx[0];
          CCTK_REALVEC y1 = x0[1] + (lbnd[1] + j) * dx[1];
          CCTK_REALVEC z1 = x0[2] + (lbnd[2] + k) * dx[2];
          CCTK_REALVEC r2 = pow2(x1) + (pow2(y1) + pow2(z1));

Replacement in Hydro/src/initialize.cxx at line 133 [2.15599]

B:BD[2.20992] → [2.20992:21184]

          CCTK_REALVEC momx1i = CCTK_REAL(0.0);
          CCTK_REALVEC momy1i = CCTK_REAL(0.0);
          CCTK_REALVEC momz1i = CCTK_REAL(0.0);
          CCTK_REALVEC etot1i = CCTK_REAL(2.0);

[2.20992]

[2.21184]

          CCTK_REALVEC velx1i = CCTK_REAL(0.0);
          CCTK_REALVEC vely1i = CCTK_REAL(0.0);
          CCTK_REALVEC velz1i = CCTK_REAL(0.0);
          CCTK_REALVEC eint1i = CCTK_REAL(2.0);

Replacement in Hydro/src/initialize.cxx at line 139 [2.15599]

B:BD[2.21256] → [2.21256:21448]

          CCTK_REALVEC momx1o = CCTK_REAL(0.0);
          CCTK_REALVEC momy1o = CCTK_REAL(0.0);
          CCTK_REALVEC momz1o = CCTK_REAL(0.0);
          CCTK_REALVEC etot1o = CCTK_REAL(1.0);

[2.21256]

[2.21448]

          CCTK_REALVEC velx1o = CCTK_REAL(0.0);
          CCTK_REALVEC vely1o = CCTK_REAL(0.0);
          CCTK_REALVEC velz1o = CCTK_REAL(0.0);
          CCTK_REALVEC eint1o = CCTK_REAL(1.0);

Replacement in Hydro/src/initialize.cxx at line 145 [2.15599]

B:BD[2.21531] → [2.21531:21795]

          CCTK_REALVEC momx1 = ifthen(r2 <= sr2, momx1i, momx1o);
          CCTK_REALVEC momy1 = ifthen(r2 <= sr2, momy1i, momy1o);
          CCTK_REALVEC momz1 = ifthen(r2 <= sr2, momz1i, momz1o);
          CCTK_REALVEC etot1 = ifthen(r2 <= sr2, etot1i, etot1o);

[2.21531]

[2.21795]

          CCTK_REALVEC velx1 = ifthen(r2 <= sr2, velx1i, velx1o);
          CCTK_REALVEC vely1 = ifthen(r2 <= sr2, vely1i, vely1o);
          CCTK_REALVEC velz1 = ifthen(r2 <= sr2, velz1i, velz1o);
          CCTK_REALVEC eint1 = ifthen(r2 <= sr2, eint1i, eint1o);

Replacement in Hydro/src/initialize.cxx at line 150 [2.15599]

B:BD[2.21857] → [2.21857:22113]

          momx1.storeu_partial(momx[ind], i, imin[0], imax[0]);
          momy1.storeu_partial(momy[ind], i, imin[0], imax[0]);
          momz1.storeu_partial(momz[ind], i, imin[0], imax[0]);
          etot1.storeu_partial(etot[ind], i, imin[0], imax[0]);

[2.21857]

[2.22113]

          velx1.storeu_partial(velx[ind], i, imin[0], imax[0]);
          vely1.storeu_partial(vely[ind], i, imin[0], imax[0]);
          velz1.storeu_partial(velz[ind], i, imin[0], imax[0]);
          eint1.storeu_partial(eint[ind], i, imin[0], imax[0]);

Insertion in Hydro/src/initialize.cxx at line 160 [2.15599]

[2.22164]

  }
  // Set an ignorable pressure
  for (int k = imin[2]; k < imax[2]; ++k) {
    for (int j = imin[1]; j < imax[1]; ++j) {
      for (int i = imin[0]; i < imax[0]; i += vsize) {
        ptrdiff_t ind = i + dj * j + dk * k;
        CCTK_REALVEC press1 = CCTK_REAL(-1.0);
        press1.storeu_partial(press[ind], i, imin[0], imax[0]);
      }
    }

Insertion in Hydro/src/make.code.defn at line 9 [2.22210]
[2.22395]
[2.22395]
```
	prim2con.cxx				\
```

File addition: prim2con.cxx (-rw-r--r--)

[2.4577]

#include "defs.hxx"
#include <loop.hxx>
#include <cctk.h>
#include <cctk_Arguments_Checked.h>
#include <cctk_Parameters.h>
#include <array>
#include <cmath>
namespace Hydro {
using namespace std;
////////////////////////////////////////////////////////////////////////////////
extern "C" void Hydro_Prim2Con(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_Hydro_Prim2Con;
  DECLARE_CCTK_PARAMETERS;
  const Loop::vect<int, dim> lsh{{cctk_lsh[0], cctk_lsh[1], cctk_lsh[2]}};
  array<CCTK_INT, dim> tmin_, tmax_;
  GetTileExtent(cctkGH, tmin_.data(), tmax_.data());
  const Loop::vect<int, dim> tmin(tmin_);
  const Loop::vect<int, dim> tmax(tmax_);
  const Loop::vect<int, dim> imin = tmin;
  const Loop::vect<int, dim> imax = tmax;
  const Loop::vect<int, dim> ash{{cctk_ash[0], cctk_ash[1], cctk_ash[2]}};
  constexpr ptrdiff_t di = 1;
  const ptrdiff_t dj = di * cctk_ash[0];
  const ptrdiff_t dk = dj * cctk_ash[1];
  for (int k = imin[2]; k < imax[2]; ++k) {
    for (int j = imin[1]; j < imax[1]; ++j) {
      for (int i = imin[0]; i < imax[0]; i += vsize) {
        ptrdiff_t ind = i + dj * j + dk * k;
        CCTK_REALVEC rho1 = vloadu(rho[ind]);
        CCTK_REALVEC velx1 = vloadu(velx[ind]);
        CCTK_REALVEC vely1 = vloadu(vely[ind]);
        CCTK_REALVEC velz1 = vloadu(velz[ind]);
        CCTK_REALVEC eint1 = vloadu(eint[ind]);
        CCTK_REALVEC ekin1 =
            CCTK_REAL(0.5) * rho1 * (pow2(velx1) + pow2(vely1) + pow2(velz1));
        CCTK_REALVEC dens1 = rho1;
        CCTK_REALVEC momx1 = rho1 * velx1;
        CCTK_REALVEC momy1 = rho1 * vely1;
        CCTK_REALVEC momz1 = rho1 * velz1;
        CCTK_REALVEC etot1 = ekin1 + eint1;
        dens1.storeu_partial(dens[ind], i, imin[0], imax[0]);
        momx1.storeu_partial(momx[ind], i, imin[0], imax[0]);
        momy1.storeu_partial(momy[ind], i, imin[0], imax[0]);
        momz1.storeu_partial(momz[ind], i, imin[0], imax[0]);
        etot1.storeu_partial(etot[ind], i, imin[0], imax[0]);
      }
    }
  }
}
} // namespace Hydro

Insertion in Hydro/src/rhs.cxx at line 51 [2.22468]

[2.24072]

  const ptrdiff_t offx =
      nghostzones[0] * dix + nghostzones[1] * djx + nghostzones[2] * dkx;

Insertion in Hydro/src/rhs.cxx at line 60 [2.22468]

[2.24294]

  const ptrdiff_t offy =
      nghostzones[0] * diy + nghostzones[1] * djy + nghostzones[2] * dky;

Insertion in Hydro/src/rhs.cxx at line 69 [2.22468]

[2.24516]

  const ptrdiff_t offz =
      nghostzones[0] * diz + nghostzones[1] * djz + nghostzones[2] * dkz;

Replacement in Hydro/src/rhs.cxx at line 79 [2.22468]

B:BD[2.24862] → [2.24862:24920]

            for (int i = imin[0]; i < imax[0]; i += VS) {

[2.24862]

[2.24920]

            for (int i = imin[0]; i < imax[0]; i += vsize) {

Replacement in Hydro/src/rhs.cxx at line 81 [2.22468]

B:BD[2.24971] → [2.24971:25133]

              ptrdiff_t indx = i + djx * j + dkx * k;
              ptrdiff_t indy = i + djy * j + dky * k;
              ptrdiff_t indz = i + djz * j + dkz * k;

[2.24971]

[2.25133]

              ptrdiff_t indx = i + djx * j + dkx * k - offx;
              ptrdiff_t indy = i + djy * j + dky * k - offy;
              ptrdiff_t indz = i + djz * j + dkz * k - offz;

Replacement in Hydro/src/rhs.cxx at line 90 [2.22468]

B:BD[2.25402] → [2.25402:25471]

              dtvar1.store_partial(dtvar[ind], i, imin[0], imax[0]);

[2.25402]

[2.25471]

              dtvar1.storeu_partial(dtvar[ind], i, imin[0], imax[0]);

Replacement in Hydro/src/rhs.cxx at line 96 [2.22468]
B:BD[2.25518] → [2.25518:25557]
```
  calcrhs(fxrho, fyrho, fzrho, dtrho);
```
[2.25518]
[2.25557]
```
  calcrhs(fxdens, fydens, fzdens, dtdens);
```