Add BC convenience interface for single plane

ORNL · Aug 8, 2024 · 111f2fa · 111f2fa
1 parent b50761c
commit 111f2fa
Show file tree

Hide file tree

Showing 2 changed files with 24 additions and 3 deletions.
diff --git a/examples/mechanics/kalthoff_winkler.cpp b/examples/mechanics/kalthoff_winkler.cpp
@@ -101,10 +101,8 @@ void kalthoffWinklerExample( const std::string filename )
     CabanaPD::RegionBoundary<CabanaPD::RectangularPrism> plane(
         x_bc - dx, x_bc + dx * 1.25, y_prenotch1 - dx * 0.25,
         y_prenotch2 + dx * 0.25, -thickness, thickness );
-    std::vector<CabanaPD::RegionBoundary<CabanaPD::Rectilinear>> planes = {
-        plane };
     auto bc = createBoundaryCondition( CabanaPD::ForceValueBCTag{}, 0.0,
-                                       exec_space{}, *particles, planes );
+                                       exec_space{}, *particles, plane );
 
     // ====================================================
     //            Custom particle initialization

diff --git a/src/CabanaPD_Boundary.hpp b/src/CabanaPD_Boundary.hpp
@@ -357,6 +357,29 @@ auto createBoundaryCondition( UserFunctor user_functor, ExecSpace exec_space,
         bc_indices, user_functor, force_update );
 }
 
+template <class BoundaryType, class BCTag, class ExecSpace, class Particles>
+auto createBoundaryCondition( BCTag tag, const double value,
+                              ExecSpace exec_space, Particles particles,
+                              BoundaryType plane,
+                              const double initial_guess = 0.5 )
+{
+    std::vector<BoundaryType> plane_vec = { plane };
+    return createBoundaryCondition( tag, value, exec_space, particles,
+                                    plane_vec, initial_guess );
+}
+
+template <class UserFunctor, class BoundaryType, class ExecSpace,
+          class Particles>
+auto createBoundaryCondition( UserFunctor user_functor, ExecSpace exec_space,
+                              Particles particles, BoundaryType plane,
+                              const bool force_update,
+                              const double initial_guess = 0.5 )
+{
+    std::vector<BoundaryType> plane_vec = { plane };
+    return createBoundaryCondition( user_functor, exec_space, particles,
+                                    plane_vec, force_update, initial_guess );
+}
+
 } // namespace CabanaPD
 
 #endif