code: def geom_check_axis(g, atwts, ax,
          nmax=_DEF.SYMM_MATCH_NMAX,
          tol=_DEF.SYMM_MATCH_TOL):
      
  
      # Imports
      import numpy as np
  
      # Store the max found rotation order of the geometry.
      order = geom_find_rotsymm(g, atwts, ax, \
                                          False, nmax, tol)[0]
  
      # Store the presence/absence of a reflection plane.
      refl = geom_symm_match(g, atwts, ax, 0, True) < tol
  
      # Return the pair of values for outside handling
      return order, refl
docs: [Get max proper order and reflection for an axis]
  
      .. todo:: Complete geom_parse_axis docstring