eschnett/cactusamrex - Change COCGSXFZAM7XB6JFRIAEAZGS5RQKFGOS7OGBLLUCSW6WT2PVWD6AC

TestProlongate: Update thorn

Created by Erik Schnetter on January 1, 2020

COCGSXFZAM7XB6JFRIAEAZGS5RQKFGOS7OGBLLUCSW6WT2PVWD6AC

Dependencies

In channels

master

main

Change contents

Replacement in TestProlongate/interface.ccl at line 2 [4.4774]
B:BD[4.4823] → [4.4823:4850]
B:BD[4.4850] → [3.16311:16360]
```
implements: TestProlongate
inherits: CarpetX # Grid # grid contains scalars
```
[4.4823]
[4.4896]
```
IMPLEMENTS: TestProlongate
INHERITS: CarpetX
```
Insertion in TestProlongate/interface.ccl at line 8 [4.4774]
[4.4927]
[4.4927]
Replacement in TestProlongate/interface.ccl at line 10 [4.4774]
B:BD[4.4928] → [4.4928:4937]
```
private:
```
[4.4928]
[4.4937]
```
PRIVATE:
```

Replacement in TestProlongate/interface.ccl at line 13 [4.4774]

B:BD[4.4938] → [4.4938:4995]

CCTK_REAL data TYPE=GF TIMELEVELS=1 "test grid function"

[4.4938]

[5.0]

CCTK_REAL data TYPE=gf TIMELEVELS=1 "test grid function"
# CCTK_REAL gf000 TYPE=gf TAGS='index={0 0 0}' "Test grid function"
# CCTK_REAL gf001 TYPE=gf TAGS='index={0 0 1}' "Test grid function"
# CCTK_REAL gf010 TYPE=gf TAGS='index={0 1 0}' "Test grid function"
# CCTK_REAL gf011 TYPE=gf TAGS='index={0 1 1}' "Test grid function"
# CCTK_REAL gf100 TYPE=gf TAGS='index={1 0 0}' "Test grid function"
# CCTK_REAL gf101 TYPE=gf TAGS='index={1 0 1}' "Test grid function"
# CCTK_REAL gf110 TYPE=gf TAGS='index={1 1 0}' "Test grid function"
# CCTK_REAL gf111 TYPE=gf TAGS='index={1 1 1}' "Test grid function"

Replacement in TestProlongate/interface.ccl at line 24 [4.4774]

B:BD[5.1] → [5.1:76]

CCTK_REAL max_diff TYPE=SCALAR TIMELEVELS=1 "difference to expected value"

[5.1]

CCTK_REAL max_diff TYPE=scalar TIMELEVELS=1 "difference to expected value"

Replacement in TestProlongate/schedule.ccl at line 3 [4.5074]

B:BD[4.5124] → [4.5124:5142]

#STORAGE: data[1]

[4.5124]

[4.5142]

STORAGE: data
# STORAGE: gf000
# STORAGE: gf001
# STORAGE: gf010
# STORAGE: gf011
# STORAGE: gf100
# STORAGE: gf101
# STORAGE: gf110
# STORAGE: gf111

Replacement in TestProlongate/schedule.ccl at line 14 [4.5074]
B:BD[4.5143] → [6.48:99]
```
SCHEDULE TestProlongate_Regrid AT CCTK_POSTINITIAL
```
[4.5143]
[6.99]
```
SCHEDULE TestProlongate_Regrid AT postinitial
```
Replacement in TestProlongate/schedule.ccl at line 18 [4.5074]
B:BD[6.111] → [3.16361:16405]
```
  WRITES: CarpetX::regrid_error(everywhere)
```
[6.111]
[6.153]
```
  WRITES: CarpetX::regrid_error(interior)
```
Replacement in TestProlongate/schedule.ccl at line 21 [4.5074]
B:BD[6.175] → [6.175:223]
```
SCHEDULE TestProlongate_Regrid AT CCTK_POSTSTEP
```
[6.175]
[4.5188]
```
SCHEDULE TestProlongate_Regrid AT poststep
```
Replacement in TestProlongate/schedule.ccl at line 24 [4.5074]
B:BD[4.5200] → [3.16406:16450]
```
  WRITES: CarpetX::regrid_error(everywhere)
```
[4.5200]
[6.266]
```
  WRITES: CarpetX::regrid_error(interior)
```
Insertion in TestProlongate/schedule.ccl at line 26 [4.5074]
[6.287]
[2.0]

Replacement in TestProlongate/schedule.ccl at line 29 [4.5074]

B:BD[6.288] → [2.2:57]

SCHEDULE GROUP TestProlongate_SetValue AT CCTK_INITIAL

[6.288]

[2.57]

SCHEDULE GROUP TestProlongate_SetValue AT initial

Replacement in TestProlongate/schedule.ccl at line 33 [4.5074]
B:BD[2.81] → [2.81:133]
```
SCHEDULE GROUP TestProlongate_SetValue AT CCTK_EVOL
```
[2.81]
[2.133]
```
SCHEDULE GROUP TestProlongate_SetValue AT evol
```

Replacement in TestProlongate/schedule.ccl at line 37 [4.5074]

B:BD[2.157] → [2.157:213]

SCHEDULE TestProlongate_Test IN TestProlongate_SetValue

[2.157]

[6.333]

SCHEDULE TestProlongate_Set IN TestProlongate_SetValue

Replacement in TestProlongate/schedule.ccl at line 40 [4.5074]

B:BD[6.345] → [5.77:121]

  WRITES: data(interior) max_diff(interior)

[6.345]

[4.5200]

  WRITES: max_diff(interior)
  WRITES: data(interior)
  # WRITES: gf000(interior)
  # WRITES: gf001(interior)
  # WRITES: gf010(interior)
  # WRITES: gf011(interior)
  # WRITES: gf100(interior)
  # WRITES: gf101(interior)
  # WRITES: gf110(interior)
  # WRITES: gf111(interior)

Replacement in TestProlongate/schedule.ccl at line 52 [4.5074]

B:BD[4.5222] → [2.214:296]

SCHEDULE TestProlongate_Sync IN TestProlongate_SetValue AFTER TestProlongate_Test

[4.5222]

[2.296]

SCHEDULE TestProlongate_Sync IN TestProlongate_SetValue AFTER TestProlongate_Set

Insertion in TestProlongate/schedule.ccl at line 55 [4.5074]
[2.308]
[2.308]
```
  WRITES: data(boundary)
```

Insertion in TestProlongate/schedule.ccl at line 57 [4.5074]

[2.321]

  # SYNC: gf000
  # SYNC: gf001
  # SYNC: gf010
  # SYNC: gf011
  # SYNC: gf100
  # SYNC: gf101
  # SYNC: gf110
  # SYNC: gf111

Replacement in TestProlongate/schedule.ccl at line 67 [4.5074]

B:BD[2.344] → [2.344:430]

# need to apply b/c in POSTRESTRICT
SCHEDULE TestProlongate_Sync IN CCTK_POSTRESTRICT

[2.344]

[4.5293]

SCHEDULE TestProlongate_Sync IN postregrid
{
  LANG: C
  WRITES: data(boundary)
  SYNC: data
  # SYNC: gf000
  # SYNC: gf001
  # SYNC: gf010
  # SYNC: gf011
  # SYNC: gf100
  # SYNC: gf101
  # SYNC: gf110
  # SYNC: gf111
} "Sync up test data"
SCHEDULE TestProlongate_Sync IN postrestrict

Insertion in TestProlongate/schedule.ccl at line 85 [4.5074]
[4.5305]
[4.5305]
```
  WRITES: data(boundary)
```

Insertion in TestProlongate/schedule.ccl at line 87 [4.5074]

[4.5318]

  # SYNC: gf000
  # SYNC: gf001
  # SYNC: gf010
  # SYNC: gf011
  # SYNC: gf100
  # SYNC: gf101
  # SYNC: gf110
  # SYNC: gf111

Replacement in TestProlongate/schedule.ccl at line 97 [4.5074]
B:BD[5.123] → [5.123:170]
```
SCHEDULE TestProlongate_Check IN CCTK_ANALYSIS
```
[5.123]
[5.170]
```
SCHEDULE TestProlongate_Check IN analysis
```
Replacement in TestProlongate/src/make.code.defn at line 4 [4.5365]
B:BD[4.5452] → [4.5452:5467]
```
SRCS = test.cc
```
[4.5452]
[4.5467]
```
SRCS = test.cxx
```
File move: test.cc → test.cxx
BF:BFD[4.5347] → [4.5521:5530]
[4.5347]
[4.5531]

Replacement in TestProlongate/src/test.cxx at line 15 [4.5531]

B:BD[7.188] → [7.188:242]

extern "C" void TestProlongate_Test(CCTK_ARGUMENTS) {

[7.188]

[4.5693]

template <typename T>
T fun1d(const T x, const T dx, const bool avgx, const int order) {
  if (!avgx)
    return pow(x, order);
  else
    return (pow(x + dx / 2, order + 1) - pow(x - dx / 2, order + 1)) /
           ((order + 1) * dx);
}
// The grid stores the average values of the underlying function (x*y*z)**n
// which amounts to storing differences of the anti-derivative 1/(n+1)**3 *
// (x*y*z)**(n+1)
template <typename T>
T fun(const T x, const T y, const T z, const T dx, const T dy, const T dz,
      const bool avgx, const bool avgy, const bool avgz, const int order) {
  return fun1d(x, dx, avgx, order) * fun1d(y, dy, avgy, order) *
         fun1d(z, dz, avgz, order);
}
extern "C" void TestProlongate_Set(CCTK_ARGUMENTS) {

Replacement in TestProlongate/src/test.cxx at line 38 [4.5531]

B:BD[4.5747] → [8.0:78]

  // get prolongation order from driver, the parmeter is private since really

[4.5747]

[8.78]

  const CCTK_REAL dx = CCTK_DELTA_SPACE(0);
  const CCTK_REAL dy = CCTK_DELTA_SPACE(1);
  const CCTK_REAL dz = CCTK_DELTA_SPACE(2);
  // Get prolongation order from driver, the parmeter is private since really

Replacement in TestProlongate/src/test.cxx at line 45 [4.5531]

B:BD[8.145] → [3.16451:16526]

∅:D[3.16526] → [8.218:274]

B:BD[8.218] → [8.218:274]

  const void *order_p = CCTK_ParameterGet("prolongation_order", "CarpetX",
                                          &order_type);

[8.145]

[8.274]

  const void *order_p =
      CCTK_ParameterGet("prolongation_order", "CarpetX", &order_type);

Replacement in TestProlongate/src/test.cxx at line 49 [4.5531]

B:BD[8.332] → [8.332:406]

  const CCTK_INT operator_order = *static_cast<const CCTK_INT*>(order_p);

[8.332]

[9.0]

  const CCTK_INT operator_order = *static_cast<const CCTK_INT *>(order_p);

Replacement in TestProlongate/src/test.cxx at line 52 [4.5531]

B:BD[10.3862] → [10.3862:3898]

B:BD[10.3898] → [3.16527:16607]

  const void *prolongation_type_p =
    CCTK_ParameterGet("prolongation_type", "CarpetX", &prolongation_type_type);

[10.3862]

[10.3976]

  const void *prolongation_type_p = CCTK_ParameterGet(
      "prolongation_type", "CarpetX", &prolongation_type_type);

Replacement in TestProlongate/src/test.cxx at line 56 [4.5531]

B:BD[10.4062] → [10.4062:4155]

  const char* prolongation_type =
    *static_cast<const char* const*>(prolongation_type_p);

[10.4062]

[10.4155]

  const char *prolongation_type =
      *static_cast<const char *const *>(prolongation_type_p);

Replacement in TestProlongate/src/test.cxx at line 59 [4.5531]

B:BD[10.4196] → [10.4196:4248]

∅:D[9.310] → [8.406:407]

∅:D[10.4248] → [8.406:407]

B:BD[8.406] → [8.406:407]

B:BD[8.407] → [9.311:744]

    CCTK_EQUALS(prolongation_type, "conservative");
  // the grid stores the average values of the underlying function
  // (x*y*z)**n which amounts to storing differences of the anti-derivative
  // 1/(n+1)**3 * (x*y*z)**(n+1)
  auto avg_fun = [&](CCTK_REAL xi, int dir) {
    return 1./(operator_order + 1) *
           (pow(xi+CCTK_DELTA_SPACE(dir)/2., operator_order + 1) -
            pow(xi-CCTK_DELTA_SPACE(dir)/2., operator_order + 1)) /
           CCTK_DELTA_SPACE(dir);
  };

[10.4196]

[9.744]

      CCTK_EQUALS(prolongation_type, "conservative");

Replacement in TestProlongate/src/test.cxx at line 62 [4.5531]

B:BD[7.413] → [9.746:954]

    CCTK_REAL good_data;
    if (conservative_prolongation) {
      good_data = avg_fun(p.x, 0) * avg_fun(p.y, 1) * avg_fun(p.z, 2);
    } else {
      good_data = pow(p.x * p.y * p.z, operator_order);
    }

[7.413]

[9.954]

    const CCTK_REAL good_data = fun(
        p.x, p.y, p.z, dx, dy, dz, conservative_prolongation,
        conservative_prolongation, conservative_prolongation, operator_order);

Replacement in TestProlongate/src/test.cxx at line 68 [4.5531]
B:BD[5.290] → [5.290:308]
```
  *max_diff = 0.;
```
[5.290]
[4.6090]
```
  *max_diff = 0;
```
Insertion in TestProlongate/src/test.cxx at line 77 [4.5531]
[6.483]
[7.469]
Deletion in TestProlongate/src/test.cxx at line 81 [4.5531]
B:BD[6.582] → [6.582:583]
B:BD[6.583] → [11.270:307]
∅:D[11.307] → [6.608:620]
∅:D[7.552] → [6.608:620]
B:BD[6.608] → [6.608:620]
```
  if (cctk_iteration < regrid_after)
    return;
```

Replacement in TestProlongate/src/test.cxx at line 84 [4.5531]

B:BD[6.716] → [11.363:472]

∅:D[11.472] → [6.883:916]

∅:D[7.716] → [6.883:916]

B:BD[6.883] → [6.883:916]

    if (fabs(p.x) <= refined_radius && fabs(p.y) <= refined_radius &&
        fabs(p.z) <= refined_radius) {
      regrid_error[p.idx] = 1e3;

[6.716]

[6.916]

    if (cctk_iteration < regrid_after) {
      regrid_error[p.idx] = 0;

Replacement in TestProlongate/src/test.cxx at line 87 [4.5531]

B:BD[6.929] → [6.929:961]

      regrid_error[p.idx] = 0.;

[6.929]

[6.961]

      if (fabs(p.x) <= refined_radius && fabs(p.y) <= refined_radius &&
          fabs(p.z) <= refined_radius) {
        regrid_error[p.idx] = 1;
      } else {
        regrid_error[p.idx] = 0;
      }

Insertion in TestProlongate/src/test.cxx at line 100 [4.5531]

[5.417]


  const CCTK_REAL dx = CCTK_DELTA_SPACE(0);
  const CCTK_REAL dy = CCTK_DELTA_SPACE(1);
  const CCTK_REAL dz = CCTK_DELTA_SPACE(2);

Replacement in TestProlongate/src/test.cxx at line 112 [4.5531]

B:BD[5.714] → [3.16608:16683]

∅:D[3.16683] → [5.787:843]

B:BD[5.787] → [5.787:843]

  const void *order_p = CCTK_ParameterGet("prolongation_order", "CarpetX",
                                          &order_type);

[5.714]

[5.843]

  const void *order_p =
      CCTK_ParameterGet("prolongation_order", "CarpetX", &order_type);

Replacement in TestProlongate/src/test.cxx at line 116 [4.5531]

B:BD[5.901] → [5.901:975]

  const CCTK_INT operator_order = *static_cast<const CCTK_INT*>(order_p);

[5.901]

[5.975]

  const CCTK_INT operator_order = *static_cast<const CCTK_INT *>(order_p);

Replacement in TestProlongate/src/test.cxx at line 119 [4.5531]

B:BD[10.4279] → [10.4279:4315]

B:BD[10.4315] → [3.16684:16764]

  const void *prolongation_type_p =
    CCTK_ParameterGet("prolongation_type", "CarpetX", &prolongation_type_type);

[10.4279]

[10.4393]

  const void *prolongation_type_p = CCTK_ParameterGet(
      "prolongation_type", "CarpetX", &prolongation_type_type);

Replacement in TestProlongate/src/test.cxx at line 123 [4.5531]

B:BD[10.4479] → [10.4479:4572]

  const char* prolongation_type =
    *static_cast<const char* const*>(prolongation_type_p);

[10.4479]

[10.4572]

  const char *prolongation_type =
      *static_cast<const char *const *>(prolongation_type_p);

Replacement in TestProlongate/src/test.cxx at line 126 [4.5531]

B:BD[10.4613] → [10.4613:4665]

∅:D[10.4665] → [9.1293:1727]

B:BD[9.1293] → [9.1293:1727]

    CCTK_EQUALS(prolongation_type, "conservative");
  // the grid stores the average values of the underlying function
  // (x*y*z)**n which amounts to storing differences of the anti-derivative
  // 1/(n+1)**3 * (x*y*z)**(n+1)
  auto avg_fun = [&](CCTK_REAL xi, int dir) {
    return 1./(operator_order + 1) *
           (pow(xi+CCTK_DELTA_SPACE(dir)/2., operator_order + 1) -
            pow(xi-CCTK_DELTA_SPACE(dir)/2., operator_order + 1)) /
           CCTK_DELTA_SPACE(dir);
  };

[10.4613]

[9.1727]

      CCTK_EQUALS(prolongation_type, "conservative");

Replacement in TestProlongate/src/test.cxx at line 130 [4.5531]

B:BD[5.1071] → [9.1729:1937]

    CCTK_REAL good_data;
    if (conservative_prolongation) {
      good_data = avg_fun(p.x, 0) * avg_fun(p.y, 1) * avg_fun(p.z, 2);
    } else {
      good_data = pow(p.x * p.y * p.z, operator_order);
    }

[5.1071]

[9.1937]

    const CCTK_REAL good_data = fun(
        p.x, p.y, p.z, dx, dy, dz, conservative_prolongation,
        conservative_prolongation, conservative_prolongation, operator_order);

Replacement in TestProlongate/src/test.cxx at line 134 [4.5531]
∅:D[9.1995] → [5.1188:1230]
B:BD[5.1188] → [5.1188:1230]
```
    my_max_diff = max(diff, my_max_diff);
```
[9.1995]
[5.1230]
```
    my_max_diff = fmax(diff, my_max_diff);
```

Replacement in TestProlongate/src/test.cxx at line 137 [4.5531]

B:BD[5.1237] → [5.1237:1313]

  #pragma omp critical
  {
    *max_diff = max(*max_diff, my_max_diff);
  }

[5.1237]

[5.1313]

#pragma omp critical
  *max_diff = max(*max_diff, my_max_diff);

Deletion in TestProlongate/test/test_cc_co_o1.par at line 3 [12.103]
∅:D[3.16777] → [12.131:145]
B:BD[12.131] → [12.131:145]
```
    Formaline
```
Deletion in TestProlongate/test/test_cc_nc_o1.par at line 3 [4.6236]
∅:D[3.17293] → [4.6264:6278]
B:BD[4.6264] → [4.6264:6278]
```
    Formaline
```
Deletion in TestProlongate/test/test_cc_nc_o3.par at line 3 [13.104]
∅:D[3.17765] → [13.132:146]
B:BD[13.132] → [13.132:146]
```
    Formaline
```

TestProlongate: Update thorn

Dependencies

In channels

Change contents

Replacement in TestProlongate/interface.ccl at line 2 [4.4774]

Insertion in TestProlongate/interface.ccl at line 8 [4.4774]

Replacement in TestProlongate/interface.ccl at line 10 [4.4774]

Replacement in TestProlongate/interface.ccl at line 13 [4.4774]

Replacement in TestProlongate/interface.ccl at line 24 [4.4774]

Replacement in TestProlongate/schedule.ccl at line 3 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 14 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 18 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 21 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 24 [4.5074]

Insertion in TestProlongate/schedule.ccl at line 26 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 29 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 33 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 37 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 40 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 52 [4.5074]

Insertion in TestProlongate/schedule.ccl at line 55 [4.5074]

Insertion in TestProlongate/schedule.ccl at line 57 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 67 [4.5074]

Insertion in TestProlongate/schedule.ccl at line 85 [4.5074]

Insertion in TestProlongate/schedule.ccl at line 87 [4.5074]

Replacement in TestProlongate/schedule.ccl at line 97 [4.5074]

Replacement in TestProlongate/src/make.code.defn at line 4 [4.5365]

File move: test.cc → test.cxx