eschnett/cactusamrex - Change FAKFUJ4VZQFKOSTNHZGLRIIMZDXFAO6QOPI7UJKMKSKPXWJSZZMQC

MaxwellToyCarpetX: Modernize thorn; calculate initial fields from potentials

Created by Erik Schnetter on November 2, 2019

FAKFUJ4VZQFKOSTNHZGLRIIMZDXFAO6QOPI7UJKMKSKPXWJSZZMQC

Dependencies

In channels

master

main

Change contents

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 2 [3.30]
B:BD[3.48] → [3.48:58]
```
    AMReX
```
[3.48]
[3.58]
```
    CarpetX
```
Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 5 [3.30]
B:BD[3.83] → [3.83:103]
```
    MaxwellToyAMReX
```
[3.83]
[4.9213]
```
    MaxwellToyCarpetX
```
Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 21 [3.30]
B:BD[3.285] → [5.0:20]
```
AMReX::verbose = no
```
[3.285]
[3.308]
```
CarpetX::verbose = no
```
Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 23 [3.30]
B:BD[3.309] → [3.309:366]
```
AMReX::xmin = -1.0
AMReX::ymin = -1.0
AMReX::zmin = -1.0
```
[3.309]
[3.366]
```
CarpetX::xmin = -1.0
CarpetX::ymin = -1.0
CarpetX::zmin = -1.0
```
Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 27 [3.30]
B:BD[3.367] → [3.367:424]
```
AMReX::xmax = +1.0
AMReX::ymax = +1.0
AMReX::zmax = +1.0
```
[3.367]
[3.424]
```
CarpetX::xmax = +1.0
CarpetX::ymax = +1.0
CarpetX::zmax = +1.0
```

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 31 [3.30]

B:BD[3.425] → [3.425:503]

AMReX::ncells_x = $ncells
AMReX::ncells_y = $ncells
AMReX::ncells_z = $ncells

[3.425]

[6.35]

CarpetX::ncells_x = $ncells
CarpetX::ncells_y = $ncells
CarpetX::ncells_z = $ncells

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 35 [3.30]
B:BD[6.36] → [6.36:69]
```
AMReX::prolongation_type = "ddf"
```
[6.36]
[3.503]
```
CarpetX::prolongation_type = "ddf"
```

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 37 [3.30]

B:BD[3.504] → [3.504:562]

B:BD[3.562] → [7.2301:2338]

AMReX::max_num_levels = $nlevels
AMReX::regrid_every = 16
AMReX::regrid_error_threshold = 0.01

[3.504]

[3.598]

CarpetX::max_num_levels = $nlevels
CarpetX::regrid_every = 16
CarpetX::regrid_error_threshold = 0.01

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 41 [3.30]
B:BD[3.599] → [3.599:618]
```
AMReX::dtfac = 0.5
```
[3.599]
[3.618]
```
CarpetX::dtfac = 0.5
```
Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 43 [3.30]
B:BD[3.619] → [4.9274:9306]
```
MaxwellToyAMReX::evolve_b = yes
```
[3.619]
[4.9306]
```
MaxwellToyCarpetX::evolve_b = yes
```
Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 45 [3.30]
∅:D[4.9307] → [3.619:672]
B:BD[3.619] → [3.619:672]
```
MaxwellToyAMReX::initial_condition = "Gaussian wave"
```
[4.9307]
[3.672]
```
MaxwellToyCarpetX::initial_condition = "Gaussian wave"
```

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 47 [3.30]

B:BD[3.673] → [3.673:802]

MaxwellToyAMReX::spatial_frequency_x = 0.5
MaxwellToyAMReX::spatial_frequency_y = 0.0
MaxwellToyAMReX::spatial_frequency_z = 0.0

[3.673]

[3.802]

MaxwellToyCarpetX::spatial_frequency_x = 0.5
MaxwellToyCarpetX::spatial_frequency_y = 0.0
MaxwellToyCarpetX::spatial_frequency_z = 0.0

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 51 [3.30]

B:BD[3.803] → [6.70:149]

MaxwellToyAMReX::analyse_every = 1   # TODO $ncells * 2 ** ($nlevels - 1) / 32

[3.803]

[4.9376]

MaxwellToyCarpetX::analyse_every = 1   # TODO $ncells * 2 ** ($nlevels - 1) / 32

Insertion in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 56 [3.30]

[3.878]

[2.13019]

CarpetX::out_plotfile_groups = "
    MaxwellToyCarpetX::potential_phi
    MaxwellToyCarpetX::constraints_dive
    MaxwellToyCarpetX::constraints_divb
"

Deletion in MaxwellToyCarpetX/schedule.ccl at line 2 [8.15053]
∅:D[2.14081] → [8.15103:15104]
B:BD[8.15103] → [8.15103:15104]
B:BD[8.15104] → [4.10010:10026]
```
if (evolve_b) {
```

Replacement in MaxwellToyCarpetX/schedule.ccl at line 3 [8.15053]

B:BD[9.1957] → [2.14082:14133]

∅:D[2.14133] → [4.10075:10091]

B:BD[4.10075] → [4.10075:10091]

B:BD[4.10091] → [10.13387:13388]

B:BD[10.13388] → [9.1958:2322]

  SCHEDULE MaxwellToyCarpetX_Initialize AT initial
  {
    LANG: C
    WRITES: potential_phi(interior)
    WRITES: potential_ax(interior) potential_ay(interior) potential_az(interior)
    WRITES: field_ex(interior) field_ey(interior) field_ez(interior)
    WRITES: field_bx(interior) field_by(interior) field_bz(interior)
    WRITES: current_rho(interior)
    WRITES: current_jx(interior) current_jy(interior) current_jz(interior)

[9.1957]

[10.13780]

SCHEDULE MaxwellToyCarpetX_Initialize AT initial
{
  LANG: C

Replacement in MaxwellToyCarpetX/schedule.ccl at line 7 [8.15053]

∅:D[10.13781] → [4.10091:10361]

B:BD[4.10091] → [4.10091:10361]

B:BD[4.10361] → [9.2323:2324]

∅:D[9.2324] → [4.10361:10370]

B:BD[4.10361] → [4.10361:10370]

B:BD[4.10370] → [9.2325:2326]

B:BD[9.2326] → [2.14134:14185]

∅:D[2.14185] → [4.10419:10435]

B:BD[4.10419] → [4.10419:10435]

B:BD[4.10435] → [10.13782:13783]

B:BD[10.13783] → [9.2327:2622]

∅:D[9.2622] → [10.14100:14102]

B:BD[10.14100] → [10.14100:14102]

∅:D[10.14102] → [4.10435:10668]

B:BD[4.10435] → [4.10435:10668]

    SYNC: potential_phi
    SYNC: potential_ax potential_ay potential_az
    SYNC: field_ex field_ey field_ez
    SYNC: field_bx field_by field_bz
    SYNC: current_rho
    SYNC: current_jx current_jy current_jz
  } "Set up initial conditions for the Maxwell equations"
} else {
  SCHEDULE MaxwellToyCarpetX_Initialize AT initial
  {
    LANG: C
    WRITES: potential_phi(interior)
    WRITES: potential_ax(interior) potential_ay(interior) potential_az(interior)
    WRITES: field_ex(interior) field_ey(interior) field_ez(interior)
    WRITES: current_rho(interior)
    WRITES: current_jx(interior) current_jy(interior) current_jz(interior)
 
    SYNC: potential_phi
    SYNC: potential_ax potential_ay potential_az
    SYNC: field_ex field_ey field_ez
    SYNC: current_rho
    SYNC: current_jx current_jy current_jz
  } "Set up initial conditions for the Maxwell equations"

[10.13781]

[4.10668]

  WRITES: potential_phi(interior)
  WRITES: potential_ax(interior) potential_ay(interior) potential_az(interior)
  WRITES: current_rho(interior)
  WRITES: current_jx(interior) current_jy(interior) current_jz(interior)

Replacement in MaxwellToyCarpetX/schedule.ccl at line 12 [8.15053]

B:BD[4.10669] → [2.14186:14273]

∅:D[2.14273] → [4.10752:10768]

B:BD[4.10752] → [4.10752:10768]

  SCHEDULE MaxwellToyCarpetX_Dependents1 AT initial AFTER MaxwellToyCarpetX_Initialize
  {
    LANG: C

[4.10669]

[10.14103]

  SYNC: potential_phi
  SYNC: potential_ax potential_ay potential_az
  SYNC: current_rho
  SYNC: current_jx current_jy current_jz
} "Set up initial conditions for the Maxwell equations"

Replacement in MaxwellToyCarpetX/schedule.ccl at line 18 [8.15053]

B:BD[10.14104] → [9.2623:2772]

    READS: potential_ax(interior) potential_ay(interior) potential_az(interior)
    WRITES: field_bx(interior) field_by(interior) field_bz(interior)

[10.14104]

[10.14179]

SCHEDULE MaxwellToyCarpetX_InitializeFields AT initial AFTER MaxwellToyCarpetX_Initialize
{
  LANG: C

Replacement in MaxwellToyCarpetX/schedule.ccl at line 22 [8.15053]

∅:D[10.14180] → [4.10768:10840]

B:BD[4.10768] → [4.10768:10840]

    SYNC: field_bx field_by field_bz
  } "Calculate the magnetic field"

[10.14180]

[9.2773]

  READS: potential_phi(everywhere)
  READS: potential_ax(everywhere) potential_ay(everywhere) potential_az(everywhere)
  READS: current_rho(everywhere)
  READS: current_jx(everywhere) current_jy(everywhere) current_jz(everywhere)
  WRITES: field_ex(interior) field_ey(interior) field_ez(interior)
  WRITES: field_bx(interior) field_by(interior) field_bz(interior)

Replacement in MaxwellToyCarpetX/schedule.ccl at line 29 [8.15053]
∅:D[9.2774] → [4.10840:10842]
B:BD[4.10840] → [4.10840:10842]
```
}
```
[9.2774]
[11.1772]
```
  SYNC: field_ex field_ey field_ez
  SYNC: field_bx field_by field_bz
} "Calculate Maxwell fields from initial conditions"
```

Deletion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 123 [8.16301]

B:BD[8.18735] → [7.3078:3210]

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

Deletion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 130 [8.16301]

B:BD[12.23235] → [12.23235:23573]

  const Loop::GF3D<CCTK_REAL, 1, 0, 0> ex_(cctkGH, ex);
  const Loop::GF3D<CCTK_REAL, 0, 1, 0> ey_(cctkGH, ey);
  const Loop::GF3D<CCTK_REAL, 0, 0, 1> ez_(cctkGH, ez);
  const Loop::GF3D<CCTK_REAL, 0, 1, 1> bx_(cctkGH, bx);
  const Loop::GF3D<CCTK_REAL, 1, 0, 1> by_(cctkGH, by);
  const Loop::GF3D<CCTK_REAL, 1, 1, 0> bz_(cctkGH, bz);

Deletion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 151 [8.16301]

B:BD[4.12972] → [4.12972:14458]


    Loop::loop_all<1, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
      ex_(p.i, p.j, p.k) =
          -xderiv(linear<CCTK_REAL>, dx)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(linear<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).ax;
    });
    Loop::loop_all<0, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
      ey_(p.i, p.j, p.k) =
          -yderiv(linear<CCTK_REAL>, dy)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(linear<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).ay;
    });
    Loop::loop_all<0, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
      ez_(p.i, p.j, p.k) =
          -zderiv(linear<CCTK_REAL>, dz)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(linear<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).az;
    });
    if (evolve_b) {
      Loop::loop_all<0, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
        bx_(p.i, p.j, p.k) =
            yderiv(linear<CCTK_REAL>, dy)(t, p.x, p.y, p.z).az -
            zderiv(linear<CCTK_REAL>, dz)(t, p.x, p.y, p.z).ay;
      });
      Loop::loop_all<1, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
        by_(p.i, p.j, p.k) =
            zderiv(linear<CCTK_REAL>, dz)(t, p.x, p.y, p.z).ax -
            xderiv(linear<CCTK_REAL>, dx)(t, p.x, p.y, p.z).az;
      });
      Loop::loop_all<1, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
        bz_(p.i, p.j, p.k) =
            xderiv(linear<CCTK_REAL>, dx)(t, p.x, p.y, p.z).ay -
            yderiv(linear<CCTK_REAL>, dy)(t, p.x, p.y, p.z).ax;
      });
    }

Deletion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 178 [8.16301]

B:BD[8.19348] → [12.24358:24453]

B:BD[12.24453] → [7.3211:3289]

∅:D[7.3289] → [12.24533:24711]

B:BD[12.24533] → [12.24533:24711]

B:BD[12.24711] → [7.3290:3368]

∅:D[7.3368] → [12.24791:24969]

B:BD[12.24791] → [12.24791:24969]

B:BD[12.24969] → [7.3369:3447]

∅:D[7.3447] → [12.25049:25132]

B:BD[12.25049] → [12.25049:25132]

∅:D[12.25132] → [8.20392:20393]

B:BD[8.20392] → [8.20392:20393]

B:BD[8.20393] → [4.14952:15716]

∅:D[4.15716] → [7.3857:3858]

B:BD[7.3857] → [7.3857:3858]

    Loop::loop_all<1, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
      ex_(p.i, p.j, p.k) =
          -xderiv(plane_wave<CCTK_REAL>, dx)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(plane_wave<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).ax;
    });
    Loop::loop_all<0, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
      ey_(p.i, p.j, p.k) =
          -yderiv(plane_wave<CCTK_REAL>, dy)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(plane_wave<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).ay;
    });
    Loop::loop_all<0, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
      ez_(p.i, p.j, p.k) =
          -zderiv(plane_wave<CCTK_REAL>, dz)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(plane_wave<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).az;
    });
    if (evolve_b) {
      Loop::loop_all<0, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
        bx_(p.i, p.j, p.k) =
            yderiv(plane_wave<CCTK_REAL>, dy)(t, p.x, p.y, p.z).az -
            zderiv(plane_wave<CCTK_REAL>, dz)(t, p.x, p.y, p.z).ay;
      });
      Loop::loop_all<1, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
        by_(p.i, p.j, p.k) =
            zderiv(plane_wave<CCTK_REAL>, dz)(t, p.x, p.y, p.z).ax -
            xderiv(plane_wave<CCTK_REAL>, dx)(t, p.x, p.y, p.z).az;
      });
      Loop::loop_all<1, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
        bz_(p.i, p.j, p.k) =
            xderiv(plane_wave<CCTK_REAL>, dx)(t, p.x, p.y, p.z).ay -
            yderiv(plane_wave<CCTK_REAL>, dy)(t, p.x, p.y, p.z).ax;
      });
    }

Deletion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 204 [8.16301]

B:BD[7.4923] → [7.4923:5710]

B:BD[7.5710] → [4.15717:16499]

∅:D[4.16499] → [12.25846:25847]

B:BD[8.20486] → [12.25846:25847]

    Loop::loop_all<1, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
      ex_(p.i, p.j, p.k) =
          -xderiv(gaussian_wave<CCTK_REAL>, dx)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(gaussian_wave<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).ax;
    });
    Loop::loop_all<0, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
      ey_(p.i, p.j, p.k) =
          -yderiv(gaussian_wave<CCTK_REAL>, dy)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(gaussian_wave<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).ay;
    });
    Loop::loop_all<0, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
      ez_(p.i, p.j, p.k) =
          -zderiv(gaussian_wave<CCTK_REAL>, dz)(t + dt / 2, p.x, p.y, p.z).phi -
          tderiv(gaussian_wave<CCTK_REAL>, dt)(t + dt / 2, p.x, p.y, p.z).az;
    });
    if (evolve_b) {
      Loop::loop_all<0, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
        bx_(p.i, p.j, p.k) =
            yderiv(gaussian_wave<CCTK_REAL>, dy)(t, p.x, p.y, p.z).az -
            zderiv(gaussian_wave<CCTK_REAL>, dz)(t, p.x, p.y, p.z).ay;
      });
      Loop::loop_all<1, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
        by_(p.i, p.j, p.k) =
            zderiv(gaussian_wave<CCTK_REAL>, dz)(t, p.x, p.y, p.z).ax -
            xderiv(gaussian_wave<CCTK_REAL>, dx)(t, p.x, p.y, p.z).az;
      });
      Loop::loop_all<1, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
        bz_(p.i, p.j, p.k) =
            xderiv(gaussian_wave<CCTK_REAL>, dx)(t, p.x, p.y, p.z).ay -
            yderiv(gaussian_wave<CCTK_REAL>, dy)(t, p.x, p.y, p.z).ax;
      });
    }

Replacement in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 219 [8.16301]

B:BD[4.16503] → [2.15341:15461]

extern "C" void MaxwellToyCarpetX_Dependents1(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_MaxwellToyCarpetX_Dependents1;

[4.16503]

[4.16591]

extern "C" void MaxwellToyCarpetX_InitializeFields(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_MaxwellToyCarpetX_InitializeFields;

Insertion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 226 [8.16301]
[4.16751]
[4.16751]
```
  const Loop::GF3D<const CCTK_REAL, 0, 0, 0> phi_(cctkGH, phi);
```

Insertion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 233 [8.16301]

[4.16939]

  const Loop::GF3D<CCTK_REAL, 1, 0, 0> ex_(cctkGH, ex);
  const Loop::GF3D<CCTK_REAL, 0, 1, 0> ey_(cctkGH, ey);
  const Loop::GF3D<CCTK_REAL, 0, 0, 1> ez_(cctkGH, ez);

Insertion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 241 [8.16301]

[4.17108]

  Loop::loop_int<1, 0, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    ex_(p.i, p.j, p.k) = (phi_(p.i + 1, p.j, p.k) - phi_(p.i, p.j, p.k)) / dx;
  });
  Loop::loop_int<0, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    ey_(p.i, p.j, p.k) = (phi_(p.i, p.j + 1, p.k) - phi_(p.i, p.j, p.k)) / dy;
  });
  Loop::loop_int<0, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    ez_(p.i, p.j, p.k) = (phi_(p.i, p.j, p.k + 1) - phi_(p.i, p.j, p.k)) / dz;
  });

Insertion in MaxwellToyCarpetX/src/maxwelltoy.cxx at line 638 [8.16301]

[4.20491]

[8.22984]

}
extern "C" void MaxwellToyCarpetX_Dependents1(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_MaxwellToyCarpetX_Dependents1;
  DECLARE_CCTK_PARAMETERS;
  const CCTK_REAL dx = CCTK_DELTA_SPACE(0);
  const CCTK_REAL dy = CCTK_DELTA_SPACE(1);
  const CCTK_REAL dz = CCTK_DELTA_SPACE(2);
  const Loop::GF3D<const CCTK_REAL, 1, 0, 0> ax_(cctkGH, ax);
  const Loop::GF3D<const CCTK_REAL, 0, 1, 0> ay_(cctkGH, ay);
  const Loop::GF3D<const CCTK_REAL, 0, 0, 1> az_(cctkGH, az);
  const Loop::GF3D<CCTK_REAL, 0, 1, 1> bx_(cctkGH, bx);
  const Loop::GF3D<CCTK_REAL, 1, 0, 1> by_(cctkGH, by);
  const Loop::GF3D<CCTK_REAL, 1, 1, 0> bz_(cctkGH, bz);
  Loop::loop_int<0, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    bx_(p.i, p.j, p.k) = (az_(p.i, p.j + 1, p.k) - az_(p.i, p.j, p.k)) / dy -
                         (ay_(p.i, p.j, p.k + 1) - ay_(p.i, p.j, p.k)) / dz;
  });
  Loop::loop_int<1, 0, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    by_(p.i, p.j, p.k) = (ax_(p.i, p.j, p.k + 1) - ax_(p.i, p.j, p.k)) / dz -
                         (az_(p.i + 1, p.j, p.k) - az_(p.i, p.j, p.k)) / dx;
  });
  Loop::loop_int<1, 1, 0>(cctkGH, [&](const Loop::PointDesc &p) {
    bz_(p.i, p.j, p.k) = (ay_(p.i + 1, p.j, p.k) - ay_(p.i, p.j, p.k)) / dx -
                         (ax_(p.i, p.j + 1, p.k) - ax_(p.i, p.j, p.k)) / dy;
  });

MaxwellToyCarpetX: Modernize thorn; calculate initial fields from potentials

Dependencies

In channels

Change contents

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 2 [3.30]

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 5 [3.30]

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 21 [3.30]

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 23 [3.30]

Replacement in MaxwellToyCarpetX/par/maxwell-gaussian.par at line 27 [3.30]