@silence_errors
    def _cos_factorization(self, alpha, omega, phase):
        r1 = np.sum(alpha*np.cos(phase),axis=1)[:,None]
        r2 = np.sum(alpha*np.sin(phase),axis=1)[:,None]
        r =  np.sqrt(r1**2 + r2**2)
        psi = np.where(r1 != 0, (np.arctan(r2/r1) + (0. < r1)*np.pi),np.arcsin(r2))
        return r,omega[:,0:1], psi
dataset/ETHPy150Open SheffieldML/GPy/GPy/kern/src/periodic.py/Periodic._cos_factorization/swapped operands of `<`