Mercurial > repos > public > sbplib_julia
comparison test/testSbpOperators.jl @ 645:c5ccab5cf164 feature/volume_and_boundary_operators
Remove out-commented testset BoundaryQuadrature since it is obsolete.
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Mon, 04 Jan 2021 18:28:52 +0100 |
| parents | e3fd4b7aa789 |
| children | 27b4bcb22154 |
comparison
equal
deleted
inserted
replaced
| 644:e3fd4b7aa789 | 645:c5ccab5cf164 |
|---|---|
| 901 # @test_broken d_e'*g_y ≈ G_e | 901 # @test_broken d_e'*g_y ≈ G_e |
| 902 # @test d_s'*g_x ≈ G_s | 902 # @test d_s'*g_x ≈ G_s |
| 903 # @test_broken d_n'*g_x ≈ G_n | 903 # @test_broken d_n'*g_x ≈ G_n |
| 904 end | 904 end |
| 905 end | 905 end |
| 906 # | 906 |
| 907 # @testset "BoundaryQuadrature" begin | 907 end |
| 908 # op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4) | |
| 909 # g = EquidistantGrid((10,11), (0.0, 0.0), (1.0,1.0)) | |
| 910 # | |
| 911 # H_w = BoundaryQuadrature(op, g, CartesianBoundary{1,Lower}()) | |
| 912 # H_e = BoundaryQuadrature(op, g, CartesianBoundary{1,Upper}()) | |
| 913 # H_s = BoundaryQuadrature(op, g, CartesianBoundary{2,Lower}()) | |
| 914 # H_n = BoundaryQuadrature(op, g, CartesianBoundary{2,Upper}()) | |
| 915 # | |
| 916 # v = evalOn(g, (x,y)-> x^2 + (y-1)^2 + x*y) | |
| 917 # | |
| 918 # function get_quadrature(N) | |
| 919 # qc = op.quadratureClosure | |
| 920 # q = (qc..., ones(N-2*closuresize(op))..., reverse(qc)...) | |
| 921 # @assert length(q) == N | |
| 922 # return q | |
| 923 # end | |
| 924 # | |
| 925 # v_w = v[1,:] | |
| 926 # v_e = v[10,:] | |
| 927 # v_s = v[:,1] | |
| 928 # v_n = v[:,11] | |
| 929 # | |
| 930 # q_x = spacing(g)[1].*get_quadrature(10) | |
| 931 # q_y = spacing(g)[2].*get_quadrature(11) | |
| 932 # | |
| 933 # @test H_w isa TensorOperator{T,1} where T | |
| 934 # | |
| 935 # @test domain_size(H_w, (3,)) == (3,) | |
| 936 # @test domain_size(H_n, (3,)) == (3,) | |
| 937 # | |
| 938 # @test range_size(H_w, (3,)) == (3,) | |
| 939 # @test range_size(H_n, (3,)) == (3,) | |
| 940 # | |
| 941 # @test size(H_w*v_w) == (11,) | |
| 942 # @test size(H_e*v_e) == (11,) | |
| 943 # @test size(H_s*v_s) == (10,) | |
| 944 # @test size(H_n*v_n) == (10,) | |
| 945 # | |
| 946 # @test H_w*v_w .≈ q_y.*v_w | |
| 947 # @test H_e*v_e .≈ q_y.*v_e | |
| 948 # @test H_s*v_s .≈ q_x.*v_s | |
| 949 # @test H_n*v_n .≈ q_x.*v_n | |
| 950 # | |
| 951 # @test H_w'*v_w == H_w'*v_w | |
| 952 # @test H_e'*v_e == H_e'*v_e | |
| 953 # @test H_s'*v_s == H_s'*v_s | |
| 954 # @test H_n'*v_n == H_n'*v_n | |
| 955 # end | |
| 956 | |
| 957 end |