Mercurial > repos > public > sbplib_julia

comparison test/testSbpOperators.jl @ 698:5ddf28ddee18 refactor/operator_naming

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Test inverse_inner_product on 0-dimensional grid
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Sun, 14 Feb 2021 13:52:13 +0100
parents 1b3b8f82349e
children 3cd582257072
comparison
equal deleted inserted replaced
697:1b3b8f82349e 698:5ddf28ddee18
490 Ly = Float64(π) 490 Ly = Float64(π)
491 g_1D = EquidistantGrid(77, 0.0, Lx) 491 g_1D = EquidistantGrid(77, 0.0, Lx)
492 g_2D = EquidistantGrid((77,66), (0.0, 0.0), (Lx,Ly)) 492 g_2D = EquidistantGrid((77,66), (0.0, 0.0), (Lx,Ly))
493 @testset "inverse_inner_product" begin 493 @testset "inverse_inner_product" begin
494 op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4) 494 op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4)
495 @testset "0D" begin
496 Hi = inverse_inner_product(EquidistantGrid{Float64}(),op.quadratureClosure)
497 @test Hi == IdentityMapping{Float64}()
498 @test Hi isa TensorMapping{T,0,0} where T
499 end
495 @testset "1D" begin 500 @testset "1D" begin
496 Hi = inverse_inner_product(g_1D, op.quadratureClosure); 501 Hi = inverse_inner_product(g_1D, op.quadratureClosure);
497 inner_stencil = CenteredStencil(1.) 502 inner_stencil = CenteredStencil(1.)
498 closures = () 503 closures = ()
499 for i = 1:length(op.quadratureClosure) 504 for i = 1:length(op.quadratureClosure)