InterpolationRestriction.h

Go to the documentation of this file.

// This file is part of the ExaHyPE2 project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
 
#include "toolbox/blockstructured/Interpolation.h"
#include "toolbox/blockstructured/Restriction.h"
 
 
namespace toolbox {
  namespace blockstructured {
    template <typename C>
    void interpolateHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      interpolateHaloLayer_AoS_tensor_product(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::NormalInterpolationMatrix1d,
        C::TangentialInterpolationMatrix1d,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
 
    template <typename C>
    void interpolateHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      interpolateHaloLayer_AoS_tensor_product(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::NormalInterpolationMatrix1d,
        C::TangentialInterpolationMatrix1d,
        coarseGridCellValues,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
 
    template <typename C>
    void interpolateHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      int normal = marker.getSelectedFaceNumber() % 3;
      int patchIndexLinear = (Dimensions == 2) ? marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) : marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) * 3 + marker.getRelativePositionWithinFatherCell()((normal + 2) % Dimensions);
 
      interpolateHaloLayer_AoS_matrix(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::InterpolationData,
        C::InterpolationColumnIndices,
        C::InterpolationRowIndices,
        coarseGridCellValues,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
    template <typename C>
    void interpolateHaloLayer_AoS_second_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      interpolateHaloLayer_AoS_second_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        coarseGridCellValues,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
    template <typename C>
    void interpolateHaloLayer_AoS_third_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      interpolateHaloLayer_AoS_third_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        coarseGridCellValues,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
 
    template <typename C>
    void interpolateCell_AoS_tensor_product(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    ) {
      interpolateCell_AoS_tensor_product(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        C::TangentialInterpolationMatrix1d,
        coarseGridCellValues,
        fineGridCellValues
      );
    }
 
 
    template <typename C>
    void restrictCell_AoS_tensor_product(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      restrictCell_AoS_tensor_product(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        C::TangentialRestrictionMatrix1d,
        fineGridValues,
        coarseGridValues
      );
    }
 
 
    template <typename C>
    void interpolateCell_AoS_matrix(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    ) {
      interpolateCell_AoS_matrix(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        C::InterpolationData,
        C::InterpolationColumnIndices,
        C::InterpolationRowIndices,
        coarseGridCellValues,
        fineGridCellValues
      );
    }
 
    template <typename C>
    void interpolateCell_AoS_second_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    ) {
      interpolateCell_AoS_second_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        coarseGridCellValues,
        fineGridCellValues
      );
    }
 
    template <typename C>
    void interpolateCell_AoS_third_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    ) {
      interpolateCell_AoS_third_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        coarseGridCellValues,
        fineGridCellValues
      );
    }
 
 
    template <typename C>
    void restrictCell_AoS_matrix(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      restrictCell_AoS_matrix(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        C::RestrictionData,
        C::RestrictionColumnIndices,
        C::RestrictionRowIndices,
        fineGridValues,
        coarseGridValues
      );
    }
 
    template <typename C>
    void restrictCell_AoS_second_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      restrictCell_AoS_second_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        fineGridValues,
        coarseGridValues
      );
    }
 
    template <typename C>
    void restrictCell_AoS_third_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      restrictCell_AoS_third_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        unknowns,
        fineGridValues,
        coarseGridValues
      );
    }
 
 
    template <typename C>
    void restrictHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      restrictHaloLayer_AoS_tensor_product(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::NormalRestrictionMatrix1d,
        C::TangentialRestrictionMatrix1d,
        fineGridValues,
        coarseGridValues
      );
    }
 
 
    template <typename C>
    void restrictInnerHalfOfHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    ) {
      restrictInnerHalfOfHaloLayer_AoS_tensor_product(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::NormalRestrictionMatrix1d,
        C::TangentialRestrictionMatrix1d,
        fineGridValues,
        coarseGridValues,
        swapInsideOutside
      );
    }
 
 
    template <typename C>
    void interpolateHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      int normal = marker.getSelectedFaceNumber() % 3;
      int patchIndexLinear = (Dimensions == 2) ? marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) : marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) * 3 + marker.getRelativePositionWithinFatherCell()((normal + 2) % Dimensions);
 
      interpolateHaloLayer_AoS_matrix(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::InterpolationData,
        C::InterpolationColumnIndices,
        C::InterpolationRowIndices,
        patchIndexLinear * C::InterpolationDataOffset,
        patchIndexLinear * C::InterpolationRowIndicesOffset,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
 
    template <typename C>
    void interpolateHaloLayer_AoS_second_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      interpolateHaloLayer_AoS_second_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
    template <typename C>
    void interpolateHaloLayer_AoS_third_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    ) {
      interpolateHaloLayer_AoS_third_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        coarseGridFaceValues,
        fineGridFaceValues
      );
    }
 
 
    template <typename C>
    void restrictHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      int normal = marker.getSelectedFaceNumber() % 3;
      int patchIndexLinear = (Dimensions == 2) ? marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) : marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) + 3 * marker.getRelativePositionWithinFatherCell()((normal + 2) % Dimensions);
 
      restrictHaloLayer_AoS_matrix(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::RestrictionData,
        C::RestrictionColumnIndices,
        C::RestrictionRowIndices,
        patchIndexLinear * C::RestrictionDataOffset,
        patchIndexLinear * C::RestrictionRowIndicesOffset,
        fineGridValues,
        coarseGridValues
      );
    }
 
    template <typename C>
    void restrictHaloLayer_AoS_second_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      restrictHaloLayer_AoS_second_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        fineGridValues,
        coarseGridValues
      );
    }
 
    template <typename C>
    void restrictHaloLayer_AoS_third_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    ) {
      restrictHaloLayer_AoS_third_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        fineGridValues,
        coarseGridValues
      );
    }
 
 
    template <typename C>
    void restrictInnerHalfOfHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    ) {
      int normal = marker.getSelectedFaceNumber() % 3;
      int patchIndexLinear = (Dimensions == 2) ? marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) : marker.getRelativePositionWithinFatherCell()((normal + 1) % Dimensions) + 3 * marker.getRelativePositionWithinFatherCell()((normal + 2) % Dimensions);
 
      restrictInnerHalfOfHaloLayer_AoS_matrix(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        C::RestrictionData,
        C::RestrictionColumnIndices,
        C::RestrictionRowIndices,
        patchIndexLinear * C::RestrictionDataOffset,
        patchIndexLinear * C::RestrictionRowIndicesOffset,
        fineGridValues,
        coarseGridValues,
        swapInsideOutside
      );
    }
 
    template <typename C>
    void restrictInnerHalfOfHaloLayer_AoS_second_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    ) {
      restrictInnerHalfOfHaloLayer_AoS_second_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        fineGridValues,
        coarseGridValues,
        swapInsideOutside
      );
    }
 
    template <typename C>
    void restrictInnerHalfOfHaloLayer_AoS_third_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    ) {
      restrictInnerHalfOfHaloLayer_AoS_third_order(
        marker,
        numberOfDoFsPerAxisInPatch,
        overlap,
        unknowns,
        fineGridValues,
        coarseGridValues,
        swapInsideOutside
      );
    }
 
  } // namespace blockstructured
} // namespace toolbox