A preconditioned conjugate gradient (PCG)-based domain decomposition method was given in  and  for the solution of linear equations arising in the finite element method applied to the elliptic Neumann problem. The novelty of the proposed algorithm was that the recommended preconditioner was constructed by using symmetric-cyclic matrix. But we could give only the definitions of the entries of this cyclic matrix. Here we give a short description of this algorithm, the method of calculation of matrix entries and the results of calculation. The numerical experiments presented show, that this construction of precondition in the practice works well.