This class stores statistical values for a given spatial position; to be used in combination with StatisticsVector. More...

#include <StatisticsPoint.h>

Public Member Functions
	StatisticsPoint ()
	Constructor sets sensible values. More...

	StatisticsPoint (const StatisticsPoint &other)
	Copy constructor; simply copies everything. More...

void	setCGShape (const char *CG_type)
	see StatisticsVector::setCGShape More...

CG_TYPE	getCGShape () const
	see StatisticsVector::getCGShape More...

void	setCGWidth2 (Mdouble new_)
	see StatisticsVector::setCGWidth2 More...

Mdouble	getCGWidthSquared () const
	see StatisticsVector::getCGWidthSquared More...

Mdouble	getCGWidth () const
	see StatisticsVector::getCGWidth More...

Mdouble	getCutoff ()
	see StatisticsVector::getCutoff More...

Mdouble	getCutoff2 ()
	see StatisticsVector::getCutoff2 More...

Mdouble	getXMaxStat ()
	see StatisticsVector::getXMaxStat More...

Mdouble	getYMaxStat ()
	see StatisticsVector::getYMaxStat More...

Mdouble	getZMaxStat ()
	see StatisticsVector::getZMaxStat More...

Mdouble	getXMinStat ()
	see StatisticsVector::getXMinStat More...

Mdouble	getYMinStat ()
	see StatisticsVector::getYMinStat More...

Mdouble	getZMinStat ()
	see StatisticsVector::getZMinStat More...

void	getN (int &nx_, int &ny_, int &nz_)
	see StatisticsVector::get_n More...

Mdouble	evaluatePolynomial (Mdouble r)
	see StatisticsVector::evaluatePolynomial More...

Mdouble	evaluatePolynomialGradient (Mdouble r)
	see StatisticsVector::evaluatePolynomialGradient More...

Mdouble	evaluateIntegral (Mdouble n1, Mdouble n2, Mdouble t)
	see StatisticsVector::evaluateIntegral More...

void	setCGInverseVolume ()
	sets CG_invvolume More...

int	nonaveragedDim ()

double	averagingVolume ()

Mdouble	getCGInverseVolume ()
	returns CG_invvolume More...

void	set_Gaussian_invvolume (int dim)
	sets CG_invvolume if CG_type=Gaussian More...

double	compute_Gaussian_invvolume (int dim)
	computes CG_invvolume if CG_type=Gaussian More...

void	set_Heaviside_invvolume ()
	sets CG_invvolume if CG_type=HeaviSideSphere More...

void	set_Polynomial_invvolume (int dim)
	sets CG_invvolume if CG_type=Polynomial More...

void	setPosition (Vec3D new_)
	sets Position More...

Vec3D	getPosition () const
	returns Position More...

void	set_zero ()
	Sets all statistical variables to zero. More...

StatisticsPoint< T >	getSquared ()
	Squares all statistical variables; needed for variance. More...

StatisticsPoint< T > &	operator= (const StatisticsPoint< T > &P)
	Defines a equal operator needed for copy constructor. More...

StatisticsPoint< T > &	operator+= (const StatisticsPoint< T > &P)
	Defines a plus operator needed to average values ( $ \bar{v} = (\sum_{i=1}^n v_i)/n $ ) More...

StatisticsPoint< T > &	operator-= (const StatisticsPoint< T > &P)
	Defines a plus operator needed to calculate variance. More...

StatisticsPoint< T > &	operator/= (const Mdouble a)
	Defines a division operator needed to average values ( $ \bar{v} = (\sum_{i=1}^n v_i)/n $ ) More...

void	firstTimeAverage (const int n)
	Defines a division operator needed to time-average values (because the displacement does not have a value at the first time step, this is slightly different than /=) More...

Mdouble	getDistanceSquaredNonAveraged (const Vec3D &P)
	returns the coarse graining distance in the coordinates that are not averaged about More...

Mdouble	dotNonAveraged (const Vec3D &P, const Vec3D &Q)
	Returns the dot product of two vectors in the coordinates that are not averaged about. More...

Vec3D	clearAveragedDirections (Vec3D P)
	Returns a vector where the averaged directions are zero. More...

Vec3D	crossNonAveraged (Vec3D P, Vec3D &Q)
	Returns the cross product of two vectors in the coordinates that are not averaged about. More...

Matrix3D	matrixCrossNonAveraged (Vec3D P, Matrix3D &Q)
	Returns the cross product of two vectors in the coordinates that are not averaged about. More...

Mdouble	CG_function (const Vec3D &PI)
	Returns the value of the course graining function phi(P,PI) More...

Mdouble	CG_function_2D (const Vec3D &PI)

Mdouble	CG_function_1D (const Vec3D &PI)

Vec3D	CG_gradient (const Vec3D &P, const Mdouble phi)

Vec3D	CG_integral_gradient (Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance)

Mdouble	CG_integral_gradient_1D (Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance)

Mdouble	CG_integral (Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance, Vec3D &rpsi)
	Returns the value of the coarse graining integral $ \psi(P,PI,PJ) = \int_0^1 \phi(P - PI + s PJ) ds $. More...

Mdouble	CG_integral_2D (Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance, Mdouble &rpsi_scalar)

Mdouble	CG_integral_1D (Vec3D &P1, Vec3D &P2, Vec3D &P1_P2_normal, Mdouble P1_P2_distance, Mdouble &rpsi_scalar)

std::string	print () const
	Outputs statistical variables in human-readable format. More...

std::string	print_sqrt () const
	Outputs root of statistical variables in human-readable format. More...

std::string	write_variable_names ()
	Outputs names of statistical variables in computer-readable format. More...

std::string	write () const
	Outputs statistical variables in computer-readable format. More...

Static Public Member Functions
static void	set_gb (StatisticsVector< T > *new_gb)
	see StatisticsVector::setCGShape More...

Public Attributes
Mdouble	Nu
	Particle volume fraction, $\sum_i m_i/\rho_i \phi_i$. More...

Mdouble	Density
	Density, $\sum_i m_i \phi_i$. More...

Vec3D	Momentum
	Momentum, $\sum_i m_i v_{ia} \phi_i$. More...

Vec3D	DisplacementMomentum
	Momentum from linear displacement, $\sum_i m_i u_{ia} \phi_i$, where $\vec{u}_i=(\vec{r}_i(t)-\vec{r}_i(t-\tau))/\tau$, with $\tau=savecount \Delta t$ the time intervall between outputs. More...

MatrixSymmetric3D	Displacement
	Linear displacement, $2/(2 \rho_i^2) \sum_{ij} m_i m_j \phi_j (v_{ija} \partial_b \phi_i + v_{ijb} \partial_a \phi_i) $. More...

MatrixSymmetric3D	MomentumFlux
	Momentum flux, $\sum_i m_i v_{ia} v_{ib} \phi_i$. More...

MatrixSymmetric3D	DisplacementMomentumFlux
	Momentum flux from linear displacement, $\sum_i m_i u_{ia} u_{ib} \phi_i$. More...

Vec3D	EnergyFlux
	Energy flux, $2/2 \sum_i m_i v_{ia} v_{ia} v_{ib} \phi_i$. More...

Matrix3D	NormalStress
	Stress from normal forces, $\sum_{ij} f^n_{ija} l_{ijb} \psi_{ij}$. More...

Matrix3D	TangentialStress
	Stress from tangential forces, $\sum_{ij} f^t_{ija} l_{ijb} \psi_{ij}$. More...

Vec3D	NormalTraction
	Traction from normal forces, $\sum_{ij} f^n_{ija} l_{ijb} \phi_i$. More...

Vec3D	TangentialTraction
	Traction from tangential forces, $\sum_{ij} f^t_{ija} l_{ijb} \phi_i$. More...

MatrixSymmetric3D	Fabric
	Fabric tensor, $\sum_{ij} m_i/\rho_i n_{ija} n_{ijb} \psi_{ij}$. More...

Vec3D	CollisionalHeatFlux

Mdouble	Dissipation

Mdouble	Potential

Vec3D	LocalAngularMomentum

Matrix3D	LocalAngularMomentumFlux

Matrix3D	ContactCoupleStress

Mdouble	CG_invvolume
	Prefactor of CG function which depends on $w. More...

int	mirrorParticle
	indicates that a position is a (fake) particles used for periodic walls More...

Private Attributes
Vec3D	Position
	Position at which evaluation occurs. More...

Static Private Attributes
static StatisticsVector< T > *	gb
	Pointer to StatisticsVector (to obtain global parameters) More...

Friends
std::ostream &	operator (std::ostream &os, const StatisticsPoint< T > &stat)
	Outputs statistical variables to ostream. More...

Detailed Description

template<StatType T>
struct StatisticsPoint< T >

This class stores statistical values for a given spatial position; to be used in combination with StatisticsVector.

Constructor & Destructor Documentation

◆ StatisticsPoint() [1/2]

template<StatType T>

StatisticsPoint< T >::StatisticsPoint ( )

inline

Constructor sets sensible values.

     {
         this->Position.setZero();
         mirrorParticle = -1;
         //~ this->set_zero(); 
     }

References StatisticsPoint< T >::mirrorParticle, StatisticsPoint< T >::Position, and Vec3D::setZero().

◆ StatisticsPoint() [2/2]

template<StatType T>

StatisticsPoint< T >::StatisticsPoint ( const StatisticsPoint< T > & other )

inline

Copy constructor; simply copies everything.

     {
         *this = other;
         mirrorParticle = -1;
     }

References StatisticsPoint< T >::mirrorParticle.

Member Function Documentation

◆ averagingVolume()

template<StatType T>

double StatisticsPoint< T >::averagingVolume ( )

◆ CG_function()

template<StatType T>

Mdouble StatisticsPoint< T >::CG_function ( const Vec3D & PI )

Returns the value of the course graining function phi(P,PI)

◆ CG_function_1D()

template<StatType T>

Mdouble StatisticsPoint< T >::CG_function_1D ( const Vec3D & PI )

Returns the value of the course graining function phi(P,PI) averaged along a plane

◆ CG_function_2D()

template<StatType T>

Mdouble StatisticsPoint< T >::CG_function_2D ( const Vec3D & PI )

returns the value of the course graining function phi(P,PI) averaged along a line

◆ CG_gradient()

template<StatType T>

Vec3D StatisticsPoint< T >::CG_gradient	(	const Vec3D &	P,
		const Mdouble	phi
	)

gradient of phi

Todo:

◆ CG_integral()

template<StatType T>

Mdouble StatisticsPoint< T >::CG_integral	(	Vec3D &	P1,
		Vec3D &	P2,
		Vec3D &	P1_P2_normal,
		Mdouble	P1_P2_distance,
		Vec3D &	rpsi
	)

Returns the value of the coarse graining integral $ \psi(P,PI,PJ) = \int_0^1 \phi(P - PI + s PJ) ds $.

◆ CG_integral_1D()

template<StatType T>

Mdouble StatisticsPoint< T >::CG_integral_1D	(	Vec3D &	P1,
		Vec3D &	P2,
		Vec3D &	P1_P2_normal,
		Mdouble	P1_P2_distance,
		Mdouble &	rpsi_scalar
	)

Returns the value of the coarse graining integral $ \psi(P,PI,PJ) = \int_0^1 \phi(P - PI + s PJ) ds $ averaged along a plane

◆ CG_integral_2D()

template<StatType T>

Mdouble StatisticsPoint< T >::CG_integral_2D	(	Vec3D &	P1,
		Vec3D &	P2,
		Vec3D &	P1_P2_normal,
		Mdouble	P1_P2_distance,
		Mdouble &	rpsi_scalar
	)

Returns the value of the coarse graining integral $ \psi(P,PI,PJ) = \int_0^1 \phi(P - PI + s PJ) ds $ averaged along a line

◆ CG_integral_gradient()

template<StatType T>

Vec3D StatisticsPoint< T >::CG_integral_gradient	(	Vec3D &	P1,
		Vec3D &	P2,
		Vec3D &	P1_P2_normal,
		Mdouble	P1_P2_distance
	)

gradient of phi

Todo:

◆ CG_integral_gradient_1D()

template<StatType T>

Mdouble StatisticsPoint< T >::CG_integral_gradient_1D	(	Vec3D &	P1,
		Vec3D &	P2,
		Vec3D &	P1_P2_normal,
		Mdouble	P1_P2_distance
	)

◆ clearAveragedDirections()

template<StatType T>

Vec3D StatisticsPoint< T >::clearAveragedDirections ( Vec3D P )

Returns a vector where the averaged directions are zero.

◆ compute_Gaussian_invvolume()

template<StatType T>

double StatisticsPoint< T >::compute_Gaussian_invvolume ( int dim )

computes CG_invvolume if CG_type=Gaussian

◆ crossNonAveraged()

template<StatType T>

Vec3D StatisticsPoint< T >::crossNonAveraged	(	Vec3D	P,
		Vec3D &	Q
	)

Returns the cross product of two vectors in the coordinates that are not averaged about.

◆ dotNonAveraged()

template<StatType T>

Mdouble StatisticsPoint< T >::dotNonAveraged	(	const Vec3D &	P,
		const Vec3D &	Q
	)

Returns the dot product of two vectors in the coordinates that are not averaged about.

◆ evaluateIntegral()

template<StatType T>

Mdouble StatisticsPoint< T >::evaluateIntegral	(	Mdouble	n1,
		Mdouble	n2,
		Mdouble	t
	)

inline

see StatisticsVector::evaluateIntegral

     {
         return gb->evaluateIntegral(n1, n2, t);
     }

References StatisticsPoint< T >::gb, and plotPSD::t.

◆ evaluatePolynomial()

template<StatType T>

Mdouble StatisticsPoint< T >::evaluatePolynomial ( Mdouble r )

inline

see StatisticsVector::evaluatePolynomial

     {
         return gb->evaluatePolynomial(r);
     }

References StatisticsPoint< T >::gb, and UniformPSDSelfTest::r.

◆ evaluatePolynomialGradient()

template<StatType T>

Mdouble StatisticsPoint< T >::evaluatePolynomialGradient ( Mdouble r )

inline

see StatisticsVector::evaluatePolynomialGradient

     {
         return gb->evaluatePolynomialGradient(r);
     }

References StatisticsPoint< T >::gb, and UniformPSDSelfTest::r.

◆ firstTimeAverage()

template<StatType T>

void StatisticsPoint< T >::firstTimeAverage ( const int n )

inline

Defines a division operator needed to time-average values (because the displacement does not have a value at the first time step, this is slightly different than /=)

◆ getCGInverseVolume()

template<StatType T>

Mdouble StatisticsPoint< T >::getCGInverseVolume ( )

inline

returns CG_invvolume

     {
         return this->CG_invvolume;
     }

References StatisticsPoint< T >::CG_invvolume.

◆ getCGShape()

template<StatType T>

CG_TYPE StatisticsPoint< T >::getCGShape ( ) const

inline

see StatisticsVector::getCGShape

     {
         return this->gb->getCGShape();
     }

References StatisticsPoint< T >::gb.

◆ getCGWidth()

template<StatType T>

Mdouble StatisticsPoint< T >::getCGWidth ( ) const

inline

see StatisticsVector::getCGWidth

     {
         return this->gb->getCGWidth();
     }

References StatisticsPoint< T >::gb.

◆ getCGWidthSquared()

template<StatType T>

Mdouble StatisticsPoint< T >::getCGWidthSquared ( ) const

inline

see StatisticsVector::getCGWidthSquared

     {
         return this->gb->getCGWidthSquared();
     }

References StatisticsPoint< T >::gb.

◆ getCutoff()

template<StatType T>

Mdouble StatisticsPoint< T >::getCutoff ( )

inline

see StatisticsVector::getCutoff

     {
         return this->gb->getCutoff();
     }

References StatisticsPoint< T >::gb.

◆ getCutoff2()

template<StatType T>

Mdouble StatisticsPoint< T >::getCutoff2 ( )

inline

see StatisticsVector::getCutoff2

     {
         return this->gb->getCutoff2();
     }

References StatisticsPoint< T >::gb.

◆ getDistanceSquaredNonAveraged()

template<StatType T>

Mdouble StatisticsPoint< T >::getDistanceSquaredNonAveraged ( const Vec3D & P )

returns the coarse graining distance in the coordinates that are not averaged about

◆ getN()

template<StatType T>

void StatisticsPoint< T >::getN	(	int &	nx_,
		int &	ny_,
		int &	nz_
	)

inline

see StatisticsVector::get_n

     {
         this->gb->getN(nx_, ny_, nz_);
     }

References StatisticsPoint< T >::gb.

◆ getPosition()

template<StatType T>

Vec3D StatisticsPoint< T >::getPosition ( ) const

inline

returns Position

     {
         return this->Position;
     }

References StatisticsPoint< T >::Position.

◆ getSquared()

template<StatType T>

StatisticsPoint<T> StatisticsPoint< T >::getSquared ( )

Squares all statistical variables; needed for variance.

◆ getXMaxStat()

template<StatType T>

Mdouble StatisticsPoint< T >::getXMaxStat ( )

inline

see StatisticsVector::getXMaxStat

     {
         return this->gb->getXMaxStat();
     }

References StatisticsPoint< T >::gb.

◆ getXMinStat()

template<StatType T>

Mdouble StatisticsPoint< T >::getXMinStat ( )

inline

see StatisticsVector::getXMinStat

     {
         return this->gb->getXMinStat();
     }

References StatisticsPoint< T >::gb.

◆ getYMaxStat()

template<StatType T>

Mdouble StatisticsPoint< T >::getYMaxStat ( )

inline

see StatisticsVector::getYMaxStat

     {
         return this->gb->getYMaxStat();
     }

References StatisticsPoint< T >::gb.

◆ getYMinStat()

template<StatType T>

Mdouble StatisticsPoint< T >::getYMinStat ( )

inline

see StatisticsVector::getYMinStat

     {
         return this->gb->getYMinStat();
     }

References StatisticsPoint< T >::gb.

◆ getZMaxStat()

template<StatType T>

Mdouble StatisticsPoint< T >::getZMaxStat ( )

inline

see StatisticsVector::getZMaxStat

     {
         return this->gb->getZMaxStat();
     }

References StatisticsPoint< T >::gb.

◆ getZMinStat()

template<StatType T>

Mdouble StatisticsPoint< T >::getZMinStat ( )

inline

see StatisticsVector::getZMinStat

     {
         return this->gb->getZMinStat();
     }

References StatisticsPoint< T >::gb.

◆ matrixCrossNonAveraged()

template<StatType T>

Matrix3D StatisticsPoint< T >::matrixCrossNonAveraged	(	Vec3D	P,
		Matrix3D &	Q
	)

Returns the cross product of two vectors in the coordinates that are not averaged about.

◆ nonaveragedDim()

template<StatType T>

int StatisticsPoint< T >::nonaveragedDim ( )

◆ operator+=()

template<StatType T>

StatisticsPoint<T>& StatisticsPoint< T >::operator+= ( const StatisticsPoint< T > & P )

inline

Defines a plus operator needed to average values ( $ \bar{v} = (\sum_{i=1}^n v_i)/n $ )

◆ operator-=()

template<StatType T>

StatisticsPoint<T>& StatisticsPoint< T >::operator-= ( const StatisticsPoint< T > & P )

inline

Defines a plus operator needed to calculate variance.

◆ operator/=()

template<StatType T>

StatisticsPoint<T>& StatisticsPoint< T >::operator/= ( const Mdouble a )

inline

Defines a division operator needed to average values ( $ \bar{v} = (\sum_{i=1}^n v_i)/n $ )

◆ operator=()

template<StatType T>

StatisticsPoint<T>& StatisticsPoint< T >::operator= ( const StatisticsPoint< T > & P )

inline

Defines a equal operator needed for copy constructor.

◆ print()

template<StatType T>

std::string StatisticsPoint< T >::print ( ) const

Outputs statistical variables in human-readable format.

◆ print_sqrt()

template<StatType T>

std::string StatisticsPoint< T >::print_sqrt ( ) const

Outputs root of statistical variables in human-readable format.

◆ set_Gaussian_invvolume()

template<StatType T>

void StatisticsPoint< T >::set_Gaussian_invvolume ( int dim )

sets CG_invvolume if CG_type=Gaussian

◆ set_gb()

template<StatType T>

static void StatisticsPoint< T >::set_gb ( StatisticsVector< T > * new_gb )

inlinestatic

see StatisticsVector::setCGShape

     {
         gb = new_gb;
     }

References StatisticsPoint< T >::gb.

◆ set_Heaviside_invvolume()

template<StatType T>

void StatisticsPoint< T >::set_Heaviside_invvolume ( )

sets CG_invvolume if CG_type=HeaviSideSphere

◆ set_Polynomial_invvolume()

template<StatType T>

void StatisticsPoint< T >::set_Polynomial_invvolume ( int dim )

sets CG_invvolume if CG_type=Polynomial

◆ set_zero()

template<StatType T>

void StatisticsPoint< T >::set_zero ( )

Sets all statistical variables to zero.

◆ setCGInverseVolume()

template<StatType T>

void StatisticsPoint< T >::setCGInverseVolume ( )

sets CG_invvolume

◆ setCGShape()

template<StatType T>

void StatisticsPoint< T >::setCGShape ( const char * CG_type )

inline

see StatisticsVector::setCGShape

     {
         this->gb->setCGShape(CG_type);
     }

References StatisticsPoint< T >::gb.

◆ setCGWidth2()

template<StatType T>

void StatisticsPoint< T >::setCGWidth2 ( Mdouble new_ )

inline

see StatisticsVector::setCGWidth2

     {
         this->gb->setCGWidth2(new_);
     }

References StatisticsPoint< T >::gb.

◆ setPosition()

template<StatType T>

void StatisticsPoint< T >::setPosition ( Vec3D new_ )

inline

sets Position

     {
         this->Position = new_;
     }

References StatisticsPoint< T >::Position.

◆ write()

template<StatType T>

std::string StatisticsPoint< T >::write ( ) const

Outputs statistical variables in computer-readable format.

◆ write_variable_names()

template<StatType T>

std::string StatisticsPoint< T >::write_variable_names ( )

Outputs names of statistical variables in computer-readable format.

Friends And Related Function Documentation

◆ operator

template<StatType T>

std::ostream& operator	(	std::ostream &	os,
		const StatisticsPoint< T > &	stat
	)

friend

Outputs statistical variables to ostream.

Member Data Documentation

◆ CG_invvolume

template<StatType T>

Mdouble StatisticsPoint< T >::CG_invvolume

Prefactor of CG function which depends on $w.

Referenced by StatisticsPoint< T >::getCGInverseVolume().

◆ CollisionalHeatFlux

template<StatType T>

Vec3D StatisticsPoint< T >::CollisionalHeatFlux

Heat flux from collisions, $\sum_{ij} (v_{ia}+v_ja)/2 f_{ija} l_{ijb} \psi_i$

Todo:

◆ ContactCoupleStress

template<StatType T>

Matrix3D StatisticsPoint< T >::ContactCoupleStress

◆ Density

template<StatType T>

Mdouble StatisticsPoint< T >::Density

Density, $\sum_i m_i \phi_i$.

◆ Displacement

template<StatType T>

MatrixSymmetric3D StatisticsPoint< T >::Displacement

Linear displacement, $2/(2 \rho_i^2) \sum_{ij} m_i m_j \phi_j (v_{ija} \partial_b \phi_i + v_{ijb} \partial_a \phi_i) $.

◆ DisplacementMomentum

template<StatType T>

Vec3D StatisticsPoint< T >::DisplacementMomentum

Momentum from linear displacement, $\sum_i m_i u_{ia} \phi_i$, where $\vec{u}_i=(\vec{r}_i(t)-\vec{r}_i(t-\tau))/\tau$, with $\tau=savecount \Delta t$ the time intervall between outputs.

◆ DisplacementMomentumFlux

template<StatType T>

MatrixSymmetric3D StatisticsPoint< T >::DisplacementMomentumFlux

Momentum flux from linear displacement, $\sum_i m_i u_{ia} u_{ib} \phi_i$.

◆ Dissipation

template<StatType T>

Mdouble StatisticsPoint< T >::Dissipation

Dissipation form collisions, $\sum_{ij} (v_{ia}+v_ja)/2 f_{ija} \psi_i$

Todo:

◆ EnergyFlux

template<StatType T>

Vec3D StatisticsPoint< T >::EnergyFlux

Energy flux, $2/2 \sum_i m_i v_{ia} v_{ia} v_{ib} \phi_i$.

◆ Fabric

template<StatType T>

MatrixSymmetric3D StatisticsPoint< T >::Fabric

Fabric tensor, $\sum_{ij} m_i/\rho_i n_{ija} n_{ijb} \psi_{ij}$.

◆ gb

template<StatType T>

StatisticsVector<T>* StatisticsPoint< T >::gb

staticprivate

Pointer to StatisticsVector (to obtain global parameters)

◆ LocalAngularMomentum

template<StatType T>

Vec3D StatisticsPoint< T >::LocalAngularMomentum

◆ LocalAngularMomentumFlux

template<StatType T>

Matrix3D StatisticsPoint< T >::LocalAngularMomentumFlux

◆ mirrorParticle

template<StatType T>

int StatisticsPoint< T >::mirrorParticle

indicates that a position is a (fake) particles used for periodic walls

Referenced by StatisticsPoint< T >::StatisticsPoint().

◆ Momentum

template<StatType T>

Vec3D StatisticsPoint< T >::Momentum

Momentum, $\sum_i m_i v_{ia} \phi_i$.

◆ MomentumFlux

template<StatType T>

MatrixSymmetric3D StatisticsPoint< T >::MomentumFlux

Momentum flux, $\sum_i m_i v_{ia} v_{ib} \phi_i$.

◆ NormalStress

template<StatType T>

Matrix3D StatisticsPoint< T >::NormalStress

Stress from normal forces, $\sum_{ij} f^n_{ija} l_{ijb} \psi_{ij}$.

◆ NormalTraction

template<StatType T>

Vec3D StatisticsPoint< T >::NormalTraction

Traction from normal forces, $\sum_{ij} f^n_{ija} l_{ijb} \phi_i$.

◆ Nu

template<StatType T>

Mdouble StatisticsPoint< T >::Nu

Particle volume fraction, $\sum_i m_i/\rho_i \phi_i$.

◆ Position

template<StatType T>

Vec3D StatisticsPoint< T >::Position

private

Position at which evaluation occurs.

Referenced by StatisticsPoint< T >::getPosition(), StatisticsPoint< T >::setPosition(), and StatisticsPoint< T >::StatisticsPoint().

◆ Potential

template<StatType T>

Mdouble StatisticsPoint< T >::Potential

Elastic energy $2/4 \sum_{ij} (k \delta_{ij}^2 +k^t {\delta^t_{ij}}^2) \phi_i$

Todo:: {Potential calculations have not been checked, only implemented; fstat file does not include the torsional and rolling spring and hence cannot account for their potential}

◆ TangentialStress

template<StatType T>

Matrix3D StatisticsPoint< T >::TangentialStress

Stress from tangential forces, $\sum_{ij} f^t_{ija} l_{ijb} \psi_{ij}$.

◆ TangentialTraction

template<StatType T>

Vec3D StatisticsPoint< T >::TangentialTraction

Traction from tangential forces, $\sum_{ij} f^t_{ija} l_{ijb} \phi_i$.

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

StatisticsPoint.h

Public Attributes
Mdouble	Nu
	Particle volume fraction, \(\sum_i m_i/\rho_i \phi_i\). More...

Mdouble	Density
	Density, \(\sum_i m_i \phi_i\). More...

Vec3D	Momentum
	Momentum, \(\sum_i m_i v_{ia} \phi_i\). More...

Vec3D	DisplacementMomentum
	Momentum from linear displacement, \(\sum_i m_i u_{ia} \phi_i\), where \(\vec{u}_i=(\vec{r}_i(t)-\vec{r}_i(t-\tau))/\tau\), with \(\tau=savecount \Delta t\) the time intervall between outputs. More...

MatrixSymmetric3D	Displacement
	Linear displacement, \(2/(2 \rho_i^2) \sum_{ij} m_i m_j \phi_j (v_{ija} \partial_b \phi_i + v_{ijb} \partial_a \phi_i) \). More...

MatrixSymmetric3D	MomentumFlux
	Momentum flux, \(\sum_i m_i v_{ia} v_{ib} \phi_i\). More...

MatrixSymmetric3D	DisplacementMomentumFlux
	Momentum flux from linear displacement, \(\sum_i m_i u_{ia} u_{ib} \phi_i\). More...

Vec3D	EnergyFlux
	Energy flux, \(2/2 \sum_i m_i v_{ia} v_{ia} v_{ib} \phi_i\). More...

Matrix3D	NormalStress
	Stress from normal forces, \(\sum_{ij} f^n_{ija} l_{ijb} \psi_{ij}\). More...

Matrix3D	TangentialStress
	Stress from tangential forces, \(\sum_{ij} f^t_{ija} l_{ijb} \psi_{ij}\). More...

Vec3D	NormalTraction
	Traction from normal forces, \(\sum_{ij} f^n_{ija} l_{ijb} \phi_i\). More...

Vec3D	TangentialTraction
	Traction from tangential forces, \(\sum_{ij} f^t_{ija} l_{ijb} \phi_i\). More...

MatrixSymmetric3D	Fabric
	Fabric tensor, \(\sum_{ij} m_i/\rho_i n_{ija} n_{ijb} \psi_{ij}\). More...

Vec3D	CollisionalHeatFlux

Mdouble	Dissipation

Mdouble	Potential

Vec3D	LocalAngularMomentum

Matrix3D	LocalAngularMomentumFlux

Matrix3D	ContactCoupleStress

Mdouble	CG_invvolume
	Prefactor of CG function which depends on $w. More...

int	mirrorParticle
	indicates that a position is a (fake) particles used for periodic walls More...

Public Member Functions

Static Public Member Functions

Public Attributes

Private Attributes

Static Private Attributes

Friends

Detailed Description

template<StatType T> struct StatisticsPoint< T >

Constructor & Destructor Documentation

◆ StatisticsPoint() [1/2]

◆ StatisticsPoint() [2/2]

Member Function Documentation

◆ averagingVolume()

◆ CG_function()

◆ CG_function_1D()

◆ CG_function_2D()

◆ CG_gradient()

◆ CG_integral()

◆ CG_integral_1D()

◆ CG_integral_2D()

◆ CG_integral_gradient()

◆ CG_integral_gradient_1D()

◆ clearAveragedDirections()

◆ compute_Gaussian_invvolume()

◆ crossNonAveraged()

◆ dotNonAveraged()

◆ evaluateIntegral()

◆ evaluatePolynomial()

◆ evaluatePolynomialGradient()

◆ firstTimeAverage()

◆ getCGInverseVolume()

◆ getCGShape()

◆ getCGWidth()

◆ getCGWidthSquared()

◆ getCutoff()

◆ getCutoff2()

◆ getDistanceSquaredNonAveraged()

◆ getN()

◆ getPosition()

◆ getSquared()

◆ getXMaxStat()

◆ getXMinStat()

◆ getYMaxStat()

◆ getYMinStat()

◆ getZMaxStat()

◆ getZMinStat()

◆ matrixCrossNonAveraged()

◆ nonaveragedDim()

◆ operator+=()

◆ operator-=()

◆ operator/=()

◆ operator=()

◆ print()

◆ print_sqrt()

◆ set_Gaussian_invvolume()

◆ set_gb()

◆ set_Heaviside_invvolume()

◆ set_Polynomial_invvolume()

◆ set_zero()

◆ setCGInverseVolume()

◆ setCGShape()

◆ setCGWidth2()

◆ setPosition()

◆ write()

◆ write_variable_names()

Friends And Related Function Documentation

◆ operator

Member Data Documentation

◆ CG_invvolume

◆ CollisionalHeatFlux

◆ ContactCoupleStress

◆ Density

◆ Displacement

◆ DisplacementMomentum

◆ DisplacementMomentumFlux

◆ Dissipation

◆ EnergyFlux

◆ Fabric

◆ gb

◆ LocalAngularMomentum

template<StatType T>
struct StatisticsPoint< T >