Use Gram-Schmidt for create_mrt_orthogonal
1 unresolved thread
1 unresolved thread
- orthogonalize automatically if no moments from literature are provided
- must specify whether weighted or unweighted orthogonality is desired
- add
lbmpy.moments.is_bulk_mode
and makelbmpy.moments.is_shear_mode
also detect the diagonal modes - use the bulk mode () in orthogonalization so we get the same physical modes as in literature
We 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)
Edited by Michael Kuron