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

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

#include <smoothIndicatorF2D.h>

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

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

Public Member Functions
	SmoothIndicatorFactoredCircle2D (Vector< S, 2 > center, S radius, 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::SmoothIndicatorFactoredCircle2D< T, S, PARTICLE >

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

Definition at line 174 of file smoothIndicatorF2D.h.

Constructor & Destructor Documentation

SmoothIndicatorFactoredCircle2D()

template<typename T , typename S , bool PARTICLE>

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

Definition at line 544 of file smoothIndicatorF2D.hh.

  : _radius(radius), _factor(factor)
{
  this->_epsilon = epsilon;
  if constexpr (!PARTICLE) {
    this->_pos = center;
  }
 
  this->_circumRadius = radius + 0.5*epsilon;
  if constexpr (!PARTICLE) {
    this->_myMin = {center[0] - this->getCircumRadius(), center[1] - this->getCircumRadius()};
    this->_myMax = {center[0] + this->getCircumRadius(), center[1] + this->getCircumRadius()};
    this->init();
  }
}

Member Function Documentation

calcMofiAndMass()

template<typename T , typename S , bool PARTICLE>

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

override

Definition at line 561 of file smoothIndicatorF2D.hh.

{
  T mass = M_PI*_radius*_radius*density;
  T mofi = 0.5 * mass * _radius * _radius;
  return Vector<S,2>(mofi, mass);
}

References M_PI.

operator()()

template<typename T , typename S , bool PARTICLE>

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

override

Definition at line 570 of file smoothIndicatorF2D.hh.

{
  double distToCenter2 = util::pow(this->getPos()[0]-input[0], 2) +
                         util::pow(this->getPos()[1]-input[1], 2);
  /*double d = util::sqrt(distToCenter2) - this->_radius + 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;
  }*/
 
  double d = util::sqrt(distToCenter2) - this->_radius;
  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