Mercurial > repos > public > sbplib_julia

changeset 928:453fd1a2e858 feature/variable_derivatives

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge bugfix/normal_derivative_sign
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 18 Feb 2022 08:03:37 +0100
parents eb054537fc63 (current diff) 35be8253de89 (diff)
children e9dd43cbd127
files src/SbpOperators/operators/standard_diagonal.toml
diffstat 2 files changed, 11 insertions(+), 11 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/operators/standard_diagonal.toml	Fri Feb 18 07:39:07 2022 +0100
+++ b/src/SbpOperators/operators/standard_diagonal.toml	Fri Feb 18 08:03:37 2022 +0100
@@ -21,7 +21,7 @@
 H.closure = ["1/2"]
 
 e.closure = ["1"]
-d1.closure = {s = ["-3/2", "2", "-1/2"], c = 1}
+d1.closure = {s = ["3/2", "-2", "1/2"], c = 1}
 
 
 D1.inner_stencil = ["-1/2", "0", "1/2"]
@@ -47,7 +47,7 @@
 H.closure = ["17/48", "59/48", "43/48", "49/48"]
 
 e.closure = ["1"]
-d1.closure = {s = ["-11/6", "3", "-3/2", "1/3"], c = 1}
+d1.closure = {s = ["11/6", "-3", "3/2", "-1/3"], c = 1}
 
 D1.inner_stencil = ["1/12","-2/3","0","2/3","-1/12"]
 D1.closure_stencils = [
--- a/test/SbpOperators/boundaryops/normal_derivative_test.jl	Fri Feb 18 07:39:07 2022 +0100
+++ b/test/SbpOperators/boundaryops/normal_derivative_test.jl	Fri Feb 18 08:03:37 2022 +0100
@@ -45,24 +45,24 @@
             d_s = normal_derivative(g_2D, d_closure, CartesianBoundary{2,Lower}())
             d_n = normal_derivative(g_2D, d_closure, CartesianBoundary{2,Upper}())
 
-            @test d_w*v ≈ v∂x[1,:] atol = 1e-13
-            @test d_e*v ≈ -v∂x[end,:] atol = 1e-13
-            @test d_s*v ≈ v∂y[:,1] atol = 1e-13
-            @test d_n*v ≈ -v∂y[:,end] atol = 1e-13
+            @test d_w*v ≈ -v∂x[1,:] atol = 1e-13
+            @test d_e*v ≈ v∂x[end,:] atol = 1e-13
+            @test d_s*v ≈ -v∂y[:,1] atol = 1e-13
+            @test d_n*v ≈ v∂y[:,end] atol = 1e-13
         end
 
         @testset "4th order" begin
-            stencil_set = read_stencil_set(sbp_operators_path()*"standard_diagonal.toml"; order=2)
+            stencil_set = read_stencil_set(sbp_operators_path()*"standard_diagonal.toml"; order=4)
         	d_closure = parse_stencil(stencil_set["d1"]["closure"])
             d_w = normal_derivative(g_2D, d_closure, CartesianBoundary{1,Lower}())
             d_e = normal_derivative(g_2D, d_closure, CartesianBoundary{1,Upper}())
             d_s = normal_derivative(g_2D, d_closure, CartesianBoundary{2,Lower}())
             d_n = normal_derivative(g_2D, d_closure, CartesianBoundary{2,Upper}())
 
-            @test d_w*v ≈ v∂x[1,:] atol = 1e-13
-            @test d_e*v ≈ -v∂x[end,:] atol = 1e-13
-            @test d_s*v ≈ v∂y[:,1] atol = 1e-13
-            @test d_n*v ≈ -v∂y[:,end] atol = 1e-13
+            @test d_w*v ≈ -v∂x[1,:] atol = 1e-13
+            @test d_e*v ≈ v∂x[end,:] atol = 1e-13
+            @test d_s*v ≈ -v∂y[:,1] atol = 1e-13
+            @test d_n*v ≈ v∂y[:,end] atol = 1e-13
         end
     end
 end