lbmpy.moments.is_bulk_mode
and make lbmpy.moments.is_shear_mode
also detect the diagonal modesx^2+y^2+z^2-1
) in orthogonalization so we get the same physical modes as in literatureWe now match the literature MRTs in all the physical moments. The only exception is D3Q27, which chooses a different linear combination for the shear moments. Differences in the higher-order moments remain, but they have no effect on the model.
Fixes #5 (closed)