eschnett/cactusamrex - Change X7X5L4QNRVTEEOSUCKE3YCJLNZHJPGKTSR4H63SKNWA3YS4ID7GQC

Z4c: Update scheduling

Created by Erik Schnetter on June 2, 2020

X7X5L4QNRVTEEOSUCKE3YCJLNZHJPGKTSR4H63SKNWA3YS4ID7GQC

Dependencies

In channels

master

main

Change contents

Replacement in Z4c/schedule.ccl at line 69 [4.5677]
B:BD[3.5571] → [3.5571:5611]
```
SCHEDULE Z4c_Boundaries IN Z4c_PostStep
```
[3.5571]
[4.6859]
```
SCHEDULE Z4c_Sync IN Z4c_PostStep
```

Deletion in Z4c/schedule.ccl at line 72 [4.5677]

B:BD[4.6871] → [4.6871:7090]

  WRITES: chi(boundary)
  WRITES: gamma_tilde(boundary)
  WRITES: K_hat(boundary)
  WRITES: A_tilde(boundary)
  WRITES: Gam_tilde(boundary)
  WRITES: Theta(boundary)
  WRITES: alphaG(boundary)
  WRITES: betaG(boundary)

Insertion in Z4c/schedule.ccl at line 80 [4.5677]

[3.5735]

[4.7090]

} "Synchronize Z4c variables"
SCHEDULE Z4c_Boundaries IN Z4c_PostStep AFTER Z4c_Sync
{
  LANG: C
  WRITES: chi(boundary)
  WRITES: gamma_tilde(boundary)
  WRITES: K_hat(boundary)
  WRITES: A_tilde(boundary)
  WRITES: Gam_tilde(boundary)
  WRITES: Theta(boundary)
  WRITES: alphaG(boundary)
  WRITES: betaG(boundary)

Insertion in Z4c/schedule.ccl at line 127 [4.5677]

[4.7368]

SCHEDULE Z4c_ADM2 IN Z4c_PostStep AFTER Z4c_Enforce
{
  LANG: C
  READS: chi(everywhere)
  READS: gamma_tilde(everywhere)
  READS: K_hat(everywhere)
  READS: A_tilde(everywhere)
  READS: Gam_tilde(everywhere)
  READS: Theta(everywhere)
  READS: alphaG(everywhere)
  READS: betaG(everywhere)
  WRITES: ADMBase::dtcurv(interior)
  WRITES: ADMBase::dt2lapse(interior)
  WRITES: ADMBase::dt2shift(interior)
} "Calculate second time derivatives of ADM variables"

Insertion in Z4c/schedule.ccl at line 143 [4.5677]

[4.7369]

SCHEDULE Z4c_ADM2Boundaries IN Z4c_PostStep
{
  LANG: C
  WRITES: ADMBase::dtcurv(boundary)
  WRITES: ADMBase::dt2lapse(boundary)
  WRITES: ADMBase::dt2shift(boundary)
  SYNC: ADMBase::dtcurv
  SYNC: ADMBase::dt2lapse
  SYNC: ADMBase::dt2shift
} "Apply boundary conditions to second time derivatives of ADM variables"

Deletion in Z4c/schedule.ccl at line 205 [4.5677]

B:BD[4.8670] → [5.40:156]

  #
  WRITES: ADMBase::dtcurv(interior)
  WRITES: ADMBase::dt2lapse(interior)
  WRITES: ADMBase::dt2shift(interior)

Deletion in Z4c/schedule.ccl at line 226 [4.5677]

B:BD[6.55019] → [5.157:349]

  #
  WRITES: ADMBase::dtcurv(boundary)
  WRITES: ADMBase::dt2lapse(boundary)
  WRITES: ADMBase::dt2shift(boundary)
  SYNC: ADMBase::dtcurv
  SYNC: ADMBase::dt2lapse
  SYNC: ADMBase::dt2shift

Replacement in Z4c/schedule.ccl at line 266 [4.5677]

B:BD[3.7325] → [3.7325:7402]

SCHEDULE GROUP Z4c_PostSTep AT postrestrict
{
} "Post-process Z4c variables"

[3.7325]

# SCHEDULE GROUP Z4c_PostStep AT postrestrict
# {
# } "Post-process Z4c variables"

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

[4.8854]

#include "derivs.hxx"
#include "physics.hxx"
#include "tensor.hxx"
#include "z4c_vars.hxx"
#include <loop.hxx>
#include <cctk.h>
#include <cctk_Arguments_Checked.h>
#include <cctk_Parameters.h>
#include <cmath>
namespace Z4c {
using namespace Loop;
using namespace std;
extern "C" void Z4c_ADM2(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_Z4c_ADM2;
  DECLARE_CCTK_PARAMETERS;
  for (int d = 0; d < 3; ++d)
    if (cctk_nghostzones[d] < deriv_order / 2 + 1)
      CCTK_VERROR("Need at least %d ghost zones", deriv_order / 2 + 1);
  const vec3<CCTK_REAL, UP> dx{
      CCTK_DELTA_SPACE(0),
      CCTK_DELTA_SPACE(1),
      CCTK_DELTA_SPACE(2),
  };
  //
  const GF3D<const CCTK_REAL, 0, 0, 0> gf_chi_(cctkGH, chi);
  const mat3<GF3D<const CCTK_REAL, 0, 0, 0>, DN, DN> gf_gammat_(
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, gammatxx),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, gammatxy),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, gammatxz),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, gammatyy),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, gammatyz),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, gammatzz));
  const GF3D<const CCTK_REAL, 0, 0, 0> gf_Kh_(cctkGH, Kh);
  const mat3<GF3D<const CCTK_REAL, 0, 0, 0>, DN, DN> gf_At_(
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Atxx),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Atxy),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Atxz),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Atyy),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Atyz),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Atzz));
  const vec3<GF3D<const CCTK_REAL, 0, 0, 0>, UP> gf_Gamt_(
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Gamtx),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Gamty),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, Gamtz));
  const GF3D<const CCTK_REAL, 0, 0, 0> gf_Theta_(cctkGH, Theta);
  const GF3D<const CCTK_REAL, 0, 0, 0> gf_alphaG_(cctkGH, alphaG);
  const vec3<GF3D<const CCTK_REAL, 0, 0, 0>, UP> gf_betaG_(
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, betaGx),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, betaGy),
      GF3D<const CCTK_REAL, 0, 0, 0>(cctkGH, betaGz));
  //
  const vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN> gf_dchi_(cctkGH, allocate());
  const mat3<GF3D<CCTK_REAL, 0, 0, 0>, DN, DN> gf_ddchi_(cctkGH, allocate());
  calc_derivs2(cctkGH, gf_chi_, gf_dchi_, gf_ddchi_);
  const mat3<vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN>, DN, DN> gf_dgammat_(
      cctkGH, allocate());
  const mat3<mat3<GF3D<CCTK_REAL, 0, 0, 0>, DN, DN>, DN, DN> gf_ddgammat_(
      cctkGH, allocate());
  calc_derivs2(cctkGH, gf_gammat_, gf_dgammat_, gf_ddgammat_);
  const vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN> gf_dKh_(cctkGH, allocate());
  calc_derivs(cctkGH, gf_Kh_, gf_dKh_);
  const mat3<vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN>, DN, DN> gf_dAt_(cctkGH,
                                                                 allocate());
  calc_derivs(cctkGH, gf_At_, gf_dAt_);
  const vec3<vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN>, UP> gf_dGamt_(cctkGH,
                                                               allocate());
  calc_derivs(cctkGH, gf_Gamt_, gf_dGamt_);
  const vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN> gf_dTheta_(cctkGH, allocate());
  calc_derivs(cctkGH, gf_Theta_, gf_dTheta_);
  const vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN> gf_dalphaG_(cctkGH, allocate());
  const mat3<GF3D<CCTK_REAL, 0, 0, 0>, DN, DN> gf_ddalphaG_(cctkGH, allocate());
  calc_derivs2(cctkGH, gf_alphaG_, gf_dalphaG_, gf_ddalphaG_);
  const vec3<vec3<GF3D<CCTK_REAL, 0, 0, 0>, DN>, UP> gf_dbetaG_(cctkGH,
                                                                allocate());
  const vec3<mat3<GF3D<CCTK_REAL, 0, 0, 0>, DN, DN>, UP> gf_ddbetaG_(
      cctkGH, allocate());
  calc_derivs2(cctkGH, gf_betaG_, gf_dbetaG_, gf_ddbetaG_);
  //
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxx_(cctkGH, dtkxx);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxy_(cctkGH, dtkxy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxz_(cctkGH, dtkxz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyy_(cctkGH, dtkyy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyz_(cctkGH, dtkyz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkzz_(cctkGH, dtkzz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2alp_(cctkGH, dt2alp);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betax_(cctkGH, dt2betax);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betay_(cctkGH, dt2betay);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betaz_(cctkGH, dt2betaz);
  //
  loop_int<0, 0, 0>(cctkGH, [&](const PointDesc &p) Z4C_INLINE {
    // Load and calculate
    const CCTK_REAL rho{0};
    const vec3<CCTK_REAL, DN> Si{0, 0, 0};
    const mat3<CCTK_REAL, DN, DN> Sij{0, 0, 0, 0, 0, 0};
    const z4c_vars<CCTK_REAL> vars(
        kappa1, kappa2, f_mu_L, f_mu_S, eta,                  //
        gf_chi_(p.I), gf_dchi_(p.I), gf_ddchi_(p.I),          //
        gf_gammat_(p.I), gf_dgammat_(p.I), gf_ddgammat_(p.I), //
        gf_Kh_(p.I), gf_dKh_(p.I),                            //
        gf_At_(p.I), gf_dAt_(p.I),                            //
        gf_Gamt_(p.I), gf_dGamt_(p.I),                        //
        gf_Theta_(p.I), gf_dTheta_(p.I),                      //
        gf_alphaG_(p.I), gf_dalphaG_(p.I), gf_ddalphaG_(p.I), //
        gf_betaG_(p.I), gf_dbetaG_(p.I), gf_ddbetaG_(p.I),    //
        rho,                                                  //
        Si,                                                   //
        Sij);
    // Store
    vars.K_rhs.store(gf_dtkxx_, gf_dtkxy_, gf_dtkxz_, gf_dtkyy_, gf_dtkyz_,
                     gf_dtkzz_, p.I);
    gf_dt2alp_(p.I) = vars.dtalpha_rhs;
    vars.dtbeta_rhs.store(gf_dt2betax_, gf_dt2betay_, gf_dt2betaz_, p.I);
  });
}
} // namespace Z4c

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

[4.8854]

#include <loop.hxx>
#include <cctk.h>
#include <cctk_Arguments_Checked.h>
namespace Z4c {
using namespace Loop;
extern "C" void Z4c_ADM2Boundaries(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_Z4c_ADM2Boundaries;
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxx_(cctkGH, dtkxx);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxy_(cctkGH, dtkxy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxz_(cctkGH, dtkxz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyy_(cctkGH, dtkyy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyz_(cctkGH, dtkyz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkzz_(cctkGH, dtkzz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2alp_(cctkGH, dt2alp);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betax_(cctkGH, dt2betax);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betay_(cctkGH, dt2betay);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betaz_(cctkGH, dt2betaz);
  //
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxx_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxy_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxz_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkyy_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkyz_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkzz_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2alp_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betax_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betay_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betaz_(p.I) = 0; });
}
} // namespace Z4c

Insertion in Z4c/src/make.code.defn at line 6 [4.30723]
[4.30832]
[4.30832]
```
	adm2.cxx				\
	adm2boundaries.cxx			\
```
Insertion in Z4c/src/make.code.defn at line 16 [4.30723]
[3.13559]
[3.13559]
```
        sync.cxx				\
```

Deletion in Z4c/src/rhs.cxx at line 135 [4.35315]

B:BD[4.38764] → [5.1276:1693]


  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxx_(cctkGH, dtkxx);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxy_(cctkGH, dtkxy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxz_(cctkGH, dtkxz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyy_(cctkGH, dtkyy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyz_(cctkGH, dtkyz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkzz_(cctkGH, dtkzz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2alp_(cctkGH, dt2alp);

Deletion in Z4c/src/rhs.cxx at line 136 [4.35315]

B:BD[5.1694] → [5.1694:1889]

∅:D[5.1889] → [4.38764:38765]

B:BD[4.38764] → [4.38764:38765]

  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betax_(cctkGH, dt2betax);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betay_(cctkGH, dt2betay);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betaz_(cctkGH, dt2betaz);

Deletion in Z4c/src/rhs.cxx at line 171 [4.35315]

B:BD[7.2838] → [2.180:409]


    vars.K_rhs.store(gf_dtkxx_, gf_dtkxy_, gf_dtkxz_, gf_dtkyy_, gf_dtkyz_,
                     gf_dtkzz_, p.I);
    gf_dt2alp_(p.I) = vars.dtalpha_rhs;
    vars.dtbeta_rhs.store(gf_dt2betax_, gf_dt2betay_, gf_dt2betaz_, p.I);

Deletion in Z4c/src/rhsboundaries.cxx at line 41 [3.20735]

B:BD[3.22451] → [5.3908:4325]


  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxx_(cctkGH, dtkxx);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxy_(cctkGH, dtkxy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkxz_(cctkGH, dtkxz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyy_(cctkGH, dtkyy);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkyz_(cctkGH, dtkyz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtkzz_(cctkGH, dtkzz);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2alp_(cctkGH, dt2alp);

Deletion in Z4c/src/rhsboundaries.cxx at line 42 [3.20735]

B:BD[3.22452] → [5.4326:4522]

  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betax_(cctkGH, dt2betax);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betay_(cctkGH, dt2betay);
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dt2betaz_(cctkGH, dt2betaz);

Deletion in Z4c/src/rhsboundaries.cxx at line 85 [3.20735]

B:BD[3.24554] → [5.4523:5330]


  // ADMBase
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxx_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxy_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxz_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkyy_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkyz_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkzz_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2alp_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betax_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betay_(p.I) = 0; });
  loop_bnd<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betaz_(p.I) = 0; });

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

[4.8854]

#include <loop.hxx>
#include <cctk.h>
#include <cctk_Arguments_Checked.h>
namespace Z4c {
using namespace Loop;
extern "C" void Z4c_Sync(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_Z4c_Sync;
}
} // namespace Z4c

Z4c: Update scheduling

Dependencies

In channels

Change contents

Replacement in Z4c/schedule.ccl at line 69 [4.5677]

Deletion in Z4c/schedule.ccl at line 72 [4.5677]

Insertion in Z4c/schedule.ccl at line 80 [4.5677]

Insertion in Z4c/schedule.ccl at line 127 [4.5677]

Insertion in Z4c/schedule.ccl at line 143 [4.5677]

Deletion in Z4c/schedule.ccl at line 205 [4.5677]

Deletion in Z4c/schedule.ccl at line 226 [4.5677]

Replacement in Z4c/schedule.ccl at line 266 [4.5677]

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

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

Insertion in Z4c/src/make.code.defn at line 6 [4.30723]

Insertion in Z4c/src/make.code.defn at line 16 [4.30723]

Deletion in Z4c/src/rhs.cxx at line 135 [4.35315]

Deletion in Z4c/src/rhs.cxx at line 136 [4.35315]

Deletion in Z4c/src/rhs.cxx at line 171 [4.35315]

Deletion in Z4c/src/rhsboundaries.cxx at line 41 [3.20735]

Deletion in Z4c/src/rhsboundaries.cxx at line 42 [3.20735]

Deletion in Z4c/src/rhsboundaries.cxx at line 85 [3.20735]

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