diff --git a/boundaries/boundaryconditions.py b/boundaries/boundaryconditions.py index e7146e3fb953159ac586fec7dd4f8ec088416077..c49e613927d13204189246619d012790e80aa6c3 100644 --- a/boundaries/boundaryconditions.py +++ b/boundaries/boundaryconditions.py @@ -232,9 +232,3 @@ class NeumannByCopy(Boundary): def __eq__(self, other): return type(other) == NeumannByCopy - -if __name__ == '__main__': - from lbmpy.scenarios import * - sc = createForceDrivenChannel(domainSize=[60, 30]) - sc.boundaryHandling.setBoundary(FixedDensity(1.01), makeSlice[0.2:0.25, :]) - sc.boundaryHandling.prepare() diff --git a/creationfunctions.py b/creationfunctions.py index 280d39d42d1ac7c979a737f7a0532d4fb4882d15..aacbe746dc765c86738bc10173d00720c3103973 100644 --- a/creationfunctions.py +++ b/creationfunctions.py @@ -385,7 +385,11 @@ def createLatticeBoltzmannUpdateRule(lbMethod=None, optimizationParams={}, **kwa else: collisionRule = addEntropyCondition(collisionRule, omegaOutputField=params['omegaOutputField']) - fieldDtype = 'float64' if optimizationParams['doublePrecision'] else 'float32' + if 'doublePrecision' in optimizationParams: + fieldDtype = 'float64' if optimizationParams['doublePrecision'] else 'float32' + else: + fieldDtype = 'float64' + if optParams['fieldSize']: fieldSize = [s + 2 for s in optParams['fieldSize']] + [len(stencil)] srcField = Field.createFixedSize(params['fieldName'], fieldSize, indexDimensions=1, diff --git a/methods/momentbasedsimplifications.py b/methods/momentbasedsimplifications.py index 07558fc1436ded4070bf9d7072b3ec998db153be..fcf3ee14f957719f7aa42488eacd90232a96c1cf 100644 --- a/methods/momentbasedsimplifications.py +++ b/methods/momentbasedsimplifications.py @@ -35,7 +35,7 @@ def factorRelaxationRates(lbmCollisionEqs): """ sh = lbmCollisionEqs.simplificationHints assert 'relaxationRates' in sh, "Needs simplification hint 'relaxationRates': Sequence of relaxation rates" - if len(sh['relaxationRates']) > 12: # heuristics - for too many relaxation rates this simplification makes no sense + if len(sh['relaxationRates']) > 19: # heuristics - for too many relaxation rates this simplification makes no sense return lbmCollisionEqs relaxationRates = sp.Matrix(sh['relaxationRates']).atoms(sp.Symbol)