Function initializeClusteredLts

src/Initializer/InitProcedure/InitModel.cpp:271–323

}

void initializeClusteredLts(LtsInfo& ltsInfo, seissol::SeisSol& seissolInstance) {
  const auto& seissolParams = seissolInstance.getSeisSolParameters();

  assert(seissolParams.timeStepping.lts.getRate() > 0);

  if (seissolParams.timeStepping.lts.getRate() == 1) {
    seissolInstance.getLtsLayout().deriveLayout(TimeClustering::Single, 1);
  } else {
    seissolInstance.getLtsLayout().deriveLayout(TimeClustering::MultiRate,
                                                seissolParams.timeStepping.lts.getRate());
  }

  seissolInstance.getLtsLayout().getMeshStructure(ltsInfo.meshStructure);
  seissolInstance.getLtsLayout().getCrossClusterTimeStepping(ltsInfo.timeStepping);

  seissolInstance.getMemoryManager().initializeFrictionLaw();

  unsigned* numberOfDRCopyFaces = nullptr;
  unsigned* numberOfDRInteriorFaces = nullptr;

  seissolInstance.getLtsLayout().getDynamicRuptureInformation(
      ltsInfo.ltsMeshToFace, numberOfDRCopyFaces, numberOfDRInteriorFaces);

  seissolInstance.getMemoryManager().fixateLtsTree(ltsInfo.timeStepping,
                                                   ltsInfo.meshStructure,
                                                   numberOfDRCopyFaces,
                                                   numberOfDRInteriorFaces,
                                                   seissolParams.model.plasticity);

  seissolInstance.getMemoryManager().setLtsToFace(ltsInfo.ltsMeshToFace);

  delete[] numberOfDRCopyFaces;
  delete[] numberOfDRInteriorFaces;

  const auto& ltsTree = seissolInstance.getMemoryManager().getLtsTree();
  const auto& lts = seissolInstance.getMemoryManager().getLts();

  unsigned* ltsToMesh = nullptr;
  unsigned numberOfMeshCells = 0;

  seissolInstance.getLtsLayout().getCellInformation(
      ltsTree->var(lts->cellInformation), ltsToMesh, numberOfMeshCells);

  // TODO(David): move all of this method to the MemoryManager
  seissolInstance.getMemoryManager().getLtsLutUnsafe().createLuts(
      ltsTree, ltsToMesh, numberOfMeshCells);

  delete[] ltsToMesh;

  seissol::initializer::time_stepping::deriveLtsSetups(ltsInfo.timeStepping.numberOfLocalClusters,
                                                       ltsInfo.meshStructure,
                                                       ltsTree->var(lts->cellInformation));
}

void initializeMemoryLayout(LtsInfo& ltsInfo, seissol::SeisSol& seissolInstance) {
  const auto& seissolParams = seissolInstance.getSeisSolParameters();