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

SuperSum3D sums all components of f over a indicated subset. More...

#include <superIntegralF3D.h>

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

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

Public Member Functions
	SuperSum3D (FunctorPtr< SuperF3D< T, W > > &&f, FunctorPtr< SuperIndicatorF3D< T > > &&indicatorF)
	Constructor for summing f on a indicated subset.
	SuperSum3D (FunctorPtr< SuperF3D< T, W > > &&f, SuperGeometry< T, 3 > &superGeometry, const int material)
	Constructor for summing f on a given material.
bool	operator() (W output[], const int input[]) override
	has to be implemented for 'every' derived class
Public Member Functions inherited from olb::SuperF3D< T, T >
SuperF3D< T, W > &	operator- (SuperF3D< T, W > &rhs)
SuperF3D< T, W > &	operator+ (SuperF3D< T, W > &rhs)
SuperF3D< T, W > &	operator* (SuperF3D< T, W > &rhs)
SuperF3D< T, W > &	operator/ (SuperF3D< T, W > &rhs)
SuperStructure< T, 3 > &	getSuperStructure ()
int	getBlockFSize () const
BlockF3D< 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::SuperF3D< 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)
Static Public Attributes inherited from olb::SuperF3D< T, T >
static constexpr bool	isSuper
static constexpr unsigned	d
Protected Member Functions inherited from olb::SuperF3D< T, T >
	SuperF3D (SuperStructure< T, 3 > &superStructure, int targetDim)
Protected Member Functions inherited from olb::GenericF< W, int >
	GenericF (int targetDim, int sourceDim)
Protected Attributes inherited from olb::SuperF3D< T, T >
SuperStructure< T, 3 > &	_superStructure
std::vector< std::unique_ptr< BlockF3D< W > > >	_blockF
	Super functors may consist of several BlockF3D<W> derived functors.

Detailed Description

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

SuperSum3D sums all components of f over a indicated subset.

Definition at line 39 of file superIntegralF3D.h.

Constructor & Destructor Documentation

SuperSum3D() [1/2]

template<typename T , typename W >

olb::SuperSum3D< T, W >::SuperSum3D	(	FunctorPtr< SuperF3D< T, W > > &&	f,
		FunctorPtr< SuperIndicatorF3D< T > > &&	indicatorF )

Constructor for summing f on a indicated subset.

Parameters

f	functor to be summed
indicatorF	indicator describing the subset on which to evaluate f

Definition at line 39 of file superIntegralF3D.hh.

  : SuperF3D<T,W>(f->getSuperStructure(), f->getTargetDim()+1),
    _f(std::move(f)),
    _indicatorF(std::move(indicatorF))
{
  this->getName() = "Sum("+_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 BlockSum3D<T,W>(_f->getBlockF(iC),
                            _indicatorF->getBlockIndicatorF(iC))
      );
    }
  }
}

References olb::GenericF< W, int >::getName().

Here is the call graph for this function:

SuperSum3D() [2/2]

template<typename T , typename W >

olb::SuperSum3D< T, W >::SuperSum3D	(	FunctorPtr< SuperF3D< T, W > > &&	f,
		SuperGeometry< T, 3 > &	superGeometry,
		const int	material )

Constructor for summing f on a given material.

Parameters

f	functor to be summed
superGeometry	super geometry for constructing material indicator
material	number of the relevant material

Definition at line 61 of file superIntegralF3D.hh.

  : SuperSum3D(
      std::forward<decltype(f)>(f),
      superGeometry.getMaterialIndicator(material))
{ }

Member Function Documentation

operator()()

template<typename T , typename W >

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

overridevirtual

has to be implemented for 'every' derived class

Reimplemented from olb::SuperF3D< T, T >.

Definition at line 70 of file superIntegralF3D.hh.

{
  _f->getSuperStructure().communicate();
  auto& geometry = _f->getSuperStructure().getCuboidDecomposition();
  LoadBalancer<T>&     load     = _f->getSuperStructure().getLoadBalancer();
 
  for (int i = 0; i < this->getTargetDim(); ++i) {
    output[i] = W(0);
  }
 
  if (this->_blockF.empty()) {
    W outputTmp[_f->getTargetDim()];
    int inputTmp[_f->getSourceDim()];
    std::size_t voxels(0);
    std::vector<util::KahanSummator<W>> summators(_f->getTargetDim(), util::KahanSummator<W>());
 
    for (int iC = 0; iC < load.size(); ++iC) {
      const Cuboid3D<T> cuboid = geometry.get(load.glob(iC));
      inputTmp[0] = load.glob(iC);
      for (inputTmp[1] = 0; inputTmp[1] < cuboid.getNx(); ++inputTmp[1]) {
        for (inputTmp[2] = 0; inputTmp[2] < cuboid.getNy(); ++inputTmp[2]) {
          for (inputTmp[3] = 0; inputTmp[3] < cuboid.getNz(); ++inputTmp[3]) {
            if (_indicatorF(inputTmp)) {
              _f(outputTmp,inputTmp);
              for (int i = 0; i < _f->getTargetDim(); ++i) {
                summators[i].add(outputTmp[i]);
              }
              voxels += 1;
            }
          }
        }
      }
    }
    for (int i = 0; i < _f->getTargetDim(); ++i) {
      output[i] = summators[i].getSum();
    }
    output[_f->getTargetDim()] = voxels;
  }
  else {
    for (int iC = 0; iC < load.size(); ++iC) {
      this->getBlockF(iC)(output, input);
    }
  }
 
#ifdef PARALLEL_MODE_MPI
  for (int i = 0; i < this->getTargetDim(); ++i) {
    singleton::mpi().reduceAndBcast(output[i], MPI_SUM);
  }
#endif
  return true;
}

References olb::SuperStructure< T, D >::communicate(), olb::Cuboid< T, D >::getNx(), olb::Cuboid< T, D >::getNy(), olb::Cuboid< T, D >::getNz(), olb::SuperF3D< T, U >::getSuperStructure(), olb::LoadBalancer< T >::glob(), olb::singleton::mpi(), olb::singleton::MpiManager::reduceAndBcast(), and olb::LoadBalancer< T >::size().

Here is the call graph for this function:

---

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

• src/functors/lattice/integral/superIntegralF3D.h
• src/functors/lattice/integral/superIntegralF3D.hh