Mercurial > repos > public > sbplib_julia

comparison test/testSbpOperators.jl @ 687:5c229c889614 feature/boundary_quads

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in test.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 09 Feb 2021 14:40:26 +0100
parents 43cf58c69f91
children 728fd5a2455a
comparison
equal deleted inserted replaced
686:27dcac8fb350 687:5c229c889614
402 g_3D = EquidistantGrid((10,10, 10), (0.0, 0.0, 0.0), (Lx,Ly,Lz)) 402 g_3D = EquidistantGrid((10,10, 10), (0.0, 0.0, 0.0), (Lx,Ly,Lz))
403 integral(H,v) = sum(H*v) 403 integral(H,v) = sum(H*v)
404 @testset "quadrature" begin 404 @testset "quadrature" begin
405 op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4) 405 op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4)
406 @testset "0D" begin 406 @testset "0D" begin
407 H = quadrature(EquidistantGrid((),(),()),op.quadratureClosure) 407 H = quadrature(EquidistantGrid{Float64}(),op.quadratureClosure)
408 @test H == IdentityMapping{Float64}() 408 @test H == IdentityMapping{Float64}()
409 @test H isa TensorMapping{T,0,0} where T 409 @test H isa TensorMapping{T,0,0} where T
410 end 410 end
411 @testset "1D" begin 411 @testset "1D" begin
412 H = quadrature(g_1D,op.quadratureClosure) 412 H = quadrature(g_1D,op.quadratureClosure)