diff --git a/test/mesh/test_coloring.jl b/test/mesh/test_coloring.jl
new file mode 100644
index 0000000..5e43966
--- /dev/null
+++ b/test/mesh/test_coloring.jl
@@ -0,0 +1,20 @@
+using JuliaFEM, LinearAlgebra, Test
+
+datadir = first(splitext(basename(@__FILE__)))
+
+@testset "Test that mesh coloring provides a good coloring"
+    fn = joinpath(datadir, "cube_tet4.inp")
+    mesh = JuliaFEM.Mesh(open(parse_abaqus, fn))
+
+    coloring = JuliaFEM.create_coloring(mesh)
+    for colors in coloring
+        for ele_i in colors
+            for ele_j in colors
+                if ele_i == ele_j
+                    continue
+                end
+                @test intersect(Set(mesh.elements[ele_i]), Set(mesh.elements[ele_j])) |> isempty
+            end
+        end
+    end
+end
\ No newline at end of file