Using two_level_mesh in example001

GalerkinToolkit · Feb 5, 2024 · fdc7408 · fdc7408
1 parent 69583d3
commit fdc7408
Show file tree

Hide file tree

Showing 2 changed files with 14 additions and 1 deletion.
diff --git a/extensions/GalerkinToolkitExamples/Project.toml b/extensions/GalerkinToolkitExamples/Project.toml
@@ -11,7 +11,6 @@ PartitionedArrays = "5a9dfac6-5c52-46f7-8278-5e2210713be9"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
-prrte_jll = "eb928a42-fffd-568d-ab9c-3f5d54fc65b9"
 
 [extras]
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

diff --git a/extensions/GalerkinToolkitExamples/test/example001_tests.jl b/extensions/GalerkinToolkitExamples/test/example001_tests.jl
@@ -42,5 +42,19 @@ results = Example001.main(params)
 @test results[:eh1] < tol
 @test results[:el2] < tol
 
+domain = (0,1,0,1)
+cells = (10,10)
+fine_mesh = gk.cartesian_mesh(domain,cells)
+domain = (0,30,0,10)
+cells = (2,2)
+coarse_mesh = gk.cartesian_mesh(domain,cells)
+mesh, = gk.two_level_mesh(coarse_mesh,fine_mesh)
+
+params = Dict{Symbol,Any}()
+params[:mesh] = mesh
+results = Example001.main(params)
+@test results[:eh1] < tol
+@test results[:el2] < tol
+
 
 end # module