from scipy import linalg
import numpy as np
a = np.random.randn(3, 2) + 1.j*np.random.randn(3, 2)
U, s, Vh = linalg.svd(a)
print(U, Vh, s)