Fix claim for number of wells in an unoptimized stellarator

In general the number of wells is upper bound by the maximum frequency
along the field line of the Fourier representation in straight field line
coordinates of M + iota N.

g : ϑ, ϕ ↦ ∑ₘₙ aₘₙ exp(j [mϑ + nϕ])
g : α, ϕ ↦ ∑ₘₙ aₘₙ exp(j [mα + (m ι + n)ϕ])

desc/integrals/bounce_integral.py

@@ -223,8 +223,9 @@ def points(self, pitch, num_well=None):
         num_well : int or None
             Specify to return the first ``num_well`` pairs of bounce points for each
             pitch along each field line. This is useful if ``num_well`` tightly
-            bounds the actual number. As a reference, there are typically at most 5
-            wells per toroidal transit for a given pitch.
+            bounds the actual number. As a reference, there are typically 20 wells
+            per toroidal transit for a given pitch. You can check this by plotting
+            the field lines with the ``check_points`` method.
 
             If not specified, then all bounce points are returned. If there were fewer
             wells detected along a field line than the size of the last axis of the
@@ -327,8 +328,9 @@ def integrate(
         num_well : int or None
             Specify to return the first ``num_well`` pairs of bounce points for each
             pitch along each field line. This is useful if ``num_well`` tightly
-            bounds the actual number. As a reference, there are typically at most 5
-            wells per toroidal transit for a given pitch.
+            bounds the actual number. As a reference, there are typically 20 wells
+            per toroidal transit for a given pitch. You can check this by plotting
+            the field lines with the ``check_points`` method.
 
             If not specified, then all bounce points are returned. If there were fewer
             wells detected along a field line than the size of the last axis of the

desc/integrals/bounce_utils.py

@@ -138,8 +138,9 @@ def bounce_points(
     num_well : int or None
         Specify to return the first ``num_well`` pairs of bounce points for each
         pitch along each field line. This is useful if ``num_well`` tightly
-        bounds the actual number. As a reference, there are typically at most 5
-        wells per toroidal transit for a given pitch.
+        bounds the actual number. As a reference, there are typically 20 wells
+        per toroidal transit for a given pitch. You can check this by plotting
+        the field lines with the ``_check_bounce_points`` method.
 
         If not specified, then all bounce points are returned. If there were fewer
         wells detected along a field line than the size of the last axis of the