Mercurial > repos > public > sbplib_julia

changeset 1697:a63278c25c40 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix test of laplace on mapped grid
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 28 Aug 2024 13:28:13 +0200
parents 29b96fc75bee
children 3eb4e584425c
files test/SbpOperators/volumeops/laplace/laplace_test.jl
diffstat 1 files changed, 3 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
diff -r 29b96fc75bee -r a63278c25c40 test/SbpOperators/volumeops/laplace/laplace_test.jl
--- a/test/SbpOperators/volumeops/laplace/laplace_test.jl	Wed Aug 28 10:50:15 2024 +0200
+++ b/test/SbpOperators/volumeops/laplace/laplace_test.jl	Wed Aug 28 13:28:13 2024 +0200
@@ -95,18 +95,18 @@
         end
         Grids.jacobian(c::typeof(c), (ξ,η)) = @SMatrix[2 1-2η; (2+η)*ξ 3+ξ^2/2]
 
-        g = equidistant_grid(c, 30,30)
+        g = equidistant_grid(c, 60,60)
 
         @test laplace(g, stencil_set) isa LazyTensor{<:Any,2,2}
 
         f((x,y)) = sin(4(x + y))
-        Δf((x,y)) = -16sin(4(x + y))
+        Δf((x,y)) = -32sin(4(x + y))
         gf = map(f,g)
 
         Δ = laplace(g, stencil_set)
 
         @test collect(Δ*gf) isa Array{<:Any,2}
-        @test Δ*gf ≈ map(Δf, g)
+        @test Δ*gf ≈ map(Δf, g) rtol=2e-2
     end
 end