Support SoA-Only Particles

Src/AmrCore/AMReX_AmrParticles.H

@@ -10,10 +10,10 @@
 
 namespace amrex {
 
-template <int NStructReal, int NStructInt, int NArrayReal, int NArrayInt,
+template <typename ParticleType, int NArrayReal, int NArrayInt,
           template<class> class Allocator>
 void
-ParticleContainer<NStructReal, NStructInt, NArrayReal, NArrayInt, Allocator>
+ParticleContainer_impl<ParticleType, NArrayReal, NArrayInt, Allocator>
 ::AssignDensity (int rho_index,
                  Vector<std::unique_ptr<MultiFab> >& mf_to_be_filled,
                  int lev_min, int ncomp, int finest_level, int ngrow) const
@@ -249,40 +249,47 @@ ParticleToMesh (PC const& pc, const Vector<MultiFab*>& mf,
     }
 }
 
-template <int NStructReal, int NStructInt=0, int NArrayReal=0, int NArrayInt=0,
+template <typename T_ParticleType, int NArrayReal=0, int NArrayInt=0,
           template<class> class Allocator=DefaultAllocator>
-class AmrParticleContainer // NOLINT
-    : public ParticleContainer<NStructReal, NStructInt, NArrayReal, NArrayInt, Allocator>
+class AmrParticleContainer_impl // NOLINT(cppcoreguidelines-virtual-class-destructor)
+    : public ParticleContainer_impl<T_ParticleType, NArrayReal, NArrayInt, Allocator>
 {
 
 public:
 
-    typedef Particle<NStructReal, NStructInt> ParticleType;
+    using ParticleType = T_ParticleType;
 
-    AmrParticleContainer () = default;
+    AmrParticleContainer_impl ()
+        : ParticleContainer_impl<ParticleType, NArrayReal, NArrayInt, Allocator>()
+    {
+    }
 
-    AmrParticleContainer (AmrCore* amr_core)
-        : ParticleContainer<NStructReal, NStructInt, NArrayReal, NArrayInt, Allocator>(amr_core->GetParGDB())
+    AmrParticleContainer_impl (AmrCore* amr_core)
+        : ParticleContainer_impl<ParticleType, NArrayReal, NArrayInt, Allocator>(amr_core->GetParGDB())
     {
     }
 
-    AmrParticleContainer (const Vector<Geometry>            & geom,
-                          const Vector<DistributionMapping> & dmap,
-                          const Vector<BoxArray>            & ba,
-                          const Vector<int>                 & rr)
-        : ParticleContainer<NStructReal, NStructInt, NArrayReal, NArrayInt, Allocator>(geom, dmap, ba, rr)
+    AmrParticleContainer_impl (const Vector<Geometry>            & geom,
+                               const Vector<DistributionMapping> & dmap,
+                               const Vector<BoxArray>            & ba,
+                               const Vector<int>                 & rr)
+        : ParticleContainer_impl<ParticleType, NArrayReal, NArrayInt, Allocator>(geom, dmap, ba, rr)
     {
     }
 
-    ~AmrParticleContainer () override = default;
+    ~AmrParticleContainer_impl () override = default;
 
-    AmrParticleContainer ( const AmrParticleContainer &) = delete;
-    AmrParticleContainer& operator= ( const AmrParticleContainer & ) = delete;
+    AmrParticleContainer_impl ( const AmrParticleContainer_impl &) = delete;
+    AmrParticleContainer_impl& operator= ( const AmrParticleContainer_impl & ) = delete;
 
-    AmrParticleContainer ( AmrParticleContainer && ) noexcept = default;
-    AmrParticleContainer& operator= ( AmrParticleContainer && ) noexcept = default;
+    AmrParticleContainer_impl ( AmrParticleContainer_impl && ) noexcept = default;
+    AmrParticleContainer_impl& operator= ( AmrParticleContainer_impl && ) noexcept = default;
 };
 
+template <int T_NStructReal, int T_NStructInt=0, int T_NArrayReal=0, int T_NArrayInt=0,
+          template<class> class Allocator=DefaultAllocator>
+using AmrParticleContainer = AmrParticleContainer_impl<Particle<T_NStructReal, T_NStructInt>, T_NArrayReal, T_NArrayInt, Allocator>;
+
 class AmrTracerParticleContainer
     : public TracerParticleContainer
 {

Src/Base/AMReX_Array.H

@@ -897,3 +897,4 @@ namespace amrex
 }
 
 #endif
+

Src/Base/AMReX_Geometry.H

@@ -194,6 +194,7 @@ public:
     [[nodiscard]] GpuArray<ParticleReal,AMREX_SPACEDIM> ProbLoArrayInParticleReal () const noexcept {
         return roundoff_lo;
     }
+
     [[nodiscard]] GpuArray<ParticleReal,AMREX_SPACEDIM> ProbHiArrayInParticleReal () const noexcept {
         return roundoff_hi;
     }

Src/Base/AMReX_TypeTraits.H

@@ -54,17 +54,17 @@ namespace amrex
     struct IsMultiFabIterator : public std::is_base_of<MFIter, T>::type {};
 
 #ifdef AMREX_PARTICLES
-    template <bool is_const, int NStructReal, int NStructInt, int NArrayReal, int NArrayInt,
-              template<class> class Allocator>
-    class ParIterBase;
+    // template <bool is_const, int NStructReal, int NStructInt, int NArrayReal, int NArrayInt,
+    //           template<class> class Allocator>
+    // class ParIterBase;
 
-    template <int NStructReal, int NStructInt, int NArrayReal, int NArrayInt,
-              template<class> class Allocator>
-    class ParIter;
+    // template <int NStructReal, int NStructInt, int NArrayReal, int NArrayInt,
+    //           template<class> class Allocator>
+    // class ParIter;
 
-    template <int NStructReal, int NStructInt, int NArrayReal, int NArrayInt,
-              template<class> class Allocator>
-    class ParConstIter;
+    // template <int NStructReal, int NStructInt, int NArrayReal, int NArrayInt,
+    //           template<class> class Allocator>
+    // class ParConstIter;
 
     class ParticleContainerBase;
 

Src/EB/AMReX_EB2.H

@@ -134,7 +134,6 @@ void Build (const Geometry& geom,
             bool extend_domain_face = ExtendDomainFace(),
             int num_coarsen_opt = NumCoarsenOpt());
 
-
 void BuildFromChkptFile (std::string const& fname,
                          const Geometry& geom,
                          int required_coarsening_level,