Method PUMLReader

src/Geometry/PUMLReader.cpp:184–206

} // namespace

seissol::geometry::PUMLReader::PUMLReader(
    const char* meshFile,
    const char* partitioningLib,
    double maximumAllowedTimeStep,
    seissol::initializer::parameters::BoundaryFormat boundaryFormat,
    initializer::time_stepping::LtsWeights* ltsWeights,
    double tpwgt)
    : seissol::geometry::MeshReader(MPI::mpi.rank()), boundaryFormat(boundaryFormat) {
  PUML::TETPUML puml;
  puml.setComm(MPI::mpi.comm());

  read(puml, meshFile);

  generatePUML(puml); // We need to call generatePUML in order to create the dual graph of the mesh
  if (ltsWeights != nullptr) {
    ltsWeights->computeWeights(puml, maximumAllowedTimeStep);
  }
  partition(puml, ltsWeights, tpwgt, meshFile, partitioningLib);

  generatePUML(puml);

  getMesh(puml);
}

void seissol::geometry::PUMLReader::read(PUML::TETPUML& puml, const char* meshFile) {
  SCOREP_USER_REGION("PUMLReader_read", SCOREP_USER_REGION_TYPE_FUNCTION);