eschnett/cactusamrex - Change FHE6M7GL6TJSS4VHGV6XY5PQSZBCU6SIXAGZG67Z76NRGKZZAGTQC

Z4c: Much work, much debugging

Created by Erik Schnetter on January 23, 2020

FHE6M7GL6TJSS4VHGV6XY5PQSZBCU6SIXAGZG67Z76NRGKZZAGTQC

Dependencies

In channels

master

main

Change contents

File addition: linear-wave-ctgamma.par (-rw-r--r--)

[3.2485]

ActiveThorns = "
    ADMBase
    CTGBase
    CTGEvolution
    CTGGauge
    Carpet
    CarpetIOASCII
    CarpetIOBasic
    CarpetIOScalar
    CarpetLib
    CartGrid3D
    CoordBase
    CartesianCoordinates
    Formaline
    GlobalDerivative
    IDLinearWaves
    IOUtil
    MoL
    PeriodicCarpet
    SphericalSurface
    StaticConformal
    SummationByParts
    SymBase
    Time
    CarpetReduce
"
Cactus::cctk_show_schedule = yes
# Cactus::terminate = "time"
# Cactus::cctk_final_time = 0.0
Cactus::cctk_itlast = 1
CartGrid3D::type = "coordbase"
CartGrid3D::domain = "full"
CartGrid3D::avoid_origin = no
CoordBase::xmin = 0.0
CoordBase::ymin = 0.0
CoordBase::zmin = 0.0
CoordBase::xmax = 1.0
CoordBase::ymax = 2.0 / 64
CoordBase::zmax = 2.0 / 64
CoordBase::dx = 1.0 / 64
CoordBase::dy = 1.0 / 64
CoordBase::dz = 1.0 / 64
CoordBase::boundary_size_x_lower = 2
CoordBase::boundary_size_y_lower = 2
CoordBase::boundary_size_z_lower = 2
CoordBase::boundary_size_x_upper = 2
CoordBase::boundary_size_y_upper = 2
CoordBase::boundary_size_z_upper = 2
CoordBase::boundary_shiftout_x_lower = 1
CoordBase::boundary_shiftout_y_lower = 1
CoordBase::boundary_shiftout_z_lower = 1
CoordBase::boundary_shiftout_x_upper = 0
CoordBase::boundary_shiftout_y_upper = 0
CoordBase::boundary_shiftout_z_upper = 0
Carpet::domain_from_coordbase = yes
PeriodicCarpet::periodic = yes
driver::ghost_size = 2
ADMBase::initial_data = "sine_planewaves"
ADMBase::initial_lapse = "one"
ADMBase::initial_shift = "zero"
IDLinearWaves::amplitude = 1.0e-8
IDLinearWaves::wavelength = 1.0
IDLinearWaves::wavetheta = 90
ADMBase::evolution_method = "ctgamma"
SummationByParts::order = 2
CTGBase::schedule_rhs_at_analysis = yes
CTGBase::conformal_factor_type = "chi"
CTGBase::evolution_system_type = "Z4c"
CTGEvolution::MaxNumEvolvedVars = 18
CTGEvolution::MaxNumConstrainedVars = 13
CTGEvolution::force_lndetg_zero = yes
CTGEvolution::kappa1 = 0.02
CTGEvolution::kappa2 = 0.0
Time::dtfac = 0.5
MoL::ODE_Method = "Euler"
MoL::MoL_Intermediate_Steps = 1
IO::out_dir = $parfile
IO::out_every = 1
IOASCII::out1d_vars = "
    ADMBase::metric
    ADMBase::curv
    ADMBase::lapse
    ADMBase::shift
    ADMBase::dtlapse
    ADMBase::dtshift
    CTGBase::conformal_factor
    CTGBase::conformal_metric
    CTGBase::curvature_scalar
    CTGBase::curvature_scalar_a
    CTGBase::curvature_scalar_b
    CTGBase::curvature_tensor
    CTGBase::Gamma
    CTGEvolution::conformal_factor_source
    CTGEvolution::conformal_metric_source
    CTGEvolution::curvature_scalar_source
    CTGEvolution::Theta_source
    CTGEvolution::curvature_tensor_source
    CTGEvolution::Gamma_source
    CTGGauge::shift_aux
    CTGGauge::shift_aux_source
"

Deletion in Z4c/par/linear-wave.par at line 4 [2.100]
B:BD[2.142] → [2.142:161]
```
    ErrorEstimator
```
Deletion in Z4c/par/linear-wave.par at line 10 [2.100]
B:BD[2.213] → [2.213:226]
```
$nlevels = 1
```

Replacement in Z4c/par/linear-wave.par at line 25 [2.100]

B:BD[2.454] → [2.454:514]

CarpetX::ymax = 5.0 / $ncells
CarpetX::zmax = 5.0 / $ncells

[2.454]

[2.514]

CarpetX::ymax = 2.0 / $ncells
CarpetX::zmax = 2.0 / $ncells

Deletion in Z4c/par/linear-wave.par at line 36 [2.100]

B:BD[2.683] → [2.683:988]

# CarpetX::periodic_x = no
# CarpetX::periodic_y = no
# CarpetX::periodic_z = no
# CarpetX::periodic = no
# CarpetX::reflection_x = yes
# CarpetX::reflection_y = yes
# CarpetX::reflection_z = yes
# CarpetX::reflection_upper_x = yes
# CarpetX::reflection_upper_y = yes
# CarpetX::reflection_upper_z = yes

Deletion in Z4c/par/linear-wave.par at line 38 [2.100]

B:BD[2.1013] → [2.1013:1269]

CarpetX::max_num_levels = $nlevels
CarpetX::regrid_every = 16
CarpetX::regrid_error_threshold = 1.0 / 4.0
ErrorEstimator::region_shape = "cube"
ErrorEstimator::scale_by_resolution = yes
CarpetX::prolongation_type = "ddf"
CarpetX::prolongation_order = 1

Replacement in Z4c/par/linear-wave.par at line 47 [2.100]
B:BD[2.1502] → [2.1502:1521]
```
Z4c::epsdiss = 0.0
```
[2.1502]
[2.1521]
```
Z4c::epsdiss = 0.32
```
Replacement in Z4c/par/linear-wave.par at line 50 [2.100]
B:BD[2.1545] → [2.1545:1605]
```
IO::out_every = 1   #TODO 2 * $ncells * 2 ** ($nlevels - 1)
```
[2.1545]
[2.1605]
```
IO::out_every = 2 * $ncells
```
Replacement in Z4c/par/linear-wave.par at line 53 [2.100]
B:BD[2.1640] → [2.1640:1674]
```
CarpetX::out_tsv = yes   #TODO no
```
[2.1640]
```
CarpetX::out_tsv = yes
```
Replacement in Z4c/src/constraints.cxx at line 25 [3.16341]
B:BD[3.16635] → [3.16635:16663]
```
  const vec3<CCTK_REAL> dx{
```
[3.16635]
[3.16663]
```
  const vec3<CCTK_REAL, UP> dx{
```

Replacement in Z4c/src/enforce.cxx at line 33 [3.22673]

B:BD[3.23693] → [3.23693:23747]

  loop_int<0, 0, 0>(cctkGH, [&](const PointDesc &p) {

[3.23693]

[3.23747]

  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) {

Replacement in Z4c/src/enforce.cxx at line 35 [3.22673]

B:BD[3.23759] → [3.23759:23828]

B:BD[3.23828] → [2.9198:9273]

∅:D[2.9273] → [3.23901:23974]

B:BD[3.23901] → [3.23901:23974]

B:BD[3.23974] → [2.9274:9313]

    mat3<CCTK_REAL> gammat(gf_gammatxx_, gf_gammatxy_, gf_gammatxz_,
                           gf_gammatyy_, gf_gammatyz_, gf_gammatzz_, p.I);
    mat3<CCTK_REAL> At(gf_Atxx_, gf_Atxy_, gf_Atxz_, gf_Atyy_, gf_Atyz_,
                       gf_Atzz_, p.I);

[3.23759]

[3.24011]

    const mat3<CCTK_REAL, DN, DN> gammat_old(gf_gammatxx_, gf_gammatxy_,
                                             gf_gammatxz_, gf_gammatyy_,
                                             gf_gammatyz_, gf_gammatzz_, p.I);
    const mat3<CCTK_REAL, DN, DN> At_old(gf_Atxx_, gf_Atxy_, gf_Atxz_, gf_Atyy_,
                                         gf_Atyz_, gf_Atzz_, p.I);

Replacement in Z4c/src/enforce.cxx at line 44 [3.22673]

B:BD[3.24086] → [2.9314:9411]

∅:D[2.9411] → [3.24182:24248]

B:BD[3.24182] → [3.24182:24248]

B:BD[3.24248] → [2.9412:9445]

    const CCTK_REAL detgammat = gammat.det();
    const CCTK_REAL chi_new = 1 / cbrt(detgammat);
    for (int a = 0; a < 3; ++a)
      for (int b = a; b < 3; ++b)
        gammat(a, b) *= chi_new;

[3.24086]

[2.9445]

    const CCTK_REAL detgammat_old = gammat_old.det();
    const CCTK_REAL chi_old = 1 / cbrt(detgammat_old);
    const mat3<CCTK_REAL, DN, DN> gammat(
        [&](int a, int b) { return chi_old * gammat_old(a, b); });

Replacement in Z4c/src/enforce.cxx at line 49 [3.22673]

B:BD[2.9463] → [2.9463:9513]

    const CCTK_REAL detgammat_new = gammat.det();

[2.9463]

[2.9513]

    const CCTK_REAL detgammat = gammat.det();

Replacement in Z4c/src/enforce.cxx at line 52 [3.22673]

B:BD[2.9612] → [2.9612:9676]

    if (!(fabs(detgammat_new - 1) <= 1.0e-12 * gammat_scale)) {

[2.9612]

[2.9676]

    if (!(fabs(detgammat - 1) <= 1.0e-12 * gammat_scale)) {

Replacement in Z4c/src/enforce.cxx at line 55 [3.22673]

B:BD[2.9757] → [2.9757:9804]

          << " det(gammat)=" << detgammat_new;

[2.9757]

[2.9804]

          << " det(gammat)=" << detgammat;

Replacement in Z4c/src/enforce.cxx at line 58 [3.22673]

B:BD[2.9854] → [2.9854:9917]

    assert(fabs(detgammat_new - 1) <= 1.0e-12 * gammat_scale);

[2.9854]

[2.9917]

    assert(fabs(detgammat - 1) <= 1.0e-12 * gammat_scale);

Replacement in Z4c/src/enforce.cxx at line 61 [3.22673]

B:BD[3.24279] → [3.24279:24330]

    const mat3<CCTK_REAL> gammatu = gammat.inv(1);

[3.24279]

[3.24330]

    const mat3<CCTK_REAL, UP, UP> gammatu = gammat.inv(1);

Replacement in Z4c/src/enforce.cxx at line 63 [3.22673]

B:BD[3.24331] → [2.9925:10025]

∅:D[2.10025] → [3.24470:24536]

B:BD[3.24470] → [3.24470:24536]

B:BD[3.24536] → [2.10026:10074]

    const CCTK_REAL traceAt =
        sum2([&](int x, int y) { return gammatu(x, y) * At(x, y); });
    for (int a = 0; a < 3; ++a)
      for (int b = a; b < 3; ++b)
        At(a, b) -= traceAt / 3 * gammat(a, b);

[3.24331]

[2.10074]

    const CCTK_REAL traceAt_old =
        sum2([&](int x, int y) { return gammatu(x, y) * At_old(x, y); });
    const mat3<CCTK_REAL, DN, DN> At([&](int a, int b) {
      return At_old(a, b) - traceAt_old / 3 * gammat(a, b);
    });

Replacement in Z4c/src/enforce.cxx at line 69 [3.22673]
B:BD[2.10092] → [2.10092:10126]
```
    const CCTK_REAL traceAt_new =
```
[2.10092]
[2.10126]
```
    const CCTK_REAL traceAt =
```

Replacement in Z4c/src/enforce.cxx at line 74 [3.22673]

B:BD[2.10371] → [2.10371:10425]

    if (!(fabs(traceAt_new) <= 1.0e-12 * At_scale)) {

[2.10371]

[2.10425]

    if (!(fabs(traceAt) <= 1.0e-12 * At_scale)) {

Replacement in Z4c/src/enforce.cxx at line 76 [3.22673]

B:BD[2.10450] → [2.10450:10512]

      buf << "tr At: At=" << At << " tr(At)=" << traceAt_new;

[2.10450]

[2.10512]

      buf << "tr At: At=" << At << " tr(At)=" << traceAt;

Replacement in Z4c/src/enforce.cxx at line 79 [3.22673]

B:BD[2.10562] → [2.10562:10615]

    assert(fabs(traceAt_new) <= 1.0e-12 * At_scale);

[2.10562]

[2.10615]

    assert(fabs(traceAt) <= 1.0e-12 * At_scale);

Replacement in Z4c/src/initial1.cxx at line 65 [3.24845]

B:BD[3.27318] → [3.27318:27391]

B:BD[3.27391] → [2.12597:12640]

∅:D[2.12640] → [3.27432:27505]

B:BD[3.27432] → [3.27432:27505]

B:BD[3.27505] → [2.12641:12684]

    const mat3<CCTK_REAL> g(gf_gxx_, gf_gxy_, gf_gxz_, gf_gyy_, gf_gyz_,
                            gf_gzz_, p.I);
    const mat3<CCTK_REAL> k(gf_kxx_, gf_kxy_, gf_kxz_, gf_kyy_, gf_kyz_,
                            gf_kzz_, p.I);

[3.27318]

[3.27546]

    const mat3<CCTK_REAL, DN, DN> g(gf_gxx_, gf_gxy_, gf_gxz_, gf_gyy_, gf_gyz_,
                                    gf_gzz_, p.I);
    const mat3<CCTK_REAL, DN, DN> k(gf_kxx_, gf_kxy_, gf_kxz_, gf_kyy_, gf_kyz_,
                                    gf_kzz_, p.I);

Replacement in Z4c/src/initial1.cxx at line 70 [3.24845]

B:BD[3.27586] → [2.12685:12755]

    const vec3<CCTK_REAL> beta(gf_betax_, gf_betay_, gf_betaz_, p.I);

[3.27586]

[3.27654]

    const vec3<CCTK_REAL, UP> beta(gf_betax_, gf_betay_, gf_betaz_, p.I);

Replacement in Z4c/src/initial1.cxx at line 74 [3.24845]
B:BD[3.27742] → [3.27742:27786]
```
    const mat3<CCTK_REAL> gu = g.inv(detg);
```
[3.27742]
[3.27786]
```
    const mat3<CCTK_REAL, UP, UP> gu = g.inv(detg);
```

Replacement in Z4c/src/initial1.cxx at line 78 [3.24845]

B:BD[3.27862] → [2.12799:12878]

    const mat3<CCTK_REAL> gammat([&](int a, int b) { return chi * g(a, b); });

[3.27862]

[3.27995]

    const mat3<CCTK_REAL, DN, DN> gammat(
        [&](int a, int b) { return chi * g(a, b); });

Replacement in Z4c/src/initial1.cxx at line 87 [3.24845]
B:BD[3.28189] → [2.12962:12992]
```
    const mat3<CCTK_REAL> At(
```
[3.28189]
[2.12992]
```
    const mat3<CCTK_REAL, DN, DN> At(
```

Replacement in Z4c/src/initial1.cxx at line 92 [3.24845]

B:BD[2.13068] → [2.13068:13133]

    const vec3<CCTK_REAL> betaG([&](int a) { return beta(a); });

[2.13068]

[3.28454]

    const vec3<CCTK_REAL, UP> betaG([&](int a) { return beta(a); });

Replacement in Z4c/src/initial2.cxx at line 19 [3.28871]
B:BD[3.29159] → [3.29159:29187]
```
  const vec3<CCTK_REAL> dx{
```
[3.29159]
[3.29187]
```
  const vec3<CCTK_REAL, UP> dx{
```

Replacement in Z4c/src/initial2.cxx at line 38 [3.28871]

B:BD[3.29945] → [3.29945:30020]

B:BD[3.30020] → [2.13156:13237]

    const mat3<CCTK_REAL> gammat(gf_gammatxx_, gf_gammatxy_, gf_gammatxz_,
                                 gf_gammatyy_, gf_gammatyz_, gf_gammatzz_, p.I);

[3.29945]

[3.30099]

    const mat3<CCTK_REAL, DN, DN> gammat(gf_gammatxx_, gf_gammatxy_,
                                         gf_gammatxz_, gf_gammatyy_,
                                         gf_gammatyz_, gf_gammatzz_, p.I);

Replacement in Z4c/src/initial2.cxx at line 43 [3.28871]

B:BD[3.30143] → [3.30143:30194]

    const mat3<CCTK_REAL> gammatu = gammat.inv(1);

[3.30143]

[3.30194]

    const mat3<CCTK_REAL, UP, UP> gammatu = gammat.inv(1);

Replacement in Z4c/src/initial2.cxx at line 45 [3.28871]
B:BD[3.30195] → [3.30195:30237]
```
    const mat3<vec3<CCTK_REAL> > dgammat{
```
[3.30195]
[3.30237]
```
    const mat3<vec3<CCTK_REAL, DN>, DN, DN> dgammat{
```

Replacement in Z4c/src/initial2.cxx at line 51 [3.28871]

B:BD[3.30449] → [2.13238:13419]

    const vec3<mat3<CCTK_REAL> > Gammatl = calc_gammal(dgammat);
    const vec3<mat3<CCTK_REAL> > Gammat = calc_gamma(gammatu, Gammatl);
    const vec3<CCTK_REAL> Gamt([&](int a) {

[3.30449]

[2.13419]

    const vec3<mat3<CCTK_REAL, DN, DN>, DN> Gammatl = calc_gammal(dgammat);
    const vec3<mat3<CCTK_REAL, DN, DN>, UP> Gammat =
        calc_gamma(gammatu, Gammatl);
    const vec3<CCTK_REAL, UP> Gamt([&](int a) {

Replacement in Z4c/src/physics.hxx at line 9 [3.31007]

B:BD[3.31111] → [2.13570:13684]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<vec3<T> >
calc_dgu(const mat3<T> &gu, const mat3<vec3<T> > &dg) {

[3.31111]

[3.31182]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<vec3<T, DN>, UP, UP>
calc_dgu(const mat3<T, UP, UP> &gu, const mat3<vec3<T, DN>, DN, DN> &dg) {

Replacement in Z4c/src/physics.hxx at line 18 [3.31007]

B:BD[3.31495] → [2.13685:13761]

  return mat3<vec3<T> >([&](int a, int b) {
    return vec3<T>([&](int c) {

[3.31495]

[2.13761]

  return mat3<vec3<T, DN>, UP, UP>([&](int a, int b) {
    return vec3<T, DN>([&](int c) {

Replacement in Z4c/src/physics.hxx at line 27 [3.31007]

B:BD[3.31831] → [2.13872:14040]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<vec3<T> >
calc_dAu(const mat3<T> &gu, const mat3<vec3<T> > &dgu, const mat3<T> &A,
         const mat3<vec3<T> > &dA) {

[3.31831]

[3.31971]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<vec3<T, DN>, UP, UP>
calc_dAu(const mat3<T, UP, UP> &gu, const mat3<vec3<T, DN>, UP, UP> &dgu,
         const mat3<T, DN, DN> &A, const mat3<vec3<T, DN>, DN, DN> &dA) {

Replacement in Z4c/src/physics.hxx at line 32 [3.31007]

B:BD[3.32073] → [2.14041:14117]

  return mat3<vec3<T> >([&](int a, int b) {
    return vec3<T>([&](int c) {

[3.32073]

[2.14117]

  return mat3<vec3<T, DN>, UP, UP>([&](int a, int b) {
    return vec3<T, DN>([&](int c) {

Replacement in Z4c/src/physics.hxx at line 44 [3.31007]

B:BD[3.32515] → [2.14340:14437]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<mat3<T> >
calc_gammal(const mat3<vec3<T> > dg) {

[3.32515]

[3.32715]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<mat3<T, DN, DN>, DN>
calc_gammal(const mat3<vec3<T, DN>, DN, DN> dg) {

Replacement in Z4c/src/physics.hxx at line 47 [3.31007]

B:BD[3.32731] → [2.14438:14514]

  return vec3<mat3<T> >([&](int a) {
    return mat3<T>([&](int b, int c) {

[3.32731]

[2.14514]

  return vec3<mat3<T, DN, DN>, DN>([&](int a) {
    return mat3<T, DN, DN>([&](int b, int c) {

Replacement in Z4c/src/physics.hxx at line 55 [3.31007]

B:BD[3.33355] → [2.14589:14709]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<mat3<T> >
calc_gamma(const mat3<T> &gu, const vec3<mat3<T> > &Gammal) {

[3.33355]

[2.14709]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<mat3<T, DN, DN>, UP>
calc_gamma(const mat3<T, UP, UP> &gu, const vec3<mat3<T, DN, DN>, DN> &Gammal) {

Replacement in Z4c/src/physics.hxx at line 58 [3.31007]

B:BD[2.14725] → [2.14725:14801]

  return vec3<mat3<T> >([&](int a) {
    return mat3<T>([&](int b, int c) {

[2.14725]

[2.14801]

  return vec3<mat3<T, DN, DN>, UP>([&](int a) {
    return mat3<T, DN, DN>([&](int b, int c) {

Replacement in Z4c/src/rhs.cxx at line 27 [3.35315]
B:BD[2.15320] → [2.15320:15348]
```
  const vec3<CCTK_REAL> dx{
```
[2.15320]
[2.15348]
```
  const vec3<CCTK_REAL, UP> dx{
```

Replacement in Z4c/src/rhs.cxx at line 34 [3.35315]

B:BD[2.15489] → [2.15489:15636]

    const vec3<CCTK_REAL> betaG(gf_betaGx_, gf_betaGy_, gf_betaGz_, p.I);
    const vec3<CCTK_REAL> dgf_upwind(deriv_upwind(gf_, p.I, betaG, dx));

[2.15489]

[2.15636]

    const vec3<CCTK_REAL, UP> betaG(gf_betaGx_, gf_betaGy_, gf_betaGz_, p.I);
    const vec3<CCTK_REAL, DN> dgf_upwind(deriv_upwind(gf_, p.I, betaG, dx));

Replacement in Z4c/src/rhs.cxx at line 47 [3.35315]
B:BD[2.15946] → [2.15946:15974]
```
  const vec3<CCTK_REAL> dx{
```
[2.15946]
[2.15974]
```
  const vec3<CCTK_REAL, UP> dx{
```
Replacement in Z4c/src/rhs.cxx at line 66 [3.35315]
B:BD[3.35686] → [3.35686:35714]
```
  const vec3<CCTK_REAL> dx{
```
[3.35686]
[3.35714]
```
  const vec3<CCTK_REAL, UP> dx{
```

Insertion in Z4c/src/tensor.hxx at line 65 [3.48124]

[3.48805]

enum class dnup_t : bool { dn, up };
constexpr dnup_t DN = dnup_t::dn;
constexpr dnup_t UP = dnup_t::up;
constexpr dnup_t operator!(const dnup_t dnup) { return dnup_t(!bool(dnup)); }

Replacement in Z4c/src/tensor.hxx at line 71 [3.48124]
B:BD[3.48817] → [3.48817:48852]
```
template <typename T> class vec3 {
```
[3.48817]
[3.48852]
```
template <typename T, dnup_t dnup> class vec3 {
```

Replacement in Z4c/src/tensor.hxx at line 114 [3.48124]

B:BD[3.49677] → [2.27100:27180]

    if (!((!CCTK_isnan(v(0))) && (!CCTK_isnan(v(1))) && (!CCTK_isnan(v(2))))) {

[3.49677]

[2.27180]

    if (!((CCTK_isfinite(v(0))) && (CCTK_isfinite(v(1))) &&
          (CCTK_isfinite(v(2))))) {

Replacement in Z4c/src/tensor.hxx at line 120 [3.48124]

∅:D[2.27290] → [3.49677:49770]

B:BD[3.49677] → [3.49677:49770]

    assert(!CCTK_isnan(v(0)));
    assert(!CCTK_isnan(v(1)));
    assert(!CCTK_isnan(v(2)));

[2.27290]

[3.49770]

    assert(CCTK_isfinite(v(0)));
    assert(CCTK_isfinite(v(1)));
    assert(CCTK_isfinite(v(2)));

Replacement in Z4c/src/tensor.hxx at line 134 [3.48124]

B:BD[2.27392] → [2.27392:27484]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T>
  operator+(const vec3<T> &x) {

[2.27392]

[2.27484]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T, dnup>
  operator+(const vec3<T, dnup> &x) {

Replacement in Z4c/src/tensor.hxx at line 138 [3.48124]

B:BD[2.27510] → [2.27510:27602]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T>
  operator-(const vec3<T> &x) {

[2.27510]

[2.27602]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T, dnup>
  operator-(const vec3<T, dnup> &x) {

Replacement in Z4c/src/tensor.hxx at line 142 [3.48124]

B:BD[2.27628] → [2.27628:27738]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T>
  operator+(const vec3<T> &x, const vec3<T> &y) {

[2.27628]

[2.27738]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T, dnup>
  operator+(const vec3<T, dnup> &x, const vec3<T, dnup> &y) {

Replacement in Z4c/src/tensor.hxx at line 146 [3.48124]

B:BD[2.27772] → [2.27772:27882]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T>
  operator-(const vec3<T> &x, const vec3<T> &y) {

[2.27772]

[2.27882]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T, dnup>
  operator-(const vec3<T, dnup> &x, const vec3<T, dnup> &y) {

Replacement in Z4c/src/tensor.hxx at line 150 [3.48124]

B:BD[2.27916] → [2.27916:28020]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T>
  operator*(const T &a, const vec3<T> &x) {

[2.27916]

[2.28020]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T, dnup>
  operator*(const T &a, const vec3<T, dnup> &x) {

Replacement in Z4c/src/tensor.hxx at line 156 [3.48124]

B:BD[2.28107] → [2.28107:28158]

  operator==(const vec3<T> &x, const vec3<T> &y) {

[2.28107]

[2.28158]

  operator==(const vec3<T, dnup> &x, const vec3<T, dnup> &y) {

Replacement in Z4c/src/tensor.hxx at line 160 [3.48124]

B:BD[2.28271] → [2.28271:28322]

  operator!=(const vec3<T> &x, const vec3<T> &y) {

[2.28271]

[2.28322]

  operator!=(const vec3<T, dnup> &x, const vec3<T, dnup> &y) {

Replacement in Z4c/src/tensor.hxx at line 170 [3.48124]
∅:D[2.28474] → [3.49963:50025]
B:BD[3.49963] → [3.49963:50025]
```
  friend ostream &operator<<(ostream &os, const vec3<T> &v) {
```
[2.28474]
[3.50025]
```
  friend ostream &operator<<(ostream &os, const vec3<T, dnup> &v) {
```

Replacement in Z4c/src/tensor.hxx at line 175 [3.48124]

∅:D[2.28476] → [3.50099:50144]

B:BD[3.50099] → [3.50099:50144]

B:BD[3.50144] → [2.28477:28573]

template <typename T> struct nan<vec3<T> > {
  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T> operator()() const {
    return vec3<T>();

[2.28476]

[2.28573]

template <typename T, dnup_t dnup> struct nan<vec3<T, dnup> > {
  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr vec3<T, dnup> operator()() const {
    return vec3<T, dnup>();

Replacement in Z4c/src/tensor.hxx at line 181 [3.48124]

B:BD[2.28581] → [2.28581:28628]

template <typename T> struct norm1<vec3<T> > {

[2.28581]

[2.28628]

template <typename T, dnup_t dnup> struct norm1<vec3<T, dnup> > {

Replacement in Z4c/src/tensor.hxx at line 184 [3.48124]
B:BD[2.28749] → [2.28749:28788]
```
  operator()(const vec3<T> &x) const {
```
[2.28749]
[2.28788]
```
  operator()(const vec3<T, dnup> &x) const {
```

Replacement in Z4c/src/tensor.hxx at line 192 [3.48124]

B:BD[3.50313] → [3.50313:50348]

template <typename T> class mat3 {

[3.50313]

[3.50348]

template <typename T, dnup_t dnup1, dnup_t dnup2> class mat3 {
  static_assert(dnup1 == dnup2, "");

Replacement in Z4c/src/tensor.hxx at line 298 [3.48124]

B:BD[3.52253] → [3.52253:52457]

    assert(!CCTK_isnan(A(0, 0)));
    assert(!CCTK_isnan(A(0, 1)));
    assert(!CCTK_isnan(A(0, 2)));
    assert(!CCTK_isnan(A(1, 1)));
    assert(!CCTK_isnan(A(1, 2)));
    assert(!CCTK_isnan(A(2, 2)));

[3.52253]

[3.52457]

    assert(CCTK_isfinite(A(0, 0)));
    assert(CCTK_isfinite(A(0, 1)));
    assert(CCTK_isfinite(A(0, 2)));
    assert(CCTK_isfinite(A(1, 1)));
    assert(CCTK_isfinite(A(1, 2)));
    assert(CCTK_isfinite(A(2, 2)));

Replacement in Z4c/src/tensor.hxx at line 322 [3.48124]

B:BD[2.32229] → [2.32229:32321]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T>
  operator+(const mat3<T> &x) {

[2.32229]

[2.32321]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, dnup1, dnup2>
  operator+(const mat3<T, dnup1, dnup2> &x) {

Replacement in Z4c/src/tensor.hxx at line 326 [3.48124]

B:BD[2.32347] → [2.32347:32439]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T>
  operator-(const mat3<T> &x) {

[2.32347]

[2.32439]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, dnup1, dnup2>
  operator-(const mat3<T, dnup1, dnup2> &x) {

Replacement in Z4c/src/tensor.hxx at line 330 [3.48124]

B:BD[2.32465] → [2.32465:32575]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T>
  operator+(const mat3<T> &x, const mat3<T> &y) {

[2.32465]

[2.32575]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, dnup1, dnup2>
  operator+(const mat3<T, dnup1, dnup2> &x, const mat3<T, dnup1, dnup2> &y) {

Replacement in Z4c/src/tensor.hxx at line 334 [3.48124]

B:BD[2.32609] → [2.32609:32719]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T>
  operator-(const mat3<T> &x, const mat3<T> &y) {

[2.32609]

[2.32719]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, dnup1, dnup2>
  operator-(const mat3<T, dnup1, dnup2> &x, const mat3<T, dnup1, dnup2> &y) {

Replacement in Z4c/src/tensor.hxx at line 338 [3.48124]

B:BD[2.32753] → [2.32753:32857]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T>
  operator*(const T &a, const mat3<T> &x) {

[2.32753]

[2.32857]

  friend /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, dnup1, dnup2>
  operator*(const T &a, const mat3<T, dnup1, dnup2> &x) {

Replacement in Z4c/src/tensor.hxx at line 344 [3.48124]

B:BD[2.32944] → [2.32944:32995]

  operator==(const mat3<T> &x, const mat3<T> &y) {

[2.32944]

[2.32995]

  operator==(const mat3<T, dnup1, dnup2> &x, const mat3<T, dnup1, dnup2> &y) {

Replacement in Z4c/src/tensor.hxx at line 348 [3.48124]

B:BD[2.33108] → [2.33108:33159]

  operator!=(const mat3<T> &x, const mat3<T> &y) {

[2.33108]

[2.33159]

  operator!=(const mat3<T, dnup1, dnup2> &x, const mat3<T, dnup1, dnup2> &y) {

Replacement in Z4c/src/tensor.hxx at line 365 [3.48124]

B:BD[3.53078] → [2.33373:33449]

  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3 inv(const T detA) const {

[3.53078]

[3.53121]

  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, !dnup1, !dnup2>
  inv(const T detA) const {

Replacement in Z4c/src/tensor.hxx at line 369 [3.48124]

B:BD[3.53178] → [2.33450:33775]

∅:D[2.33775] → [3.53503:53569]

B:BD[3.53503] → [3.53503:53569]

    return mat3{detA1 * (A(1, 1) * A(2, 2) - A(1, 2) * A(2, 1)),
                detA1 * (A(1, 2) * A(2, 0) - A(1, 0) * A(2, 2)),
                detA1 * (A(1, 0) * A(2, 1) - A(1, 1) * A(2, 0)),
                detA1 * (A(2, 2) * A(0, 0) - A(2, 0) * A(0, 2)),
                detA1 * (A(2, 0) * A(0, 1) - A(2, 1) * A(0, 0)),
                detA1 * (A(0, 0) * A(1, 1) - A(0, 1) * A(1, 0))};

[3.53178]

[3.53569]

    return mat3<T, !dnup1, !dnup2>{
        detA1 * (A(1, 1) * A(2, 2) - A(1, 2) * A(2, 1)),
        detA1 * (A(1, 2) * A(2, 0) - A(1, 0) * A(2, 2)),
        detA1 * (A(1, 0) * A(2, 1) - A(1, 1) * A(2, 0)),
        detA1 * (A(2, 2) * A(0, 0) - A(2, 0) * A(0, 2)),
        detA1 * (A(2, 0) * A(0, 1) - A(2, 1) * A(0, 0)),
        detA1 * (A(0, 0) * A(1, 1) - A(0, 1) * A(1, 0))};

Replacement in Z4c/src/tensor.hxx at line 378 [3.48124]

B:BD[3.53574] → [2.33776:33853]

  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr T trace(const mat3 &gu) const {

[3.53574]

[2.33853]

  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr T
  trace(const mat3<T, !dnup1, !dnup2> &gu) const {

Replacement in Z4c/src/tensor.hxx at line 384 [3.48124]

B:BD[2.33952] → [2.33952:34054]

  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3
  trace_free(const mat3 &g, const mat3 &gu) const {

[2.33952]

[3.53604]

  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3 trace_free(
      const mat3<T, dnup1, dnup2> &g, const mat3<T, !dnup1, !dnup2> &gu) const {

Replacement in Z4c/src/tensor.hxx at line 391 [3.48124]

B:BD[3.53676] → [3.53676:53738]

  friend ostream &operator<<(ostream &os, const mat3<T> &A) {

[3.53676]

[3.53738]

  friend ostream &operator<<(ostream &os, const mat3<T, dnup1, dnup2> &A) {

Replacement in Z4c/src/tensor.hxx at line 397 [3.48124]

B:BD[3.53963] → [3.53963:54008]

B:BD[3.54008] → [2.34163:34259]

template <typename T> struct nan<mat3<T> > {
  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T> operator()() const {
    return mat3<T>();

[3.53963]

[2.34259]

template <typename T, dnup_t dnup1, dnup_t dnup2>
struct nan<mat3<T, dnup1, dnup2> > {
  /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, dnup1, dnup2>
  operator()() const {
    return mat3<T, dnup1, dnup2>();

Replacement in Z4c/src/tensor.hxx at line 405 [3.48124]

B:BD[2.34267] → [2.34267:34314]

template <typename T> struct norm1<mat3<T> > {

[2.34267]

[2.34314]

template <typename T, dnup_t dnup1, dnup_t dnup2>
struct norm1<mat3<T, dnup1, dnup2> > {

Replacement in Z4c/src/tensor.hxx at line 409 [3.48124]
B:BD[2.34435] → [2.34435:34474]
```
  operator()(const mat3<T> &x) const {
```
[2.34435]
[2.34474]
```
  operator()(const mat3<T, dnup1, dnup2> &x) const {
```

Replacement in Z4c/src/tensor.hxx at line 414 [3.48124]

B:BD[3.54073] → [3.54073:54095]

B:BD[3.54095] → [2.34520:34668]

template <typename T>
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T> mul(const mat3<T> &A,
                                                       const mat3<T> &B) {

[3.54073]

[3.54168]

template <typename T, dnup_t dnup1, dnup_t dnup2, dnup_t dnup3>
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr mat3<T, dnup1, dnup2>
mul(const mat3<T, dnup1, dnup3> &A, const mat3<T, !dnup3, dnup2> &B) {

Replacement in Z4c/src/tensor.hxx at line 418 [3.48124]
B:BD[3.54196] → [2.34669:34706]
```
  return mat3<T>([&](int a, int b) {
```
[3.54196]
[2.34706]
```
  return mat3<T, dnup1, dnup2>([&](int a, int b) {
```

Replacement in Z4c/src/tensor.hxx at line 423 [3.48124]

B:BD[2.34774] → [2.34774:34859]

template <typename T>
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr T add(const T &x) {

[2.34774]

[2.34859]

template <typename F, typename T>
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr T fold(const F &f, const T &x) {

Insertion in Z4c/src/tensor.hxx at line 426 [3.48124]

[2.34871]

}
template <typename F, typename T, typename... Ts>
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr T fold(const F &f, const T &x0,
                                                  const T &x1,
                                                  const Ts &... xs) {
  return fold(f, fold(f, x0, x1), xs...);

Insertion in Z4c/src/tensor.hxx at line 433 [3.48124]

[2.34873]


template <typename T>
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr T add() {
  return T(0);
}
// template <typename T>
// /*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ constexpr T add(const T &x) {
//   return x;
// }

Replacement in Z4c/src/tensor.hxx at line 482 [3.48124]

B:BD[2.35879] → [2.35879:35956]

                                         const vec3<T> &dx, const int dir) {

[2.35879]

[3.54597]

                                         const vec3<T, UP> &dx, const int dir) {

Replacement in Z4c/src/tensor.hxx at line 490 [3.48124]

B:BD[2.36065] → [2.36065:36131]

             const bool sign, const vec3<T> &dx, const int dir) {

[2.36065]

[2.36131]

             const bool sign, const vec3<T, UP> &dx, const int dir) {

Replacement in Z4c/src/tensor.hxx at line 514 [3.48124]

B:BD[2.36893] → [2.36893:36953]

       const vec3<T> &dx, const int dir1, const int dir2) {

[2.36893]

[3.55126]

       const vec3<T, UP> &dx, const int dir1, const int dir2) {

Replacement in Z4c/src/tensor.hxx at line 518 [3.48124]
∅:D[2.37022] → [3.55287:55321]
B:BD[3.55287] → [3.55287:55321]
```
           (dx(dir1) * dx(dir2));
```
[2.37022]
[3.55321]
```
           pow2(dx(dir1));
```

Replacement in Z4c/src/tensor.hxx at line 520 [3.48124]

B:BD[3.55328] → [3.55328:55454]

    return (deriv(gf_, I + DI[dir2], dx, dir1) -
            deriv(gf_, I - DI[dir2], dx, dir1)) /
           (2 * dx(dir2));

[3.55328]

[3.55454]

    // return (deriv(gf_, I + DI[dir2], dx, dir1) -
    //         deriv(gf_, I - DI[dir2], dx, dir1)) /
    //        (2 * dx(dir2));
    return ((gf_(I + DI[dir1] + DI[dir2])       //
             - gf_(I - DI[dir1] + DI[dir2]))    //
            - (gf_(I + DI[dir1] - DI[dir2])     //
               - gf_(I - DI[dir1] - DI[dir2]))) //
           / (4 * dx(dir1) * dx(dir2));

Replacement in Z4c/src/tensor.hxx at line 531 [3.48124]

B:BD[3.55479] → [2.37023:37245]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ T deriv4(const GF3D<const T, 0, 0, 0> &gf_,
                                          const vect<int, dim> &I,
                                          const vec3<T> &dx, const int dir) {

[3.55479]

[2.37245]

/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ T
deriv4(const GF3D<const T, 0, 0, 0> &gf_, const vect<int, dim> &I,
       const vec3<T, UP> &dx, const int dir) {

Replacement in Z4c/src/tensor.hxx at line 542 [3.48124]
B:BD[2.37513] → [2.37513:37554]
```
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ vec3<T>
```
[2.37513]
[2.37554]
```
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ vec3<T, DN>
```
Replacement in Z4c/src/tensor.hxx at line 544 [3.48124]
B:BD[2.37620] → [2.37620:37647]
```
      const vec3<T> &dx) {
```
[2.37620]
[2.37647]
```
      const vec3<T, UP> &dx) {
```
Replacement in Z4c/src/tensor.hxx at line 553 [3.48124]
B:BD[2.37772] → [2.37772:37813]
```
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ vec3<T>
```
[2.37772]
[2.37813]
```
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ vec3<T, DN>
```

Replacement in Z4c/src/tensor.hxx at line 555 [3.48124]

B:BD[2.37886] → [2.37886:37940]

             const vec3<T> &dir, const vec3<T> &dx) {

[2.37886]

[2.37940]

             const vec3<T, UP> &dir, const vec3<T, UP> &dx) {

Replacement in Z4c/src/tensor.hxx at line 564 [3.48124]
B:BD[2.38137] → [2.38137:38178]
```
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ mat3<T>
```
[2.38137]
[2.38178]
```
/*CCTK_ATTRIBUTE_ALWAYS_INLINE*/ mat3<T, DN, DN>
```
Replacement in Z4c/src/tensor.hxx at line 566 [3.48124]
B:BD[2.38245] → [2.38245:38273]
```
       const vec3<T> &dx) {
```
[2.38245]
[3.55590]
```
       const vec3<T, UP> &dx) {
```

Replacement in Z4c/src/tensor.hxx at line 577 [3.48124]

B:BD[2.38436] → [2.38436:38497]

                                        const vec3<T> &dx) {

[2.38436]

[2.38497]

                                        const vec3<T, UP> &dx) {

Deletion in Z4c/src/test.cxx at line 1 [2.38726]
B:BD[2.38726] → [2.38727:38764]
```
#warning "TODO"
#include <iostream>
```
Replacement in Z4c/src/test.cxx at line 27 [2.38726]
B:BD[2.39349] → [2.39349:39374]
```
    return mat3<double>(
```
[2.39349]
[2.39374]
```
    return mat3<double, DN, DN>(
```

Replacement in Z4c/src/test.cxx at line 31 [2.38726]

B:BD[2.39447] → [2.39447:39582]

  const mat3<double> Z([&](int a, int b) { return double(0); });
  const mat3<double> I([&](int a, int b) { return double(a == b); });

[2.39447]

[2.39582]

  const mat3<double, DN, DN> Z([&](int a, int b) { return double(0); });
  const mat3<double, DN, DN> I([&](int a, int b) { return double(a == b); });

Insertion in Z4c/src/test.cxx at line 34 [2.38726]

[2.39600]

  const mat3<double, UP, UP> Zup([&](int a, int b) { return double(0); });
  const mat3<double, UP, UP> Iup([&](int a, int b) { return double(a == b); });
  assert(Iup != Zup);

Replacement in Z4c/src/test.cxx at line 39 [2.38726]

B:BD[2.39635] → [2.39635:39755]

    const mat3<double> A = randmat10();
    const mat3<double> B = randmat10();
    const mat3<double> C = randmat10();

[2.39635]

[2.39755]

    const mat3<double, DN, DN> A = randmat10();
    const mat3<double, DN, DN> B = randmat10();
    const mat3<double, DN, DN> C = randmat10();

Replacement in Z4c/src/test.cxx at line 62 [2.38726]

B:BD[2.40321] → [2.40321:40433]

    assert(mul(I, A) == A);
    assert(mul(A, I) == A);
    assert(mul(Z, A) == Z);
    assert(mul(A, Z) == Z);

[2.40321]

[2.40433]

    // DNUP  assert(mul(I, A) == A);
    // DNUP  assert(mul(A, I) == A);
    // DNUP  assert(mul(Z, A) == Z);
    // DNUP  assert(mul(A, Z) == Z);

Replacement in Z4c/src/test.cxx at line 71 [2.38726]

B:BD[2.40535] → [2.40535:40589]

    assert(Z.inv(1) == Z);
    assert(I.inv(1) == I);

[2.40535]

[2.40589]

    assert(Z.inv(1) == Zup);
    assert(I.inv(1) == Iup);

Replacement in Z4c/src/test.cxx at line 74 [2.38726]

B:BD[2.40590] → [2.40590:40680]

    assert(mul(A.inv(1), A) == A.det() * I);
    assert(mul(A, A.inv(1)) == A.det() * I);

[2.40590]

[2.40680]

    // DNUP assert(mul(A.inv(1), A) == A.det() * Iup);
    // DNUP assert(mul(A, A.inv(1)) == A.det() * Iup);

Replacement in Z4c/src/z4c_vars.hxx at line 27 [2.40773]

B:BD[2.41225] → [2.41225:41501]

  const T chi;          // W = -2/3
  const mat3<T> gammat; // W = -2/3
  const T Kh;           // W = 0
  const mat3<T> At;     // W = -2/3
  const vec3<T> Gamt;   // W = +2/3
  const T Theta;        // W = 0
  const T alphaG;       // W = 0
  const vec3<T> betaG;  // W = 0

[2.41225]

[2.41501]

  const T chi;                  // W = -2/3
  const mat3<T, DN, DN> gammat; // W = -2/3
  const T Kh;                   // W = 0
  const mat3<T, DN, DN> At;     // W = -2/3
  const vec3<T, UP> Gamt;       // W = +2/3
  const T Theta;                // W = 0
  const T alphaG;               // W = 0
  const vec3<T, UP> betaG;      // W = 0

Replacement in Z4c/src/z4c_vars.hxx at line 38 [2.40773]
B:BD[2.41538] → [2.41538:41579]
```
  const vec3<T> Si;
  const mat3<T> Sij;
```
[2.41538]
[2.41579]
```
  const vec3<T, DN> Si;
  const mat3<T, DN, DN> Sij;
```

Replacement in Z4c/src/z4c_vars.hxx at line 42 [2.40773]

B:BD[2.41599] → [2.41599:41723]

  const mat3<T> g;    // W = 0
  const mat3<T> k;    // W = 0
  const T alp;        // W = 0
  const vec3<T> beta; // W = 0

[2.41599]

[2.41723]

  const mat3<T, DN, DN> g; // W = 0
  const mat3<T, DN, DN> k; // W = 0
  const T alp;             // W = 0
  const vec3<T, UP> beta;  // W = 0

Replacement in Z4c/src/z4c_vars.hxx at line 51 [2.40773]

B:BD[2.41904] → [2.41904:42109]

                    const T &chi, const mat3<T> &gammat, const T &Kh,
                    const mat3<T> &At, const vec3<T> &Gamt, const T &Theta,
                    const T &alphaG, const vec3<T> &betaG,

[2.41904]

[2.42109]

                    const T &chi, const mat3<T, DN, DN> &gammat, const T &Kh,
                    const mat3<T, DN, DN> &At, const vec3<T, UP> &Gamt,
                    const T &Theta, const T &alphaG, const vec3<T, UP> &betaG,

Replacement in Z4c/src/z4c_vars.hxx at line 55 [2.40773]

B:BD[2.42132] → [2.42132:42205]

                    const T &rho, const vec3<T> &Si, const mat3<T> &Sij)

[2.42132]

[2.42205]

                    const T &rho, const vec3<T, DN> &Si,
                    const mat3<T, DN, DN> &Sij)

Replacement in Z4c/src/z4c_vars.hxx at line 113 [2.40773]

B:BD[2.44704] → [2.44704:44903]

                             mat3<T>(gf_gammatxx_, gf_gammatxy_, gf_gammatxz_,
                                     gf_gammatyy_, gf_gammatyz_, gf_gammatzz_,
                                     I),

[2.44704]

[2.44903]

                             mat3<T, DN, DN>(gf_gammatxx_, gf_gammatxy_,
                                             gf_gammatxz_, gf_gammatyy_,
                                             gf_gammatyz_, gf_gammatzz_, I),

Replacement in Z4c/src/z4c_vars.hxx at line 119 [2.40773]

B:BD[2.45007] → [2.45007:45145]

                             mat3<T>(gf_Atxx_, gf_Atxy_, gf_Atxz_, gf_Atyy_,
                                     gf_Atyz_, gf_Atzz_, I),

[2.45007]

[2.45145]

                             mat3<T, DN, DN>(gf_Atxx_, gf_Atxy_, gf_Atxz_,
                                             gf_Atyy_, gf_Atyz_, gf_Atzz_, I),

Replacement in Z4c/src/z4c_vars.hxx at line 122 [2.40773]

B:BD[2.45177] → [2.45177:45251]

                             vec3<T>(gf_Gamtx_, gf_Gamty_, gf_Gamtz_, I),

[2.45177]

[2.45251]

                             vec3<T, UP>(gf_Gamtx_, gf_Gamty_, gf_Gamtz_, I),

Replacement in Z4c/src/z4c_vars.hxx at line 128 [2.40773]

B:BD[2.45434] → [2.45434:45511]

                             vec3<T>(gf_betaGx_, gf_betaGy_, gf_betaGz_, I),

[2.45434]

[2.45511]

                             vec3<T, UP>(gf_betaGx_, gf_betaGy_, gf_betaGz_, I),

Replacement in Z4c/src/z4c_vars.hxx at line 132 [2.40773]

B:BD[2.45607] → [2.45607:45654]

                             vec3<T>{0, 0, 0},

[2.45607]

[2.45654]

                             vec3<T, DN>{0, 0, 0},

Replacement in Z4c/src/z4c_vars.hxx at line 134 [2.40773]

B:BD[2.45686] → [2.45686:45742]

                             mat3<T>{0, 0, 0, 0, 0, 0})

[2.45686]

[2.45742]

                             mat3<T, DN, DN>{0, 0, 0, 0, 0, 0})

Replacement in Z4c/src/z4c_vars.hxx at line 172 [2.40773]

B:BD[2.46672] → [2.46672:46999]

  const vec3<T> dchi;
  const mat3<T> ddchi;
  const mat3<vec3<T> > dgammat;
  const mat3<mat3<T> > ddgammat;
  const vec3<T> dKh;
  const mat3<vec3<T> > dAt;
  const vec3<vec3<T> > dGamt;
  const vec3<T> dTheta;
  const vec3<T> dalphaG;
  const mat3<T> ddalphaG;
  const vec3<vec3<T> > dbetaG;
  const vec3<mat3<T> > ddbetaG;

[2.46672]

[2.46999]

  const vec3<T, DN> dchi;
  const mat3<T, DN, DN> ddchi;
  const mat3<vec3<T, DN>, DN, DN> dgammat;
  const mat3<mat3<T, DN, DN>, DN, DN> ddgammat;
  const vec3<T, DN> dKh;
  const mat3<vec3<T, DN>, DN, DN> dAt;
  const vec3<vec3<T, DN>, UP> dGamt;
  const vec3<T, DN> dTheta;
  const vec3<T, DN> dalphaG;
  const mat3<T, DN, DN> ddalphaG;
  const vec3<vec3<T, DN>, UP> dbetaG;
  const vec3<mat3<T, DN, DN>, UP> ddbetaG;

Replacement in Z4c/src/z4c_vars.hxx at line 186 [2.40773]

B:BD[2.47028] → [2.47028:47390]

  const mat3<T> gammatu;
  const mat3<vec3<T> > dgammatu;
  const vec3<mat3<T> > Gammatl;
  const vec3<mat3<T> > Gammat;
  const vec3<T> Gamtd;
  const mat3<T> DDchi;
  const mat3<T> gu;
  const mat3<vec3<T> > dg;
  const vec3<mat3<T> > Gammal;
  const vec3<mat3<T> > Gamma;
  const mat3<T> DDalphaG;
  const mat3<T> Rchi;
  const mat3<T> Rt;
  const mat3<T> R;

[2.47028]

[2.47390]

  const mat3<T, UP, UP> gammatu;
  const mat3<vec3<T, DN>, UP, UP> dgammatu;
  const vec3<mat3<T, DN, DN>, DN> Gammatl;
  const vec3<mat3<T, DN, DN>, UP> Gammat;
  const vec3<T, UP> Gamtd;
  const mat3<T, DN, DN> DDchi;
  const mat3<T, UP, UP> gu;
  const mat3<vec3<T, DN>, DN, DN> dg;
  const vec3<mat3<T, DN, DN>, DN> Gammal;
  const vec3<mat3<T, DN, DN>, UP> Gamma;
  const mat3<T, DN, DN> DDalphaG;
  const mat3<T, DN, DN> Rchi;
  const mat3<T, DN, DN> Rt;
  const mat3<T, DN, DN> R;

Replacement in Z4c/src/z4c_vars.hxx at line 201 [2.40773]

B:BD[2.47405] → [2.47405:47455]

  const mat3<T> Atu;
  const mat3<vec3<T> > dAtu;

[2.47405]

[2.47455]

  const mat3<T, UP, UP> Atu;
  const mat3<vec3<T, DN>, UP, UP> dAtu;

Replacement in Z4c/src/z4c_vars.hxx at line 206 [2.40773]
B:BD[2.47493] → [2.47493:47514]
```
  const vec3<T> ZtC;
```
[2.47493]
[2.47514]
```
  const vec3<T, UP> ZtC;
```
Replacement in Z4c/src/z4c_vars.hxx at line 208 [2.40773]
B:BD[2.47528] → [2.47528:47549]
```
  const vec3<T> MtC;
```
[2.47528]
[2.47549]
```
  const vec3<T, UP> MtC;
```
Replacement in Z4c/src/z4c_vars.hxx at line 213 [2.40773]
B:BD[2.47604] → [2.47604:47632]
```
  const mat3<T> gammat_rhs;
```
[2.47604]
[2.47632]
```
  const mat3<T, DN, DN> gammat_rhs;
```

Replacement in Z4c/src/z4c_vars.hxx at line 215 [2.40773]

B:BD[2.47650] → [2.47650:47700]

  const mat3<T> At_rhs;
  const vec3<T> Gamt_rhs;

[2.47650]

[2.47700]

  const mat3<T, DN, DN> At_rhs;
  const vec3<T, UP> Gamt_rhs;

Replacement in Z4c/src/z4c_vars.hxx at line 219 [2.40773]
B:BD[2.47743] → [2.47743:47770]
```
  const vec3<T> betaG_rhs;
```
[2.47743]
[2.47770]
```
  const vec3<T, UP> betaG_rhs;
```

Replacement in Z4c/src/z4c_vars.hxx at line 226 [2.40773]

B:BD[2.47957] → [2.47957:48678]

           const T &chi, const vec3<T> &dchi, const mat3<T> &ddchi, //
           const mat3<T> &gammat, const mat3<vec3<T> > &dgammat,
           const mat3<mat3<T> > &ddgammat,                                   //
           const T &Kh, const vec3<T> &dKh,                                  //
           const mat3<T> &At, const mat3<vec3<T> > &dAt,                     //
           const vec3<T> &Gamt, const vec3<vec3<T> > &dGamt,                 //
           const T &Theta, const vec3<T> &dTheta,                            //
           const T &alphaG, const vec3<T> &dalphaG, const mat3<T> &ddalphaG, //
           const vec3<T> &betaG, const vec3<vec3<T> > &dbetaG,
           const vec3<mat3<T> > &ddbetaG,

[2.47957]

[2.48678]

           const T &chi, const vec3<T, DN> &dchi,
           const mat3<T, DN, DN> &ddchi, //
           const mat3<T, DN, DN> &gammat,
           const mat3<vec3<T, DN>, DN, DN> &dgammat,
           const mat3<mat3<T, DN, DN>, DN, DN> &ddgammat,                   //
           const T &Kh, const vec3<T, DN> &dKh,                             //
           const mat3<T, DN, DN> &At, const mat3<vec3<T, DN>, DN, DN> &dAt, //
           const vec3<T, UP> &Gamt, const vec3<vec3<T, DN>, UP> &dGamt,     //
           const T &Theta, const vec3<T, DN> &dTheta,                       //
           const T &alphaG, const vec3<T, DN> &dalphaG,
           const mat3<T, DN, DN> &ddalphaG, //
           const vec3<T, UP> &betaG, const vec3<vec3<T, DN>, UP> &dbetaG,
           const vec3<mat3<T, DN, DN>, UP> &ddbetaG,

Replacement in Z4c/src/z4c_vars.hxx at line 240 [2.40773]

B:BD[2.48692] → [2.48692:48756]

           const T &rho, const vec3<T> &Si, const mat3<T> &Sij)

[2.48692]

[2.48756]

           const T &rho, const vec3<T, DN> &Si, const mat3<T, DN, DN> &Sij)

Replacement in Z4c/src/z4c_vars.hxx at line 268 [2.40773]
B:BD[2.50049] → [2.50049:50087]
```
          return vec3<T>([&](int c) {
```
[2.50049]
[2.50087]
```
          return vec3<T, DN>([&](int c) {
```

Insertion in Z4c/src/z4c_vars.hxx at line 278 [2.40773]

[2.50483]

          // return ddalphaG(a, b) //
          //        - sum1([&](int x) { return Gammat(x)(a, b) * dalphaG(x); })
          //        //
          //        // chipsipower = -4
          //        // oochipsipower = -1
          //        // df = oochipsipower dchi(a) / chi
          //        // ddf == oochipsipower ddchi / chi - (-4) df df
          //        + 1 / (2 * chi) *
          //              (dchi(a) * dalphaG(b)    //
          //               + dchi(b) * dalphaG(a)) //
          //        + 1 / (4 * chi) * sum2([&](int x, int y) {
          //            return gammat(a, b) * gammatu(x, y) * dchi(x) *
          //            dalphaG(y);
          //          });

Replacement in Z4c/src/z4c_vars.hxx at line 393 [2.40773]

B:BD[2.54846] → [2.54846:54927]

        Kh_rhs(-sum2([&](int x, int y) { return gu(x, y) * DDalphaG(x, y); }) //

[2.54846]

[2.54927]

        Kh_rhs(sum2([&](int x, int y) { return -gu(x, y) * DDalphaG(x, y); }) //

Replacement in Z4c/src/z4c_vars.hxx at line 425 [2.40773]

B:BD[2.56496] → [2.56496:56573]

          return -sum1([&](int x) { return 2 * Atu(a, x) * dalphaG(x); }) //

[2.56496]

[2.56573]

          return sum1([&](int x) { return -2 * Atu(a, x) * dalphaG(x); }) //

Replacement in Z4c/src/z4c_vars.hxx at line 508 [2.40773]

B:BD[2.60044] → [2.60044:60097]

           const vect<int, 3> &I, const vec3<T> &dx)

[2.60044]

[2.60097]

           const vect<int, 3> &I, const vec3<T, UP> &dx)

Replacement in Z4c/src/z4c_vars.hxx at line 513 [2.40773]

B:BD[2.60267] → [2.60267:60438]

                 mat3<T>(gf_gammatxx_, gf_gammatxy_, gf_gammatxz_, gf_gammatyy_,
                         gf_gammatyz_, gf_gammatzz_, I),
                 mat3<vec3<T> >{

[2.60267]

[2.60438]

                 mat3<T, DN, DN>(gf_gammatxx_, gf_gammatxy_, gf_gammatxz_,
                                 gf_gammatyy_, gf_gammatyz_, gf_gammatzz_, I),
                 mat3<vec3<T, DN>, DN, DN>{

Replacement in Z4c/src/z4c_vars.hxx at line 523 [2.40773]

B:BD[2.60752] → [2.60752:60785]

                 mat3<mat3<T> >{

[2.60752]

[2.60785]

                 mat3<mat3<T, DN, DN>, DN, DN>{

Replacement in Z4c/src/z4c_vars.hxx at line 534 [2.40773]

B:BD[2.61195] → [2.61195:61342]

                 mat3<T>(gf_Atxx_, gf_Atxy_, gf_Atxz_, gf_Atyy_, gf_Atyz_,
                         gf_Atzz_, I),
                 mat3<vec3<T> >{

[2.61195]

[2.61342]

                 mat3<T, DN, DN>(gf_Atxx_, gf_Atxy_, gf_Atxz_, gf_Atyy_,
                                 gf_Atyz_, gf_Atzz_, I),
                 mat3<vec3<T, DN>, DN, DN>{

Replacement in Z4c/src/z4c_vars.hxx at line 545 [2.40773]

B:BD[2.61652] → [2.61652:61747]

                 vec3<T>(gf_Gamtx_, gf_Gamty_, gf_Gamtz_, I),
                 vec3<vec3<T> >{

[2.61652]

[2.61747]

                 vec3<T, UP>(gf_Gamtx_, gf_Gamty_, gf_Gamtz_, I),
                 vec3<vec3<T, DN>, UP>{

Replacement in Z4c/src/z4c_vars.hxx at line 557 [2.40773]

B:BD[2.62123] → [2.62123:62221]

                 vec3<T>(gf_betaGx_, gf_betaGy_, gf_betaGz_, I),
                 vec3<vec3<T> >{

[2.62123]

[2.62221]

                 vec3<T, UP>(gf_betaGx_, gf_betaGy_, gf_betaGz_, I),
                 vec3<vec3<T, DN>, UP>{

Replacement in Z4c/src/z4c_vars.hxx at line 563 [2.40773]

B:BD[2.62382] → [2.62382:62415]

                 vec3<mat3<T> >{

[2.62382]

[2.62415]

                 vec3<mat3<T, DN, DN>, UP>{

Replacement in Z4c/src/z4c_vars.hxx at line 571 [2.40773]
B:BD[2.62639] → [2.62639:62674]
```
                 vec3<T>{0, 0, 0},
```
[2.62639]
[2.62674]
```
                 vec3<T, DN>{0, 0, 0},
```

Replacement in Z4c/src/z4c_vars.hxx at line 573 [2.40773]

B:BD[2.62694] → [2.62694:62738]

                 mat3<T>{0, 0, 0, 0, 0, 0})

[2.62694]

[2.62738]

                 mat3<T, DN, DN>{0, 0, 0, 0, 0, 0})