Method initializeFaceDisplacements

src/Initializer/MemoryManager.cpp:737–769

#endif

void seissol::initializer::MemoryManager::initializeFaceDisplacements()
{
  for (auto& layer : m_ltsTree.leaves(m_lts.faceDisplacements.mask)) {
    if (layer.getBucketSize(m_lts.faceDisplacementsBuffer) == 0) {
      continue;
    }
    real* (*displacements)[4] = layer.var(m_lts.faceDisplacements);
    real* bucket = static_cast<real*>(layer.bucket(m_lts.faceDisplacementsBuffer));
    real* (*displacementsDevice)[4] = layer.var(m_lts.faceDisplacementsDevice);
    real* bucketDevice = static_cast<real*>(layer.bucket(m_lts.faceDisplacementsBuffer, seissol::initializer::AllocationPlace::Device));

#ifdef _OPENMP
#pragma omp parallel for schedule(static) default(none) shared(layer, displacements, bucket, displacementsDevice, bucketDevice)
#endif // _OPENMP
    for (unsigned cell = 0; cell < layer.getNumberOfCells(); ++cell) {
      for (unsigned face = 0; face < 4; ++face) {
        if (displacements[cell][face] != nullptr) {
          // Remove constant part that was added in deriveDisplacementsBucket.
          // We then have the pointer offset that needs to be added to the bucket.
          // The final value of this pointer then points to a valid memory address
          // somewhere in the bucket.
          auto offset = (reinterpret_cast<std::intptr_t>(displacements[cell][face]) - 1);
          displacements[cell][face] = bucket + offset;
          displacementsDevice[cell][face] = bucketDevice + offset;
          for (unsigned dof = 0; dof < tensor::faceDisplacement::size(); ++dof) {
            // zero displacements
            displacements[cell][face][dof] = static_cast<real>(0.0);
          }
        }
      }
    }
  }
}

void seissol::initializer::MemoryManager::initializeMemoryLayout()
{