Mercurial > repos > public > sbplib_julia

comparison test/SbpOperators/volumeops/derivatives/dissipation_test.jl @ 1291:356ec6a72974 refactor/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Implement changes in SbpOperators
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 07 Mar 2023 09:48:00 +0100
parents 7d52c4835d15
children 760a4a1ec4b7
comparison
equal deleted inserted replaced
1290:31d0b7638304 1291:356ec6a72974
24 return zero(x) 24 return zero(x)
25 end 25 end
26 x^k/factorial(k) 26 x^k/factorial(k)
27 end 27 end
28 28
29 @test_skip @testset "undivided_skewed04" begin 29 @testset "undivided_skewed04" begin
30 g = EquidistantGrid(20, 0., 11.) 30 g = equidistant_grid(20, 0., 11.)
31 D,Dᵀ = undivided_skewed04(g, 1) 31 D,Dᵀ = undivided_skewed04(g, 1)
32 32
33 @test D isa LazyTensor{Float64,1,1} 33 @test D isa LazyTensor{Float64,1,1}
34 @test Dᵀ isa LazyTensor{Float64,1,1} 34 @test Dᵀ isa LazyTensor{Float64,1,1}
35 35
36 @testset "Accuracy conditions" begin 36 @testset "Accuracy conditions" begin
37 N = 20 37 N = 20
38 g = EquidistantGrid(N, 0//1,2//1) 38 g = equidistant_grid(N, 0//1,2//1)
39 h = only(spacing(g)) 39 h = only(spacing(g))
40 @testset "D_$p" for p ∈ [1,2,3,4] 40 @testset "D_$p" for p ∈ [1,2,3,4]
41 D,Dᵀ = undivided_skewed04(g, p) 41 D,Dᵀ = undivided_skewed04(g, p)
42 42
43 @testset "x^$k" for k ∈ 0:p 43 @testset "x^$k" for k ∈ 0:p
44 v = evalOn(g, x->monomial(x,k)) 44 v = eval_on(g, x->monomial(x,k))
45 vₚₓ = evalOn(g, x->monomial(x,k-p)) 45 vₚₓ = eval_on(g, x->monomial(x,k-p))
46 46
47 @test D*v == h^p * vₚₓ 47 @test D*v == h^p * vₚₓ
48 end 48 end
49 end 49 end
50
51 # TODO: Add 2D tests
50 end 52 end
51 53
52 @testset "transpose equality" begin 54 @testset "transpose equality" begin
53 function get_matrix(D) 55 function get_matrix(D)
54 N = only(range_size(D)) 56 N = only(range_size(D))
65 end 67 end
66 68
67 return Dmat 69 return Dmat
68 end 70 end
69 71
70 g = EquidistantGrid(11, 0., 1.) 72 g = equidistant_grid(11, 0., 1.)
71 @testset "D_$p" for p ∈ [1,2,3,4] 73 @testset "D_$p" for p ∈ [1,2,3,4]
72 D,Dᵀ = undivided_skewed04(g, p) 74 D,Dᵀ = undivided_skewed04(g, p)
73 75
74 D̄ = get_matrix(D) 76 D̄ = get_matrix(D)
75 D̄ᵀ = get_matrix(Dᵀ) 77 D̄ᵀ = get_matrix(Dᵀ)