Restriction.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 "peano4/datamanagement/CellMarker.h"
#include "peano4/datamanagement/FaceMarker.h"
#include "tarch/la/DynamicMatrix.h"
 
namespace toolbox {
  namespace blockstructured {
    void clearHaloLayerAoS(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   values
    );
 
    void clearCell(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   values
    );
 
    void restrictInnerHalfOfHaloLayer_AoS_inject(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    );
 
    void restrictHaloLayer_AoS_inject(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
    void restrictCell_AoS_inject(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
    void restrictCellIntoOverlappingCell_inject(
      int     numberOfDoFsPerAxisInSourcePatch,
      int     numberOfDoFsPerAxisInDestinationPatch,
      int     unknowns,
      double* sourceValues,
      double* destinationValues
    );
 
    void restrictCellIntoOverlappingCell_inject_and_average(
      int     numberOfDoFsPerAxisInSourcePatch,
      int     numberOfDoFsPerAxisInDestinationPatch,
      int     unknowns,
      double* sourceValues,
      double* destinationValues,
      double  weightOfInjectedValue = 0.5
    );
 
 
    void restrictCell_AoS_averaging(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
    void restrictHaloLayer_AoS_averaging(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
    void restrictInnerHalfOfHaloLayer_AoS_averaging(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    );
 
    void restrictCell_AoS_tensor_product(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ tangentialRestrictionMatrix1d,
      double* fineGridValues,
      double* coarseGridValues
    );
 
    void restrictCell_AoS_matrix(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      const double* __restrict__ restrictionData,
      const int* __restrict__ columnIndices,
      const int* __restrict__ rowIndices,
      double* fineGridValues,
      double* coarseGridValues
    );
 
    void restrictCell_AoS_second_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
    void restrictCell_AoS_third_order(
      const peano4::datamanagement::CellMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
 
    void restrictHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ normalRestrictionMatrix1d,
      const double* __restrict__ tangentialRestrictionMatrix1d,
      double* fineGridValues,
      double* coarseGridValues
    );
 
    void restrictHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ restrictionData,
      const int* __restrict__ columnIndices,
      const int* __restrict__ rowIndices,
      int     dataOffset,
      int     rowIndicesOffset,
      double* fineGridValues,
      double* coarseGridValues
    );
 
    void restrictHaloLayer_AoS_second_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
    void restrictHaloLayer_AoS_third_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues
    );
 
    void restrictInnerHalfOfHaloLayer_AoS_tensor_product(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ normalRestrictionMatrix1d,
      const double* __restrict__ tangentialRestrictionMatrix1d,
      double* fineGridValues,
      double* coarseGridValues,
      bool    swapInsideOutside = false
    );
 
    void restrictInnerHalfOfHaloLayer_AoS_matrix(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      const double* __restrict__ restrictionMatrixData,
      const int* __restrict__ columnIndices,
      const int* __restrict__ rowIndices,
      int     dataOffset,
      int     rowIndicesOffset,
      double* fineGridValues,
      double* coarseGridValues,
      bool    swapInsideOutside = false
    );
 
    void restrictInnerHalfOfHaloLayer_AoS_second_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    );
 
    void restrictInnerHalfOfHaloLayer_AoS_third_order(
      const peano4::datamanagement::FaceMarker& marker,
      int                                       numberOfDoFsPerAxisInPatch,
      int                                       overlap,
      int                                       unknowns,
      double*                                   fineGridValues,
      double*                                   coarseGridValues,
      bool                                      swapInsideOutside = false
    );
 
 
    namespace internal {
 
      void clearHalfOfHaloLayerAoS(
        const peano4::datamanagement::FaceMarker& marker,
        int                                       numberOfDoFsPerAxisInPatch,
        int                                       overlap,
        int                                       unknowns,
        bool                                      clearInnerPart,
        double*                                   values
      );
 
      void projectCells_AoS(
        const peano4::datamanagement::CellMarker& fineGridCellMarker,
        int                                       numberOfDoFsPerAxisInPatch,
        std::function<void(
          tarch::la::Vector<Dimensions, int>    coarseVolume,
          tarch::la::Vector<Dimensions, int>    fineVolume,
          tarch::la::Vector<Dimensions, double> coarseVolumeCentre,
          tarch::la::Vector<Dimensions, double> fineVolumeCentre,
          double                                coarseVolumeH,
          double                                fineVolumeH
        )>                                        update
      );
    } // namespace internal
  }   // namespace blockstructured
} // namespace toolbox