olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, true > Struct Template Reference

#include <finiteDifference3D.h>

Collaboration diagram for olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, true >:

Static Public Member Functions
static bool	canInterpolateVector (BlockLattice< T, DESCRIPTOR > const &blockLattice, int iX, int iY, int iZ)
static void	interpolateVector (T velDeriv[DESCRIPTOR::d], BlockLattice< T, DESCRIPTOR > const &blockLattice, int iX, int iY, int iZ)
static void	interpolateScalar (T &rhoDeriv, BlockLattice< T, DESCRIPTOR > const &blockLattice, int iX, int iY, int iZ)
template<typename CELL >
static void	interpolateVector (T velDeriv[DESCRIPTOR::d], CELL &cell) any_platform

Detailed Description

template<typename T, typename DESCRIPTOR, int direction, int orientation, int deriveDirection>
struct olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, true >

Definition at line 75 of file finiteDifference3D.h.

Member Function Documentation

canInterpolateVector()

template<typename T , typename DESCRIPTOR , int direction, int orientation, int deriveDirection>

static bool olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, true >::canInterpolateVector ( BlockLattice< T, DESCRIPTOR > const & blockLattice, int iX, int iY, int iZ )

inlinestatic

Definition at line 77 of file finiteDifference3D.h.

  {
    return blockLattice.isInside(iX,iY,iZ)
           && blockLattice.isInside(iX+(direction==0 ? (-orientation):0),
                                    iY+(direction==1 ? (-orientation):0),
                                    iZ+(direction==2 ? (-orientation):0))
           && blockLattice.isInside(iX+(direction==0 ? (-2*orientation):0),
                                    iY+(direction==1 ? (-2*orientation):0),
                                    iZ+(direction==2 ? (-2*orientation):0));
  }

References olb::BlockStructureD< D >::isInside().

Here is the call graph for this function:

interpolateScalar()

template<typename T , typename DESCRIPTOR , int direction, int orientation, int deriveDirection>

static void olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, true >::interpolateScalar ( T & rhoDeriv, BlockLattice< T, DESCRIPTOR > const & blockLattice, int iX, int iY, int iZ )

inlinestatic

Definition at line 112 of file finiteDifference3D.h.

  {
    using namespace fd;
 
    // note that the derivative runs along direction.
    T rho0 = blockLattice.get(iX,iY,iZ).computeRho();
    T rho1 = blockLattice.get (
               iX+(direction==0 ? (-orientation):0),
               iY+(direction==1 ? (-orientation):0),
               iZ+(direction==2 ? (-orientation):0)  ).computeRho();
    T rho2 = blockLattice.get (
               iX+(direction==0 ? (-2*orientation):0),
               iY+(direction==1 ? (-2*orientation):0),
               iZ+(direction==2 ? (-2*orientation):0) ).computeRho();
 
    rhoDeriv = -orientation * boundaryGradient(rho0, rho1, rho2);
  }

References olb::fd::boundaryGradient(), and olb::BlockLattice< T, DESCRIPTOR >::get().

Here is the call graph for this function:

interpolateVector() [1/2]

template<typename T , typename DESCRIPTOR , int direction, int orientation, int deriveDirection>

static void olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, true >::interpolateVector ( T velDeriv[DESCRIPTOR::d], BlockLattice< T, DESCRIPTOR > const & blockLattice, int iX, int iY, int iZ )

inlinestatic

Definition at line 89 of file finiteDifference3D.h.

  {
    using namespace fd;
 
    T u0[DESCRIPTOR::d], u1[DESCRIPTOR::d], u2[DESCRIPTOR::d];
 
    blockLattice.get(iX,iY,iZ).computeU(u0);
    blockLattice.get (
      iX+(direction==0 ? (-orientation):0),
      iY+(direction==1 ? (-orientation):0),
      iZ+(direction==2 ? (-orientation):0)  ).computeU(u1);
    blockLattice.get (
      iX+(direction==0 ? (-2*orientation):0),
      iY+(direction==1 ? (-2*orientation):0),
      iZ+(direction==2 ? (-2*orientation):0) ).computeU(u2);
 
    for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
      velDeriv[iD] = -orientation * boundaryGradient(u0[iD], u1[iD], u2[iD]);
    }
  }

References olb::fd::boundaryGradient(), and olb::BlockLattice< T, DESCRIPTOR >::get().

Here is the call graph for this function:

interpolateVector() [2/2]

template<typename T , typename DESCRIPTOR , int direction, int orientation, int deriveDirection>

template<typename CELL >

static void olb::fd::DirectedGradients3D< T, DESCRIPTOR, direction, orientation, deriveDirection, true >::interpolateVector ( T velDeriv[DESCRIPTOR::d], CELL & cell )

inlinestatic

Definition at line 133 of file finiteDifference3D.h.

  {
    using namespace fd;
 
    T u0[DESCRIPTOR::d], u1[DESCRIPTOR::d], u2[DESCRIPTOR::d];
 
    cell.computeU(u0);
    cell.neighbor({(direction==0 ? (-orientation):0),
                   (direction==1 ? (-orientation):0),
                   (direction==2 ? (-orientation):0)}).computeU(u1);
    cell.neighbor({(direction==0 ? (-2*orientation):0),
                   (direction==1 ? (-2*orientation):0),
                   (direction==2 ? (-2*orientation):0)}).computeU(u2);
 
    for (int iD=0; iD<DESCRIPTOR::d; ++iD) {
      velDeriv[iD] = -orientation * boundaryGradient(u0[iD], u1[iD], u2[iD]);
    }
  }

References olb::fd::boundaryGradient().

Here is the call graph for this function:

---

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

• src/utilities/finiteDifference3D.h