Mercurial > repos > public > sbplib_julia

comparison test/testDiffOps.jl @ 594:cc86b920531a refactor/toml_operator_format

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change the readoperator function to use the .toml format
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 02 Dec 2020 15:26:13 +0100
parents ffddaf053085
children
comparison
equal deleted inserted replaced
593:fa03dae0ff0b 594:cc86b920531a
6 using Sbplib.LazyTensors 6 using Sbplib.LazyTensors
7 7
8 @testset "DiffOps" begin 8 @testset "DiffOps" begin
9 # 9 #
10 # @testset "BoundaryValue" begin 10 # @testset "BoundaryValue" begin
11 # op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt") 11 # op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4)
12 # g = EquidistantGrid((4,5), (0.0, 0.0), (1.0,1.0)) 12 # g = EquidistantGrid((4,5), (0.0, 0.0), (1.0,1.0))
13 # 13 #
14 # e_w = BoundaryValue(op, g, CartesianBoundary{1,Lower}()) 14 # e_w = BoundaryValue(op, g, CartesianBoundary{1,Lower}())
15 # e_e = BoundaryValue(op, g, CartesianBoundary{1,Upper}()) 15 # e_e = BoundaryValue(op, g, CartesianBoundary{1,Upper}())
16 # e_s = BoundaryValue(op, g, CartesianBoundary{2,Lower}()) 16 # e_s = BoundaryValue(op, g, CartesianBoundary{2,Lower}())
67 # @test_broken collect(e_s*g_x) == G_s 67 # @test_broken collect(e_s*g_x) == G_s
68 # @test_broken collect(e_n*g_x) == G_n 68 # @test_broken collect(e_n*g_x) == G_n
69 # end 69 # end
70 # 70 #
71 # @testset "NormalDerivative" begin 71 # @testset "NormalDerivative" begin
72 # op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt") 72 # op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4)
73 # g = EquidistantGrid((5,6), (0.0, 0.0), (4.0,5.0)) 73 # g = EquidistantGrid((5,6), (0.0, 0.0), (4.0,5.0))
74 # 74 #
75 # d_w = NormalDerivative(op, g, CartesianBoundary{1,Lower}()) 75 # d_w = NormalDerivative(op, g, CartesianBoundary{1,Lower}())
76 # d_e = NormalDerivative(op, g, CartesianBoundary{1,Upper}()) 76 # d_e = NormalDerivative(op, g, CartesianBoundary{1,Upper}())
77 # d_s = NormalDerivative(op, g, CartesianBoundary{2,Lower}()) 77 # d_s = NormalDerivative(op, g, CartesianBoundary{2,Lower}())
144 # @test_broken collect(d_s*g_x) ≈ G_s 144 # @test_broken collect(d_s*g_x) ≈ G_s
145 # @test_broken collect(d_n*g_x) ≈ G_n 145 # @test_broken collect(d_n*g_x) ≈ G_n
146 # end 146 # end
147 # 147 #
148 # @testset "BoundaryQuadrature" begin 148 # @testset "BoundaryQuadrature" begin
149 # op = readOperator(sbp_operators_path()*"d2_4th.txt",sbp_operators_path()*"h_4th.txt") 149 # op = read_D2_operator(sbp_operators_path()*"standard_diagonal.toml"; order=4)
150 # g = EquidistantGrid((10,11), (0.0, 0.0), (1.0,1.0)) 150 # g = EquidistantGrid((10,11), (0.0, 0.0), (1.0,1.0))
151 # 151 #
152 # H_w = BoundaryQuadrature(op, g, CartesianBoundary{1,Lower}()) 152 # H_w = BoundaryQuadrature(op, g, CartesianBoundary{1,Lower}())
153 # H_e = BoundaryQuadrature(op, g, CartesianBoundary{1,Upper}()) 153 # H_e = BoundaryQuadrature(op, g, CartesianBoundary{1,Upper}())
154 # H_s = BoundaryQuadrature(op, g, CartesianBoundary{2,Lower}()) 154 # H_s = BoundaryQuadrature(op, g, CartesianBoundary{2,Lower}())