Mercurial > repos > public > sbplib_julia

changeset 2020:7f7207b9bd6c feature/grids/multiblock_grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add constructor tests for 2d grids and using matrix for the grid collection
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 12 Mar 2025 10:48:58 +0100
parents fb2dc185c197
children d6618d628515
files test/Grids/multiblockgrids_test.jl
diffstat 1 files changed, 32 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/test/Grids/multiblockgrids_test.jl	Wed Mar 12 10:40:28 2025 +0100
+++ b/test/Grids/multiblockgrids_test.jl	Wed Mar 12 10:48:58 2025 +0100
@@ -1,15 +1,15 @@
 using Diffinitive.Grids
+using StaticArrays
 
 @testset "MultiBlockGrid" begin
     @test MultiBlockGrid <: Grid
 
     @testset "Constructors" begin
+        # Vector
         g₁ = equidistant_grid(0,1,5)
         g₂ = equidistant_grid(0,1,5)
         g₃ = equidistant_grid(0,1,5)
 
-
-        # Vector
         C = connection.([
             (1, UpperBoundary(), 2, LowerBoundary()),
             (2, UpperBoundary(), 3, LowerBoundary()),
@@ -18,6 +18,36 @@
         @test MultiBlockGrid([g₁,g₂,g₃], C) isa Grid{Float64,1}
         @test MultiBlockGrid([g₁,g₂,g₃], C) isa MultiBlockGrid{Float64,1}
 
+
+        g₁₁ = equidistant_grid((0,0),(1,1),5,5)
+        g₁₂ = equidistant_grid((0,0),(1,1),5,5)
+        g₂₁ = equidistant_grid((0,0),(1,1),5,5)
+        g₂₂ = equidistant_grid((0,0),(1,1),5,5)
+
+        C = map(((i1,d1,b1,i2,d2,b2),)->connection(i1,CartesianBoundary{d1,b1}(), i2, CartesianBoundary{d2,b2}()),[
+            (1, 1, UpperBoundary, 2, 1, LowerBoundary),
+            (3, 1, UpperBoundary, 4, 1, LowerBoundary),
+            (1, 2, UpperBoundary, 3, 2, LowerBoundary),
+            (2, 2, UpperBoundary, 4, 2, LowerBoundary),
+        ])
+
+        @test MultiBlockGrid([g₁₁, g₁₂, g₂₁, g₂₂], C) isa MultiBlockGrid{SVector{2,Float64},2}
+
+        # Matrix
+        g₁₁ = equidistant_grid((0,0),(1,1),5,5)
+        g₁₂ = equidistant_grid((0,0),(1,1),5,5)
+        g₂₁ = equidistant_grid((0,0),(1,1),5,5)
+        g₂₂ = equidistant_grid((0,0),(1,1),5,5)
+
+        C = map(((i1,d1,b1,i2,d2,b2),)->connection(CartesianIndex(i1),CartesianBoundary{d1,b1}(), CartesianIndex(i2), CartesianBoundary{d2,b2}()),[
+            ((1,1), 1, UpperBoundary, (1,2), 1, LowerBoundary),
+            ((2,1), 1, UpperBoundary, (2,2), 1, LowerBoundary),
+            ((1,1), 2, UpperBoundary, (1,2), 2, LowerBoundary),
+            ((2,1), 2, UpperBoundary, (2,2), 2, LowerBoundary),
+        ])
+
+        @test MultiBlockGrid([g₁₁ g₁₂; g₂₁ g₂₂], C) isa MultiBlockGrid{SVector{2,Float64},2}
+
         # Dictionary
         @test MultiBlockGrid(Dict(:a=>g₁, :b=>g₂), [connection(:a, UpperBoundary(), :b, LowerBoundary())]) isa MultiBlockGrid{Float64,1}
     end