From 3befbd5d92e020d62fe160df559ae5e00d463827 Mon Sep 17 00:00:00 2001 From: Isuru Fernando Date: Sun, 1 Dec 2019 00:02:48 -0600 Subject: [PATCH] Add Biharmonic 3D --- sumpy/kernel.py | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/sumpy/kernel.py b/sumpy/kernel.py index b0e245e7..6383d8d5 100644 --- a/sumpy/kernel.py +++ b/sumpy/kernel.py @@ -455,11 +455,12 @@ class BiharmonicKernel(ExpressionKernel): expr = r**2 * var("log")(r) scaling = 1/(8*var("pi")) elif dim == 3: + # Ref: Jiang, Shidong, Bo Ren, Paul Tsuji, and Lexing Ying. + # "Second kind integral equations for the first kind Dirichlet problem + # of the biharmonic equation in three dimensions." + # Journal of Computational Physics 230, no. 19 (2011): 7488-7501. expr = r - scaling = 1 # FIXME: Unknown - from warnings import warn - warn("scaling factor for Biharmonic 3D is unknown.", - stacklevel=2) + scaling = -1/(8*var("pi")) else: raise RuntimeError("unsupported dimensionality") -- GitLab