Converts super functions to analytical functions. More...

#include <interpolationF2D.h>

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

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

Public Member Functions
	AnalyticalFfromSuperF2D (SuperF2D< T > &f, bool communicateToAll=false, bool communicateOverlap=true)

bool	operator() (T output[], const T physC[]) override

int	getBlockFSize () const

AnalyticalFfromBlockF2D< T, W > &	getBlockF (int iCloc)

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
const bool	_communicateToAll

const bool	_communicateOverlap

SuperF2D< T > &	_f

CuboidDecomposition< T, 2 > &	_cuboidDecomposition

std::vector< std::unique_ptr< AnalyticalFfromBlockF2D< T, W > > >	_blockF

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::AnalyticalFfromSuperF2D< T, W >

Converts super functions to analytical functions.

Definition at line 63 of file interpolationF2D.h.

Constructor & Destructor Documentation

◆ AnalyticalFfromSuperF2D()

template<typename T , typename W >

olb::AnalyticalFfromSuperF2D< T, W >::AnalyticalFfromSuperF2D	(	SuperF2D< T > &	f,
		bool	communicateToAll = false,
		bool	communicateOverlap = true )

Definition at line 198 of file interpolationF2D.hh.

  : AnalyticalF2D<T,W>(f.getTargetDim()),
    _communicateToAll(communicateToAll),
    _communicateOverlap(communicateOverlap),
    _f(f),
    _cuboidDecomposition(_f.getSuperStructure().getCuboidDecomposition())
{
  this->getName() = "fromSuperF("+ f.getName()+")";
 
  LoadBalancer<T>& load = _f.getSuperStructure().getLoadBalancer();
  for (int iC = 0; iC < load.size(); ++iC) {
    this->_blockF.emplace_back(
      new AnalyticalFfromBlockF2D<T>(_f.getBlockF(iC),
                                     _cuboidDecomposition.get(load.glob(iC)))
    );
  }
}

References olb::AnalyticalFfromSuperF2D< T, W >::_blockF, olb::AnalyticalFfromSuperF2D< T, W >::_cuboidDecomposition, olb::AnalyticalFfromSuperF2D< T, W >::_f, olb::CuboidDecomposition< T, D >::get(), olb::SuperF2D< T, U >::getBlockF(), olb::SuperStructure< T, D >::getLoadBalancer(), olb::GenericF< T, S >::getName(), olb::SuperF2D< T, U >::getSuperStructure(), olb::LoadBalancer< T >::glob(), and olb::LoadBalancer< T >::size().

Here is the call graph for this function:

Member Function Documentation

◆ getBlockF()

template<typename T , typename W >

AnalyticalFfromBlockF2D< T, W > & olb::AnalyticalFfromSuperF2D< T, W >::getBlockF ( int iCloc )

Returns: _blockF[iCloc]

Definition at line 277 of file interpolationF2D.hh.

{
  OLB_ASSERT(size_t(iCloc) < _blockF.size() && iCloc >= 0,
             "block functor index within bounds");
  return *(_blockF[iCloc]);
}

References OLB_ASSERT.

◆ getBlockFSize()

template<typename T , typename W >

int olb::AnalyticalFfromSuperF2D< T, W >::getBlockFSize ( ) const

Returns: Size of _blockF vector

Definition at line 269 of file interpolationF2D.hh.

{
  OLB_ASSERT(_blockF.size() < UINT32_MAX,
             "it is safe to cast std::size_t to int");
  return _blockF.size();
}

References OLB_ASSERT.

◆ operator()()

template<typename T , typename W >

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

override

Definition at line 218 of file interpolationF2D.hh.

{
  const auto targetDim = _f.getTargetDim();
  for (int iD = 0; iD < targetDim; ++iD) {
    output[iD] = W();
  }
  Vector<T,2> physCV(physC);
  auto latticeR = _cuboidDecomposition.getLatticeR(physCV);
  if (!latticeR) {
    return false;
  }
 
  if (_communicateOverlap) {
    _f.getSuperStructure().communicate();
  }
 
  int dataFound = 0;
 
  const LoadBalancer<T>& load = _f.getSuperStructure().getLoadBalancer();
  for (int iC = 0; iC < load.size(); ++iC) {
    if (_blockF[iC]->operator()(output, physC)) {
      ++dataFound;
    }
  }
 
  if (_communicateToAll) {
#ifdef PARALLEL_MODE_MPI
    singleton::mpi().reduceAndBcast(dataFound, MPI_SUM);
    for (int iD = 0; iD < targetDim; ++iD) {
      singleton::mpi().reduceAndBcast(output[iD], MPI_SUM);
    }
#endif
    for (int iD = 0; iD < targetDim; ++iD) {
      output[iD]/=dataFound;
    }
  }
  else {
    if (dataFound!=0) {
      for (int iD = 0; iD < targetDim; ++iD) {
        output[iD]/=dataFound;
      }
    }
  }
 
  if (dataFound>0) {
    return true;
  }
  return false;
}

References olb::singleton::mpi(), olb::singleton::MpiManager::reduceAndBcast(), and olb::LoadBalancer< T >::size().

Here is the call graph for this function: