Mercurial > repos > public > sbplib_julia

comparison test/testSbpOperators.jl @ 646:27b4bcb22154 feature/volume_and_boundary_operators

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove outcommented obsolete code.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 04 Jan 2021 18:34:38 +0100
parents c5ccab5cf164
children 351937390162
comparison
equal deleted inserted replaced
645:c5ccab5cf164 646:27b4bcb22154
111 @test get_tuple(parsed_toml, "order2", "d1", "closure") == (-3/2, 2, -1/2) 111 @test get_tuple(parsed_toml, "order2", "d1", "closure") == (-3/2, 2, -1/2)
112 112
113 @test_throws AssertionError get_tuple(parsed_toml, "meta", "type") 113 @test_throws AssertionError get_tuple(parsed_toml, "meta", "type")
114 end 114 end
115 end 115 end
116
117 # @testset "apply_quadrature" begin
118 # op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4)
119 # h = 0.5
120 #
121 # @test apply_quadrature(op, h, 1.0, 10, 100) == h
122 #
123 # N = 10
124 # qc = op.quadratureClosure
125 # q = h.*(qc..., ones(N-2*closuresize(op))..., reverse(qc)...)
126 # @assert length(q) == N
127 #
128 # for i ∈ 1:N
129 # @test apply_quadrature(op, h, 1.0, i, N) == q[i]
130 # end
131 #
132 # v = [2.,3.,2.,4.,5.,4.,3.,4.,5.,4.5]
133 # for i ∈ 1:N
134 # @test apply_quadrature(op, h, v[i], i, N) == q[i]*v[i]
135 # end
136 # end
137 116
138 @testset "VolumeOperator" begin 117 @testset "VolumeOperator" begin
139 inner_stencil = Stencil(1/4 .* (1.,2.,1.),center=2) 118 inner_stencil = Stencil(1/4 .* (1.,2.,1.),center=2)
140 closure_stencils = (Stencil(1/2 .* (1.,1.),center=1),Stencil((0.,1.),center=2)) 119 closure_stencils = (Stencil(1/2 .* (1.,1.),center=1),Stencil((0.,1.),center=2))
141 g_1D = EquidistantGrid(11,0.,1.) 120 g_1D = EquidistantGrid(11,0.,1.)