getEquilibriumIncomp Function

public pure function getEquilibriumIncomp(dens, vel, layout, rho0) result(equil)

Calculate the equilibrium distribution function in all directions This is the incompressible formulation with reference density rho0

The equilibrium distribution function is:\n $$f^{eq}_i = w_i ( \rho + \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 = \sum_i f_i$ is the macroscopic density,\n $c_s$ is the speed of sound,\n $\vec c_i$ is the lattice unit velocity in each direction,\n $\vec u = \sum_i c_i f_i$ 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
real(kind=rk),	intent(in)			::	rho0

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