diff --git a/Manifest.toml b/Manifest.toml
index 29e0f1c..d92b83f 100644
--- a/Manifest.toml
+++ b/Manifest.toml
@@ -2,7 +2,7 @@
 
 julia_version = "1.12.6"
 manifest_format = "2.0"
-project_hash = "9a67e6793f96d4456a3ac13f063b9d5ad1df8c6a"
+project_hash = "5d2e867a4299eb5619f304e6a153e51462cc7b54"
 
 [[deps.Adapt]]
 deps = ["LinearAlgebra", "Requires"]
@@ -192,10 +192,14 @@ uuid = "f80590ac-b429-510a-8a99-e7c46989f22d"
 version = "0.5.2"
 
     [deps.JuliaFEM.extensions]
+    JuliaFEMGmshExt = ["Gmsh"]
     JuliaFEMMPIExt = ["MPI"]
+    JuliaFEMWriteVTKExt = ["WriteVTK"]
 
     [deps.JuliaFEM.weakdeps]
+    Gmsh = "705231aa-382f-11e9-3f0c-b7cb4346fdeb"
     MPI = "da04e1cc-30fd-572f-bb4f-1f8673147195"
+    WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 
 [[deps.JuliaSyntaxHighlighting]]
 deps = ["StyledStrings"]
diff --git a/Project.toml b/Project.toml
index a60d930..16f6dd8 100644
--- a/Project.toml
+++ b/Project.toml
@@ -4,8 +4,8 @@ version = "0.5.2"
 author = ["Jukka Aho <ahojukka5@gmail.com>", "Tero Frondelius <tero.frondelius@gmail.com>", "Olli Väinölä"]
 
 [deps]
-CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
+CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
 IterativeSolvers = "42fd0dbc-a981-5370-80f2-aaf504508153"
 KernelAbstractions = "63c18a36-062a-441e-b654-da1e3ab1ce7c"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
@@ -15,17 +15,29 @@ SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Tensors = "48a634ad-e948-5137-8d70-aa71f2a747f4"
 
+[weakdeps]
+Gmsh = "705231aa-382f-11e9-3f0c-b7cb4346fdeb"
+MPI = "da04e1cc-30fd-572f-bb4f-1f8673147195"
+WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
+
+[extensions]
+JuliaFEMGmshExt = ["Gmsh"]
+JuliaFEMMPIExt = ["MPI"]
+JuliaFEMWriteVTKExt = ["WriteVTK"]
+
 [compat]
-Aqua = "0.8"
-CSV = "0.10"
 Adapt = "4.5.2"
+Aqua = "0.8"
 BenchmarkTools = "1.6.3"
+CSV = "0.10"
+Gmsh = "0.3"
 IterativeSolvers = "0.9.4"
 KernelAbstractions = "0.9.41"
 LinearOperators = "2.11.0"
 MPI = "0.20"
 StaticArrays = "1.9.15"
 Tensors = "≥ 1.0.0"
+WriteVTK = "1.21.2"
 julia = "1.9"
 
 [extras]
@@ -36,12 +48,7 @@ MPI = "da04e1cc-30fd-572f-bb4f-1f8673147195"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
+WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 
 [targets]
-test = ["Test", "Statistics", "BenchmarkTools", "InteractiveUtils", "Random", "MPI", "Aqua"]
-
-[weakdeps]
-MPI = "da04e1cc-30fd-572f-bb4f-1f8673147195"
-
-[extensions]
-JuliaFEMMPIExt = ["MPI"]
+test = ["Test", "Statistics", "BenchmarkTools", "InteractiveUtils", "Random", "MPI", "Aqua", "WriteVTK"]