deriveEquilMSGas Subroutine

private recursive subroutine deriveEquilMSGas(fun, varSys, elempos, time, tree, nElems, nDofs, res)

Calculate the equlibrium of a given element number with the given input state vector.

The equilibrium distribution function is:\n $$f^{eq}_rest = \rho w_i( (3-2 mRatio_iField ) + 3(cx . u) + 9/2 (cx . u )^2 - 3/2 u^2 )$$ $$f^{eq}_i = \rho w_i( mRatio_iField + 3(cx . u) + 9/2 (cx . u )^2 - 3/2 u^2 )$$ 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.

The interface has to comply to the abstract interface tem_varSys_proc_element.

Arguments