From 2fc42b2cd26d4e7dc165c6d63bce313d28bd82cf Mon Sep 17 00:00:00 2001 From: Natalie Beams Date: Tue, 7 Mar 2017 13:21:29 -0600 Subject: [PATCH 1/2] change sign on stresslet to be consistent between 2d/3d --- sumpy/kernel.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/sumpy/kernel.py b/sumpy/kernel.py index e29ac516..646df2a5 100644 --- a/sumpy/kernel.py +++ b/sumpy/kernel.py @@ -518,7 +518,7 @@ class StressletKernel(ExpressionKernel): expr = ( d[icomp]*d[jcomp]*d[kcomp]/r**5 ) - scaling = -3/(4*var("pi")) + scaling = 3/(4*var("pi")) elif dim is None: expr = None -- GitLab From 8d9a18f9bd0bc3ac6848272ecd7cf954395d3c7a Mon Sep 17 00:00:00 2001 From: Natalie Beams Date: Tue, 7 Mar 2017 13:48:41 -0600 Subject: [PATCH 2/2] flip sign on 2d Stokeslet to use same convention in 2d/3d --- sumpy/kernel.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/sumpy/kernel.py b/sumpy/kernel.py index 646df2a5..83f23d2f 100644 --- a/sumpy/kernel.py +++ b/sumpy/kernel.py @@ -443,7 +443,7 @@ class StokesletKernel(ExpressionKernel): + d[icomp]*d[jcomp]/r**2 ) - scaling = 1/(4*var("pi")*mu) + scaling = -1/(4*var("pi")*mu) elif dim == 3: d = make_sym_vector("d", dim) -- GitLab