Support Blueprint mesh in shaping

src/axom/core/execution/execution_space.hpp

@@ -8,6 +8,7 @@
 
 #include "axom/config.hpp"
 #include "axom/core/memory_management.hpp"
+#include "axom/core/execution/runtime_policy.hpp"
 
 /*!
  * \file execution_space.hpp
@@ -88,6 +89,18 @@ struct execution_space
   static constexpr bool onDevice() noexcept { return false; }
   static constexpr char* name() noexcept { return (char*)"[UNDEFINED]"; }
   static int allocatorID() noexcept { return axom::INVALID_ALLOCATOR_ID; }
+  static constexpr runtime_policy::Policy runtimePolicy() noexcept
+  {
+    return runtime_policy::Policy::seq;
+  }
+  static bool usesMemorySpace(axom::MemorySpace m) noexcept
+  {
+    return m == memory_space;
+  }
+  static bool usesAllocId(int allocId) noexcept
+  {
+    return usesMemorySpace(axom::detail::getAllocatorSpace(allocId));
+  }
 };
 
 }  // namespace axom

src/axom/core/execution/internal/cuda_exec.hpp

@@ -65,6 +65,18 @@ struct execution_space<CUDA_EXEC<BLOCK_SIZE, SYNCHRONOUS>>
   {
     return axom::getUmpireResourceAllocatorID(umpire::resource::Device);
   }
+  static constexpr runtime_policy::Policy runtimePolicy() noexcept
+  {
+    return runtime_policy::Policy::cuda;
+  }
+  static bool usesMemorySpace(axom::MemorySpace m) noexcept
+  {
+    return m == memory_space || m == MemorySpace::Unified;
+  }
+  static bool usesAllocId(int allocId) noexcept
+  {
+    return usesMemorySpace(axom::detail::getAllocatorSpace(allocId));
+  }
 };
 
 /*!

src/axom/core/execution/internal/hip_exec.hpp

@@ -63,6 +63,18 @@ struct execution_space<HIP_EXEC<BLOCK_SIZE, SYNCHRONOUS>>
   {
     return axom::getUmpireResourceAllocatorID(umpire::resource::Device);
   }
+  static constexpr runtime_policy::Policy runtimePolicy() noexcept
+  {
+    return runtime_policy::Policy::hip;
+  }
+  static bool usesMemorySpace(axom::MemorySpace m) noexcept
+  {
+    return m == memory_space || m == MemorySpace::Unified;
+  }
+  static bool usesAllocId(int allocId) noexcept
+  {
+    return usesMemorySpace(axom::detail::getAllocatorSpace(allocId));
+  }
 };
 
 /*!

src/axom/core/execution/internal/omp_exec.hpp

@@ -60,6 +60,22 @@ struct execution_space<OMP_EXEC>
     return axom::getDefaultAllocatorID();
 #endif
   }
+  static constexpr runtime_policy::Policy runtimePolicy() noexcept
+  {
+    return runtime_policy::Policy::omp;
+  }
+  static bool usesMemorySpace(axom::MemorySpace m) noexcept
+  {
+    return m == memory_space
+#ifdef AXOM_USE_UMPIRE
+      || m == MemorySpace::Unified
+#endif
+      ;
+  }
+  static bool usesAllocId(int allocId) noexcept
+  {
+    return usesMemorySpace(axom::detail::getAllocatorSpace(allocId));
+  }
 };
 
 }  // namespace axom

src/axom/core/execution/internal/seq_exec.hpp

@@ -69,6 +69,22 @@ struct execution_space<SEQ_EXEC>
     return axom::getDefaultAllocatorID();
 #endif
   }
+  static constexpr runtime_policy::Policy runtimePolicy() noexcept
+  {
+    return runtime_policy::Policy::seq;
+  }
+  static bool usesMemorySpace(axom::MemorySpace m) noexcept
+  {
+    return m == memory_space
+#ifdef AXOM_USE_UMPIRE
+      || m == MemorySpace::Unified
+#endif
+      ;
+  }
+  static bool usesAllocId(int allocId) noexcept
+  {
+    return usesMemorySpace(axom::detail::getAllocatorSpace(allocId));
+  }
 };
 
 }  // namespace axom

src/axom/core/memory_management.hpp

@@ -110,6 +110,27 @@ inline int getDefaultAllocatorID()
 #endif
 }
 
+/*!
+ * \brief Get the allocator id from which data has been allocated.
+ * \return Allocator id.  If Umpire doesn't have an allocator for
+ * the pointer, or Axom wasn't configured with Umpire, return 0.
+ *
+ * \pre ptr has a valid pointer value.
+ */
+inline int getAllocatorIDFromPointer(const void* ptr)
+{
+#ifdef AXOM_USE_UMPIRE
+  umpire::ResourceManager& rm = umpire::ResourceManager::getInstance();
+  if(rm.hasAllocator(const_cast<void*>(ptr)))
+  {
+    umpire::Allocator allocator = rm.getAllocator(const_cast<void*>(ptr));
+    return allocator.getId();
+  }
+#endif
+  AXOM_UNUSED_VAR(ptr);
+  return 0;
+}
+
 /*!
  * \brief Allocates a chunk of memory of type T.
  *

src/axom/mint/mesh/Mesh.cpp

@@ -123,22 +123,19 @@ Mesh::Mesh(sidre::Group* group, const std::string& topo)
   m_coordset =
     blueprint::getCoordsetGroup(m_group, getTopologyGroup())->getName();
 
-  SLIC_ERROR_IF(!m_group->hasChildGroup("state"),
-                "root group does not have a state group.");
-
   sidre::Group* state_group = m_group->getGroup("state");
-  SLIC_ERROR_IF(!state_group->hasChildGroup(m_topology),
-                "state group has no " << m_topology << " child group.");
-
-  state_group = state_group->getGroup(m_topology);
-  if(state_group->hasChildView("block_id"))
+  if(state_group != nullptr && state_group->hasChildGroup(m_topology))
   {
-    m_block_idx = state_group->getView("block_id")->getScalar();
-  }
+    state_group = state_group->getGroup(m_topology);
+    if(state_group->hasChildView("block_id"))
+    {
+      m_block_idx = state_group->getView("block_id")->getScalar();
+    }
 
-  if(state_group->hasChildView("partition_id"))
-  {
-    m_part_idx = state_group->getView("partition_id")->getScalar();
+    if(state_group->hasChildView("partition_id"))
+    {
+      m_part_idx = state_group->getView("partition_id")->getScalar();
+    }
   }
 
   SLIC_ERROR_IF(!validMeshType(), "invalid mesh type=" << m_type);

src/axom/multimat/multimat.cpp

@@ -49,6 +49,7 @@ bool AllocatorOnDevice(int allocatorId)
     return true;
   }
 #endif
+  AXOM_UNUSED_VAR(allocatorId);
   return false;
 }
 
@@ -509,6 +510,9 @@ void ScanRelationOffsetsRAJA(const axom::ArrayView<const IndexType> counts,
       }
     });
 #else
+  AXOM_UNUSED_VAR(counts);
+  AXOM_UNUSED_VAR(begins);
+  AXOM_UNUSED_VAR(firstIndices);
   SLIC_ASSERT_MSG(
     false,
     "Calling ScanRelationOffsetsRAJA requires support for RAJA and Umpire.");

src/axom/quest/CMakeLists.txt

@@ -130,16 +130,20 @@ blt_list_append(
 
 blt_list_append( TO quest_depends_on ELEMENTS conduit::conduit IF CONDUIT_FOUND )
 
-if(MFEM_FOUND AND AXOM_ENABLE_KLEE AND AXOM_ENABLE_SIDRE AND AXOM_ENABLE_MFEM_SIDRE_DATACOLLECTION)
+if(AXOM_ENABLE_KLEE AND AXOM_ENABLE_SIDRE AND ((MFEM_FOUND AND AXOM_ENABLE_MFEM_SIDRE_DATACOLLECTION) OR CONDUIT_FOUND))
     list(APPEND quest_headers Shaper.hpp
-                              SamplingShaper.hpp
-                              IntersectionShaper.hpp
-                              DiscreteShape.hpp
-                              detail/shaping/shaping_helpers.hpp)
+                              DiscreteShape.hpp)
     list(APPEND quest_sources Shaper.cpp
-                              DiscreteShape.cpp
-                              detail/shaping/shaping_helpers.cpp)
+                              DiscreteShape.cpp)
     list(APPEND quest_depends_on klee)
+    if(RAJA_FOUND)
+        list(APPEND quest_headers IntersectionShaper.hpp)
+    endif()
+endif()
+
+if(AXOM_ENABLE_KLEE AND AXOM_ENABLE_SIDRE AND (MFEM_FOUND AND AXOM_ENABLE_MFEM_SIDRE_DATACOLLECTION))
+    list(APPEND quest_headers SamplingShaper.hpp detail/shaping/shaping_helpers.hpp)
+    list(APPEND quest_sources detail/shaping/shaping_helpers.cpp)
 endif()
 
 if(C2C_FOUND)

src/axom/quest/DiscreteShape.cpp

@@ -374,7 +374,7 @@ void DiscreteShape::createPlaneRepresentation()
   const auto& plane = geometry.getPlane();
   // Generate a big bounding hex on the positive side of the plane.
   axom::primal::Hexahedron<double, 3> boundingHex;
-  const double len = 1e6;  // Big enough to contain anticipated mesh.
+  const double len = 1e4;  // Big enough to contain anticipated mesh.
   // We should compute based on the mesh.
   // clang-format off
   boundingHex[0] = Point3D{0.0, -len, -len};