olb::SuperLocalAverage2D< T, W > Class Template Reference

Averages given functor inside the local sphere. More...

#include <superLocalAverage2D.h>

Inheritance diagram for olb::SuperLocalAverage2D< T, W >:

Collaboration diagram for olb::SuperLocalAverage2D< T, W >:

Public Member Functions
	SuperLocalAverage2D (FunctorPtr< SuperF2D< T >> &&f, FunctorPtr< SuperIndicatorF2D< T >> &&indicatorF, T radius)
	Primary constructor. More...
bool	operator() (W output[], const int input[]) override
	Returns average of functor _f evaluated on all cells both inside a sphere of _radius around input and indicated by _indicatorF. More...
Public Member Functions inherited from olb::SuperF2D< T, W >
SuperF2D< T, W > &	operator- (SuperF2D< T, W > &rhs)
SuperF2D< T, W > &	operator+ (SuperF2D< T, W > &rhs)
SuperF2D< T, W > &	operator* (SuperF2D< T, W > &rhs)
SuperF2D< T, W > &	operator/ (SuperF2D< T, W > &rhs)
SuperStructure2D< T > &	getSuperStructure ()
int	getBlockFSize () const
BlockF2D< W > &	getBlockF (int iCloc)
Public Member Functions inherited from olb::GenericF< W, int >
virtual	~GenericF ()
int	getSourceDim () const
	read only access to member variable _m More...
int	getTargetDim () const
	read only access to member variable _n More...
std::string &	getName ()
	read and write access to name More...
std::string const &	getName () const
	read only access to name More...
bool	operator() (Woutput[])
	wrapper that call the pure virtual operator() (T output[], const S input[]) from above More...
bool	operator() (Woutput[], intinput0)
bool	operator() (Woutput[], intinput0, intinput1)
bool	operator() (Woutput[], intinput0, intinput1, intinput2)
bool	operator() (Woutput[], intinput0, intinput1, intinput2, intinput3)

Additional Inherited Members
Public Types inherited from olb::SuperF2D< T, W >
using	identity_functor_type = SuperIdentity2D< T, W >
Public Attributes inherited from olb::GenericF< W, int >
std::shared_ptr< GenericF< W, int > >	_ptrCalcC
	memory management, frees resouces (calcClass) More...
Protected Member Functions inherited from olb::SuperF2D< T, W >
	SuperF2D (SuperStructure2D< T > &superStructure, int targetDim)
Protected Member Functions inherited from olb::GenericF< W, int >
	GenericF (int targetDim, int sourceDim)
Protected Attributes inherited from olb::SuperF2D< T, W >
SuperStructure2D< T > &	_superStructure
std::vector< std::unique_ptr < BlockF2D< W > > >	_blockF
	Super functors may consist of several BlockF2D<W> derived functors. More...

Detailed Description

template<typename T, typename W = T>
class olb::SuperLocalAverage2D< T, W >

Averages given functor inside the local sphere.

Definition at line 36 of file superLocalAverage2D.h.

Constructor & Destructor Documentation

template<typename T , typename W >

olb::SuperLocalAverage2D< T, W >::SuperLocalAverage2D	(	FunctorPtr< SuperF2D< T >> &&	f,
		FunctorPtr< SuperIndicatorF2D< T >> &&	indicatorF,
		T	radius
	)

Primary constructor.

Parameters

f	Functor to be locally averaged
indicatorF	Indicator describing relevant cells
radius	Radius of the locality sphere

Definition at line 35 of file superLocalAverage2D.hh.

References olb::SuperF2D< T, W >::_blockF, olb::GenericF< W, int >::getName(), and olb::LoadBalancer< T >::size().

  : SuperF2D<T,W>(f->getSuperStructure(), f->getTargetDim()),
    _f(std::move(f)),
    _indicatorF(std::move(indicatorF)),
    _radius(radius)
{
  this->getName() = "LocalAverage(" + _f->getName() + ")";

  LoadBalancer<T>& load = _f->getSuperStructure().getLoadBalancer();

  if ( _f->getBlockFSize()          == load.size() &&
       _indicatorF->getBlockFSize() == load.size() ) {
    for (int iC = 0; iC < load.size(); ++iC) {
      this->_blockF.emplace_back(
        new BlockLocalAverage2D<T,W>(_f->getBlockF(iC),
                                     _indicatorF->getBlockIndicatorF(iC),
                                     _radius)
      );
    }
  }
}

Here is the call graph for this function:

Member Function Documentation

template<typename T , typename W >

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

overridevirtual

Returns average of functor _f evaluated on all cells both inside a sphere of _radius around input and indicated by _indicatorF.

Reimplemented from olb::SuperF2D< T, W >.

Definition at line 61 of file superLocalAverage2D.hh.

References olb::singleton::mpi().

{
  const auto& geometry = this->getSuperStructure().getCuboidGeometry();
  const auto& load     = this->getSuperStructure().getLoadBalancer();

  for (int i = 0; i < this->getTargetDim(); ++i) {
    output[i] = 0.;
  }

  if (!_indicatorF(input)) {
    return true;
  }

  T centerOfCircle[2];
  geometry.getPhysR(centerOfCircle, input);
  IndicatorCircle2D<T> analyticalCircle(centerOfCircle, _radius);
  SuperIndicatorFfromIndicatorF2D<T> latticeCircle(
    analyticalCircle,
    _indicatorF->getSuperGeometry());

  std::size_t voxels(0);
  int inputTmp[3];

  for (int iC = 0; iC < load.size(); ++iC) {
    inputTmp[0] = load.glob(iC);
    const auto& cuboid = geometry.get(inputTmp[0]);

    for (inputTmp[1] = 0; inputTmp[1] < cuboid.getNx(); ++inputTmp[1]) {
      for (inputTmp[2] = 0; inputTmp[2] < cuboid.getNy(); ++inputTmp[2]) {
        if (latticeCircle(inputTmp) && _indicatorF(inputTmp)) {
          T outputTmp[_f->getTargetDim()];
          _f(outputTmp, inputTmp);
          for (int i = 0; i < this->getTargetDim(); ++i) {
            output[i] += outputTmp[i];
          }
          voxels += 1;
        }
      }
    }
  }

#ifdef PARALLEL_MODE_MPI
  singleton::mpi().reduceAndBcast(voxels, MPI_SUM);
#endif

  if (voxels > 0) {
    for (int i = 0; i < this->getTargetDim(); ++i) {
#ifdef PARALLEL_MODE_MPI
      singleton::mpi().reduceAndBcast(output[i], MPI_SUM);
#endif
      output[i] /= voxels;
    }
  }

  return true;
}

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• src/functors/lattice/superLocalAverage2D.h
• src/functors/lattice/superLocalAverage2D.hh