eschnett/cactusamrex - Change W62UFG4ZAE6JRINNKSA3EIH35S4NSKAJSAVSCCNN4ME75ZXCSVTQC

TestProlongate: Test ddf operators

Created by Erik Schnetter on January 3, 2020

W62UFG4ZAE6JRINNKSA3EIH35S4NSKAJSAVSCCNN4ME75ZXCSVTQC

Dependencies

In channels

master

main

Change contents

Replacement in TestProlongate/interface.ccl at line 15 [3.4774]

B:BD[2.119] → [2.119:663]

# 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"

[2.119]

[4.0]

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/schedule.ccl at line 4 [3.5074]

B:BD[2.754] → [2.754:890]

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

[2.754]

[2.890]

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

Replacement in TestProlongate/schedule.ccl at line 42 [3.5074]

B:BD[2.1281] → [2.1281:1505]

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

[2.1281]

[3.5200]

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

Insertion in TestProlongate/schedule.ccl at line 56 [3.5074]

[2.1613]

[5.308]

  WRITES: gf000(boundary)
  WRITES: gf001(boundary)
  WRITES: gf010(boundary)
  WRITES: gf011(boundary)
  WRITES: gf100(boundary)
  WRITES: gf101(boundary)
  WRITES: gf110(boundary)
  WRITES: gf111(boundary)

Replacement in TestProlongate/schedule.ccl at line 65 [3.5074]

B:BD[5.321] → [2.1614:1742]

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

[5.321]

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

Insertion in TestProlongate/schedule.ccl at line 79 [3.5074]

[2.1823]

  WRITES: gf000(boundary)
  WRITES: gf001(boundary)
  WRITES: gf010(boundary)
  WRITES: gf011(boundary)
  WRITES: gf100(boundary)
  WRITES: gf101(boundary)
  WRITES: gf110(boundary)
  WRITES: gf111(boundary)

Replacement in TestProlongate/schedule.ccl at line 88 [3.5074]

B:BD[2.1836] → [2.1836:1964]

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

[2.1836]

[2.1964]

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

Insertion in TestProlongate/schedule.ccl at line 102 [3.5074]

[2.2058]

[3.5305]

  WRITES: gf000(boundary)
  WRITES: gf001(boundary)
  WRITES: gf010(boundary)
  WRITES: gf011(boundary)
  WRITES: gf100(boundary)
  WRITES: gf101(boundary)
  WRITES: gf110(boundary)
  WRITES: gf111(boundary)

Replacement in TestProlongate/schedule.ccl at line 111 [3.5074]

B:BD[3.5318] → [2.2059:2187]

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

[3.5318]

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

Insertion in TestProlongate/src/test.cxx at line 66 [3.5531]

[6.983]

[3.6084]

  });
  Loop::loop_int<0, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 0, 0, operator_order);
    gf000[p.idx] = good_data;
  });
  Loop::loop_int<0, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 0, 1, operator_order);
    gf001[p.idx] = good_data;
  });
  Loop::loop_int<0, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 1, 0, operator_order);
    gf010[p.idx] = good_data;
  });
  Loop::loop_int<0, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 1, 1, operator_order);
    gf011[p.idx] = good_data;
  });
  Loop::loop_int<1, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 0, 0, operator_order);
    gf100[p.idx] = good_data;
  });
  Loop::loop_int<1, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 0, 1, operator_order);
    gf101[p.idx] = good_data;

Insertion in TestProlongate/src/test.cxx at line 98 [3.5531]

[3.6090]

[4.289]

  Loop::loop_int<1, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 1, 0, operator_order);
    gf110[p.idx] = good_data;
  });
  Loop::loop_int<1, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 1, 1, operator_order);
    gf111[p.idx] = good_data;
  });

Insertion in TestProlongate/src/test.cxx at line 170 [3.5531]
[4.1005]
[4.1005]

Insertion in TestProlongate/src/test.cxx at line 179 [3.5531]

[4.1237]

[2.4934]

  Loop::loop_int<0, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 0, 0, operator_order);
    const CCTK_REAL diff = fabs(gf000[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });
  Loop::loop_int<0, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 0, 1, operator_order);
    const CCTK_REAL diff = fabs(gf001[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });
  Loop::loop_int<0, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 1, 0, operator_order);
    const CCTK_REAL diff = fabs(gf010[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });
  Loop::loop_int<0, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 0, 1, 1, operator_order);
    const CCTK_REAL diff = fabs(gf011[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });
  Loop::loop_int<1, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 0, 0, operator_order);
    const CCTK_REAL diff = fabs(gf100[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });
  Loop::loop_int<1, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 0, 1, operator_order);
    const CCTK_REAL diff = fabs(gf101[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });
  Loop::loop_int<1, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 1, 0, operator_order);
    const CCTK_REAL diff = fabs(gf110[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });
  Loop::loop_int<1, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    const CCTK_REAL good_data =
        fun(p.x, p.y, p.z, dx, dy, dz, 1, 1, 1, operator_order);
    const CCTK_REAL diff = fabs(gf111[p.idx] - good_data);
    my_max_diff = fmax(diff, my_max_diff);
  });

File addition: test_ddf_o1 (drwxr-xr-x)
[3.6224]
File addition: testprolongate-max_diff.tsv (-rw-r--r--)
[0.5606]
```
0	0	0
1	0.025	0
2	0.05	0
```

File addition: test_ddf_o1.par (-rw-r--r--)

[3.6224]

ActiveThorns = "
    CarpetX
    IOUtil
    TestProlongate
"
 
$nlevels = 2
$ncells = 20
$blocking_factor = 2
Cactus::cctk_show_schedule = no
Cactus::cctk_itlast = 2
CarpetX::ghost_size = 3
CarpetX::prolongation_type = "ddf"
CarpetX::prolongation_order = 1
CarpetX::verbose = yes
CarpetX::ncells_x = $ncells
CarpetX::ncells_y = $ncells
CarpetX::ncells_z = $ncells
CarpetX::blocking_factor_x = $blocking_factor
CarpetX::blocking_factor_y = 1 # $blocking_factor
CarpetX::blocking_factor_z = 1 # $blocking_factor
CarpetX::max_num_levels = $nlevels
CarpetX::regrid_every = 1
CarpetX::regrid_error_threshold = 0.01
CarpetX::dtfac = 0.5
CarpetX::out_plotfile_groups = ""
CarpetX::out_tsv = yes
IO::out_dir = $parfile
IO::out_every = 1

File addition: test_ddf_o3 (drwxr-xr-x)
[3.6224]
File addition: testprolongate-max_diff.tsv (-rw-r--r--)
[0.6436]
```
0	0	0
1	0.025	0
2	0.05	0
```

File addition: test_ddf_o3.par (-rw-r--r--)