comparison test/Grids/curvilinear_grid_test.jl @ 1431:6adf31ba6cfd feature/grids/curvilinear

Add `jacobian` and `logicalgrid`
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 25 Aug 2023 15:41:19 +0200
parents 9fc3c1af33e5
children 64b60b42d367
comparison
equal deleted inserted replaced
1430:9fc3c1af33e5 1431:6adf31ba6cfd
7 x̄ = map(ξ̄ -> 2ξ̄, g) 7 x̄ = map(ξ̄ -> 2ξ̄, g)
8 J = map(ξ̄ -> @SArray(fill(2., 2, 2)), g) 8 J = map(ξ̄ -> @SArray(fill(2., 2, 2)), g)
9 9
10 @test CurvilinearGrid(g, x̄, J) isa Grid{SVector{2, Float64},2} 10 @test CurvilinearGrid(g, x̄, J) isa Grid{SVector{2, Float64},2}
11 11
12 12 cg = CurvilinearGrid(g, x̄, J)
13 @test_broken jacobian(cg) isa Array{<:AbstractVector} 13 @test jacobian(cg) isa Array{<:AbstractMatrix}
14 @test_broken logicalgrid(cg) isa Grid 14 @test logicalgrid(cg) isa Grid
15
16 15
17 16
18 @testset "Indexing Interface" begin 17 @testset "Indexing Interface" begin
19 # cg = CurvilinearGrid(g, x̄, J) 18 # cg = CurvilinearGrid(g, x̄, J)
20 # @test cg[1,1] == [0.0, 0.0] 19 # @test cg[1,1] == [0.0, 0.0]