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

factorizable output smooth sphere in 3D with a tangiant or ramp epsilon sector More...

#include <smoothIndicatorF3D.h>

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

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

Public Member Functions
	SmoothIndicatorFactoredCircle3D (Vector< S, 3 > center, S radius, S epsilon, S density=0, Vector< S, 3 > vel=Vector< S, 3 >(0., 0., 0.), S omega=0, S factor=1.)
bool	operator() (T output[], const S input[]) override
Vector< S, 4 >	calcMofiAndMass (const S density) override

Detailed Description

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

factorizable output smooth sphere in 3D with a tangiant or ramp epsilon sector

Definition at line 196 of file smoothIndicatorF3D.h.

Constructor & Destructor Documentation

SmoothIndicatorFactoredCircle3D()

template<typename T , typename S , bool PARTICLE>

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

Definition at line 951 of file smoothIndicatorF3D.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(), center[2] - this->getCircumRadius()};
    this->_myMax = {center[0] + this->getCircumRadius(), center[1] + this->getCircumRadius(), center[2] + this->getCircumRadius()};
    this->init();
  }
}

Member Function Documentation

calcMofiAndMass()

template<typename T , typename S , bool PARTICLE>

Vector< S, 4 > olb::SmoothIndicatorFactoredCircle3D< T, S, PARTICLE >::calcMofiAndMass ( const S density )

override

Definition at line 968 of file smoothIndicatorF3D.hh.

{
  T mass = 4/3*M_PI*_radius*_radius*_radius*density;
  Vector<S,3> mofi;
  mofi[0] = 2./5.*mass*_radius*_radius;
  mofi[1] = 2./5.*mass*_radius*_radius;
  mofi[2] = 2./5.*mass*_radius*_radius;
  return Vector<S,4>(mofi[0], mofi[1], mofi[2], mass);
}

References M_PI.

operator()()

template<typename T , typename S , bool PARTICLE>

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

override

Definition at line 980 of file smoothIndicatorF3D.hh.

{
  double distToCenter2 = util::pow(this->getPos()[0]-input[0], 2) +
                         util::pow(this->getPos()[1]-input[1], 2) +
                         util::pow(this->getPos()[2]-input[2], 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/smoothIndicatorF3D.h
• src/functors/analytical/indicator/smoothIndicatorF3D.hh