eschnett/cactusamrex - Change YRL5XGH3TJKZKX5L6TBQF7FL3OTVIEBMC23ISQROZ3HQ3MCZQ7HQC

Maxwell: Implement standing wave with fixed (non-periodic) boundaries

Created by Erik Schnetter on December 18, 2019

YRL5XGH3TJKZKX5L6TBQF7FL3OTVIEBMC23ISQROZ3HQ3MCZQ7HQC

Dependencies

In channels

master

main

Change contents

Insertion in Maxwell/param.ccl at line 6 [3.5230]
[3.5331]
[3.5331]
```
  "standing wave" :: ""
```

Insertion in Maxwell/schedule.ccl at line 24 [3.5933]

[3.6413]

[2.2402]

} "Set up initial conditions"
SCHEDULE Maxwell_Boundaries AT initial AFTER Maxwell_Initial
{
  LANG: C
  READS: dyz(interior) dzx(interior) dxy(interior)
  READS: byz(interior) bzx(interior) bxy(interior)
  WRITES: dyz(boundary) dzx(boundary) dxy(boundary)
  WRITES: byz(boundary) bzx(boundary) bxy(boundary)

Replacement in Maxwell/schedule.ccl at line 35 [3.5933]
B:BD[3.6479] → [2.2423:2453]
```
} "Set up initial conditions"
```
[3.6479]
[3.6516]
```
} "Apply boundary conditions"
```

Insertion in Maxwell/schedule.ccl at line 44 [3.5933]

[3.6586]

[2.2454]

  READS: dyz(interior) dzx(interior) dxy(interior)
  READS: byz(interior) bzx(interior) bxy(interior)
  WRITES: dyz(boundary) dzx(boundary) dxy(boundary)
  WRITES: byz(boundary) bzx(boundary) bxy(boundary)

Insertion in Maxwell/schedule.ccl at line 63 [3.5933]

[3.7130]

[2.2592]

} "Calculate RHS"
SCHEDULE Maxwell_RHSBoundaries IN ODESolvers_RHS AFTER Maxwell_RHS
{
  LANG: C
  READS: dtdyz(interior) dtdzx(interior) dtdxy(interior)
  READS: dtbyz(interior) dtbzx(interior) dtbxy(interior)
  WRITES: dtdyz(boundary) dtdzx(boundary) dtdxy(boundary)
  WRITES: dtbyz(boundary) dtbzx(boundary) dtbxy(boundary)

Replacement in Maxwell/schedule.ccl at line 74 [3.5933]
∅:D[2.2644] → [3.7256:7274]
B:BD[3.7256] → [3.7256:7274]
```
} "Calculate RHS"
```
[2.2644]
[3.7274]
```
} "Apply RHS boundary conditions"
```
Replacement in Maxwell/schedule.ccl at line 85 [3.5933]
B:BD[2.2759] → [2.2759:2784]
∅:D[2.2784] → [3.7581:7606]
B:BD[3.7581] → [3.7581:7606]
B:BD[3.7606] → [2.2785:2811]
```
  WRITES: divd(interior)
  WRITES: divb(interior)
  SYNC: divd
  SYNC: divb
```
[2.2759]
[3.7606]
```
  WRITES: divd(everywhere)
  WRITES: divb(everywhere)
```

Replacement in Maxwell/schedule.ccl at line 94 [3.5933]

B:BD[2.2964] → [2.2964:3099]

∅:D[2.3099] → [4.6789:6802]

B:BD[4.6789] → [4.6789:6802]

  WRITES: avgdyz(interior) avgdzx(interior) avgdxy(interior)
  WRITES: avgbyz(interior) avgbzx(interior) avgbxy(interior)
  SYNC: avgd
  SYNC: avgb

[2.2964]

[2.3100]

  WRITES: avgdyz(everywhere) avgdzx(everywhere) avgdxy(everywhere)
  WRITES: avgbyz(everywhere) avgbzx(everywhere) avgbxy(everywhere)

Replacement in Maxwell/src/average.cxx at line 35 [4.346014]

B:BD[4.349542] → [2.3994:4048]

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

[4.349542]

[2.4048]

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

Replacement in Maxwell/src/average.cxx at line 38 [4.346014]

B:BD[4.349651] → [2.4102:4156]

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

[4.349651]

[2.4156]

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

Replacement in Maxwell/src/average.cxx at line 41 [4.346014]

B:BD[4.349760] → [2.4210:4264]

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

[4.349760]

[2.4264]

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

Replacement in Maxwell/src/average.cxx at line 45 [4.346014]

B:BD[4.349870] → [2.4318:4372]

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

[4.349870]

[2.4372]

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

Replacement in Maxwell/src/average.cxx at line 48 [4.346014]

B:BD[4.350104] → [2.4426:4480]

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

[4.350104]

[2.4480]

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

Replacement in Maxwell/src/average.cxx at line 51 [4.346014]

B:BD[4.350221] → [2.4534:4588]

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

[4.350221]

[2.4588]

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

Insertion in Maxwell/src/boundaries.cxx at line 1 [3.7657]
[3.7657]
[3.7658]
```
#include <loop.hxx>
```
Insertion in Maxwell/src/boundaries.cxx at line 8 [3.7657]
[3.7762]
[3.7762]
```
using namespace Loop;
```

Insertion in Maxwell/src/boundaries.cxx at line 13 [3.7657]

[3.7888]


  for (int d = 0; d < dim; ++d)
    assert(cctk_nghostzones[d] <= 1);
  const GF3D<CCTK_REAL, 0, 1, 1> dyz_(cctkGH, dyz);
  const GF3D<CCTK_REAL, 1, 0, 1> dzx_(cctkGH, dzx);
  const GF3D<CCTK_REAL, 1, 1, 0> dxy_(cctkGH, dxy);
  const GF3D<CCTK_REAL, 0, 1, 1> byz_(cctkGH, byz);
  const GF3D<CCTK_REAL, 1, 0, 1> bzx_(cctkGH, bzx);
  const GF3D<CCTK_REAL, 1, 1, 0> bxy_(cctkGH, bxy);

Replacement in Maxwell/src/boundaries.cxx at line 25 [3.7657]

B:BD[3.7889] → [3.7889:7905]

  // do nothing

[3.7889]

[3.7905]

  loop_bnd<0, 1, 1>(cctkGH, [&](const PointDesc &p) {
    int f = p.NI[0] != 0 ? -1 : +1;
    dyz_(p.I) = f * dyz_(p.I - p.NI);
  });
  loop_bnd<1, 0, 1>(cctkGH, [&](const PointDesc &p) {
    int f = p.NI[1] != 0 ? -1 : +1;
    dzx_(p.I) = f * dzx_(p.I - p.NI);
  });
  loop_bnd<1, 1, 0>(cctkGH, [&](const PointDesc &p) {
    int f = p.NI[2] != 0 ? -1 : +1;
    dxy_(p.I) = f * dxy_(p.I - p.NI);
  });
  loop_bnd<0, 1, 1>(cctkGH, [&](const PointDesc &p) { byz_(p.I) = 0; });
  loop_bnd<1, 0, 1>(cctkGH, [&](const PointDesc &p) { bzx_(p.I) = 0; });
  loop_bnd<1, 1, 0>(cctkGH, [&](const PointDesc &p) { bxy_(p.I) = 0; });

Replacement in Maxwell/src/constraints.cxx at line 42 [3.7950]

B:BD[3.9649] → [2.5565:5619]

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

[3.9649]

[2.5619]

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

Replacement in Maxwell/src/constraints.cxx at line 46 [3.7950]

B:BD[3.10361] → [2.5687:5741]

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

[3.10361]

[2.5741]

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

Insertion in Maxwell/src/initial.cxx at line 64 [3.15681]

[2.7371]

[3.24749]

    });
  } else if (CCTK_EQUALS(setup, "standing wave")) {
    // wave number
    const CCTK_REAL kx = CCTK_REAL(M_PI) * spatial_frequency_x;
    const CCTK_REAL ky = CCTK_REAL(M_PI) * spatial_frequency_y;
    const CCTK_REAL kz = CCTK_REAL(M_PI) * spatial_frequency_z;
    const CCTK_REAL omega = sqrt(pow2(kx) + pow2(ky) + pow2(kz));
    // amplitude
    const CCTK_REAL hx = amplitude_x;
    const CCTK_REAL hy = amplitude_y;
    const CCTK_REAL hz = amplitude_z;
    loop_int<0, 1, 1>(cctkGH, [&](const PointDesc &p) {
      dyz_(p.I) =
          hx * cos(omega * t) * sin(kx * p.x) * sin(ky * p.y) * sin(kz * p.z);
    });
    loop_int<1, 0, 1>(cctkGH, [&](const PointDesc &p) {
      dzx_(p.I) =
          hy * cos(omega * t) * sin(kx * p.x) * sin(ky * p.y) * sin(kz * p.z);
    });
    loop_int<1, 1, 0>(cctkGH, [&](const PointDesc &p) {
      dxy_(p.I) =
          hz * cos(omega * t) * sin(kx * p.x) * sin(ky * p.y) * sin(kz * p.z);

Insertion in Maxwell/src/initial.cxx at line 91 [3.15681]

[3.24764]

[3.25127]

    loop_int<0, 1, 1>(cctkGH, [&](const PointDesc &p) {
      byz_(p.I) =
          hz * cos(omega * t) * cos(kx * p.x) * cos(ky * p.y) * cos(kz * p.z);
    });
    loop_int<1, 0, 1>(cctkGH, [&](const PointDesc &p) {
      bzx_(p.I) =
          hx * cos(omega * t) * cos(kx * p.x) * cos(ky * p.y) * cos(kz * p.z);
    });
    loop_int<1, 1, 0>(cctkGH, [&](const PointDesc &p) {
      bxy_(p.I) =
          hy * cos(omega * t) * cos(kx * p.x) * cos(ky * p.y) * cos(kz * p.z);
    });

Insertion in Maxwell/src/make.code.defn at line 8 [3.25201]
[4.367812]
[2.7397]
```
	rhsboundaries.cxx			\
```

File addition: rhsboundaries.cxx (-rw-r--r--)

[3.7639]

#include <loop.hxx>
#include <cctk.h>
#include <cctk_Arguments_Checked.h>
#include <cctk_Parameters.h>
namespace Maxwell {
using namespace Loop;
extern "C" void Maxwell_RHSBoundaries(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_Maxwell_RHSBoundaries;
  DECLARE_CCTK_PARAMETERS;
  const GF3D<CCTK_REAL, 0, 1, 1> dtdyz_(cctkGH, dtdyz);
  const GF3D<CCTK_REAL, 1, 0, 1> dtdzx_(cctkGH, dtdzx);
  const GF3D<CCTK_REAL, 1, 1, 0> dtdxy_(cctkGH, dtdxy);
  const GF3D<CCTK_REAL, 0, 1, 1> dtbyz_(cctkGH, dtbyz);
  const GF3D<CCTK_REAL, 1, 0, 1> dtbzx_(cctkGH, dtbzx);
  const GF3D<CCTK_REAL, 1, 1, 0> dtbxy_(cctkGH, dtbxy);
  loop_bnd<0, 1, 1>(cctkGH, [&](const PointDesc &p) { dtdyz_(p.I) = 0; });
  loop_bnd<1, 0, 1>(cctkGH, [&](const PointDesc &p) { dtdzx_(p.I) = 0; });
  loop_bnd<1, 1, 0>(cctkGH, [&](const PointDesc &p) { dtdxy_(p.I) = 0; });
  loop_bnd<0, 1, 1>(cctkGH, [&](const PointDesc &p) { dtbyz_(p.I) = 0; });
  loop_bnd<1, 0, 1>(cctkGH, [&](const PointDesc &p) { dtbzx_(p.I) = 0; });
  loop_bnd<1, 1, 0>(cctkGH, [&](const PointDesc &p) { dtbxy_(p.I) = 0; });
}
} // namespace Maxwell