olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE > Class Template Reference

#include <fdBoundaryPostProcessors3D.h>

Inheritance diagram for olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >:

Collaboration diagram for olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >:

Public Types
using	parameters = PARAMS

Public Member Functions
	FdBoundaryPostProcessor3D ()
int	getPriority () const
template<typename CELL , typename PARAMETERS >
void	apply (CELL &cell, PARAMETERS &vars) any_platform

Static Public Attributes
static constexpr OperatorScope	scope = OperatorScope::PerCellWithParameters

Additional Inherited Members
Protected Member Functions inherited from olb::FdBasePostProcessor3D< T, DESCRIPTOR, FIELD, SOURCE >
	FdBasePostProcessor3D ()
template<typename CELL >
void	applySourceTerm (T *fNew, CELL &cell) any_platform

Detailed Description

template<typename T, typename DESCRIPTOR, typename MODEL, typename SCHEME_BOUND, typename PARAMS, typename FIELD = descriptors::AD_FIELD, typename SOURCE = void>
class olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >

Definition at line 37 of file fdBoundaryPostProcessors3D.h.

Member Typedef Documentation

parameters

template<typename T , typename DESCRIPTOR , typename MODEL , typename SCHEME_BOUND , typename PARAMS , typename FIELD = descriptors::AD_FIELD, typename SOURCE = void>

using olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >::parameters = PARAMS

Definition at line 39 of file fdBoundaryPostProcessors3D.h.

Constructor & Destructor Documentation

FdBoundaryPostProcessor3D()

template<typename T , typename DESCRIPTOR , typename MODEL , typename SCHEME_BOUND , typename PARAMS , typename FIELD , typename SOURCE >

olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >::FdBoundaryPostProcessor3D

(

)

Definition at line 33 of file fdBoundaryPostProcessors3D.hh.

  :  FdBasePostProcessor3D<T,DESCRIPTOR,FIELD,SOURCE>()
{
  OLB_PRECONDITION( DESCRIPTOR::template provides<descriptors::NORMAL_X>() );
  OLB_PRECONDITION( DESCRIPTOR::template provides<descriptors::NORMAL_Y>() );
  OLB_PRECONDITION( DESCRIPTOR::template provides<descriptors::NORMAL_Z>() );
}

References OLB_PRECONDITION.

Member Function Documentation

apply()

template<typename T , typename DESCRIPTOR , typename MODEL , typename SCHEME_BOUND , typename PARAMS , typename FIELD , typename SOURCE >

template<typename CELL , typename PARAMETERS >

void olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >::apply	(	CELL &	cell,
		PARAMETERS &	vars )

Definition at line 49 of file fdBoundaryPostProcessors3D.hh.

{
  int normalX = (int)( cell.template getField<descriptors::NORMAL_X>() );
  int normalY = (int)( cell.template getField<descriptors::NORMAL_Y>() );
  int normalZ = (int)( cell.template getField<descriptors::NORMAL_Z>() );
  if (normalX!=0 || normalY!=0 || normalZ!=0) { // does absolutely nothing if the cell does not belong to a border
    constexpr int extraExtent = MODEL::extent() + SCHEME_BOUND::getExtraExtent();
    std::size_t iT = vars.template get<fd::fdParams::Timestep>();
    T f[MODEL::extent()*DESCRIPTOR::d], F[MODEL::extent()*DESCRIPTOR::d];
    T fGhost[MODEL::extent()*DESCRIPTOR::d], fNormal[extraExtent*DESCRIPTOR::d];
    T u[DESCRIPTOR::d];
    int normal[] {normalX, normalY, normalZ};
    std::fill(std::begin(fGhost), std::end(fGhost), 0);
    T* fNew = fd::accessNew<T,FIELD>(cell, iT);
    T* f0   = fd::accessOld<T,FIELD>(cell, iT);
    for (int iN=1; iN<=extraExtent; ++iN) {
      fNormal[fd::getArrayPos<extraExtent>(iN-1,0)] = *fd::accessOld<T,FIELD>( cell.neighbor({normalX*iN,           0,           0}), iT );
      fNormal[fd::getArrayPos<extraExtent>(iN-1,1)] = *fd::accessOld<T,FIELD>( cell.neighbor({          0, normalY*iN,           0}), iT );
      fNormal[fd::getArrayPos<extraExtent>(iN-1,2)] = *fd::accessOld<T,FIELD>( cell.neighbor({          0,           0, normalZ*iN}), iT );
    }
    cell.computeU(u);
    SCHEME_BOUND::apply(fGhost, *f0, fNormal, normal, u, vars);
    for (int iN=1; iN<=MODEL::extent(); ++iN) {
      f[fd::getArrayPos<MODEL::extent()>(iN-1,0)] = normalX<=0 ? *fd::accessOld<T,FIELD>( cell.neighbor({-iN,   0,   0}), iT )
                                                               : fGhost[fd::getArrayPos<MODEL::extent()+SCHEME_BOUND::getExtraExtent()>(iN-1,0)];
      F[fd::getArrayPos<MODEL::extent()>(iN-1,0)] = normalX>=0 ? *fd::accessOld<T,FIELD>( cell.neighbor({+iN,   0,   0}), iT )
                                                               : fGhost[fd::getArrayPos<MODEL::extent()+SCHEME_BOUND::getExtraExtent()>(iN-1,0)];
      f[fd::getArrayPos<MODEL::extent()>(iN-1,1)] = normalY<=0 ? *fd::accessOld<T,FIELD>( cell.neighbor({  0, -iN,   0}), iT )
                                                               : fGhost[fd::getArrayPos<MODEL::extent()+SCHEME_BOUND::getExtraExtent()>(iN-1,1)];
      F[fd::getArrayPos<MODEL::extent()>(iN-1,1)] = normalY>=0 ? *fd::accessOld<T,FIELD>( cell.neighbor({  0, +iN,   0}), iT )
                                                               : fGhost[fd::getArrayPos<MODEL::extent()+SCHEME_BOUND::getExtraExtent()>(iN-1,1)];
      f[fd::getArrayPos<MODEL::extent()>(iN-1,2)] = normalZ<=0 ? *fd::accessOld<T,FIELD>( cell.neighbor({  0,   0, -iN}), iT )
                                                               : fGhost[fd::getArrayPos<MODEL::extent()+SCHEME_BOUND::getExtraExtent()>(iN-1,2)];
      F[fd::getArrayPos<MODEL::extent()>(iN-1,2)] = normalZ>=0 ? *fd::accessOld<T,FIELD>( cell.neighbor({  0,   0, +iN}), iT )
                                                               : fGhost[fd::getArrayPos<MODEL::extent()+SCHEME_BOUND::getExtraExtent()>(iN-1,2)];
    }
    MODEL::apply(fNew, f0, f, F, u, vars);
    this->applySourceTerm(fNew, cell);
  }
}

getPriority()

template<typename T , typename DESCRIPTOR , typename MODEL , typename SCHEME_BOUND , typename PARAMS , typename FIELD , typename SOURCE >

int olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >::getPriority

(

)

const

Definition at line 42 of file fdBoundaryPostProcessors3D.hh.

{
  return 1;
}

Member Data Documentation

scope

template<typename T , typename DESCRIPTOR , typename MODEL , typename SCHEME_BOUND , typename PARAMS , typename FIELD = descriptors::AD_FIELD, typename SOURCE = void>

constexpr OperatorScope olb::FdBoundaryPostProcessor3D< T, DESCRIPTOR, MODEL, SCHEME_BOUND, PARAMS, FIELD, SOURCE >::scope = OperatorScope::PerCellWithParameters

staticconstexpr

Definition at line 40 of file fdBoundaryPostProcessors3D.h.

---

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

• src/reaction/explicitFiniteDifference/boundary/fdBoundaryPostProcessors3D.h
• src/reaction/explicitFiniteDifference/boundary/fdBoundaryPostProcessors3D.hh