However, Z cannot be calculated without knowing the results of other cluster nodes.
if we use MPI for parallelization, the parallelized version of this routine should resemble the algorithm shown in Figure 1.