diff --git a/jaxpm/pm.py b/jaxpm/pm.py
index 7aac781..d9870f7 100644
--- a/jaxpm/pm.py
+++ b/jaxpm/pm.py
@@ -75,6 +75,9 @@ def make_ode_fn(mesh_shape):
 def pgd_correction(pos, params):
     """
     improve the short-range interactions of PM-Nbody simulations with potential gradient descent method, based on https://arxiv.org/abs/1804.00671
+    args:
+      pos: particle positions [npart, 3]
+      params: [alpha, kl, ks] pgd parameters
     """
     kvec = fftk(mesh_shape)