StencilFactory.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 "Stencil.h"
 
#include <bitset>
#include <functional>
 
 
namespace toolbox {
  namespace finiteelements {
    void preprocessBoundaryStencil( Stencil& stencil, const std::bitset<Dimensions*2>& boundaryFaceNormals );
 
     Stencil exchangeCoordinates( const Stencil& stencil, int coord0, int coord1 );
 
     tarch::la::Vector<3,double> get1DMassStencil();
    
 
     tarch::la::Vector<3,double> get1DLaplaceStencil();
 
     tarch::la::Vector<3,double> get1DIdentity();
     tarch::la::Vector<3,double> get1DLeftIdentity();
     tarch::la::Vector<3,double> get1DRightIdentity();
 
     tarch::la::Vector<5,double> get1DLinearInterpolationStencil();
 
     tarch::la::Vector<3,double> get1DMeanValueStencil();
 
     tarch::la::Vector<3,double> get1DDownwindStencil();
     tarch::la::Vector<3,double> get1DUpwindStencil();
     tarch::la::Vector<3,double> get1DCentralDifferencesStencil();
 
     tarch::la::Vector<3*3,double> stencilProduct(
            const tarch::la::Vector<3,double>& a,
            const tarch::la::Vector<3,double>& b
            );
     tarch::la::Vector<5*5,double> stencilProduct(
            const tarch::la::Vector<5,double>& a,
            const tarch::la::Vector<5,double>& b
            );
 
     tarch::la::Vector<3*3*3,double> stencilProduct(
            const tarch::la::Vector<3,double>& a,
            const tarch::la::Vector<3,double>& b,
            const tarch::la::Vector<3,double>& c
            );
     tarch::la::Vector<5*5*5,double> stencilProduct(
            const tarch::la::Vector<5,double>& a,
            const tarch::la::Vector<5,double>& b,
            const tarch::la::Vector<5,double>& c
            );
 
     tarch::la::Vector<3*3*3*3,double> stencilProduct(
            const tarch::la::Vector<3,double>& a,
            const tarch::la::Vector<3,double>& b,
            const tarch::la::Vector<3,double>& c,
            const tarch::la::Vector<3,double>& d
            );
     tarch::la::Vector<5*5*5*5,double> stencilProduct(
            const tarch::la::Vector<5,double>& a,
            const tarch::la::Vector<5,double>& b,
            const tarch::la::Vector<5,double>& c,
            const tarch::la::Vector<5,double>& d
            );
 
     tarch::la::Vector<3*3*3*3*3,double> stencilProduct(
            const tarch::la::Vector<3,double>& a,
            const tarch::la::Vector<3,double>& b,
            const tarch::la::Vector<3,double>& c,
            const tarch::la::Vector<3,double>& d,
            const tarch::la::Vector<3,double>& e
            );
     tarch::la::Vector<5*5*5*5*5,double> stencilProduct(
            const tarch::la::Vector<5,double>& a,
            const tarch::la::Vector<5,double>& b,
            const tarch::la::Vector<5,double>& c,
            const tarch::la::Vector<5,double>& d,
            const tarch::la::Vector<5,double>& e
            );
 
     tarch::la::Vector<ThreePowerD,double>     getLaplacian(
            const tarch::la::Vector<Dimensions,double>& h,
            const tarch::la::Vector<Dimensions,double>& scaling = 1.0
            );
     tarch::la::Vector<ThreePowerD,std::complex<double> >  getLaplacian(
            const tarch::la::Vector<Dimensions,std::complex<double> >& h,
            const tarch::la::Vector<Dimensions,std::complex<double> >& scaling = std::complex<double>(1.0)
            );
     tarch::la::Vector<ThreePowerD,std::complex<double> >  getLaplacian(
            const tarch::la::Vector<Dimensions, double >& h,
            const tarch::la::Vector<Dimensions,std::complex<double> >& scaling
            );
 
     tarch::la::Vector<ThreePowerD,double>  getUpwindDiscretisedConvection(
             const tarch::la::Vector<Dimensions,double >& h,
             const tarch::la::Vector<Dimensions,double>& convCoeff
             );
 
     void applyNeumannBC(
       tarch::la::Vector<ThreePowerD,double>&      stencil,
       double&                                       rhs,
       const std::bitset<Dimensions*2>&              boundaryFaceNormals,
       const tarch::la::Vector<Dimensions,double>&   h,
       double                                        derivative
     );
 
     tarch::la::Vector<ThreePowerD,double>                 getMassMatrix(const tarch::la::Vector<Dimensions,double>& h);
     tarch::la::Vector<ThreePowerD,std::complex<double> >  getMassMatrix(const tarch::la::Vector<Dimensions,std::complex<double> >& h);
 
 
     tarch::la::Vector<ThreePowerD,std::complex<double> >  getHelmholtzShiftMatrix(
            const tarch::la::Vector<Dimensions,double>& h,
            const std::complex<double>& phi
            );
 
     tarch::la::Vector<ThreePowerD,std::complex<double> >  getHelmholtzShiftMatrix(
            const tarch::la::Vector<Dimensions,std::complex<double> >& h,
            const std::complex<double>& phi
            );
 
 
     tarch::la::Vector<ThreePowerD,double>                 getIdentity(const tarch::la::Vector<Dimensions,double>& h);
     tarch::la::Vector<ThreePowerD,std::complex<double> >  getIdentity(const tarch::la::Vector<Dimensions,std::complex<double> >& h);
 
     tarch::la::Vector<FivePowerD,double> getDLinearInterpolation();
 
     Stencil extractElementStencil( const Stencil& stencil, const tarch::la::Vector<Dimensions,int>&   adjacentCell );
 
 
     ElementWiseAssemblyMatrix getPoissonMatrixWithJumpingCoefficient(
       const tarch::la::Vector<Dimensions,double>&   cellCentre,
       const tarch::la::Vector<Dimensions,double>&   h,
       int                                           integrationPointsPerAxis,
       std::function<double(const tarch::la::Vector<Dimensions,double>&)>
     );
  }
}