diff --git a/meshmode/mesh/processing.py b/meshmode/mesh/processing.py
index 57300fea9dd06325b6b194a9b91902cc027cfd42..8cb8c13f3b16d4c0d3fbc9862b45ac944e226082 100644
--- a/meshmode/mesh/processing.py
+++ b/meshmode/mesh/processing.py
@@ -572,7 +572,6 @@ def merge_disjoint_meshes(meshes, skip_tests=False, single_group=False):
 
 # {{{ split meshes
 
-
 def split_mesh_groups(mesh, element_flags, return_subgroup_mapping=False):
     """Split all the groups in *mesh* according to the values of
     *element_flags*. The element flags are expected to be integers
diff --git a/test/test_meshmode.py b/test/test_meshmode.py
index eb69cf7b5d85e004f534b45f6589bb7a5b19a6a9..0031eb96474b887d0d39c36b99be9933661ff4d1 100644
--- a/test/test_meshmode.py
+++ b/test/test_meshmode.py
@@ -1272,27 +1272,6 @@ def test_is_affine_group_check(mesh_name):
     assert all(grp.is_affine for grp in mesh.groups) == is_affine
 
 
-def split_mesh(mesh, axis, cutoff):
-    groups = []
-    for grp in mesh.groups:
-        mask = np.any(mesh.vertices[axis, grp.vertex_indices] < cutoff, axis=1)
-
-        groups.append(grp.copy(
-                vertex_indices=grp.vertex_indices[mask, :].copy(),
-                nodes=grp.nodes[:, mask, :].copy())
-                )
-        groups.append(grp.copy(
-                vertex_indices=grp.vertex_indices[~mask, :].copy(),
-                nodes=grp.nodes[:, ~mask, :].copy())
-                )
-
-    from meshmode.mesh import Mesh
-    return Mesh(
-            vertices=mesh.vertices,
-            groups=groups,
-            is_conforming=mesh.is_conforming)
-
-
 @pytest.mark.parametrize("ambient_dim", [1, 2, 3])
 def test_mesh_multiple_groups(ctx_factory, ambient_dim, visualize=False):
     ctx = ctx_factory()