Zig-Zag Deposit + Cartesian Coordinate Push Without Explicit Coordinate Conversion

cmake/report.cmake

@@ -155,7 +155,6 @@ if(${Kokkos_DEVICES} MATCHES "CUDA")
     COMMAND bash -c ${cmd}
     OUTPUT_VARIABLE CUDACOMP_VERSION
     OUTPUT_STRIP_TRAILING_WHITESPACE)
-  message(STATUS "CUDACOMP: ${CUDACOMP_VERSION}")
   string(
     APPEND
     REPORT_TEXT

setups/srpic/shock/pgen.hpp

@@ -92,7 +92,6 @@ namespace user {
     // magnetic field properties
     real_t        Btheta, Bphi, Bmag;
     InitFields<D> init_flds;
-    BCFields<D>   bc_flds;
 
     inline PGen(const SimulationParams& p, const Metadomain<S, M>& global_domain)
       : arch::ProblemGenerator<S, M> { p }
@@ -105,7 +104,6 @@ namespace user {
       , Btheta { p.template get<real_t>("setup.Btheta", ZERO) }
       , Bphi { p.template get<real_t>("setup.Bphi", ZERO) }
       , init_flds { Bmag, Btheta, Bphi, drift_ux }
-      , bc_flds { Bmag, Btheta, Bphi, drift_ux }
       , filling_fraction { p.template get<real_t>("setup.filling_fraction", 1.0) }
       , injector_velocity { p.template get<real_t>("setup.injector_velocity", 1.0) }
       , injection_start { p.template get<real_t>("setup.injection_start", 0.0) }
@@ -114,7 +112,7 @@ namespace user {
 
     inline PGen() {}
 
-    auto MatchFields(real_t time) const -> BCFields<D> {
+    auto MatchFields(real_t time) const -> InitFields<D> {
       return init_flds;
     }
 
@@ -178,8 +176,7 @@ namespace user {
       const auto energy_dist = arch::TwoTemperatureMaxwellian<S, M>(
         local_domain.mesh.metric,
         local_domain.random_pool,
-        { temperature_ratio * temperature *
-            (local_domain.species[2].mass() / local_domain.species[1].mass()),
+        { temperature_ratio * temperature * local_domain.species[1].mass() ,
           temperature },
         { 1, 2 },
         -drift_ux,
@@ -321,8 +318,7 @@ namespace user {
       const auto energy_dist = arch::TwoTemperatureMaxwellian<S, M>(
         domain.mesh.metric,
         domain.random_pool,
-        { temperature_ratio * temperature *
-            (domain.species[2].mass() / domain.species[1].mass()),
+        { temperature_ratio * temperature * domain.species[1].mass(),
           temperature },
         { 1, 2 },
         -drift_ux,

setups/tests/deposit/deposit-mink.toml

@@ -0,0 +1,71 @@
+[simulation]
+  name    = "deposit-test-mink"
+  engine  = "srpic"
+  runtime = 5.0
+
+[grid]
+  resolution = [32, 32]
+  extent     = [[0.0, 1.0], [0.0, 1.0]]
+
+  [grid.metric]
+    metric = "minkowski"
+
+  [grid.boundaries]
+    fields    = [["PERIODIC"], ["PERIODIC"]]
+    particles = [["PERIODIC"], ["PERIODIC"]]
+
+[scales]
+  larmor0    = 0.1
+  skindepth0 = 0.1
+
+[algorithms]
+  current_filters = 4
+
+  [algorithms.timestep]
+    CFL = 0.5
+
+[particles]
+  ppc0 = 10.0
+
+  [[particles.species]]
+    label    = "e-"
+    mass     = 1.0
+    charge   = -1.0
+    maxnpart = 1e2
+
+  [[particles.species]]
+    label    = "e+"
+    mass     = 1.0
+    charge   = 1.0
+    maxnpart = 1e2
+
+[setup]
+  x1_e  = [0.25]
+  x2_e  = [0.85]
+  x3_e  = [0.33]
+  ux1_e = [0.6]
+  ux2_e = [-0.3]
+  ux3_e = [-0.2]
+
+  x1_i  = [0.25]
+  x2_i  = [0.85]
+  x3_i  = [0.33]
+  ux1_i = [-0.6]
+  ux2_i = [0.3]
+  ux3_i = [0.2]
+
+[output]
+  format   = "hdf5"
+  interval = 5
+
+  [output.fields]
+    quantities = ["N_1", "N_2", "E", "B", "J"]
+
+  [output.particles]
+    enable = false
+
+  [output.spectra]
+    enable = false
+
+[checkpoint]
+  keep = 0

setups/tests/deposit/deposit-sr.toml

@@ -0,0 +1,71 @@
+[simulation]
+  name    = "deposit-sr"
+  engine  = "srpic"
+  runtime = 10.0
+
+[grid]
+  resolution = [64, 64]
+  extent     = [[1.0, 5.0]]
+
+  [grid.metric]
+    metric = "qspherical"
+
+  [grid.boundaries]
+    fields    = [["FIXED", "FIXED"]]
+    particles = [["REFLECT", "REFLECT"]]
+
+[scales]
+  larmor0    = 0.1
+  skindepth0 = 0.1
+
+[algorithms]
+  current_filters = 4
+
+  [algorithms.timestep]
+    CFL = 0.5
+
+[particles]
+  ppc0 = 10.0
+
+  [[particles.species]]
+    label    = "e-"
+    mass     = 1.0
+    charge   = -1.0
+    maxnpart = 1e2
+
+  [[particles.species]]
+    label    = "e+"
+    mass     = 1.0
+    charge   = 1.0
+    maxnpart = 1e2
+
+[setup]
+  x1_e  = [2.25]
+  x2_e  = [1.25]
+  phi_e = [0.0]
+  ux1_e = [0.6]
+  ux2_e = [-0.3]
+  ux3_e = [-0.2]
+
+  x1_i  = [2.25]
+  x2_i  = [1.25]
+  phi_i = [0.0]
+  ux1_i = [-0.6]
+  ux2_i = [0.3]
+  ux3_i = [0.2]
+
+[output]
+  format   = "hdf5"
+  interval = 5
+
+  [output.fields]
+    quantities = ["N_1", "N_2", "E", "B", "J"]
+
+  [output.particles]
+    enable = false
+
+  [output.spectra]
+    enable = false
+
+[checkpoint]
+  keep = 0