Bugfixes

methods/momentbased.py

@@ -104,7 +104,7 @@ class MomentBasedLbMethod(AbstractLbMethod):
                 raise NotImplementedError("Shear moments seem to be not relaxed separately - "
                                           "Can not determine their relaxation rate automatically")
 
-    def getEquilibrium(self, conservedQuantityEquations=None, includeForceTerms=True):
+    def getEquilibrium(self, conservedQuantityEquations=None, includeForceTerms=False):
         D = sp.eye(len(self.relaxationRates))
         return self._getCollisionRuleWithRelaxationMatrix(D, conservedQuantityEquations=conservedQuantityEquations,
                                                           includeForceTerms=includeForceTerms)

moments.py

@@ -184,8 +184,8 @@ def polynomialToExponentRepresentation(polynomial, dim=3):
 
     Example:
         >>> x , y, z = MOMENT_SYMBOLS
-        >>> polynomialToExponentRepresentation(1 + (42 * x**2 * y**2 * z) )
-        [(1, (0, 0, 0)), (42, (2, 2, 1))]
+        >>> set(polynomialToExponentRepresentation(1 + (42 * x**2 * y**2 * z) )) == {(42, (2, 2, 1)), (1, (0, 0, 0))}
+        True
     """
     assert dim <= 3
     x, y, z = MOMENT_SYMBOLS