Commit c79890fd authored by Frederik Hennig's avatar Frederik Hennig
Browse files

Fixed sympy 1.2 incompatibility (I hope)

parent d5ad2a9e
Pipeline #33871 failed with stages
in 13 minutes and 14 seconds
......@@ -460,7 +460,7 @@ def set_up_shift_matrix(moments, stencil, velocity_symbols=sp.symbols("u_:3")):
N = sp.simplify(MN * M.inv())
assert N.is_lower, "Calculating the shift matrix gave not a lower diagonal matrix. Thus it failed"
assert sum(N.diagonal()) == len(stencil), "Calculating the shift matrix failed. " \
assert sum(N[i,i] for i in range(len(stencil))) == len(stencil), "Calculating the shift matrix failed. " \
"There are entries on the diagonal which are not equal to one"
return N
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment