BSplinesStencilFactory.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 "tarch/la/Vector.h"
#include "tarch/la/Matrix.h"
#include "peano4/utils/Globals.h"
 
 
namespace toolbox {
  namespace finiteelements {
    class BSplinesStencilFactory;
  }
}
 
 
class toolbox::finiteelements::BSplinesStencilFactory {
  public:
    static tarch::la::Vector<3,double>  get1DMassStencilP1();
    static tarch::la::Vector<5,double>  get1DMassStencilP2();
    static tarch::la::Vector<7,double>  get1DMassStencilP3();
    static tarch::la::Vector<9,double>  get1DMassStencilP4();
    static tarch::la::Vector<11,double> get1DMassStencilP5();
    static tarch::la::Vector<13,double> get1DMassStencilP6();
    static tarch::la::Vector<15,double> get1DMassStencilP7();
    static tarch::la::Vector<17,double> get1DMassStencilP8();
 
    static tarch::la::Vector<3,double>  get1DMassStencilP1(int elementOfSupport);
 
    static tarch::la::Vector<5,double>  get1DMassStencilP2(int elementOfSupport);
 
    static tarch::la::Vector<7,double>  get1DMassStencilP3(int elementOfSupport);
 
    static tarch::la::Vector<9,double>  get1DMassStencilP4(int elementOfSupport);
 
    static tarch::la::Vector<3,double>  get1DLaplaceStencilP1();
    static tarch::la::Vector<5,double>  get1DLaplaceStencilP2();
    static tarch::la::Vector<7,double>  get1DLaplaceStencilP3();
    static tarch::la::Vector<9,double>  get1DLaplaceStencilP4();
    static tarch::la::Vector<11,double> get1DLaplaceStencilP5();
    static tarch::la::Vector<13,double> get1DLaplaceStencilP6();
    static tarch::la::Vector<15,double> get1DLaplaceStencilP7();
    static tarch::la::Vector<17,double> get1DLaplaceStencilP8();
 
    static tarch::la::Vector<3,double>  get1DLaplaceStencilP1(int elementOfSupport);
    static tarch::la::Vector<5,double>  get1DLaplaceStencilP2(int elementOfSupport);
    static tarch::la::Vector<7,double>  get1DLaplaceStencilP3(int elementOfSupport);
    static tarch::la::Vector<9,double>  get1DLaplaceStencilP4(int elementOfSupport);
 
    static tarch::la::Vector<ThreePowerD,double>   getLaplacianStencilP1( const tarch::la::Vector<Dimensions,double>& scaling = 1.0 );
    static tarch::la::Vector<FivePowerD,double>    getLaplacianStencilP2( const tarch::la::Vector<Dimensions,double>& scaling = 1.0 );
    static tarch::la::Vector<SevenPowerD,double>   getLaplacianStencilP3( const tarch::la::Vector<Dimensions,double>& scaling = 1.0 );
    static tarch::la::Vector<NinePowerD,double>    getLaplacianStencilP4( const tarch::la::Vector<Dimensions,double>& scaling = 1.0 );
 
    static tarch::la::Vector<ThreePowerD,double>   getMassStencilP1();
    static tarch::la::Vector<FivePowerD,double>    getMassStencilP2();
    static tarch::la::Vector<SevenPowerD,double>   getMassStencilP3();
    static tarch::la::Vector<NinePowerD,double>    getMassStencilP4();
 
    template <int p, int N>
    static tarch::la::Matrix<N*N,N*N,double>  getElementWiseAssemblyMatrix(
      tarch::la::Vector<2*p+1,double> stencilX[N],
      tarch::la::Vector<2*p+1,double> stencilY[N]
    );
 
 
    template <int p, int N>
    static tarch::la::Matrix<N*N*N,N*N*N,double>  getElementWiseAssemblyMatrix(
      tarch::la::Vector<2*p+1,double> stencilX[N],
      tarch::la::Vector<2*p+1,double> stencilY[N],
      tarch::la::Vector<2*p+1,double> stencilZ[N]
    );
};
 
 
#include "BSplinesStencilFactory.cpph"