Skip to content
Snippets Groups Projects
Commit 8b59e483 authored by Andreas Klöckner's avatar Andreas Klöckner
Browse files

Merge branch 'biharmonic' into 'master'

Fix biharmonic kernel expressions

See merge request inducer/sumpy!117
parents 343c2309 7cf3197b
No related branches found
No related tags found
1 merge request!117Fix biharmonic kernel expressions
Pipeline #19185 passed with warnings
......@@ -431,13 +431,14 @@ class BiharmonicKernel(ExpressionKernel):
init_arg_names = ("dim",)
def __init__(self, dim=None):
# See https://arxiv.org/abs/1202.1811
r = pymbolic_real_norm_2(make_sym_vector("d", dim))
if dim == 2:
expr = r**2 * var("log")(r)
scaling = 1/(8*var("pi"))
expr = r**2 * (var("log")(r) - 1)
scaling = -1/(8*var("pi"))
elif dim == 3:
expr = r
scaling = 1 # FIXME: Unknown
scaling = 1/(8*var("pi"))
else:
raise RuntimeError("unsupported dimensionality")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment