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

#include <blockLpNorm2D.h>

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

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

Public Member Functions
	BlockLpNorm2D (BlockF2D< W > &f, BlockIndicatorF2D< T > &indicatorF)

bool	operator() (W output[], const int input[]) override
	has to be implemented for 'every' derived class

Public Member Functions inherited from olb::BlockF2D< W >
virtual BlockStructureD< 2 > &	getBlockStructure ()
	virtual destructor for defined behaviour

void	setBlockStructure (BlockStructureD< 2 > *blockStructure)

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

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

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

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

Public Member Functions inherited from olb::GenericF< W, int >
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

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

bool	operator() (W output[], int input0)

bool	operator() (W output[], int input0, int input1)

bool	operator() (W output[], int input0, int input1, int input2)

bool	operator() (W output[], int input0, int input1, int input2, int input3)

Protected Attributes
BlockF2D< W > &	_f

BlockIndicatorF2D< T > &	_indicatorF

Protected Attributes inherited from olb::BlockF2D< W >
BlockStructureD< 2 > *	_blockStructure

Additional Inherited Members
Public Types inherited from olb::GenericF< W, int >
using	targetType

using	sourceType

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

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

	BlockF2D (int targetDim)

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

Detailed Description

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

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

Instances are contained in SuperLpNorm2D::_blockF.

Definition at line 37 of file blockLpNorm2D.h.

Constructor & Destructor Documentation

◆ BlockLpNorm2D()

template<typename T , typename W , int P>

olb::BlockLpNorm2D< T, W, P >::BlockLpNorm2D	(	BlockF2D< W > &	f,
		BlockIndicatorF2D< T > &	indicatorF )

Parameters

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

Definition at line 35 of file blockLpNorm2D.hh.

  : BlockF2D<W>(f.getBlockStructure(), 1),
    _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::BlockLpNorm2D< T, W, P >::_f, olb::BlockLpNorm2D< T, W, P >::_indicatorF, olb::GenericF< T, S >::getName(), olb::GenericF< W, int >::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::BlockLpNorm2D< T, W, P >::operator()	(	W	output[],
		const int	input[] )

overridevirtual

has to be implemented for 'every' derived class

Implements olb::GenericF< W, int >.

Definition at line 48 of file blockLpNorm2D.hh.

{
  const auto& blockGeometry = _indicatorF.getBlockGeometry();
  const int nX = blockGeometry.getNx();
  const int nY = blockGeometry.getNy();
  const T weight = util::pow(blockGeometry.getDeltaR(), 2);
 
  W sum[2] = {0,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]) {
      if (_indicatorF(inputTmp)) {
        _f(outputTmp, inputTmp);
        for (int iDim = 0; iDim < _f.getTargetDim(); ++iDim) {
          LpNormImpl<T,W,P>()(sum, outputTmp[iDim], weight);
        }
      }
    }
  }
 
  output[0] = LpNormImpl<T,W,P>().enclose(sum[0]);
 
  return true;
}