def qd(A, B, output='real', lwork=None, sort=None, overwrite_a=False, overwrite_b=False, check_finite=True):
 (result, _) = _qd(A, B, output=output, lwork=lwork, sort=sort, overwrite_a=overwrite_a, overwrite_b=overwrite_b, check_finite=check_finite)
 return (result[0], result[1], result[(-4)], result[(-3)])

def qd ( A , B , output = 'real' , lwork = None , sort = None , overwrite_a = False , overwrite_b = False , check_finite = True ) : ( result , _ ) = _qd ( A , B , output = output , lwork = lwork , sort = sort , overwrite_a = overwrite_a , overwrite_b = overwrite_b , check_finite = check_finite ) return ( result [ 0 ] , result [ 1 ] , result [ ( - 4 ) ] , result [ ( - 3 ) ] )
qz decomposition for generalized eigenvalues of a pair of matrices .