eschnett/cactusamrex - Change 2HBMXCNDT3QJWGTNHCJX3IBQ7H3QEGFGMELYYO63P3WDPBTP3JWQC

Z4c: Set ADM RHS variables (but not yet correctly)

Created by Erik Schnetter on April 3, 2020

2HBMXCNDT3QJWGTNHCJX3IBQ7H3QEGFGMELYYO63P3WDPBTP3JWQC

Dependencies

In channels

master

main

Change contents

Insertion in Z4c/schedule.ccl at line 116 [6.5677]
[5.6166]
[5.6166]
```
  WRITES: ADMBase::dtlapse(everywhere)
```

Insertion in Z4c/schedule.ccl at line 172 [6.5677]

[6.8670]

[5.6375]

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

Insertion in Z4c/schedule.ccl at line 197 [6.5677]

[7.55019]

[5.6713]

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

Insertion in Z4c/src/adm.cxx at line 61 [6.8865]
[6.11101]
[6.11101]
```
  const GF3D<CCTK_REAL, 0, 0, 0> gf_dtalp_(cctkGH, dtalp);
```

Replacement in Z4c/src/adm.cxx at line 73 [6.8865]

∅:D[3.65] → [7.55110:55647]

B:BD[7.55110] → [7.55110:55647]

        // Load and calculate
        const z4c_vars_noderivs<CCTK_REAL> vars(
            kappa1, kappa2, f_mu_L, f_mu_S, eta, //
            gf_chi_, gf_gammatxx_, gf_gammatxy_, gf_gammatxz_, gf_gammatyy_,
            gf_gammatyz_, gf_gammatzz_, gf_Kh_, gf_Atxx_, gf_Atxy_, gf_Atxz_,
            gf_Atyy_, gf_Atyz_, gf_Atzz_,
            //
            // gf_Gamtx_, gf_Gamty_, gf_Gamtz_,
            gf_chi_, gf_chi_, gf_chi_,
            //
            gf_Theta_, gf_alphaG_, gf_betaGx_, gf_betaGy_, gf_betaGz_, //
            p.I);

[3.65]

[6.12514]

    // Load and calculate
    const z4c_vars_noderivs<CCTK_REAL> vars(
        kappa1, kappa2, f_mu_L, f_mu_S, eta, //
        gf_chi_, gf_gammatxx_, gf_gammatxy_, gf_gammatxz_, gf_gammatyy_,
        gf_gammatyz_, gf_gammatzz_, gf_Kh_, gf_Atxx_, gf_Atxy_, gf_Atxz_,
        gf_Atyy_, gf_Atyz_, gf_Atzz_,
        //
        // gf_Gamtx_, gf_Gamty_, gf_Gamtz_,
        gf_chi_, gf_chi_, gf_chi_,
        //
        gf_Theta_, gf_alphaG_, gf_betaGx_, gf_betaGy_, gf_betaGz_, //
        p.I);
    // Store
    vars.g.store(gf_gxx_, gf_gxy_, gf_gxz_, gf_gyy_, gf_gyz_, gf_gzz_, p.I);
    vars.k.store(gf_kxx_, gf_kxy_, gf_kxz_, gf_kyy_, gf_kyz_, gf_kzz_, p.I);
    gf_alp_(p.I) = vars.alp;
    vars.beta.store(gf_betax_, gf_betay_, gf_betaz_, p.I);
  });

Replacement in Z4c/src/adm.cxx at line 93 [6.8865]

B:BD[6.12515] → [7.55648:55933]

        // Store
        vars.g.store(gf_gxx_, gf_gxy_, gf_gxz_, gf_gyy_, gf_gyz_, gf_gzz_, p.I);
        vars.k.store(gf_kxx_, gf_kxy_, gf_kxz_, gf_kyy_, gf_kyz_, gf_kzz_, p.I);
        gf_alp_(p.I) = vars.alp;
        vars.beta.store(gf_betax_, gf_betay_, gf_betaz_, p.I);
      });

[6.12515]

[6.12750]

#warning "TODO: These are wrong!"
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtalp_(p.I) = 0; });

Insertion in Z4c/src/rhs.cxx at line 135 [6.35315]

[6.38764]


  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);

Replacement in Z4c/src/rhs.cxx at line 152 [6.35315]
∅:D[3.131] → [7.80013:80043]
B:BD[7.80013] → [7.80013:80043]
```
        // Load and calculate
```
[3.131]
[7.80043]
```
    // Load and calculate
```

Replacement in Z4c/src/rhs.cxx at line 154 [6.35315]

B:BD[7.80044] → [7.80044:80184]

        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};

[7.80044]

[7.80184]

    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};

Replacement in Z4c/src/rhs.cxx at line 158 [6.35315]

B:BD[7.80185] → [7.80185:81002]

        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);

[7.80185]

[7.81002]

    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);

Replacement in Z4c/src/rhs.cxx at line 172 [6.35315]

B:BD[7.81003] → [7.81003:81061]

        // Store
        gf_chi_rhs_(p.I) = vars.chi_rhs;

[7.81003]

[4.2431]

    // Store
    gf_chi_rhs_(p.I) = vars.chi_rhs;

Replacement in Z4c/src/rhs.cxx at line 177 [6.35315]
∅:D[4.2623] → [7.81265:81304]
B:BD[7.81265] → [7.81265:81304]
```
        gf_Kh_rhs_(p.I) = vars.Kh_rhs;
```
[4.2623]
[4.2624]
```
    gf_Kh_rhs_(p.I) = vars.Kh_rhs;
```

Replacement in Z4c/src/rhs.cxx at line 180 [6.35315]

∅:D[4.2758] → [7.81446:81617]

B:BD[7.81446] → [7.81446:81617]

        vars.Gamt_rhs.store(gf_Gamtx_rhs_, gf_Gamty_rhs_, gf_Gamtz_rhs_, p.I);
        gf_Theta_rhs_(p.I) = vars.Theta_rhs;
        gf_alphaG_rhs_(p.I) = vars.alphaG_rhs;

[4.2758]

[4.2759]

    vars.Gamt_rhs.store(gf_Gamtx_rhs_, gf_Gamty_rhs_, gf_Gamtz_rhs_, p.I);
    gf_Theta_rhs_(p.I) = vars.Theta_rhs;
    gf_alphaG_rhs_(p.I) = vars.alphaG_rhs;

Replacement in Z4c/src/rhs.cxx at line 184 [6.35315]
∅:D[4.2838] → [7.81729:81739]
B:BD[7.81729] → [7.81729:81739]
```
      });
```
[4.2838]
[6.39344]
```
  });
```

Insertion in Z4c/src/rhs.cxx at line 249 [6.35315]

[7.84129]

[2.0]


  // ADMBase
#warning "TODO: These are wrong!"
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxx_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxy_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkxz_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkyy_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkyz_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dtkzz_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2alp_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betax_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betay_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { gf_dt2betaz_(p.I) = 0; });

Insertion in Z4c/src/rhsboundaries.cxx at line 41 [5.20735]

[5.22451]


  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);

Insertion in Z4c/src/rhsboundaries.cxx at line 51 [5.20735]

[5.22452]

  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);

Insertion in Z4c/src/rhsboundaries.cxx at line 98 [5.20735]

[5.24554]


  // 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; });