Enabled surface model for the async mode. (apt-sim#361)

This just enables the surface model replacing the solid one when
ADEPT_USE_SURF is ON

Author: agheata

include/AdePT/core/AsyncAdePTTransport.cuh

@@ -28,6 +28,9 @@
 #ifdef VECGEOM_ENABLE_CUDA
 #include <VecGeom/backend/cuda/Interface.h>
 #endif
+#ifdef ADEPT_USE_SURF
+#include <VecGeom/surfaces/cuda/BrepCudaManager.h>
+#endif
 
 #include <G4HepEmState.hh>
 #include <G4HepEmData.hh>
@@ -392,6 +395,22 @@ __global__ void InitSlotManagers(SlotManager *mgr, std::size_t N)
 /// WIP: Free functions implementing the CUDA parts
 namespace async_adept_impl {
 
+void CopySurfaceModelToGPU()
+{
+// Copy surface data to GPU
+#ifdef ADEPT_USE_SURF
+#ifdef ADEPT_USE_SURF_SINGLE
+  using SurfData        = vgbrep::SurfData<float>;
+  using BrepCudaManager = vgbrep::BrepCudaManager<float>;
+#else
+  using SurfData        = vgbrep::SurfData<double>;
+  using BrepCudaManager = vgbrep::BrepCudaManager<double>;
+#endif
+  BrepCudaManager::Instance().TransferSurfData(SurfData::Instance());
+  printf("== Surface data transferred to GPU\n");
+#endif
+}
+
 G4HepEmState *InitG4HepEm()
 {
   auto state = new G4HepEmState;

include/AdePT/core/AsyncAdePTTransport.icc

@@ -5,8 +5,12 @@
 
 #include <AdePT/integration/AdePTGeant4Integration.hh>
 
+#include <VecGeom/base/Stopwatch.h>
 #include <VecGeom/management/BVHManager.h>
-#include "VecGeom/management/GeoManager.h"
+#include <VecGeom/management/GeoManager.h>
+#ifdef ADEPT_USE_SURF
+#include <VecGeom/surfaces/BrepHelper.h>
+#endif
 
 #include <G4RunManager.hh>
 #include <G4Threading.hh>
@@ -39,6 +43,7 @@
 
 /// Forward declarations
 namespace async_adept_impl {
+void CopySurfaceModelToGPU();
 bool InitializeField(double);
 bool InitializeGeneralField(GeneralMagneticField &);
 void FreeGeneralField(GeneralMagneticField &);
@@ -200,12 +205,32 @@ bool AsyncAdePTTransport<IntegrationLayer>::InitializeGeometry(const vecgeom::cx
     std::cout << "CUDA Device heap limit: " << fCUDAHeapLimit << "\n";
     cudaDeviceSetLimit(cudaLimitMallocHeapSize, fCUDAHeapLimit);
   }
+  bool success = true;
+#ifdef ADEPT_USE_SURF
+#ifdef ADEPT_USE_SURF_SINGLE
+  using SurfData   = vgbrep::SurfData<float>;
+  using BrepHelper = vgbrep::BrepHelper<float>;
+#else
+  using SurfData   = vgbrep::SurfData<double>;
+  using BrepHelper = vgbrep::BrepHelper<double>;
+#endif
+  vecgeom::Stopwatch timer;
+  timer.Start();
+  if (!BrepHelper::Instance().Convert()) return 1;
+  BrepHelper::Instance().PrintSurfData();
+  std::cout << "== Conversion to surface model done in " << timer.Stop() << " [s]\n";
+  // Upload only navigation table to the GPU
+  cudaManager.SynchronizeNavigationTable();
+  async_adept_impl::CopySurfaceModelToGPU();
+#else
+  // Upload solid geometry to GPU.
   cudaManager.LoadGeometry(world);
   auto world_dev = cudaManager.Synchronize();
+  success        = world_dev != nullptr;
   // Initialize BVH
   InitBVH();
-
-  return (world_dev != nullptr);
+#endif
+  return success;
 }
 
 template <typename IntegrationLayer>