Method initMesh

src/Initializer/InitProcedure/InitMesh.cpp:252–334

} // namespace

void seissol::initializer::initprocedure::initMesh(seissol::SeisSol& seissolInstance) {
  SCOREP_USER_REGION("init_mesh", SCOREP_USER_REGION_TYPE_FUNCTION);

  const auto& seissolParams = seissolInstance.getSeisSolParameters();
  const auto commRank = seissol::MPI::mpi.rank();
  const auto commSize = seissol::MPI::mpi.size();

  logInfo(commRank) << "Begin init mesh.";

  // Call the pre mesh initialization hook
  seissol::Modules::callHook<ModuleHook::PreMesh>();

  const auto meshFormat = seissolParams.mesh.meshFormat;

  logInfo(commRank) << "Mesh file:" << seissolParams.mesh.meshFileName;

  seissol::Stopwatch watch;
  watch.start();

  const std::string realMeshFileName = seissolParams.mesh.meshFileName;
  [[maybe_unused]] bool addNC = true;
  if (realMeshFileName.size() >= 3) {
    const auto lastCharacters = realMeshFileName.substr(realMeshFileName.size() - 3);
    addNC = lastCharacters != ".nc";
  }

  switch (meshFormat) {
  case seissol::initializer::parameters::MeshFormat::Netcdf: {
#if USE_NETCDF
    const auto realMeshFileNameNetcdf = [&]() {
      if (addNC) {
        const auto newRealMeshFileName = realMeshFileName + ".nc";
        logInfo(commRank)
            << "The Netcdf file extension \".nc\" has been appended. Updated mesh file name:"
            << newRealMeshFileName;
        return newRealMeshFileName;
      } else {
        // (suppress preference for return move)
        // NOLINTNEXTLINE
        return realMeshFileName;
      }
    }();
    seissolInstance.setMeshReader(
        new seissol::geometry::NetcdfReader(commRank, commSize, realMeshFileNameNetcdf.c_str()));
#else
    logError()
        << "Tried to load a Netcdf mesh, however this build of SeisSol is not linked to Netcdf.";
#endif
    break;
  }
  case seissol::initializer::parameters::MeshFormat::PUML: {
    readMeshPUML(seissolParams, seissolInstance);
    break;
  }
  case seissol::initializer::parameters::MeshFormat::CubeGenerator: {
    readCubeGenerator(seissolParams, seissolInstance);
    break;
  }