Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Markus Holzer
lbmpy
Commits
f6b0f841
Commit
f6b0f841
authored
Aug 20, 2021
by
Frederik Hennig
Browse files
flake8
parent
c79890fd
Pipeline
#33872
failed with stages
in 37 minutes and 53 seconds
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
lbmpy/moments.py
View file @
f6b0f841
...
...
@@ -460,8 +460,8 @@ 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
[
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"
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
...
...
@@ -604,7 +604,7 @@ def central_moment_reduced_monomial_to_polynomial_matrix(polynomials, stencil, v
dim
=
len
(
stencil
[
0
])
if
velocity_symbols
is
None
:
velocity_symbols
=
sp
.
symbols
(
f
"u_:
{
dim
}
"
)
reduced_polynomials
=
non_aliased_polynomial_raw_moments
(
polynomials
,
stencil
)
reduced_monomials
=
sorted
(
extract_monomials
(
reduced_polynomials
),
key
=
exponent_tuple_sort_key
)
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment