Interpolation.h

Go to the documentation of this file.

// This file is part of the Peano project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
#include <map>
 
#include "peano4/datamanagement/CellMarker.h"
#include "peano4/datamanagement/FaceMarker.h"
#include "peano4/utils/Loop.h"
#include "tarch/la/DynamicMatrix.h"
 
namespace toolbox {
  namespace blockstructured {
    void interpolateHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ normalInterpolationMatrix1d,
      const double* __restrict__ tangentialInterpolationMatrix1d,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
 
    void interpolateHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ normalInterpolationMatrix1d,
      const double* __restrict__ tangentialInterpolationMatrix1d,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ interpolationData,
      const int* __restrict__ columnIndices,
      const int* __restrict__ rowIndices,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    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
    );
 
 
    void interpolateHaloLayer_AoS_second_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    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
    );
 
 
    void interpolateHaloLayer_AoS_third_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
 
    void interpolateHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ interpolationData,
      const int* __restrict__ columnIndices,
      const int* __restrict__ rowIndices,
      int dataStart,
      int rowIndicesStart,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateCell_AoS_tensor_product(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ interpolationMatrix1d,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
    void interpolateCell_AoS_matrix(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ interpolationData,
      const int* __restrict__ columnIndices,
      const int* __restrict__ rowIndices,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
    void interpolateCell_AoS_second_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
    void interpolateCell_AoS_third_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
 
    void interpolateHaloLayer_AoS_piecewise_constant(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateHaloLayer_AoS_piecewise_constant(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
 
    void interpolateCell_AoS_piecewise_constant(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
 
    void interpolateHaloLayer_AoS_linear_with_constant_extrapolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateHaloLayer_AoS_linear_with_constant_extrapolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateCell_AoS_linear_with_constant_extrapolation(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__                coarseGridCellValues,
      double* __restrict__                      fineGridCellValues
    );
 
 
    void interpolateHaloLayer_AoS_linear_with_linear_extrapolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateHaloLayer_AoS_linear_with_linear_extrapolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateCell_AoS_linear_with_linear_extrapolation(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
 
    void interpolateHaloLayer_AoS_linear_with_constant_extrapolation_and_linear_normal_interpolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateHaloLayer_AoS_linear_with_constant_extrapolation_and_linear_normal_interpolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateCell_AoS_linear_with_constant_extrapolation_and_linear_normal_interpolation(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
    void interpolateHaloLayer_AoS_linear_with_linear_extrapolation_and_linear_normal_interpolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateHaloLayer_AoS_linear_with_linear_extrapolation_and_linear_normal_interpolation(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      const double* __restrict__ coarseGridFaceValues,
      double* __restrict__ fineGridFaceValues
    );
 
    void interpolateCell_AoS_linear_with_linear_extrapolation_and_linear_normal_interpolation(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ coarseGridCellValues,
      double* __restrict__ fineGridCellValues
    );
 
    void interpolateCellDataAssociatedToVolumesIntoOverlappingCell_matrix(
      int numberOfDoFsPerAxisInSourcePatch,
      int numberOfDoFsPerAxisInDestinationPatch,
      int haloSourcePatch,
      int haloDestinationPatch,
      int unknowns,
      const double* __restrict__ interpolationData,
      const int* __restrict__ columnIndices,
      const int* __restrict__ rowIndices,
      const double* __restrict__ sourceValues,
      double* __restrict__ destinationValues,
      ::peano4::utils::LoopPlacement parallelisation
    );
 
    void interpolateCellDataAssociatedToVolumesIntoOverlappingCell_secondOrder(
      int numberOfDoFsPerAxisInSourcePatch,
      int numberOfDoFsPerAxisInDestinationPatch,
      int haloSourcePatch,
      int haloDestinationPatch,
      int unknowns,
      const double* __restrict__ sourceValues,
      double* __restrict__ destinationValues,
      ::peano4::utils::LoopPlacement parallelisation
    );
 
    void interpolateCellDataAssociatedToVolumesIntoOverlappingCell_linear(
      int numberOfDoFsPerAxisInSourcePatch,
      int numberOfDoFsPerAxisInDestinationPatch,
      int haloSourcePatch,
      int haloDestinationPatch,
      int unknowns,
      const double* __restrict__ sourceValues,
      double* __restrict__ destinationValues,
      ::peano4::utils::LoopPlacement parallelisation
    );
 
    void interpolateCellDataAssociatedToVolumesIntoOverlappingCell_fourthOrder(
      int numberOfDoFsPerAxisInSourcePatch,
      int numberOfDoFsPerAxisInDestinationPatch,
      int haloSourcePatch,
      int haloDestinationPatch,
      int unknowns,
      const double* __restrict__ sourceValues,
      double* __restrict__ destinationValues,
      ::peano4::utils::LoopPlacement parallelisation
    );
 
    namespace internal {
      int getCoarseLayerOfFineLayerInOverlap(int fineLayer, int overlap);
 
      tarch::la::DynamicMatrix* createLinearInterpolationMatrix(
        int  numberOfDoFsPerAxisInPatch,
        int  normal,
        bool extrapolateLinearly,
        bool interpolateLinearlyAlongNormal
      );
 
      tarch::la::DynamicMatrix* createLinearInterpolationMatrix(
        int  numberOfDoFsPerAxisInPatch,
        bool extrapolateLinearly
      );
 
      tarch::la::DynamicMatrix* createLinearInterpolationMatrix(
        int  numberOfDoFsPerAxisInInputPatch,
        int  numberOfDoFsPerAxisInOutputPatch,
        bool extrapolateLinearly
      );
 
      tarch::la::DynamicMatrix* createPiecewiseConstantInterpolationMatrix(
        int numberOfDoFsPerAxisInPatch,
        int normal,
        int overlap
      );
 
      tarch::la::DynamicMatrix* createPiecewiseConstantInterpolationMatrix(
        int numberOfDoFsPerAxisInPatch
      );
 
      tarch::la::DynamicMatrix*
        createInterpolationMatrixFrom1dTemplateByInsertingZeroColsAndRows(
          const tarch::la::DynamicMatrix& P1d,
          int                             numberOfDoFsPerAxisInPatch,
          int                             normal
        );
 
      tarch::la::DynamicMatrix*
        createInterpolationMatrixFrom1dTemplateByLinearInterpolationAlongNormal(
          const tarch::la::DynamicMatrix& P1d,
          int                             numberOfDoFsPerAxisInPatch,
          int                             normal
        );
 
      void projectInterpolatedFineCellsOnHaloLayer_AoS(
        const peano4::datamanagement::FaceMarker& marker,
        int                                       numberOfDoFsPerAxisInPatch,
        int                                       overlap,
        int                                       unknowns,
        const double* __restrict__ fineGridCellValuesLeft,
        const double* __restrict__ fineGridCellValuesRight,
        double* fineGridFaceValues
      );
 
      void interpolateCell_AoS_piecewise_constant(
        const tarch::la::Vector<Dimensions, int>&
            relativePositionWithinFatherCell,
        int                          numberOfDoFsPerAxisInPatch,
        int                          unknowns,
        const double* __restrict__   fineGridValues,
        double* __restrict__         coarseGridValues
      );
 
      void interpolateCell_AoS_linear(
        const tarch::la::Vector<Dimensions, int>&
            relativePositionWithinFatherCell,
        int                          numberOfDoFsPerAxisInPatch,
        int                          unknowns,
        const double* __restrict__   fineGridValues,
        double* __restrict__         coarseGridValues,
        bool                         extrapolateLinearly
      );
 
      struct FaceInterpolationOperatorKey {
        int numberOfDoFsPerAxisInPatch;
        int overlap;
        int normal;
 
        inline FaceInterpolationOperatorKey(
          int numberOfDoFsPerAxisInPatch_,
          int overlap_,
          int normal_
        ):
          numberOfDoFsPerAxisInPatch(numberOfDoFsPerAxisInPatch_),
          overlap(overlap_),
          normal(normal_) {}
 
        inline bool operator<(const FaceInterpolationOperatorKey& otherKey
        ) const {
          return (this->numberOfDoFsPerAxisInPatch < otherKey.numberOfDoFsPerAxisInPatch)
              or (this->numberOfDoFsPerAxisInPatch == otherKey.numberOfDoFsPerAxisInPatch and this->overlap < otherKey.overlap)
              or (this->numberOfDoFsPerAxisInPatch == otherKey.numberOfDoFsPerAxisInPatch and this->overlap == otherKey.overlap and this->normal < otherKey.normal);
        }
      };
 
      struct CellInterpolationOperatorKey {
        int numberOfDoFsPerAxisInPatch;
 
        inline CellInterpolationOperatorKey(int numberOfDoFsPerAxisInPatch_):
          numberOfDoFsPerAxisInPatch(numberOfDoFsPerAxisInPatch_) {}
 
        inline bool operator<(const CellInterpolationOperatorKey& otherKey
        ) const {
          return this->numberOfDoFsPerAxisInPatch
                 < otherKey.numberOfDoFsPerAxisInPatch;
        }
      };
 
      typedef std::map<FaceInterpolationOperatorKey, tarch::la::DynamicMatrix*>
        FaceInterpolationMap;
 
      typedef std::map<CellInterpolationOperatorKey, tarch::la::DynamicMatrix*>
        CellInterpolationMap;
    } // namespace internal
  }   // namespace blockstructured
} // namespace toolbox