The iterative IGD method is endowed with scalability for eigen-analysis of large DCPPS. Blocks of infinitesimal generator's approximant are reformulated as Kronecker products. The dominant block is factorized as sum of Kronecker products of system state matrices. Sparsities in infinitesimal generator's approximant and system matrices are exploited. Reformulation, preconditioning and IDR(s) ensure efficiency and scalability of IIGD.