olb::SmoothBlockIndicator3D< T, DESCRIPTOR > Class Template Reference

#include <reductionF3D.h>

Inheritance diagram for olb::SmoothBlockIndicator3D< T, DESCRIPTOR >:

Collaboration diagram for olb::SmoothBlockIndicator3D< T, DESCRIPTOR >:

Public Member Functions
	SmoothBlockIndicator3D (IndicatorF3D< T > &f, T h, T eps, T sigma)
Public Member Functions inherited from olb::BlockDataF3D< T, T >
	BlockDataF3D (BlockData< 3, T, T > &blockData)
	Constructor.
	BlockDataF3D (BlockF3D< T > &f)
	to store functor data, constuctor creates _blockData with functor data
BlockData< 3, T, T > &	getBlockData ()
	returns _blockData
bool	operator() (T output[], const int input[]) override
	access to _blockData via its get()
Public Member Functions inherited from olb::BlockF3D< T >
	~BlockF3D () override
	virtual destructor for defined behaviour
virtual BlockStructureD< 3 > &	getBlockStructure () const
BlockF3D< T > &	operator- (BlockF3D< T > &rhs)
BlockF3D< T > &	operator+ (BlockF3D< T > &rhs)
BlockF3D< T > &	operator* (BlockF3D< T > &rhs)
BlockF3D< T > &	operator/ (BlockF3D< T > &rhs)
Public Member Functions inherited from olb::GenericF< T, 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
virtual bool	operator() (T output[], const int 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[], int input0)
bool	operator() (T output[], int input0, int input1)
bool	operator() (T output[], int input0, int input1, int input2)
bool	operator() (T output[], int input0, int input1, int input2, int input3)

Protected Attributes
const T	_h
	Lattice spacing of the particle grid.
const T	_sigma
	Important parameter for the Gaussian point spread Function (standard deviations)
const int	_eps
	Size (always a multiple of 2) of the epsilon layer eps_phys = _eps * _h.
const int	_wa
	size of the matrix of weight coefficients (from 3D Gaussian Function) _wa x _wa x _wa
IndicatorF3D< T > &	_f
	_f holds the geometry
Protected Attributes inherited from olb::BlockDataF3D< T, T >
std::unique_ptr< BlockData< 3, T, T > >	_blockDataStorage
BlockData< 3, T, T > &	_blockData
Protected Attributes inherited from olb::BlockF3D< T >
BlockStructureD< 3 > &	_blockStructure

Additional Inherited Members
Public Types inherited from olb::GenericF< T, int >
using	targetType
using	sourceType
Public Attributes inherited from olb::GenericF< T, int >
std::shared_ptr< GenericF< T, int > >	_ptrCalcC
	memory management, frees resouces (calcClass)
Protected Member Functions inherited from olb::BlockDataF3D< T, T >
	BlockDataF3D (int nx, int ny, int nz, int size=1)
Protected Member Functions inherited from olb::BlockF3D< T >
	BlockF3D (BlockStructureD< 3 > &blockStructure, int targetDim)
Protected Member Functions inherited from olb::GenericF< T, int >
	GenericF (int targetDim, int sourceDim)

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::SmoothBlockIndicator3D< T, DESCRIPTOR >

Definition at line 82 of file reductionF3D.h.

Constructor & Destructor Documentation

SmoothBlockIndicator3D()

template<typename T , typename DESCRIPTOR >

olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::SmoothBlockIndicator3D	(	IndicatorF3D< T > &	f,
		T	h,
		T	eps,
		T	sigma )

Definition at line 96 of file reductionF3D.hh.

  : BlockDataF3D<T, T>((int)((f.getMax()[0] - f.getMin()[0]) / h + ( util::round(eps*0.5)*2+2 )),
                       (int)((f.getMax()[1] - f.getMin()[1]) / h + ( util::round(eps*0.5)*2+2 )),
                       (int)((f.getMax()[2] - f.getMin()[2]) / h + ( util::round(eps*0.5)*2+2 ))),
    _h(h),
    _sigma(sigma),
    _eps(util::round(eps*0.5)*2),
    _wa(_eps+1),
    _f(f)
{
  this->getName() = "SmoothBlockIndicator3D";
  T value, dx, dy, dz;
  T weights[this->_wa][this->_wa][this->_wa];
  T sum = 0;
  const int iStart = util::floor(this->_wa*0.5);
  const int iEnd = util::ceil(this->_wa*0.5);
 
  // calculate weights: they are constants, but calculation here is less error-prone than hardcoding these parameters
  for (int x = -iStart; x < iEnd; ++x) {
    for (int y = -iStart; y < iEnd; ++y) {
      for (int z = -iStart; z < iEnd; ++z) {
        weights[x+iStart][y+iStart][z+iStart] = util::exp(-(x*x+y*y+z*z)/(2*this->_sigma*this->_sigma)) / (util::pow(this->_sigma,3)*util::sqrt(util::pow(2,3)*util::pow(M_PI,3)));
        // important because sum of all weigths only equals 1 for this->_wa -> infinity
        sum += weights[x+iStart][y+iStart][z+iStart];
      }
    }
  }
  const T invSum = 1./sum;
 
  for (int iX=0; iX<this->getBlockData().getNx(); ++iX) {
    for (int iY=0; iY<this->getBlockData().getNy(); ++iY) {
      for (int iZ=0; iZ<this->getBlockData().getNz(); ++iZ) {
        bool output[1];
        value = 0;
 
        // input: regarded point (center)
        T input[] = {
          _f.getMin()[0] + iX*_h,
          _f.getMin()[1] + iY*_h,
          _f.getMin()[2] + iZ*_h
        };
 
        /*
         * three loops to look at every point (which weight is not 0) around the regarded point
         * sum all weighted porosities
         */
        for (int x = -iStart; x < iEnd; ++x) {
          for (int y = -iStart; y < iEnd; ++y) {
            for (int z = -iStart; z < iEnd; ++z) {
              dx = x*_h;
              dy = y*_h;
              dz = z*_h;
 
              // move from regarded point to point in neighborhood
              input[0] += dx;
              input[1] += dy;
              input[2] += dz;
 
              // get porosity
              _f(output,input);
 
              // sum porosity
              value += output[0] * weights[x+iStart][y+iStart][z+iStart];
 
              // move back to center
              input[0] -= dx;
              input[1] -= dy;
              input[2] -= dz;
            }
          }
        }
        /*
         * Round to 3 decimals
         * See above sum != 1.0, that's the reason for devision, otherwise porosity will never reach 0
         */
        this->getBlockData().get({iX,iY,iZ},0) = value*invSum;//nearbyint(1000*value/sum)/1000.0;
      }
    }
  }
}

References olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_f, olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_h, olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_sigma, olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_wa, olb::util::ceil(), olb::util::exp(), olb::util::floor(), olb::BlockData< D, T, U >::get(), olb::BlockDataF3D< T, T >::getBlockData(), olb::IndicatorF3D< S >::getMin(), olb::GenericF< T, int >::getName(), olb::BlockStructureD< D >::getNx(), olb::BlockStructureD< D >::getNy(), olb::BlockStructureD< D >::getNz(), M_PI, olb::util::pow(), and olb::util::sqrt().

Here is the call graph for this function:

Member Data Documentation

_eps

template<typename T , typename DESCRIPTOR >

const int olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_eps

protected

Size (always a multiple of 2) of the epsilon layer eps_phys = _eps * _h.

Definition at line 97 of file reductionF3D.h.

_f

template<typename T , typename DESCRIPTOR >

IndicatorF3D<T>& olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_f

protected

_f holds the geometry

Definition at line 101 of file reductionF3D.h.

_h

template<typename T , typename DESCRIPTOR >

const T olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_h

protected

Lattice spacing of the particle grid.

Definition at line 93 of file reductionF3D.h.

_sigma

template<typename T , typename DESCRIPTOR >

const T olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_sigma

protected

Important parameter for the Gaussian point spread Function (standard deviations)

Definition at line 95 of file reductionF3D.h.

_wa

template<typename T , typename DESCRIPTOR >

const int olb::SmoothBlockIndicator3D< T, DESCRIPTOR >::_wa

protected

size of the matrix of weight coefficients (from 3D Gaussian Function) _wa x _wa x _wa

Definition at line 99 of file reductionF3D.h.

---

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

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