Python, Numpy, Julia coding

Please convert the following block of code to Julia

a = np.random.randn(9, 6) + 1j*np.random.randn(9, 6)

b = np.random.randn(2, 7, 8, 3) + 1j*np.random.randn(2, 7, 8, 3)

u, s, vh = np.linalg.svd(a, full_matrices=True) u.shape, s.shape, vh.shape ((9, 9), (6,), (6, 6)) np.allclose(a, np.dot(u[:, :6] * s, vh)) True smat = np.zeros((9, 6), dtype=complex) smat[:6, :6] = np.diag(s) np.allclose(a, np.dot(u, np.dot(smat, vh))) True