diff --git a/test/python/quantum_info/xx_decompose/test_polytopes.py b/test/python/quantum_info/xx_decompose/test_polytopes.py
index 20f3bc6ea65f..b3d383389270 100644
--- a/test/python/quantum_info/xx_decompose/test_polytopes.py
+++ b/test/python/quantum_info/xx_decompose/test_polytopes.py
@@ -46,7 +46,6 @@ def test_nearest(self, offbody, expected):
         """Check that the nearest point calculator recovers some known cases."""
         polytope = XXPolytope.from_strengths(pi / 6, pi / 8, pi / 10)
         result = polytope.nearest(np.array(offbody))
-        print(offbody, result)
         self.assertTrue(np.all(np.abs(np.array(expected) - result) < EPSILON))
 
     def test_add_strengths(self):
diff --git a/test/python/quantum_info/xx_decompose/test_weyl.py b/test/python/quantum_info/xx_decompose/test_weyl.py
index 3d8a915d82f4..ebe942d6f083 100644
--- a/test/python/quantum_info/xx_decompose/test_weyl.py
+++ b/test/python/quantum_info/xx_decompose/test_weyl.py
@@ -52,10 +52,6 @@ def test_reflections(self):
             original_matrix = canonical_matrix(*coordinate)
             reflected_matrix = canonical_matrix(*reflected_coordinate)
             reflect_matrix = Operator(reflection_circuit).data
-            with np.printoptions(precision=3, suppress=True):
-                print(name)
-                print(reflect_matrix.conjugate().transpose(1, 0) @ original_matrix @ reflect_matrix)
-                print(reflected_matrix * reflection_phase)
             self.assertTrue(
                 np.all(
                     np.abs(
@@ -76,10 +72,6 @@ def test_shifts(self):
             original_matrix = canonical_matrix(*coordinate)
             shifted_matrix = canonical_matrix(*shifted_coordinate)
             shift_matrix = Operator(shift_circuit).data
-            with np.printoptions(precision=3, suppress=True):
-                print(name)
-                print(original_matrix @ shift_matrix)
-                print(shifted_matrix * shift_phase)
             self.assertTrue(
                 np.all(
                     np.abs(original_matrix @ shift_matrix - shifted_matrix * shift_phase) < EPSILON