Block level functor that returns the Lp norm over omega of the euklid norm of the input block functor. More...

#include <blockLpNorm3D.h>

Inheritance diagram for olb::BlockLpNorm3D< T, W, P >:

Collaboration diagram for olb::BlockLpNorm3D< T, W, P >:

Public Member Functions
	BlockLpNorm3D (BlockF3D< W > &f, BlockIndicatorF3D< T > &indicatorF)

bool	operator() (W output[], const int input[]) override

Public Member Functions inherited from olb::BlockF3D< W >
	~BlockF3D () override
	virtual destructor for defined behaviour

virtual BlockStructureD< 3 > &	getBlockStructure () const

BlockF3D< W > &	operator- (BlockF3D< W > &rhs)

BlockF3D< W > &	operator+ (BlockF3D< W > &rhs)

BlockF3D< W > &	operator* (BlockF3D< W > &rhs)

BlockF3D< W > &	operator/ (BlockF3D< W > &rhs)

Public Member Functions inherited from olb::GenericF< T, S >
virtual	~GenericF ()=default

int	getSourceDim () const
	read only access to member variable _m

int	getTargetDim () const
	read only access to member variable _n

std::string &	getName ()
	read and write access to name

std::string const &	getName () const
	read only access to name

virtual bool	operator() (T output[], const S input[])=0
	has to be implemented for 'every' derived class

bool	operator() (T output[])
	wrapper that call the pure virtual operator() (T output[], const S input[]) from above

bool	operator() (T output[], S input0)

bool	operator() (T output[], S input0, S input1)

bool	operator() (T output[], S input0, S input1, S input2)

bool	operator() (T output[], S input0, S input1, S input2, S input3)

Protected Attributes
BlockF3D< W > &	_f

BlockIndicatorF3D< T > &	_indicatorF

Protected Attributes inherited from olb::BlockF3D< W >
BlockStructureD< 3 > &	_blockStructure

Additional Inherited Members
Public Types inherited from olb::GenericF< T, S >
using	targetType = T

using	sourceType = S

Public Attributes inherited from olb::GenericF< T, S >
std::shared_ptr< GenericF< T, S > >	_ptrCalcC
	memory management, frees resouces (calcClass)

Protected Member Functions inherited from olb::BlockF3D< W >
	BlockF3D (BlockStructureD< 3 > &blockStructure, int targetDim)

Protected Member Functions inherited from olb::GenericF< T, S >
	GenericF (int targetDim, int sourceDim)

Detailed Description

template<typename T, typename W, int P>
class olb::BlockLpNorm3D< T, W, P >

Block level functor that returns the Lp norm over omega of the euklid norm of the input block functor.

Definition at line 35 of file blockLpNorm3D.h.

Constructor & Destructor Documentation

◆ BlockLpNorm3D()

template<typename T , typename W , int P>

olb::BlockLpNorm3D< T, W, P >::BlockLpNorm3D	(	BlockF3D< W > &	f,
		BlockIndicatorF3D< T > &	indicatorF )

Parameters

f	data functor
indicatorF	indicator functor describing the subset to be integrated

Definition at line 35 of file blockLpNorm3D.hh.

  : BlockF3D<W>(f.getBlockStructure(), f.getTargetDim()),
    _f(f),
    _indicatorF(indicatorF)
{
  OLB_ASSERT(_f.getSourceDim() == _indicatorF.getSourceDim(),
             "functor source dimension equals indicator source dimension");
 
  this->getName() = "BlockL" + std::to_string(P) + "Norm(" + _f.getName() + ")";
}

References olb::BlockLpNorm3D< T, W, P >::_f, olb::BlockLpNorm3D< T, W, P >::_indicatorF, olb::GenericF< T, S >::getName(), olb::GenericF< T, S >::getSourceDim(), and OLB_ASSERT.

Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename W , int P>

bool olb::BlockLpNorm3D< T, W, P >::operator()	(	W	output[],
		const int	input[] )

override

Definition at line 48 of file blockLpNorm3D.hh.

{
  const auto& blockGeometry = _indicatorF.getBlockGeometry();
  const int nX = blockGeometry.getNx();
  const int nY = blockGeometry.getNy();
  const int nZ = blockGeometry.getNz();
  const T weight = util::pow(blockGeometry.getDeltaR(), 3);
 
  output[0] = W(0);
  W outputTmp[_f.getTargetDim()];
  int inputTmp[_f.getSourceDim()];
 
  for (inputTmp[0] = 0; inputTmp[0] < nX; ++inputTmp[0]) {
    for (inputTmp[1] = 0; inputTmp[1] < nY; ++inputTmp[1]) {
      for (inputTmp[2] = 0; inputTmp[2] < nZ; ++inputTmp[2]) {
        if (_indicatorF(inputTmp)) {
          _f(outputTmp, inputTmp);
          for (int iDim = 0; iDim < _f.getTargetDim(); ++iDim) {
            output[0] = LpNormImpl<T,W,P>()(output[0], outputTmp[iDim], weight);
          }
        }
      }
    }
  }
 
  output[0] = LpNormImpl<T,W,P>().enclose(output[0]);
 
  return true;
}