diff --git a/subtract/subtract_with_wsclean.py b/subtract/subtract_with_wsclean.py
index 57ad494..5fc915b 100644
--- a/subtract/subtract_with_wsclean.py
+++ b/subtract/subtract_with_wsclean.py
@@ -582,6 +582,7 @@ def run_DP3(self, phaseshift: str = None, freqavg: str = None,
             command += ['ps.type=phaseshifter',
                         'ps.phasecenter=' + phasecenter]
 
+        # For facet subtraction we can take some shortcuts and get major speed improvements
         if speedup_facet_subtract:
             # 2) AVERAGING
             if freqavg is not None or timeres is not None:
@@ -628,6 +629,7 @@ def run_DP3(self, phaseshift: str = None, freqavg: str = None,
                             ac_count += 1
                     T.close()
 
+        # Most accurate version, but more expensive than the --speedup_facet_subtract approach
         else:
             # 2) APPLY BEAM
             if applybeam: