yz_mat1 = (np.matrix(np.diag(Zeros_G)) *
                       np.matrix(np.ones([len(Zeros_G), len(Zeros_G)])))