diff --git a/sumpy/expansion/__init__.py b/sumpy/expansion/__init__.py
index 480f53fcdd1efc7dba7562f96ebcd7d5dcbb804e..647ebb01589a502c128a30aa597377c8fed1d0e4 100644
--- a/sumpy/expansion/__init__.py
+++ b/sumpy/expansion/__init__.py
@@ -458,6 +458,12 @@ class LinearRecurrenceBasedExpansionTermsWrangler(ExpansionTermsWrangler):
     def get_reduced_coeffs(self, nullspace):
         """
         Returns the indices of the reduced set of derivatives which are stored.
+        Override this method if the reduced set is known analytically.
+
+        This method does elementary row operations to figure out which rows are
+        linearly dependent on the previous rows. Partial pivoting is not done
+        to preserve the order so that a row is not linearly dependent on a row
+        that came after in the original row order.
         """
         mat = nullspace.copy()
         nrows = mat.shape[0]