code: def symm_op(g, ax, theta, do_refl):
      
  
      # Imports
      import numpy as np
  
      # Depend on lower functions' geometry vector coercion. Just
      #  do the rotation and, if indicated, the reflection.
      gx = geom_rotate(g, ax, theta)
      if do_refl:
          gx = geom_reflect(gx, ax)
      ## end if
  
      # Should be good to go
      return gx
docs: Perform general point symmetry operation on a geometry.
  
      .. todo:: Complete symm_op docstring