olb::boundary::FullSlipO2D< NX, NY > Struct Template Reference

#include <slip2D.hh>

Collaboration diagram for olb::boundary::FullSlipO2D< NX, NY >:

Public Member Functions
int	getPriority () const
template<typename CELL , typename V = typename CELL::value_t>
void	apply (CELL &x_b) any_platform

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

Detailed Description

template<int NX, int NY>
struct olb::boundary::FullSlipO2D< NX, NY >

Definition at line 38 of file slip2D.hh.

Member Function Documentation

apply()

template<int NX, int NY>

template<typename CELL , typename V = typename CELL::value_t>

void olb::boundary::FullSlipO2D< NX, NY >::apply ( CELL & x_b )

inline

Definition at line 47 of file slip2D.hh.

                                     {
    using DESCRIPTOR = typename CELL::descriptor_t;
    int reflectionPop[DESCRIPTOR::q];
    int mirrorDirection0;
    int mirrorDirection1;
    int mult = 2 / (NX*NX + NY*NY);
    reflectionPop[0] =0;
    for (int iPop = 1; iPop < DESCRIPTOR::q; iPop++) {
      reflectionPop[iPop] = 0;
      // iPop are the directions which pointing into the fluid, discreteNormal is pointing outwarts
      int scalarProduct = descriptors::c<DESCRIPTOR>(iPop,0)*NX + descriptors::c<DESCRIPTOR>(iPop,1)*NY;
      if ( scalarProduct < 0) {
        // bounce back for the case discreteNormalX = discreteNormalY = 1, that is mult=1
        if (mult == 1) {
          mirrorDirection0 = -descriptors::c<DESCRIPTOR>(iPop,0);
          mirrorDirection1 = -descriptors::c<DESCRIPTOR>(iPop,1);
        }
        else {
          mirrorDirection0 = descriptors::c<DESCRIPTOR>(iPop,0) - mult*scalarProduct*NX;
          mirrorDirection1 = descriptors::c<DESCRIPTOR>(iPop,1) - mult*scalarProduct*NY;
        }
 
        // run through all lattice directions and look for match of direction
        for (int i = 1; i < DESCRIPTOR::q; i++) {
          if (descriptors::c<DESCRIPTOR>(i,0)==mirrorDirection0
              && descriptors::c<DESCRIPTOR>(i,1)==mirrorDirection1) {
            reflectionPop[iPop] = i;
            break;
          }
        }
      }
    }
    for (int iPop = 1; iPop < DESCRIPTOR::q ; ++iPop) {
      if (reflectionPop[iPop]!=0) {
        //do reflection
        x_b[iPop] = x_b[reflectionPop[iPop]];
      }
    }
  }

References olb::descriptors::c().

Here is the call graph for this function:

getPriority()

template<int NX, int NY>

int olb::boundary::FullSlipO2D< NX, NY >::getPriority ( ) const

inline

Definition at line 41 of file slip2D.hh.

                          {
    return -1;
  }

Member Data Documentation

scope

template<int NX, int NY>

OperatorScope olb::boundary::FullSlipO2D< NX, NY >::scope = OperatorScope::PerCell

staticconstexpr

Definition at line 39 of file slip2D.hh.

---

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

• src/boundary/slip2D.h
• src/boundary/slip2D.hh