def METHOD_NAME(m, n, k, dtype=complex):
    s = np.sort(qu.randn(k)**2)[::-1]
    U = qu.gen.rand.rand_iso(m, k, dtype=dtype)
    VH = qu.gen.rand.rand_iso(n, k, dtype=dtype).conj().T
    if U.dtype in ('float32', 'complex64'):
        s = s.astype('float32')
    return usv2dense(U, s, VH)