olb::BlockLaplacian3D< T > Class Template Reference

functor to get pointwise finite difference Laplacian operator More...

#include <latticeDerivatives3D.h>

Inheritance diagram for olb::BlockLaplacian3D< T >:

Collaboration diagram for olb::BlockLaplacian3D< T >:

Public Member Functions
	BlockLaplacian3D (BlockGeometry< T, 3 > &blockGeometry, BlockF3D< T > &blockFunctor, bool forthOrder)
bool	operator() (T output[], const int input[]) override
	has to be implemented for 'every' derived class
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)

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::BlockF3D< T >
	BlockF3D (BlockStructureD< 3 > &blockStructure, int targetDim)
Protected Member Functions inherited from olb::GenericF< T, int >
	GenericF (int targetDim, int sourceDim)
Protected Attributes inherited from olb::BlockF3D< T >
BlockStructureD< 3 > &	_blockStructure

Detailed Description

template<typename T>
class olb::BlockLaplacian3D< T >

functor to get pointwise finite difference Laplacian operator

Definition at line 99 of file latticeDerivatives3D.h.

Constructor & Destructor Documentation

BlockLaplacian3D()

template<typename T >

olb::BlockLaplacian3D< T >::BlockLaplacian3D	(	BlockGeometry< T, 3 > &	blockGeometry,
		BlockF3D< T > &	blockFunctor,
		bool	forthOrder )

Definition at line 283 of file latticeDerivatives3D.hh.

  : BlockF3D<T>(blockFunctor.getBlockStructure(), blockFunctor.getTargetDim()),
    _blockGeometry(blockGeometry), _blockFunctor(blockFunctor),
    _forthOrder(forthOrder)
{
  this->getName() = "Laplacian(" + blockFunctor.getName() + ")";
  _n[0] = this-> _blockGeometry.getNx()-1;
  _n[1] = this-> _blockGeometry.getNy()-1;
  _n[2] = this-> _blockGeometry.getNz()-1;
}

References olb::GenericF< T, S >::getName(), olb::GenericF< T, int >::getName(), olb::BlockStructureD< D >::getNx(), olb::BlockStructureD< D >::getNy(), and olb::BlockStructureD< D >::getNz().

Here is the call graph for this function:

Member Function Documentation

operator()()

template<typename T >

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

overridevirtual

has to be implemented for 'every' derived class

Implements olb::GenericF< T, int >.

Definition at line 297 of file latticeDerivatives3D.hh.

{
  for (int i=0; i<this->getTargetDim(); ++i) {
    output[i] = 0;
  }
  int inputMod[3];
  for (int j=0; j<3; ++j) {
    inputMod[j] = input[j];
  }
  T u_0[this->getTargetDim()];
  T u_m[this->getTargetDim()];
  T u_p[this->getTargetDim()];
 
  // only used for forth order dq
  T u_2m[this->getTargetDim()];
  T u_2p[this->getTargetDim()];
 
  // compute discrete second derivatives for each direction and add them
  for (int j=0; j<3; ++j) {
    if (input[j] < 1) {
      // use forward difference quotient at the boundary
      _blockFunctor.operator()(u_m, inputMod);
      inputMod[j] += 1;
      _blockFunctor.operator()(u_0, inputMod);
      inputMod[j] += 1;
      _blockFunctor.operator()(u_p, inputMod);
    } else if (input[j] > _n[j]-1) {
      // use backward difference quotient at the boundary
      _blockFunctor.operator()(u_p, inputMod);
      inputMod[j] -= 1;
      _blockFunctor.operator()(u_0, inputMod);
      inputMod[j] -= 1;
      _blockFunctor.operator()(u_m, inputMod);
    } else {
      // use central difference quotient
      _blockFunctor.operator()(u_0, inputMod);
      inputMod[j] -= 1;
      _blockFunctor.operator()(u_m, inputMod);
      inputMod[j] += 2;
      _blockFunctor.operator()(u_p, inputMod);
    }
    inputMod[j] = input[j];  // reset
 
    if ((_forthOrder) && (input[j] > 1) && (input[j] < _n[j]-1)) {
      // additional evaluations
      inputMod[j] -= 2;
      _blockFunctor.operator()(u_2m, inputMod);
      inputMod[j] += 4;
      _blockFunctor.operator()(u_2p, inputMod);
      inputMod[j] = input[j];  // reset
 
      // forth order dq
      for (int i=0; i<this->getTargetDim(); ++i) {
        output[i] += fd::centralSecondDeriv(u_2m[i], u_m[i], u_0[i], u_p[i], u_2p[i]);
      }
    } else {
      // second order dq
      for (int i=0; i<this->getTargetDim(); ++i) {
        output[i] += fd::centralSecondDeriv(u_m[i], u_0[i], u_p[i]);
      }
    }
  }
  return true;
/*
  // todo: enable treating all spatial directions at once?
  _blockFunctor.operator()(u_0, inputMod);
 
  inputMod[0] -= 1;
  _blockFunctor.operator()(u_xm1, inputMod);
  inputMod[0] += 2;
  _blockFunctor.operator()(u_xp1, inputMod);
  inputMod[0] -= 1;
 
  inputMod[1] -= 1;
  _blockFunctor.operator()(u_ym1, inputMod);
  inputMod[1] += 2;
  _blockFunctor.operator()(u_yp1, inputMod);
  inputMod[1] -= 1;
 
  inputMod[2] -= 1;
  _blockFunctor.operator()(u_zm1, inputMod);
  inputMod[2] += 2;
  _blockFunctor.operator()(u_zp1, inputMod);
  inputMod[2] -= 1;
 
  for (int i=0; i<_targetDim; ++i) {
    output[i] = fd::laplacian3D(u_xm1[i], u_ym1[i], u_zm1[i], u_0[i], u_xp1[i], u_yp1[i], u_zp1[i]);
  }
*/
}

References olb::fd::centralSecondDeriv().

Here is the call graph for this function:

---

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

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