diff --git a/test/combinations.jl b/test/combinations.jl
index ae2fd299..4e015114 100644
--- a/test/combinations.jl
+++ b/test/combinations.jl
@@ -151,10 +151,6 @@ end
     # Package and run tests
     testable = TestableGeometry(integrand, curve, solution)
     runtests(testable, SupportStatus(:line); rtol=0.5e-2)
-
-    # TODO move this to Differentiation testitem
-    # Check Bezier-specific jacobian bounds
-    @test_throws DomainError jacobian(curve, [1.1])
 end
 
 @testitem "Meshes.Box 1D" setup=[Combinations] begin
@@ -188,10 +184,6 @@ end
     # Package and run tests
     testable = TestableGeometry(integrand, box, solution)
     runtests(testable, SupportStatus(:surface); rtol=1e-6)
-
-    # TODO move this to Differentiation testitem
-    # Test jacobian with wrong number of parametric coordinates
-    @test_throws ArgumentError jacobian(box, zeros(3))
 end
 
 @testitem "Meshes.Box 3D" setup=[Combinations] begin
diff --git a/test/utils.jl b/test/utils.jl
index 50cb4ea1..87ef390f 100644
--- a/test/utils.jl
+++ b/test/utils.jl
@@ -19,7 +19,7 @@
     @test _ones(Float32, 2) == (1.0f0, 1.0f0)
 end
 
-@testitem "DifferentiationMethod" setup=[Setup] begin
+@testitem "Differentiation" setup=[Setup] begin
     using MeshIntegrals: _default_diff_method
 
     # _default_diff_method
@@ -29,6 +29,10 @@ end
 
     # FiniteDifference
     @test FiniteDifference().ε ≈ 1e-6
+
+    # Test jacobian with wrong number of parametric coordinates
+    box = Box(Point(0, 0), Point(a, a))    
+    @test_throws ArgumentError jacobian(box, zeros(3))
 end
 
 @testitem "_ParametricGeometry" setup=[Setup] begin