Mercurial > repos > public > sbplib_julia
diff test/SbpOperators/boundaryops/normal_derivative_test.jl @ 927:d360fc2d9620 feature/laplace_opset
Merge with default
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Mon, 21 Feb 2022 23:36:41 +0100 |
| parents | 47425442bbc5 35be8253de89 |
| children | af670581b464 |
line wrap: on
line diff
--- a/test/SbpOperators/boundaryops/normal_derivative_test.jl Mon Feb 21 23:33:29 2022 +0100 +++ b/test/SbpOperators/boundaryops/normal_derivative_test.jl Mon Feb 21 23:36:41 2022 +0100 @@ -41,10 +41,10 @@ (d_w, d_e, d_s, d_n) = map(id -> normal_derivative(g_2D, d_closure, id), boundary_identifiers(g_2D)) - @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 @@ -53,10 +53,10 @@ (d_w, d_e, d_s, d_n) = map(id -> normal_derivative(g_2D, d_closure, id), boundary_identifiers(g_2D)) - @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