sbplib_julia: test/testSbpOperators.jl comparison

comparison test/testSbpOperators.jl @ 383:aaf8e331cb80 refactor/sbp_operators_tests/collect_and_compare

Remove collects around TensorMappingApplications

author	Jonatan Werpers <jonatan@werpers.com>
date	Thu, 01 Oct 2020 06:11:40 +0200
parents	5c10cd0ed1fe
children	3a779c31e59e

comparison

equal deleted inserted replaced

-:5c10cd0ed1fe
+:aaf8e331cb80
 accuracytol = 0.5*1e-3
 # 4th order interior stencil, 2nd order boundary stencil,
 # implies that L*v should be exact for v - monomial up to order 3.
 # Exact differentiation is measured point-wise. For other grid functions
 # the error is measured in the l2-norm.
-@test all(abs.(collect(Dₓₓ*v0)) .<= equalitytol)
+@test all(abs.(Dₓₓ*v0) .<= equalitytol)
-@test all(abs.(collect(Dₓₓ*v1)) .<= equalitytol)
+@test all(abs.(Dₓₓ*v1) .<= equalitytol)
-@test all(abs.((collect(Dₓₓ*v2) - v0)) .<= equalitytol)
+@test all(abs.((Dₓₓ*v2 - v0)) .<= equalitytol)
-@test all(abs.((collect(Dₓₓ*v3) - v1)) .<= equalitytol)
+@test all(abs.((Dₓₓ*v3 - v1)) .<= equalitytol)
-e4 = collect(Dₓₓ*v4) - v2
+e4 = Dₓₓ*v4 - v2
-e5 = collect(Dₓₓ*v5) + v5
+e5 = Dₓₓ*v5 + v5
-@test sqrt(h*sum(collect(e4.^2))) <= accuracytol
+@test sqrt(h*sum(e4.^2)) <= accuracytol
-@test sqrt(h*sum(collect(e5.^2))) <= accuracytol
+@test sqrt(h*sum(e5.^2)) <= accuracytol
 end
 @testset "Laplace2D" begin
 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
 Lx = 1.5
 accuracytol = 0.5*1e-3
 # 4th order interior stencil, 2nd order boundary stencil,
 # implies that L*v should be exact for v - monomial up to order 3.
 # Exact differentiation is measured point-wise. For other grid functions
 # the error is measured in the H-norm.
-@test all(abs.(collect(L*v0)) .<= equalitytol)
+@test all(abs.(L*v0) .<= equalitytol)
-@test all(abs.(collect(L*v1)) .<= equalitytol)
+@test all(abs.(L*v1) .<= equalitytol)
-@test all(collect(L*v2) .≈ v0) # Seems to be more accurate
+@test all(L*v2 .≈ v0) # Seems to be more accurate
-@test all(abs.((collect(L*v3) - v1)) .<= equalitytol)
+@test all(abs.((L*v3 - v1)) .<= equalitytol)
-e4 = collect(L*v4) - v2
+e4 = L*v4 - v2
-e5 = collect(L*v5) - v5ₓₓ
+e5 = L*v5 - v5ₓₓ
 h = spacing(g)
-@test sqrt(prod(h)*sum(collect(e4.^2))) <= accuracytol
+@test sqrt(prod(h)*sum(e4.^2)) <= accuracytol
-@test sqrt(prod(h)*sum(collect(e5.^2))) <= accuracytol
+@test sqrt(prod(h)*sum(e5.^2)) <= accuracytol
 end
 @testset "DiagonalInnerProduct" begin
 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
 L = 2.3
 H = DiagonalInnerProduct(g,op.quadratureClosure)
 v = ones(Float64, size(g))
 @test H isa TensorMapping{T,1,1} where T
 @test H' isa TensorMapping{T,1,1} where T
-@test sum(collect(H*v)) ≈ L
+@test sum(H*v) ≈ L
-@test collect(H*v) == collect(H'*v)
+@test H*v == H'*v
 end
 @testset "Quadrature" begin
 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
 Lx = 2.3
 Hi = InverseDiagonalInnerProduct(g,op.quadratureClosure)
 v = evalOn(g, x->sin(x))
 @test Hi isa TensorMapping{T,1,1} where T
 @test Hi' isa TensorMapping{T,1,1} where T
-@test collect(Hi*H*v)  ≈ v
+@test Hi*H*v  ≈ v
-@test collect(Hi*v) == collect(Hi'*v)
+@test Hi*v == Hi'*v
 end
 @testset "InverseQuadrature" begin
 op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
 Lx = 7.3
 Qinv = InverseQuadrature(g, op.quadratureClosure)
 v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y)
 @test Qinv isa TensorMapping{T,2,2} where T
 @test Qinv' isa TensorMapping{T,2,2} where T
-@test_broken collect(Qinv*(Q*v)) ≈ v
+@test_broken Qinv*(Q*v) ≈ v
-@test collect(Qinv*v) == collect(Qinv'*v)
+@test Qinv*v == Qinv'*v
 end
 #
 # @testset "BoundaryValue" begin
 #     op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
 #     g = EquidistantGrid((4,5), (0.0, 0.0), (1.0,1.0))
 #     @test size(e_w'*v) == (5,)
 #     @test size(e_e'*v) == (5,)
 #     @test size(e_s'*v) == (4,)
 #     @test size(e_n'*v) == (4,)
 #
-#     @test collect(e_w'*v) == [10,7,4,1.0,1]
+#     @test e_w'*v == [10,7,4,1.0,1]
-#     @test collect(e_e'*v) == [13,10,7,4,4.0]
+#     @test e_e'*v == [13,10,7,4,4.0]
-#     @test collect(e_s'*v) == [10,11,12,13.0]
+#     @test e_s'*v == [10,11,12,13.0]
-#     @test collect(e_n'*v) == [1,2,3,4.0]
+#     @test e_n'*v == [1,2,3,4.0]
 #
 #     g_x = [1,2,3,4.0]
 #     g_y = [5,4,3,2,1.0]
 #
 #     G_w = zeros(Float64, (4,5))
 #     @test size(e_e*g_y) == (UnknownDim,5)
 #     @test size(e_s*g_x) == (4,UnknownDim)
 #     @test size(e_n*g_x) == (4,UnknownDim)
 #
 #     # These tests should be moved to where they are possible (i.e we know what the grid should be)
-#     @test_broken collect(e_w*g_y) == G_w
+#     @test_broken e_w*g_y == G_w
-#     @test_broken collect(e_e*g_y) == G_e
+#     @test_broken e_e*g_y == G_e
-#     @test_broken collect(e_s*g_x) == G_s
+#     @test_broken e_s*g_x == G_s
-#     @test_broken collect(e_n*g_x) == G_n
+#     @test_broken e_n*g_x == G_n
 # end
 #
 # @testset "NormalDerivative" begin
 #     op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
 #     g = EquidistantGrid((5,6), (0.0, 0.0), (4.0,5.0))
 #     @test size(d_w'*v) == (6,)
 #     @test size(d_e'*v) == (6,)
 #     @test size(d_s'*v) == (5,)
 #     @test size(d_n'*v) == (5,)
 #
-#     @test collect(d_w'*v) ≈ v∂x[1,:]
+#     @test d_w'*v ≈ v∂x[1,:]
-#     @test collect(d_e'*v) ≈ v∂x[5,:]
+#     @test d_e'*v ≈ v∂x[5,:]
-#     @test collect(d_s'*v) ≈ v∂y[:,1]
+#     @test d_s'*v ≈ v∂y[:,1]
-#     @test collect(d_n'*v) ≈ v∂y[:,6]
+#     @test d_n'*v ≈ v∂y[:,6]
 #
 #
 #     d_x_l = zeros(Float64, 5)
 #     d_x_u = zeros(Float64, 5)
 #     for i ∈ eachindex(d_x_l)
 #     @test size(d_e*g_y) == (UnknownDim,6)
 #     @test size(d_s*g_x) == (5,UnknownDim)
 #     @test size(d_n*g_x) == (5,UnknownDim)
 #
 #     # These tests should be moved to where they are possible (i.e we know what the grid should be)
-#     @test_broken collect(d_w*g_y) ≈ G_w
+#     @test_broken d_w*g_y ≈ G_w
-#     @test_broken collect(d_e*g_y) ≈ G_e
+#     @test_broken d_e*g_y ≈ G_e
-#     @test_broken collect(d_s*g_x) ≈ G_s
+#     @test_broken d_s*g_x ≈ G_s
-#     @test_broken collect(d_n*g_x) ≈ G_n
+#     @test_broken d_n*g_x ≈ G_n
 # end
 #
 # @testset "BoundaryQuadrature" begin
 #     op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt")
 #     g = EquidistantGrid((10,11), (0.0, 0.0), (1.0,1.0))
 #     @test size(H_w*v_w) == (11,)
 #     @test size(H_e*v_e) == (11,)
 #     @test size(H_s*v_s) == (10,)
 #     @test size(H_n*v_n) == (10,)
 #
-#     @test collect(H_w*v_w) ≈ q_y.*v_w
+#     @test H_w*v_w ≈ q_y.*v_w
-#     @test collect(H_e*v_e) ≈ q_y.*v_e
+#     @test H_e*v_e ≈ q_y.*v_e
-#     @test collect(H_s*v_s) ≈ q_x.*v_s
+#     @test H_s*v_s ≈ q_x.*v_s
-#     @test collect(H_n*v_n) ≈ q_x.*v_n
+#     @test H_n*v_n ≈ q_x.*v_n
 #
-#     @test collect(H_w'*v_w) == collect(H_w'*v_w)
+#     @test H_w'*v_w == H_w'*v_w
-#     @test collect(H_e'*v_e) == collect(H_e'*v_e)
+#     @test H_e'*v_e == H_e'*v_e
-#     @test collect(H_s'*v_s) == collect(H_s'*v_s)
+#     @test H_s'*v_s == H_s'*v_s
-#     @test collect(H_n'*v_n) == collect(H_n'*v_n)
+#     @test H_n'*v_n == H_n'*v_n
 # end
 end

Mercurial > repos > public > sbplib_julia

comparison test/testSbpOperators.jl @ 383:aaf8e331cb80 refactor/sbp_operators_tests/collect_and_compare