Converts block functors to analytical functors. More...

#include <interpolationF3D.h>

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

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

Public Member Functions
	AnalyticalFfromBlockF3D (BlockF3D< W > &f, Cuboid3D< T > &cuboid)

bool	operator() (W output[], const T physC[]) override
	trilinear interpolation on cubic lattice

Public Member Functions inherited from olb::AnalyticalF< D, T, S >
AnalyticalF< D, T, S > &	operator- (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator+ (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator* (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator/ (AnalyticalF< D, T, S > &rhs)

Public Member Functions inherited from olb::GenericF< T, S >
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 S 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[], S input0)

bool	operator() (T output[], S input0, S input1)

bool	operator() (T output[], S input0, S input1, S input2)

bool	operator() (T output[], S input0, S input1, S input2, S input3)

Protected Attributes
BlockF3D< W > &	_f

Cuboid3D< T > &	_cuboid

Additional Inherited Members
Public Types inherited from olb::AnalyticalF< D, T, S >
using	identity_functor_type = AnalyticalIdentity<D,T,S>

Public Types inherited from olb::GenericF< T, S >
using	targetType = T

using	sourceType = S

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

Static Public Attributes inherited from olb::AnalyticalF< D, T, S >
static constexpr unsigned	dim = D

Protected Member Functions inherited from olb::AnalyticalF< D, T, S >
	AnalyticalF (int n)

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

Detailed Description

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

Converts block functors to analytical functors.

Definition at line 53 of file interpolationF3D.h.

Constructor & Destructor Documentation

◆ AnalyticalFfromBlockF3D()

template<typename T , typename W >

olb::AnalyticalFfromBlockF3D< T, W >::AnalyticalFfromBlockF3D	(	BlockF3D< W > &	f,
		Cuboid3D< T > &	cuboid )

Definition at line 176 of file interpolationF3D.hh.

  : AnalyticalF3D<T,W>(f.getTargetDim()),
    _f(f), _cuboid(cuboid)
{
  this->getName() = "fromBlockF";
}

References olb::GenericF< T, S >::getName().

Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename W >

bool olb::AnalyticalFfromBlockF3D< T, W >::operator()	(	W	output[],
		const T	physC[] )

override

trilinear interpolation on cubic lattice

Definition at line 186 of file interpolationF3D.hh.

{
  int latticeC[3];
  auto latticeR = _cuboid.getFloorLatticeR(physC);
 
  auto& block = _f.getBlockStructure();
  auto padding = std::min(1, block.getPadding());
 
  if (latticeR >= -padding && latticeR < block.getExtent()+padding-1) {
    const int& locX = latticeR[0];
    const int& locY = latticeR[1];
    const int& locZ = latticeR[2];
 
    Vector<T,3> physCv(physC);
    Vector<T,3> physRiC = _cuboid.getPhysR({locX, locY, locZ});
 
    // compute weights
    Vector<W,3> d = (physCv - physRiC) * (1. / _cuboid.getDeltaR());
    Vector<W,3> e = 1. - d;
 
    W output_tmp[_f.getTargetDim()];
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX;
    latticeC[1] = locY;
    latticeC[2] = locZ;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * e[0] * e[1] * e[2];
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX;
    latticeC[1] = locY + 1;
    latticeC[2] = locZ;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * e[0] * d[1] * e[2];
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX + 1;
    latticeC[1] = locY;
    latticeC[2] = locZ;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * d[0] * e[1] * e[2];
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX + 1;
    latticeC[1] = locY + 1;
    latticeC[2] = locZ;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * d[0] * d[1] * e[2];
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX;
    latticeC[1] = locY;
    latticeC[2] = locZ + 1;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * e[0] * e[1] * d[2];
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX;
    latticeC[1] = locY + 1;
    latticeC[2] = locZ + 1;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * e[0] * d[1] * d[2];
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX + 1;
    latticeC[1] = locY;
    latticeC[2] = locZ + 1;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * d[0] * e[1] * d[2];
      output_tmp[iD] = W();
    }
 
    latticeC[0] = locX + 1;
    latticeC[1] = locY + 1;
    latticeC[2] = locZ + 1;
    _f(output_tmp,latticeC);
    for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
      output[iD] += output_tmp[iD] * d[0] * d[1] * d[2];
      output_tmp[iD] = W();
    }
 
    return true;
  }
  else {
    return false;
  }
}

Member Data Documentation

◆ _cuboid

template<typename T , typename W = T>

Cuboid3D<T>& olb::AnalyticalFfromBlockF3D< T, W >::_cuboid

protected

Definition at line 56 of file interpolationF3D.h.

◆ _f

template<typename T , typename W = T>

BlockF3D<W>& olb::AnalyticalFfromBlockF3D< T, W >::_f

protected

Definition at line 55 of file interpolationF3D.h.

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

src/functors/analytical/indicator/smoothIndicatorF3D.h
src/functors/analytical/interpolationF3D.h
src/functors/analytical/interpolationF3D.hh

Public Member Functions

Protected Attributes

Additional Inherited Members