# This assumes that you have a 2-dimensional array zeros = np.zeros(shape, dtype=np.int32) zeros[:sigma.shape[0], :sigma.shape[1]] = sigma