Mercurial > repos > public > sbplib_julia

comparison test/SbpOperators/boundaryops/boundary_restriction_test.jl @ 827:beae513ef8b3 operator_storage_array_of_table

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in boundary_restriction when grids and stencil does not have matching type
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 12 Jan 2022 07:45:38 +0100
parents bea2feebbeca
children 47425442bbc5
comparison
equal deleted inserted replaced
826:4433be383840 827:beae513ef8b3
15 15
16 @testset "boundary_restriction" begin 16 @testset "boundary_restriction" begin
17 @testset "1D" begin 17 @testset "1D" begin
18 e_l = boundary_restriction(g_1D,e_closure,Lower()) 18 e_l = boundary_restriction(g_1D,e_closure,Lower())
19 @test e_l == boundary_restriction(g_1D,e_closure,CartesianBoundary{1,Lower}()) 19 @test e_l == boundary_restriction(g_1D,e_closure,CartesianBoundary{1,Lower}())
20 @test e_l == BoundaryOperator(g_1D,e_closure,Lower()) 20 @test e_l == BoundaryOperator(g_1D,Stencil{Float64}(e_closure),Lower())
21 @test e_l isa BoundaryOperator{T,Lower} where T 21 @test e_l isa BoundaryOperator{T,Lower} where T
22 @test e_l isa TensorMapping{T,0,1} where T 22 @test e_l isa TensorMapping{T,0,1} where T
23 23
24 e_r = boundary_restriction(g_1D,e_closure,Upper()) 24 e_r = boundary_restriction(g_1D,e_closure,Upper())
25 @test e_r == boundary_restriction(g_1D,e_closure,CartesianBoundary{1,Upper}()) 25 @test e_r == boundary_restriction(g_1D,e_closure,CartesianBoundary{1,Upper}())
26 @test e_r == BoundaryOperator(g_1D,e_closure,Upper()) 26 @test e_r == BoundaryOperator(g_1D,Stencil{Float64}(e_closure),Upper())
27 @test e_r isa BoundaryOperator{T,Upper} where T 27 @test e_r isa BoundaryOperator{T,Upper} where T
28 @test e_r isa TensorMapping{T,0,1} where T 28 @test e_r isa TensorMapping{T,0,1} where T
29 end 29 end
30 30
31 @testset "2D" begin 31 @testset "2D" begin