getEquilibrium_forElemfromState Function

private pure function getEquilibrium_forElemfromState(state, elem, layout, varPos, nScalars, neigh) result(equil)

Calculate the equilibrium distribution function in all directions

The equilibrim distribution function is:\n $$f^{eq}_i = w_i \rho ( 1 + \frac{\vec c_i \cdot \vec u}{c^2_s} + \frac{ {(\vec c_i \cdot \vec u)}^2}{2c^4_s} - \frac{\vec u \cdot \vec u}{2c^2_s})$$ \n where $w_i$ is the weight in each direction,\n $\rho$ is the macroscopic value of density,\n $c_s$ is the speed of sound,\n $\vec c_i$ is the lattice unit velocity in each direction,\n $\vec u$ is the macroscopic value of velocity.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in)			::	state(:)
integer,	intent(in)			::	elem
type(mus_scheme_layout_type),	intent(in)			::	layout
integer,	intent(in)			::	varPos(:)
integer,	intent(in)			::	nScalars
integer,	intent(in)			::	neigh(:)

Return Value real(kind=rk), (layout%fStencil%QQ)