olb::SmoothIndicatorFactoredCuboid2D< T, S, PARTICLE > Class Template Reference

factorizable output smooth cuboid in 2D with a tangiant or ramp epsilon sector More...

#include <smoothIndicatorF2D.h>

Inheritance diagram for olb::SmoothIndicatorFactoredCuboid2D< T, S, PARTICLE >:

Collaboration diagram for olb::SmoothIndicatorFactoredCuboid2D< T, S, PARTICLE >:

Public Member Functions
	SmoothIndicatorFactoredCuboid2D (Vector< S, 2 > center, S xLength, S yLength, S epsilon, S density=0, Vector< S, 2 > vel=Vector< S, 2 >(0., 0.), S omega=0, S factor=1.)
bool	operator() (T output[], const S input[]) override
Vector< S, 2 >	calcMofiAndMass (const S density) override

Detailed Description

template<typename T, typename S, bool PARTICLE = false>
class olb::SmoothIndicatorFactoredCuboid2D< T, S, PARTICLE >

factorizable output smooth cuboid in 2D with a tangiant or ramp epsilon sector

Definition at line 186 of file smoothIndicatorF2D.h.

Constructor & Destructor Documentation

SmoothIndicatorFactoredCuboid2D()

template<typename T , typename S , bool PARTICLE>

olb::SmoothIndicatorFactoredCuboid2D< T, S, PARTICLE >::SmoothIndicatorFactoredCuboid2D	(	Vector< S, 2 >	center,
		S	xLength,
		S	yLength,
		S	epsilon,
		S	density = 0,
		Vector< S, 2 >	vel = Vector<S,2> (0.,0.),
		S	omega = 0,
		S	factor = 1. )

Definition at line 596 of file smoothIndicatorF2D.hh.

  : _xLength(xLength), _yLength(yLength), _factor(factor)
{
  this->_epsilon = epsilon;
  if constexpr (!PARTICLE) {
    this->_pos = center;
  }
 
  this->_circumRadius = .5*(util::sqrt(util::pow( xLength, 2)+util::pow( yLength, 2))) + 0.5*epsilon;
  if constexpr (!PARTICLE) {
    this->_myMin = {
      this->_pos[0] - this->getCircumRadius(),
      this->_pos[1] - this->getCircumRadius()
    };
    this->_myMax = {
      this->_pos[0] + this->getCircumRadius(),
      this->_pos[1] + this->getCircumRadius()
    };
    this->init();
  }
}

References olb::util::pow(), and olb::util::sqrt().

Here is the call graph for this function:

Member Function Documentation

calcMofiAndMass()

template<typename T , typename S , bool PARTICLE>

Vector< S, 2 > olb::SmoothIndicatorFactoredCuboid2D< T, S, PARTICLE >::calcMofiAndMass ( const S density )

override

Definition at line 619 of file smoothIndicatorF2D.hh.

{
  T mass = _xLength*_yLength*density;
  T mofi = 1./12.*mass*(_xLength*_xLength+_yLength*_yLength);
  return Vector<S,2>(mofi, mass);
}

operator()()

template<typename T , typename S , bool PARTICLE>

bool olb::SmoothIndicatorFactoredCuboid2D< T, S, PARTICLE >::operator() ( T output[], const S input[] )

override

Definition at line 628 of file smoothIndicatorF2D.hh.

{
  T d;
  if(_xLength == 0){
    T distToCenter2 = util::pow(this->getPos()[1]-input[1], 2);
    d = util::sqrt(distToCenter2) - this->_yLength/2. /*+ this->getEpsilon()/2*/;
  }
  else if(_yLength == 0){
    T distToCenter2 = util::pow(this->getPos()[0]-input[0], 2);
    d = util::sqrt(distToCenter2) - this->_xLength/2. + this->getEpsilon()/2;
  }
  else{
    T distToCenter2x = util::pow(this->getPos()[0]-input[0], 2);
    T distToCenter2y = util::pow(this->getPos()[1]-input[1], 2);
    T ratio = util::sqrt(distToCenter2y)/util::sqrt(distToCenter2x);
    if(ratio >= _yLength/_xLength){
      d = util::sqrt(distToCenter2y) - this->_yLength/2. + this->getEpsilon()/2;
    }
    else{
      d = util::sqrt(distToCenter2x) - this->_xLength/2. + this->getEpsilon()/2;
    }
  }
 
  /*if(d <= 0){
    output[0] = T(this->_factor);
    return true;
  }
  else if(d > 0 && d <= this->getEpsilon()){
    output[0] = T(- this->_factor*(d-this->getEpsilon())/this->getEpsilon());
    return true;
  }
  else{
    output[0] = 0;
    return false;
  }*/
  output[0] = T(this->_factor*(1.-tanh(d/this->getEpsilon()))/2.);
  return true;
}

References olb::util::pow(), and olb::util::sqrt().

Here is the call graph for this function:

---

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

• src/functors/analytical/indicator/smoothIndicatorF2D.h
• src/functors/analytical/indicator/smoothIndicatorF2D.hh