FTInterpExtrap.hpp

Go to the documentation of this file.

//
// This file is subject to the terms and conditions defined in
// file 'LICENSE', which is part of this source code package.
//
 
#pragma once
 
#include "finite_elements/utils/ArborXCallbacks.hpp"
#include "finite_elements/utils/ComputeLinearSkin.hpp"
#include "finite_elements/utils/FTUtils.hpp"
#include "finite_elements/utils/SkinProjection.hpp"
 
namespace PACMAN {
 
namespace FiniteElements {
 
template <typename ExecSpace, int_t Dim>
void FTInterExtrap(Transfer<ExecSpace, Dim> &transfer) {
  Kokkos::Profiling::pushRegion("FiniteElement::FTInterpExtrap");
 
  Kokkos::Profiling::pushRegion("Compute target point FE intersection");
  ComputeBoxTargetPointIntersection<ExecSpace, Dim>(transfer);
 
  Kokkos::Profiling::popRegion(); // Compute target point FE intersection
 
  Kokkos::Profiling::pushRegion(
      "Compute target point projection on triangulated skin mesh");
  ComputeLinearSkin<ExecSpace, Dim>(transfer);
  if constexpr (Dim == 1) {
    ComputeProjectionOn1DSkin<ExecSpace>(transfer, true);
  } else if constexpr (Dim == 2) {
    ComputeProjectionOn2DSkin<ExecSpace>(transfer, true);
  } else if constexpr (Dim == 3) {
    ComputeProjectionOn3DSkin<ExecSpace>(transfer, true);
  } else {
    throw std::runtime_error("Unexpected dimension for projection");
  }
  Kokkos::Profiling::popRegion(); // Compute target point projection on
                                  // triangulated skin mesh
 
  Kokkos::Profiling::popRegion(); // FTInterpExtrap
}
} // namespace FiniteElements
} // namespace PACMAN