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.
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.
Public Member Functions inherited from olb::SuperF2D< T, T >
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)
SuperStructure< T, 2 > &	getSuperStructure ()
int	getBlockFSize () const
BlockF2D< W > &	getBlockF (int iCloc)
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)

Additional Inherited Members
Public Types inherited from olb::SuperF2D< T, T >
using	identity_functor_type
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::SuperF2D< T, T >
	SuperF2D (SuperStructure< T, 2 > &superStructure, int targetDim)
Protected Member Functions inherited from olb::GenericF< W, int >
	GenericF (int targetDim, int sourceDim)
Protected Attributes inherited from olb::SuperF2D< T, T >
SuperStructure< T, 2 > &	_superStructure
std::vector< std::unique_ptr< BlockF2D< W > > >	_blockF
	Super functors may consist of several BlockF2D<W> derived functors.

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

SuperLocalAverage2D()

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.

  : 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)
      );
    }
  }
}

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

Here is the call graph for this function:

Member Function Documentation

operator()()

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, T >.

Definition at line 61 of file superLocalAverage2D.hh.

{
  const auto& geometry = this->getSuperStructure().getCuboidDecomposition();
  const auto& load     = this->getSuperStructure().getLoadBalancer();
 
  for (int i = 0; i < this->getTargetDim(); ++i) {
    output[i] = 0.;
  }
 
  if (!_indicatorF(input)) {
    return true;
  }
  LatticeR<3> latticeR(input);
  auto centerOfCircle = geometry.getPhysR(latticeR);
  IndicatorCircle2D<T> analyticalCircle(centerOfCircle.data(), _radius);
  SuperIndicatorFfromIndicatorF2D<T> latticeCircle(
    analyticalCircle,
    _indicatorF->getSuperGeometry());
 
  std::size_t voxels(0);
  int inputTmp[3];
  std::vector<util::KahanSummator<W>> summators(_f->getTargetDim(), util::KahanSummator<W>());
 
  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) {
            summators[i].add(outputTmp[i]);
          }
          voxels += 1;
        }
      }
    }
  }
 
  for (int i = 0; i < _f->getTargetDim(); ++i) {
    output[i] = summators[i].getSum();
  }
 
#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;
}

References olb::singleton::mpi(), and olb::singleton::MpiManager::reduceAndBcast().

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