eschnett/cactusamrex - Change 67ABKXMOSVKKVRSSD6QVF4BLBQLZCEAST2UMLLKOBV7WSVTQRSEQC

Z4c: Some source reformatting

Created by Erik Schnetter on April 3, 2020

67ABKXMOSVKKVRSSD6QVF4BLBQLZCEAST2UMLLKOBV7WSVTQRSEQC

Dependencies

In channels

master

main

Change contents

Replacement in Z4c/src/enforce.cxx at line 40 [5.22673]

∅:D[3.2844] → [4.72937:72953]

B:BD[4.72937] → [4.72937:72953]

B:BD[4.72953] → [2.455:557]

        // Load
        const CCTK_REAL chi_old = gf_chi_(p.I);
        const CCTK_REAL alphaG_old = gf_alphaG_(p.I);

[3.2844]

[2.557]

    // Load
    const CCTK_REAL chi_old = gf_chi_(p.I);
    const CCTK_REAL alphaG_old = gf_alphaG_(p.I);

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

∅:D[2.558] → [4.72953:73273]

B:BD[4.72953] → [4.72953:73273]

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

[2.558]

[2.559]

    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 50 [5.22673]
B:BD[2.560] → [2.560:586]
```
        // Enforce floors
```
[2.560]
[5.24011]
```
    // Enforce floors
```

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

B:BD[5.24012] → [2.587:708]

        const CCTK_REAL chi = fmax(chi_floor, chi_old);
        const CCTK_REAL alphaG = fmax(alphaG_floor, alphaG_old);

[5.24012]

[2.708]

    const CCTK_REAL chi = fmax(chi_floor, chi_old);
    const CCTK_REAL alphaG = fmax(alphaG_floor, alphaG_old);

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

∅:D[2.709] → [4.73274:73355]

B:BD[5.24012] → [4.73274:73355]

        // Enforce algebraic constraints
        // See arXiv:1212.2901 [gr-qc].

[2.709]

[5.24085]

    // Enforce algebraic constraints
    // See arXiv:1212.2901 [gr-qc].

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

B:BD[5.24086] → [4.73356:73414]

B:BD[4.73414] → [2.710:770]

∅:D[2.770] → [4.73473:73519]

B:BD[4.73473] → [4.73473:73519]

B:BD[4.73519] → [2.771:843]

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

[5.24086]

[6.9445]

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

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

B:BD[6.9463] → [4.73591:74069]

        const CCTK_REAL detgammat = gammat.det();
        const CCTK_REAL gammat_norm = gammat.maxabs();
        const CCTK_REAL gammat_scale = gammat_norm;
        if (!(fabs(detgammat - 1) <= 1.0e-12 * gammat_scale)) {
          ostringstream buf;
          buf << "det gammat is not one: gammat=" << gammat
              << " det(gammat)=" << detgammat;
          CCTK_VERROR("%s", buf.str().c_str());
        }
        assert(fabs(detgammat - 1) <= 1.0e-12 * gammat_scale);

[6.9463]

[6.9917]

    const CCTK_REAL detgammat = gammat.det();
    const CCTK_REAL gammat_norm = gammat.maxabs();
    const CCTK_REAL gammat_scale = gammat_norm;
    if (!(fabs(detgammat - 1) <= 1.0e-12 * gammat_scale)) {
      ostringstream buf;
      buf << "det gammat is not one: gammat=" << gammat
          << " det(gammat)=" << detgammat;
      CCTK_VERROR("%s", buf.str().c_str());
    }
    assert(fabs(detgammat - 1) <= 1.0e-12 * gammat_scale);

Replacement in Z4c/src/enforce.cxx at line 75 [5.22673]

B:BD[5.24279] → [4.74070:74133]

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

[5.24279]

[5.24330]

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

Replacement in Z4c/src/enforce.cxx at line 77 [5.22673]

B:BD[5.24331] → [4.74134:74387]

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

[5.24331]

[6.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 83 [5.22673]

B:BD[6.10092] → [4.74388:74951]

        const CCTK_REAL traceAt =
            sum2([&](int x, int y) { return gammatu(x, y) * At(x, y); });
        const CCTK_REAL gammatu_norm = gammatu.maxabs();
        const CCTK_REAL At_norm = At.maxabs();
        const CCTK_REAL At_scale =
            fmax(fmax(gammat_norm, gammatu_norm), At_norm);
        if (!(fabs(traceAt) <= 1.0e-12 * At_scale)) {
          ostringstream buf;
          buf << "tr At: At=" << At << " tr(At)=" << traceAt;
          CCTK_VERROR("%s", buf.str().c_str());
        }
        assert(fabs(traceAt) <= 1.0e-12 * At_scale);

[6.10092]

[6.10615]

    const CCTK_REAL traceAt =
        sum2([&](int x, int y) { return gammatu(x, y) * At(x, y); });
    const CCTK_REAL gammatu_norm = gammatu.maxabs();
    const CCTK_REAL At_norm = At.maxabs();
    const CCTK_REAL At_scale = fmax(fmax(gammat_norm, gammatu_norm), At_norm);
    if (!(fabs(traceAt) <= 1.0e-12 * At_scale)) {
      ostringstream buf;
      buf << "tr At: At=" << At << " tr(At)=" << traceAt;
      CCTK_VERROR("%s", buf.str().c_str());
    }
    assert(fabs(traceAt) <= 1.0e-12 * At_scale);

Replacement in Z4c/src/enforce.cxx at line 96 [5.22673]

B:BD[5.24589] → [4.74952:74969]

B:BD[4.74969] → [2.844:872]

∅:D[2.872] → [4.74969:75201]

B:BD[4.74969] → [4.74969:75201]

B:BD[4.75201] → [2.873:907]

∅:D[2.907] → [4.75201:75211]

B:BD[4.75201] → [4.75201:75211]

        // Store
        gf_chi_(p.I) = chi;
        gammat.store(gf_gammatxx_, gf_gammatxy_, gf_gammatxz_, gf_gammatyy_,
                     gf_gammatyz_, gf_gammatzz_, p.I);
        At.store(gf_Atxx_, gf_Atxy_, gf_Atxz_, gf_Atyy_, gf_Atyz_, gf_Atzz_,
                 p.I);
        gf_alphaG_(p.I) = alphaG;
      });

[5.24589]

[5.24807]

    // Store
    gf_chi_(p.I) = chi;
    gammat.store(gf_gammatxx_, gf_gammatxy_, gf_gammatxz_, gf_gammatyy_,
                 gf_gammatyz_, gf_gammatzz_, p.I);
    At.store(gf_Atxx_, gf_Atxy_, gf_Atxz_, gf_Atyy_, gf_Atyz_, gf_Atzz_, p.I);
    gf_alphaG_(p.I) = alphaG;
  });

Replacement in Z4c/src/rhs.cxx at line 161 [5.35315]

B:BD[4.81061] → [4.81061:81265]

        vars.gammat_rhs.store(gf_gammatxx_rhs_, gf_gammatxy_rhs_,
                              gf_gammatxz_rhs_, gf_gammatyy_rhs_,
                              gf_gammatyz_rhs_, gf_gammatzz_rhs_, p.I);

[4.81061]

[4.81265]

    vars.gammat_rhs.store(gf_gammatxx_rhs_, gf_gammatxy_rhs_, gf_gammatxz_rhs_,
                          gf_gammatyy_rhs_, gf_gammatyz_rhs_, gf_gammatzz_rhs_,
                          p.I);

Replacement in Z4c/src/rhs.cxx at line 165 [5.35315]

B:BD[4.81304] → [4.81304:81446]

        vars.At_rhs.store(gf_Atxx_rhs_, gf_Atxy_rhs_, gf_Atxz_rhs_,
                          gf_Atyy_rhs_, gf_Atyz_rhs_, gf_Atzz_rhs_, p.I);

[4.81304]

[4.81446]

    vars.At_rhs.store(gf_Atxx_rhs_, gf_Atxy_rhs_, gf_Atxz_rhs_, gf_Atyy_rhs_,
                      gf_Atyz_rhs_, gf_Atzz_rhs_, p.I);

Replacement in Z4c/src/rhs.cxx at line 170 [5.35315]

B:BD[4.81617] → [4.81617:81729]

        vars.betaG_rhs.store(gf_betaGx_rhs_, gf_betaGy_rhs_, gf_betaGz_rhs_,
                             p.I);

[4.81617]

[4.81729]

    vars.betaG_rhs.store(gf_betaGx_rhs_, gf_betaGy_rhs_, gf_betaGz_rhs_, p.I);