Mercurial > repos > public > sbplib_julia

changeset 1029:129262c8e897 feature/dissipation_operators

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change signatures for interior stencil methods
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 22 Mar 2022 09:26:05 +0100
parents 62e80e8cc743
children 8fe630a20b64
files src/SbpOperators/volumeops/derivatives/dissipation.jl test/SbpOperators/volumeops/derivatives/dissipation_test.jl
diffstat 2 files changed, 11 insertions(+), 19 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/volumeops/derivatives/dissipation.jl	Tue Mar 22 08:53:45 2022 +0100
+++ b/src/SbpOperators/volumeops/derivatives/dissipation.jl	Tue Mar 22 09:26:05 2022 +0100
@@ -15,20 +15,12 @@
    return (0, dissipation_interior_weights(p-1)...) .- (dissipation_interior_weights(p-1)..., 0)
 end
 
-function dissipation_interior_stencil(p)
-    w = dissipation_interior_weights(p)
-    Stencil(w..., center=midpoint(w))
-end
-
-function dissipation_transpose_interior_stencil(p)
-    w = dissipation_interior_weights(p)
-    Stencil(w..., center=midpoint_transpose(w))
-end
-
-
 midpoint(weights) = length(weights)รท2 + 1
 midpoint_transpose(weights) = length(weights)+1 - midpoint(weights)
 
+dissipation_interior_stencil(weights) =           Stencil(weights..., center=midpoint(weights))
+dissipation_transpose_interior_stencil(weights) = Stencil(weights..., center=midpoint_transpose(weights))
+
 dissipation_lower_closure_size(weights) = midpoint(weights) - 1
 dissipation_upper_closure_size(weights) = length(weights) - midpoint(weights)
 
--- a/test/SbpOperators/volumeops/derivatives/dissipation_test.jl	Tue Mar 22 08:53:45 2022 +0100
+++ b/test/SbpOperators/volumeops/derivatives/dissipation_test.jl	Tue Mar 22 09:26:05 2022 +0100
@@ -59,17 +59,17 @@
 end
 
 @testset "dissipation_interior_stencil" begin
-    @test dissipation_interior_stencil(1) == Stencil(-1,1, center=2)
-    @test dissipation_interior_stencil(2) == Stencil(1,-2,1, center=2)
-    @test dissipation_interior_stencil(3) == Stencil(-1,3,-3,1, center=3)
-    @test dissipation_interior_stencil(4) == Stencil(1, -4, 6, -4, 1, center=3)
+    @test dissipation_interior_stencil(dissipation_interior_weights(1)) == Stencil(-1,1, center=2)
+    @test dissipation_interior_stencil(dissipation_interior_weights(2)) == Stencil(1,-2,1, center=2)
+    @test dissipation_interior_stencil(dissipation_interior_weights(3)) == Stencil(-1,3,-3,1, center=3)
+    @test dissipation_interior_stencil(dissipation_interior_weights(4)) == Stencil(1, -4, 6, -4, 1, center=3)
 end
 
 @testset "dissipation_transpose_interior_stencil" begin
-    @test dissipation_transpose_interior_stencil(1) == Stencil(-1,1, center=1)
-    @test dissipation_transpose_interior_stencil(2) == Stencil(1,-2,1, center=2)
-    @test dissipation_transpose_interior_stencil(3) == Stencil(-1,3,-3,1, center=2)
-    @test dissipation_transpose_interior_stencil(4) == Stencil(1, -4, 6, -4, 1, center=3)
+    @test dissipation_transpose_interior_stencil(dissipation_interior_weights(1)) == Stencil(-1,1, center=1)
+    @test dissipation_transpose_interior_stencil(dissipation_interior_weights(2)) == Stencil(1,-2,1, center=2)
+    @test dissipation_transpose_interior_stencil(dissipation_interior_weights(3)) == Stencil(-1,3,-3,1, center=2)
+    @test dissipation_transpose_interior_stencil(dissipation_interior_weights(4)) == Stencil(1, -4, 6, -4, 1, center=3)
 end
 
 @testset "midpoint" begin