diff --git a/sumpy/expansion/__init__.py b/sumpy/expansion/__init__.py
index a761d1be40da04d1a88143de203e2feb2289e3e8..d72918eb5695b9b2bf7d5c8793e440aa11a7f289 100644
--- a/sumpy/expansion/__init__.py
+++ b/sumpy/expansion/__init__.py
@@ -687,6 +687,9 @@ class LinearPDEBasedExpansionTermsWrangler(ExpansionTermsWrangler):
                 shape = (len(mis), len(mis))
                 op = CSEMatVecOperator(from_input_coeffs_by_row,
                                        from_output_coeffs_by_row, shape)
+
+                plog.done()
+
                 return mis, op
 
         ordering_key, _ = self._get_mi_ordering_key_and_axis_permutation()