def METHOD_NAME(Lmax,m,s,cos_theta,**kwargs):
    """
        Gives spin-wieghted spherical harmonic functions on the Gauss quaduature grid.
        Returns an array with
            shape = ( Lmax - Lmin(m,s) + 1, len(z) )
                 or (Lmax - Lmin(m,s) + 1,) if z is a single point.
        Parameters
        ----------
        Lmax: int >=0; spherical-harmonic degree.
        m,s : int
            spherical harmonic parameters.
        cos_theta: np.ndarray or float.
        dtype: output dtype. internal dtype = 'longdouble'.
        """
    n,a,b = spin2Jacobi(Lmax,m,s)
    init = np.exp(0.5*Jacobi.measure(a,b,cos_theta,log=True))
    init *= ((-1.)**max(m,-s))
    return Jacobi.polynomials(n,a,b,cos_theta,init,**kwargs)