Merge pull request #123 from ajnonaka/random_fix

fix random number signature

Author: ajnonaka

Source/Initialization/PlasmaInjector.cpp

@@ -691,7 +691,11 @@ amrex::XDim3 PlasmaInjector::getMomentum (amrex::Real x,
                                           amrex::Real y,
                                           amrex::Real z) const noexcept
 {
+#ifdef AMREX_USE_GPU
+    return h_inj_mom->getMomentum(x, y, z, amrex::RandomEngine{nullptr}); // gamma*beta
+#else
     return h_inj_mom->getMomentum(x, y, z, amrex::RandomEngine{}); // gamma*beta
+#endif
 }
 
 bool PlasmaInjector::insideBounds (amrex::Real x, amrex::Real y, amrex::Real z) const noexcept

Source/Particles/ParticleBoundaryBuffer.cpp

@@ -251,7 +251,11 @@ void ParticleBoundaryBuffer::gatherParticles (MultiParticleContainer& mypc,
                         amrex::ReduceData<int> reduce_data(reduce_op);
                         {
                           WARPX_PROFILE("ParticleBoundaryBuffer::gatherParticles::count_out_of_bounds");
-                          amrex::RandomEngine rng{};
+#ifdef AMREX_USE_GPU
+                          const amrex::RandomEngine rng{nullptr};
+#else
+                          const amrex::RandomEngine rng{};
+#endif
                           reduce_op.eval(np, reduce_data, [=] AMREX_GPU_HOST_DEVICE (int ip)
                                          { return predicate(ptile_data, ip, rng) ? 1 : 0; });
                         }

Source/Particles/WarpXParticleContainer.cpp

@@ -280,8 +280,13 @@ WarpXParticleContainer::AddNParticles (int /*lev*/,
         }
 
         // Default initialize the other real and integer runtime attributes
+#ifdef AMREX_USE_GPU
+        DefaultInitializeRuntimeAttributes(pinned_tile, nattr_real - 1, nattr_int,
+                                           amrex::RandomEngine{nullptr});
+#else
         DefaultInitializeRuntimeAttributes(pinned_tile, nattr_real - 1, nattr_int,
                                            amrex::RandomEngine{});
+#endif
 
         auto old_np = particle_tile.numParticles();
         auto new_np = old_np + pinned_tile.numParticles();