sbplib_julia: test/Grids/mapped_grid

Add min_spacing for 2D mapped grids

author	Jonatan Werpers <jonatan@werpers.com>
date	Thu, 22 Aug 2024 08:14:04 +0200
parents	22a9992471be
children	e11b5b6940a2

comparison

equal deleted inserted replaced

-:3ac94e8f28b3
+:72776d3d5fd6
 using Sbplib.Grids
 using Sbplib.RegionIndices
 using Test
 using StaticArrays
+using LinearAlgebra
 @testset "MappedGrid" begin
 lg = equidistant_grid((0,0), (1,1), 11, 11) # TODO: Change dims of the grid to be different
 x̄ = map(ξ̄ -> 2ξ̄, lg)
 J = map(ξ̄ -> @SArray(fill(2., 2, 2)), lg)
 end
 let g = mapped_grid(x->x + x.*(1 .- x)/2, x->@SMatrix[1.5 .- x], 11)
 @test min_spacing(g) ≈ 0.055
 end
+let g = mapped_grid(identity, x->@SMatrix[1 0; 0 1], 11,11)
+@test min_spacing(g) ≈ 0.1
+end
+let g = mapped_grid(identity, x->@SMatrix[1 0; 0 1], 11,21)
+@test min_spacing(g) ≈ 0.05
+end
+skew_grid(a,b, sz...) = mapped_grid(ξ̄->ξ̄[1]*a + ξ̄[2]*b, ξ̄->[a  b], sz...)
+@testset let a = @SVector[1,0], b = @SVector[1,1]/√2
+g = skew_grid(a,b,11,11)
+@test min_spacing(g) ≈ 0.1*norm(b-a)
+end
+@testset let a = @SVector[1,0], b = @SVector[-1,1]/√2
+g = skew_grid(a,b,11,11)
+@test min_spacing(g) ≈ 0.1*norm(a+b)
+end
+# Skevt nät
 end
 end
 @testset "mapped_grid" begin

Mercurial > repos > public > sbplib_julia