Function makePointSourceCluster

src/SourceTerm/Manager.cpp:314–354

}

auto makePointSourceCluster(const ClusterMapping& mapping,
                            const PointSources& sources,
                            const unsigned* meshIds,
                            seissol::initializer::LTSTree* ltsTree,
                            seissol::initializer::LTS* lts,
                            seissol::initializer::Lut* ltsLut)
    -> seissol::kernels::PointSourceClusterPair {
  auto hostData = std::pair<std::shared_ptr<ClusterMapping>, std::shared_ptr<PointSources>>(
      std::make_shared<ClusterMapping>(mapping), std::make_shared<PointSources>(sources));

#if defined(ACL_DEVICE) && !defined(MULTIPLE_SIMULATIONS)
  using GpuImpl = seissol::kernels::PointSourceClusterOnDevice;

  auto deviceData =
      [&]() -> std::pair<std::shared_ptr<ClusterMapping>, std::shared_ptr<PointSources>> {
    if (useUSM()) {
      return hostData;
    } else {
      constexpr auto GpuMemkind = seissol::memory::Memkind::DeviceGlobalMemory;
      auto predeviceClusterMapping = mapping;
      mapClusterToMesh(predeviceClusterMapping,
                       meshIds,
                       ltsTree,
                       lts,
                       ltsLut,
                       seissol::initializer::AllocationPlace::Device);
      auto deviceClusterMapping =
          std::make_shared<ClusterMapping>(predeviceClusterMapping, GpuMemkind);
      auto devicePointSources = std::make_shared<PointSources>(sources, GpuMemkind);
      return {deviceClusterMapping, devicePointSources};
    }
  }();
#else
  using GpuImpl = seissol::kernels::PointSourceClusterOnHost;
  auto deviceData = hostData;
#endif

  return seissol::kernels::PointSourceClusterPair{
      std::make_unique<kernels::PointSourceClusterOnHost>(hostData.first, hostData.second),
      std::make_unique<GpuImpl>(deviceData.first, deviceData.second)};
}

auto loadSourcesFromFSRM(const char* fileName,
                         const seissol::geometry::MeshReader& mesh,