Fourier-space weights needed to evaluate I_{mm'}

    This is mostly an internal function, included here for backwards compatibility.  See map2salm
    and salm2map for more useful functions.

def quadrature_weights(Ntheta):
    """Fourier-space weights needed to evaluate I_{mm'}

    This is mostly an internal function, included here for backwards compatibility.  See map2salm
    and salm2map for more useful functions.

    """
    import numpy as np
    weights = np.empty(2*(Ntheta-1), dtype=np.complex128)
    _quadrature_weights(Ntheta, weights)
    return weights