getEqByDensVel Function

public pure function getEqByDensVel(dens, vel, layout) result(equil)

Calculate the equilibrium distribution function in all directions

The equilibrium 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})$$

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)			::	dens
real(kind=rk),	intent(in)			::	vel(3)
type(mus_scheme_layout_type),	intent(in)			::	layout

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