olb::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR > Class Template Reference

functor to get the pointwise dual dissipation density on local lattices, if globIC is not on the local processor, the returned vector is empty More...

#include <dualFunctors3D.h>

Inheritance diagram for olb::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR >:

Collaboration diagram for olb::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR >:

Public Member Functions
	BlockLatticeDphysDissipationDf3D (BlockLattice< T, DESCRIPTOR > &blockLattice, int overlap, const UnitConverter< T, DESCRIPTOR > &converter)
bool	operator() (T output[], const int input[])
	has to be implemented for 'every' derived class
Public Member Functions inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
BlockLattice< T, DESCRIPTOR > &	getBlock ()
	Copy Constructor.
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
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 int	_overlap
const UnitConverter< T, DESCRIPTOR > &	_converter
Protected Attributes inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
BlockLattice< T, DESCRIPTOR > &	_blockLattice
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::BlockLatticeF3D< T, DESCRIPTOR >
	BlockLatticeF3D (BlockLattice< T, DESCRIPTOR > &blockLattice, int targetDim)
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::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR >

functor to get the pointwise dual dissipation density on local lattices, if globIC is not on the local processor, the returned vector is empty

Definition at line 38 of file dualFunctors3D.h.

Constructor & Destructor Documentation

BlockLatticeDphysDissipationDf3D()

template<typename T , typename DESCRIPTOR >

olb::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR >::BlockLatticeDphysDissipationDf3D	(	BlockLattice< T, DESCRIPTOR > &	blockLattice,
		int	overlap,
		const UnitConverter< T, DESCRIPTOR > &	converter )

Definition at line 36 of file dualFunctors3D.hh.

  : BlockLatticeF3D<T,DESCRIPTOR>(blockLattice,DESCRIPTOR::q),
    _overlap(overlap),
    _converter(converter)
{
  this->getName() = "dPhysDissipationDf";
}

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

Here is the call graph for this function:

Member Function Documentation

operator()()

template<typename T , typename DESCRIPTOR >

bool olb::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR >::operator() ( T output[], const int input[] )

virtual

has to be implemented for 'every' derived class

Implements olb::GenericF< T, int >.

Definition at line 48 of file dualFunctors3D.hh.

{
  T omega = _converter.getLatticeRelaxationFrequency();
  T dt = _converter.getConversionFactorTime();
 
  T rho;
  T u[3];
  T pi[6];
 
  /*const int localLatticeR[3] = {
    latticeR[0] + _overlap,
    latticeR[1] + _overlap,
    latticeR[2] + _overlap
  };
  this->_blockLattice.get(localLatticeR).computeAllMomenta(rho, u, pi);*/
  this->_blockLattice.get(latticeR).computeAllMomenta(rho, u, pi);
 
  for (int i=0; i<DESCRIPTOR::q; i++) { //output default 0
    dDdf[i] = T();
  }
 
  T pi2[DESCRIPTOR::d][DESCRIPTOR::d];
 
  pi2[0][0]=pi[0];
  pi2[0][1]=pi[1];
  if (util::TensorVal<DESCRIPTOR >::n == 6) {
    pi2[0][2]=pi[2];
  }
  else {
    pi2[1][1]=pi[2];
  }
  pi2[1][0]=pi[1];
  if (util::TensorVal<DESCRIPTOR >::n == 6) {
    pi2[1][1]=pi[3];
    pi2[1][2]=pi[4];
    pi2[2][0]=pi[2];
    pi2[2][1]=pi[4];
    pi2[2][2]=pi[5];
  }
 
  for (int jPop=0; jPop < DESCRIPTOR::q; ++jPop) {
    for (int iDim=0; iDim < DESCRIPTOR::d; ++iDim) {
      for (int jDim=0; jDim < DESCRIPTOR::d; ++jDim) {
        T dPidf2ndTerm = T();
        for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
          dPidf2ndTerm += descriptors::c<DESCRIPTOR>(iPop)[iDim] * descriptors::c<DESCRIPTOR>(iPop)[jDim]
                          * dualLbHelpers<T,DESCRIPTOR>::equilibrium(iPop, jPop, rho, u);
        }
        T dpidf = descriptors::c<DESCRIPTOR>(jPop,iDim)*descriptors::c<DESCRIPTOR>(jPop,jDim) - dPidf2ndTerm;
        // rho2 instead of rho vs simon!
        dDdf[jPop] += 2.*(pi2[iDim][jDim]*(dpidf - pi2[iDim][jDim]/2./rho*2.));
      }
    }
    dDdf[jPop] *= util::pow(omega*descriptors::invCs2<T,DESCRIPTOR>()/rho,2)/2.*_converter.getPhysViscosity()/_converter.getLatticeViscosity()/dt/dt;
  }
  return true;
}

References olb::opti::dualLbHelpers< T, DESCRIPTOR >::equilibrium(), and olb::util::pow().

Here is the call graph for this function:

Member Data Documentation

_converter

template<typename T , typename DESCRIPTOR >

const UnitConverter<T,DESCRIPTOR>& olb::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR >::_converter

protected

Definition at line 41 of file dualFunctors3D.h.

_overlap

template<typename T , typename DESCRIPTOR >

const int olb::opti::BlockLatticeDphysDissipationDf3D< T, DESCRIPTOR >::_overlap

protected

Definition at line 40 of file dualFunctors3D.h.

---

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

• src/optimization/dualFunctors3D.h
• src/optimization/dualFunctors3D.hh