fourier bounce

desc/grid.py

@@ -1037,8 +1037,10 @@ def _create_nodes(  # noqa: C901
         """
         self._NFP = check_posint(NFP, "NFP", False)
         self._period = (np.inf, 2 * np.pi, 2 * np.pi / self._NFP)
-        # TODO:
+        # FIXME:
         #  https://github.com/PlasmaControl/DESC/pull/1204#pullrequestreview-2246771337
+        #  Quantities like alpha, grad(alpha), etc. are computed incorrectly at
+        #  phi > 2pi / NFP and theta > 2pi / NFP.
         axis = bool(axis)
         endpoint = bool(endpoint)
         theta_period = self.period[1]

desc/integrals/__init__.py

@@ -1,6 +1,6 @@
 """Classes for function integration."""
 
-from .bounce_integral import Bounce1D
+from .bounce_integral import Bounce1D, Bounce2D
 from .singularities import (
     DFTInterpolator,
     FFTInterpolator,