Add test_bnd_prt_inflow with nascent Inflow class

src/include/particle.h

@@ -37,3 +37,11 @@ struct Inject
 
 } // namespace particle
 } // namespace psc
+
+inline std::ostream& operator<<(std::ostream& os,
+                                const psc::particle::Inject& inj)
+{
+  os << "Inject{x=" << inj.x << ", u=" << inj.u << ", w=" << inj.w
+     << ", kind=" << inj.kind << ", tag=" << inj.tag << "}";
+  return os;
+}

src/include/setup_particles.hxx

@@ -9,10 +9,10 @@
 
 struct psc_particle_npt
 {
-  int kind;    ///< particle kind
-  double n;    ///< density
-  double p[3]; ///< momentum
-  double T[3]; ///< temperature
+  int kind;  ///< particle kind
+  double n;  ///< density
+  Double3 p; ///< momentum
+  Double3 T; ///< temperature
   psc::particle::Tag tag;
 };
 
@@ -241,6 +241,18 @@ struct SetupParticles
       });
   }
 
+  // ----------------------------------------------------------------------
+  // getWeight
+
+  real_t getWeight(const psc_particle_np& np, int n_in_cell)
+  {
+    if (fractional_n_particles_per_cell) {
+      return 1.;
+    } else {
+      return np.n / (n_in_cell * norm_.cori);
+    }
+  }
+
   // ----------------------------------------------------------------------
   // setupParticles
 
@@ -269,12 +281,7 @@ struct SetupParticles
       op_cellwise(grid, p, init_np,
                   [&](int n_in_cell, psc_particle_np& np, Double3& pos) {
                     for (int cnt = 0; cnt < n_in_cell; cnt++) {
-                      real_t wni;
-                      if (fractional_n_particles_per_cell) {
-                        wni = 1.;
-                      } else {
-                        wni = np.n / (n_in_cell * norm_.cori);
-                      }
+                      real_t wni = getWeight(np, n_in_cell);
                       auto prt = setupParticle(np, pos, wni);
                       injector(prt);
                     }

src/libpsc/tests/CMakeLists.txt

@@ -24,6 +24,7 @@ add_psc_test(test_output_particles)
 add_psc_test(test_mfields)
 add_psc_test(test_mfields_cuda)
 add_psc_test(test_bnd)
+add_psc_test(test_bnd_prt_inflow)
 add_psc_test(test_deposit)
 add_psc_test(test_current_deposition)
 add_psc_test(test_push_particles)