Perform general point symmetry operation on a geometry.

    .. todo:: Complete symm_op docstring

def symm_op(g, ax, theta, do_refl):
    """ Perform general point symmetry operation on a geometry.

    .. todo:: Complete symm_op docstring

    """

    # 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