A new preconditioned conjugate gradient (PCG)-based domain decomposition method is given for the solution of linear equations
arising in the finite element method applied to the elliptic Neumann problem. The novelty of the proposed method is in the
recommended preconditioner which is constructed by using cyclic matrix. The resulting preconditioned algorithms are well suited
to parallel computation.