Abstract. Christopher in 2006 proved that under some assumptions the linear parts of the Lyapunov constants with respect to the parameters give the cyclicity of an elementary center. This paper is devote to establish a new approach, namely parallelization, to compute the linear parts of the Lyapunov constants. More concretely, it is showed that parallelization computes these linear parts in a shorter quantity of time than other traditional mechanisms.To show the power of this approach, we study the cyclicity of the holomorphic centerż = iz + z 2 + z 3 + · · · + z n under general polynomial perturbations of degree n, for n ≤ 13. We also exhibit that, from the point of view of computation, among the Hamiltonian, time-reversible, and Darboux centers, the holomorphic center is the best candidate to obtain high cyclicity examples of any degree. For n = 4, 5, . . . , 13, we prove that the cyclicity of the holomorphic center is at least n 2 + n − 2. This result give the highest lower bound for M (6), M (7), . . . , M (13) among the existing results, where M (n) is the maximum number of limit cycles bifurcating from an elementary monodromic singularity of polynomial systems of degree n. As a direct corollary we also obtain the highest lower bound for the Hilbert numbers H(6) ≥ 40, H(8) ≥ 70, and H