diff --git a/skglm/solvers/prox_newton.py b/skglm/solvers/prox_newton.py
index 8224938cc..4b8e0aaf7 100644
--- a/skglm/solvers/prox_newton.py
+++ b/skglm/solvers/prox_newton.py
@@ -140,8 +140,6 @@ def solve(self, X, y, datafit, penalty, w_init=None, Xw_init=None):
             # similar to np.argsort()[-ws_size:] but without sorting
             ws = np.argpartition(opt, -ws_size)[-ws_size:]
 
-            print(f"Iter {t:<3}: {ws}")
-
             grad_ws = grad[ws]
             tol_in = EPS_TOL * stop_crit