diff --git a/test/test_meshmode.py b/test/test_meshmode.py
index d8cede25eb59b89f5f51354a8669825057dad180..7dab556460aa459251c59c1ab7a0eb382dd31bce 100644
--- a/test/test_meshmode.py
+++ b/test/test_meshmode.py
@@ -1308,6 +1308,13 @@ def test_mesh_multiple_groups(ctx_factory, ambient_dim, visualize=False):
     assert mesh.facial_adjacency_groups
     assert mesh.nodal_adjacency
 
+    if visualize:
+        from meshmode.mesh.visualization import draw_2d_mesh
+        draw_2d_mesh(mesh, draw_element_numbers=True, draw_face_numbers=True,
+                set_bounding_box=True)
+        import matplotlib.pyplot as plt
+        plt.show()
+
     from meshmode.discretization import Discretization
     from meshmode.discretization.poly_element import \
             PolynomialWarpAndBlendGroupFactory as GroupFactory