diff --git a/extensions/GalerkinToolkitExamples/test/example002_defs.jl b/extensions/GalerkinToolkitExamples/test/example002_defs.jl
new file mode 100644
index 00000000..4e6eb1c5
--- /dev/null
+++ b/extensions/GalerkinToolkitExamples/test/example002_defs.jl
@@ -0,0 +1,39 @@
+
+import GalerkinToolkit as gk
+using GalerkinToolkitExamples: Example002
+using Test
+using PartitionedArrays
+using PetscCall
+
+function example002_tests_np_4(distribute)
+    tol = 1.0e-8
+    params = Dict{Symbol,Any}()
+    domain = (0,10,0,10)
+    cells_per_dir = (20,20)
+    parts_per_dir = (2,2)
+    np = prod(parts_per_dir)
+    parts = distribute(LinearIndices((np,)))
+    ghost_layers = 0
+    mesh = gk.cartesian_mesh(domain,cells_per_dir,parts_per_dir;parts,ghost_layers)
+    params[:mesh] = mesh
+    results = Example002.main(params)
+    results = Example002.main(params)
+    @test results[:eh1] < tol
+    @test results[:el2] < tol
+
+    options = "-pc_type gamg -ksp_type cg -ksp_error_if_not_converged true -ksp_converged_reason -ksp_rtol 1.0e-6"
+    PetscCall.init(args=split(options))
+    params = Dict{Symbol,Any}()
+    domain = (0,10,0,10,0,10)
+    cells_per_dir = (40,40,40)
+    parts_per_dir = (2,2,1)
+    np = prod(parts_per_dir)
+    parts = distribute(LinearIndices((np,)))
+    ghost_layers = 0
+    mesh = gk.cartesian_mesh(domain,cells_per_dir,parts_per_dir;parts,ghost_layers)
+    params[:mesh] = mesh
+    params[:export_vtu] = false
+    params[:solver] = Example002.ksp_solver()
+    results = Example002.main(params)
+    results = Example002.main(params)
+end