ElementMatrix.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>
 
 
namespace toolbox {
  namespace finiteelements {
     ElementWiseAssemblyMatrix getElementWiseAssemblyMatrix( const Stencil& stencil );
 
     ComplexElementWiseAssemblyMatrix getElementWiseAssemblyMatrix( const ComplexStencil& complexStencil );
 
 
     Stencil reconstructUniformStencilFragments(const ElementWiseAssemblyMatrix& matrix );
 
     tarch::la::Vector<TwoPowerD * ThreePowerD, double> reconstructStencilFragments(const ElementWiseAssemblyMatrix& matrix );
 
    ElementWiseAssemblyMatrix getElementWiseAssemblyMatrix( const VectorOfStencils& vectorOfStencils );
 
    ComplexElementWiseAssemblyMatrix getElementWiseAssemblyMatrix( const VectorOfComplexStencils& vectorOfComplexStencils );
 
     ElementWiseAssemblyMatrix getElementWiseAssemblyMatrix( const VectorOfStencils& vectorOfStencils, const std::bitset<ThreePowerD>& cellIsInside );
 
    ComplexElementWiseAssemblyMatrix getElementWiseAssemblyMatrix( const VectorOfComplexStencils& vectorOfStencils, const std::bitset<ThreePowerD>& cellIsInside );
 
    int mapElementMatrixEntryOntoStencilEntry(int row, int col);
 
    int mapElementMatrixEntryOntoStencilEntry(const tarch::la::Vector<Dimensions,int>& row, const tarch::la::Vector<Dimensions,int>&  col);
 
    double getDiagonalElement( const ElementWiseAssemblyMatrix& matrix );
    double getDiagonalElement( const Stencil& stencil );
 
    template<int StencilSize>
    tarch::la::Vector<StencilSize*StencilSize,double> stencilProduct(
      const tarch::la::Vector<StencilSize,double>& a,
      const tarch::la::Vector<StencilSize,double>& b
    );
 
    template<int StencilSize>
    tarch::la::Vector<StencilSize*StencilSize*StencilSize,double> stencilProduct(
      const tarch::la::Vector<StencilSize,double>& a,
      const tarch::la::Vector<StencilSize,double>& b,
      const tarch::la::Vector<StencilSize,double>& c
    );
 
    ElementWiseAssemblyMatrix hierarchicalTransform(
      const ElementWiseAssemblyMatrix&             matrix,
      const tarch::la::Vector<Dimensions,double>&  h,
      double                                       weight = 1.0
    );
    ElementWiseAssemblyMatrix inverseHierarchicalTransform(
      const ElementWiseAssemblyMatrix&             matrix,
      const tarch::la::Vector<Dimensions,double>&  h,
      double                                       weight = 1.0
    );
  }
}
 
#include "ElementMatrix.cpph"