SmoothIndicatorF3D is an application from $\Omega \subset R^3 \to [0,1]$ . More...

#include <smoothIndicatorBaseF3D.h>

Inheritance diagram for olb::SmoothIndicatorF3D< T, S, false >:

Collaboration diagram for olb::SmoothIndicatorF3D< T, S, false >:

Public Member Functions
void	init ()

const Vector< S, 3 > &	getMin () const

const Vector< S, 3 > &	getMax () const

const Vector< S, 3 > &	getPos () const

const Vector< S, 9 > &	getRotationMatrix () const

const Vector< S, 3 > &	getTheta () const

const S &	getCircumRadius () const

const S &	getEpsilon () const

std::string	name ()

void	setPos (Vector< S, 3 > pos)

void	setTheta (Vector< S, 3 > theta)

void	setEpsilon (S epsilon)

virtual S	getVolume ()

virtual Vector< S, 4 >	calcMofiAndMass (S density)

virtual Vector< S, 3 >	calcCenterOfMass ()

virtual Vector< S, 3 >	surfaceNormal (const Vector< S, 3 > &pos, const S meshSize)

virtual Vector< S, 3 >	surfaceNormal (const Vector< S, 3 > &pos, const S meshSize, std::function< Vector< S, 3 >(const Vector< S, 3 > &)> transformPos)

virtual const S	signedDistance (const PhysR< T, 3 > input)

virtual bool	distance (S &distance, const Vector< S, 3 > &origin, const Vector< S, 3 > &direction, S precision, S pitch)

virtual bool	operator() (T output[], const S input[])
	has to be implemented for 'every' derived class

bool	isInsideCircumRadius (const PhysR< S, 3 > &input)

SmoothIndicatorF3D< T, S, false > &	operator+ (SmoothIndicatorF3D< T, S, false > &rhs)

Public Member Functions inherited from olb::AnalyticalF< D, T, S >
AnalyticalF< D, T, S > &	operator- (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator+ (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator* (AnalyticalF< D, T, S > &rhs)

AnalyticalF< D, T, S > &	operator/ (AnalyticalF< D, T, S > &rhs)

Public Member Functions inherited from olb::GenericF< T, S >
virtual	~GenericF ()=default

int	getSourceDim () const
	read only access to member variable _m

int	getTargetDim () const
	read only access to member variable _n

std::string &	getName ()
	read and write access to name

std::string const &	getName () const
	read only access to name

bool	operator() (T output[])
	wrapper that call the pure virtual operator() (T output[], const S input[]) from above

bool	operator() (T output[], S input0)

bool	operator() (T output[], S input0, S input1)

bool	operator() (T output[], S input0, S input1, S input2)

bool	operator() (T output[], S input0, S input1, S input2, S input3)

Protected Member Functions
	SmoothIndicatorF3D ()

Protected Member Functions inherited from olb::AnalyticalF< D, T, S >
	AnalyticalF (int n)

Protected Member Functions inherited from olb::GenericF< T, S >
	GenericF (int targetDim, int sourceDim)

Protected Attributes
Vector< S, 3 >	_myMin

Vector< S, 3 >	_myMax

Vector< S, 3 >	_pos

Vector< S, 9 >	_rotMat

S	_circumRadius

Vector< S, 3 >	_theta

S	_epsilon

std::string	_name = "smoothIndicator3D"

Additional Inherited Members
Public Types inherited from olb::AnalyticalF< D, T, S >
using	identity_functor_type = AnalyticalIdentity<D,T,S>

Public Types inherited from olb::GenericF< T, S >
using	targetType = T

using	sourceType = S

Public Attributes inherited from olb::GenericF< T, S >
std::shared_ptr< GenericF< T, S > >	_ptrCalcC
	memory management, frees resouces (calcClass)

Static Public Attributes inherited from olb::AnalyticalF< D, T, S >
static constexpr unsigned	dim = D

Detailed Description

template<typename T, typename S>
class olb::SmoothIndicatorF3D< T, S, false >

SmoothIndicatorF3D is an application from $\Omega \subset R^3 \to [0,1]$ .

Parameters

_myMin	holds minimal(component wise) vector of the domain $\Omega$ .
_myMax	holds maximal(component wise) vector of the domain $\Omega$ .
_center
_diam_

Definition at line 47 of file smoothIndicatorBaseF3D.h.

Constructor & Destructor Documentation

◆ SmoothIndicatorF3D()

template<typename T , typename S >

olb::SmoothIndicatorF3D< T, S, false >::SmoothIndicatorF3D ( )

protected

Definition at line 36 of file smoothIndicatorBaseF3D.hh.

  : AnalyticalF3D<T,S>(1),
    _myMin(S()), _myMax(S()), _pos(S()), _rotMat(S()), _circumRadius(S()),
    _theta(S()), _epsilon(S())
{ }

Member Function Documentation

◆ calcCenterOfMass()

template<typename T , typename S >

Vector< S, 3 > olb::SmoothIndicatorF3D< T, S, false >::calcCenterOfMass ( )

virtual

Definition at line 132 of file smoothIndicatorBaseF3D.hh.

{
  // TODO: Fallback
  assert(false);
  return Vector<S,3>(std::numeric_limits<BaseType<S>>::quiet_NaN());
}

◆ calcMofiAndMass()

template<typename T , typename S >

Vector< S, 4 > olb::SmoothIndicatorF3D< T, S, false >::calcMofiAndMass ( S density )

virtual

Definition at line 124 of file smoothIndicatorBaseF3D.hh.

{
  // TODO: Fallback
  assert(false);
  return Vector<S,4>(std::numeric_limits<BaseType<S>>::quiet_NaN());
}

◆ distance()

template<typename T , typename S >

bool olb::SmoothIndicatorF3D< T, S, false >::distance	(	S &	distance,
		const Vector< S, 3 > &	origin,
		const Vector< S, 3 > &	direction,
		S	precision,
		S	pitch )

virtual

Definition at line 178 of file smoothIndicatorBaseF3D.hh.

{
  S const halfEps = this->getEpsilon() * 0.5;
  bool originValue;
  bool currentValue;
 
  // start at origin and move into given direction
  PhysR<S,3> currentPoint(origin);
 
  originValue = this->signedDistance(origin) <= halfEps;
  currentValue = this->signedDistance(currentPoint) <= halfEps;
 
  while (currentValue == originValue && this->isInsideCircumRadius(currentPoint)) {
    currentPoint += direction;
    // update currentValue until the first point on the other side (inside/outside) is found
    currentValue = this->signedDistance(currentPoint) <= halfEps;
  }
 
  // return false if no point was found in given direction
  if (!this->isInsideCircumRadius(currentPoint) && !originValue) {
    return false;
  }
 
 
  while (pitch >= precision) {
    if (!this->isInsideCircumRadius(currentPoint) && originValue) {
      currentPoint -= pitch * direction;
      pitch /= 2.;
    }
    else {
      currentValue = this->signedDistance(currentPoint) <= halfEps;
      if (currentValue == originValue) {
        currentPoint += pitch * direction;
        pitch /= 2.;
      }
      else {
        currentPoint -= pitch * direction;
        pitch /= 2.;
      }
    }
  }
 
  distance = norm(currentPoint - origin);
  return true;
}

References olb::norm().

Here is the call graph for this function:

◆ getCircumRadius()

template<typename T , typename S >

const S & olb::SmoothIndicatorF3D< T, S, false >::getCircumRadius ( ) const

Definition at line 67 of file smoothIndicatorBaseF3D.hh.

{
  return _circumRadius;
}

◆ getEpsilon()

template<typename T , typename S >

const S & olb::SmoothIndicatorF3D< T, S, false >::getEpsilon ( ) const

Definition at line 85 of file smoothIndicatorBaseF3D.hh.

{
  return _epsilon;
}

◆ getMax()

template<typename T , typename S >

const Vector< S, 3 > & olb::SmoothIndicatorF3D< T, S, false >::getMax ( ) const

Definition at line 55 of file smoothIndicatorBaseF3D.hh.

{
  return _myMax;
}

Here is the caller graph for this function:

◆ getMin()

template<typename T , typename S >

const Vector< S, 3 > & olb::SmoothIndicatorF3D< T, S, false >::getMin ( ) const

Definition at line 49 of file smoothIndicatorBaseF3D.hh.

{
  return _myMin;
}

Here is the caller graph for this function:

◆ getPos()

template<typename T , typename S >

const Vector< S, 3 > & olb::SmoothIndicatorF3D< T, S, false >::getPos ( ) const

Definition at line 61 of file smoothIndicatorBaseF3D.hh.

{
  return _pos;
}

◆ getRotationMatrix()

template<typename T , typename S >

const Vector< S, 9 > & olb::SmoothIndicatorF3D< T, S, false >::getRotationMatrix ( ) const

Definition at line 73 of file smoothIndicatorBaseF3D.hh.

{
  return _rotMat;
}

◆ getTheta()

template<typename T , typename S >

const Vector< S, 3 > & olb::SmoothIndicatorF3D< T, S, false >::getTheta ( ) const

Definition at line 79 of file smoothIndicatorBaseF3D.hh.

{
  return _theta;
}

◆ getVolume()

template<typename T , typename S >

S olb::SmoothIndicatorF3D< T, S, false >::getVolume ( )

virtual

Definition at line 116 of file smoothIndicatorBaseF3D.hh.

{
  // TODO: Fallback
  assert(false);
  return S(std::numeric_limits<BaseType<S>>::quiet_NaN());
}

◆ init()

template<typename T , typename S >

void olb::SmoothIndicatorF3D< T, S, false >::init ( )

Definition at line 43 of file smoothIndicatorBaseF3D.hh.

{
  _rotMat = util::calculateInverseRotationMatrix<T,3>( this->_theta );
}

◆ isInsideCircumRadius()

template<typename T , typename S >

bool olb::SmoothIndicatorF3D< T, S, false >::isInsideCircumRadius ( const PhysR< S, 3 > & input )

Definition at line 147 of file smoothIndicatorBaseF3D.hh.

{
  return norm(input-this->getPos()) <= this->getCircumRadius();
}

References olb::norm().

Here is the call graph for this function:

◆ name()

template<typename T , typename S >

std::string olb::SmoothIndicatorF3D< T, S, false >::name ( )

Definition at line 91 of file smoothIndicatorBaseF3D.hh.

{
  return _name;
}

◆ operator()()

template<typename T , typename S >

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

virtual

has to be implemented for 'every' derived class

Implements olb::GenericF< T, S >.

Reimplemented in olb::SmoothIndicatorIdentity3D< T, S >.

Definition at line 140 of file smoothIndicatorBaseF3D.hh.

{
  T const signedDist = this->signedDistance(input);
  return sdf::evalSolidVolumeFraction(output, signedDist, this->getEpsilon());
}

References olb::sdf::evalSolidVolumeFraction().

Here is the call graph for this function:

◆ operator+()

template<typename T , typename S >

SmoothIndicatorF3D< T, S, false > & olb::SmoothIndicatorF3D< T, S, false >::operator+ ( SmoothIndicatorF3D< T, S, false > & rhs )

◆ setEpsilon()

template<typename T , typename S >

void olb::SmoothIndicatorF3D< T, S, false >::setEpsilon ( S epsilon )

Definition at line 110 of file smoothIndicatorBaseF3D.hh.

{
  _epsilon = epsilon;
}

◆ setPos()

template<typename T , typename S >

void olb::SmoothIndicatorF3D< T, S, false >::setPos ( Vector< S, 3 > pos )

Definition at line 97 of file smoothIndicatorBaseF3D.hh.

{
  _pos = pos;
}

◆ setTheta()

template<typename T , typename S >

void olb::SmoothIndicatorF3D< T, S, false >::setTheta ( Vector< S, 3 > theta )

Definition at line 103 of file smoothIndicatorBaseF3D.hh.

{
  _theta = theta;
  init();
}

◆ signedDistance()

template<typename T , typename S >

const S olb::SmoothIndicatorF3D< T, S, false >::signedDistance ( const PhysR< T, 3 > input )

virtual

Definition at line 170 of file smoothIndicatorBaseF3D.hh.

{
  // TODO: Raymarching as fallback
  assert(false);
  return std::numeric_limits<double>::quiet_NaN();
}

◆ surfaceNormal() [1/2]

template<typename T , typename S >

Vector< S, 3 > olb::SmoothIndicatorF3D< T, S, false >::surfaceNormal	(	const Vector< S, 3 > &	pos,
		const S	meshSize )

virtual

Definition at line 153 of file smoothIndicatorBaseF3D.hh.

{
  return surfaceNormal(pos, meshSize, [&](const Vector<S,3>& pos) {
    return pos;
  });
}

◆ surfaceNormal() [2/2]

template<typename T , typename S >

Vector< S, 3 > olb::SmoothIndicatorF3D< T, S, false >::surfaceNormal	(	const Vector< S, 3 > &	pos,
		const S	meshSize,
		std::function< Vector< S, 3 >(const Vector< S, 3 > &)>	transformPos )

virtual

Definition at line 161 of file smoothIndicatorBaseF3D.hh.

{
  return util::surfaceNormal(pos, meshSize, [&](const Vector<S,3>& pos) {
    return this->signedDistance( transformPos(pos) );
  });
}

References olb::util::surfaceNormal().

Here is the call graph for this function:

Member Data Documentation

◆ _circumRadius

template<typename T , typename S >

S olb::SmoothIndicatorF3D< T, S, false >::_circumRadius

protected

Definition at line 55 of file smoothIndicatorBaseF3D.h.

◆ _epsilon

template<typename T , typename S >

S olb::SmoothIndicatorF3D< T, S, false >::_epsilon

protected

Definition at line 58 of file smoothIndicatorBaseF3D.h.

◆ _myMax

template<typename T , typename S >

Vector<S,3> olb::SmoothIndicatorF3D< T, S, false >::_myMax

protected

Definition at line 52 of file smoothIndicatorBaseF3D.h.

◆ _myMin

template<typename T , typename S >

Vector<S,3> olb::SmoothIndicatorF3D< T, S, false >::_myMin

protected

Definition at line 51 of file smoothIndicatorBaseF3D.h.

◆ _name

template<typename T , typename S >

std::string olb::SmoothIndicatorF3D< T, S, false >::_name = "smoothIndicator3D"

protected

Definition at line 59 of file smoothIndicatorBaseF3D.h.

◆ _pos

template<typename T , typename S >

Vector<S,3> olb::SmoothIndicatorF3D< T, S, false >::_pos

protected

Definition at line 53 of file smoothIndicatorBaseF3D.h.

◆ _rotMat

template<typename T , typename S >

Vector<S,9> olb::SmoothIndicatorF3D< T, S, false >::_rotMat

protected

Definition at line 54 of file smoothIndicatorBaseF3D.h.

◆ _theta

template<typename T , typename S >

Vector<S,3> olb::SmoothIndicatorF3D< T, S, false >::_theta

protected

Definition at line 56 of file smoothIndicatorBaseF3D.h.

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

src/functors/analytical/indicator/smoothIndicatorBaseF3D.h
src/functors/analytical/indicator/smoothIndicatorBaseF3D.hh

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SmoothIndicatorF3D()

Member Function Documentation

◆ calcCenterOfMass()

◆ calcMofiAndMass()

◆ distance()

◆ getCircumRadius()

◆ getEpsilon()

◆ getMax()

◆ getMin()

◆ getPos()

◆ getRotationMatrix()

◆ getTheta()

◆ getVolume()

◆ init()

◆ isInsideCircumRadius()

◆ name()

◆ operator()()

◆ operator+()

◆ setEpsilon()

◆ setPos()

◆ setTheta()

◆ signedDistance()

◆ surfaceNormal() [1/2]

◆ surfaceNormal() [2/2]

Member Data Documentation

◆ _circumRadius

◆ _epsilon

◆ _myMax

◆ _myMin

◆ _name

◆ _pos

◆ _rotMat

◆ _theta