For the Cellular Yule-Nielsen Spectral Neugebauer (CYNSN) model, the inverse nonlinear Gauss-Seidel iterative algorithm is studied. It is proposed to use the XYZ method of multi-illuminant dimensionality reduction to achieve the accelerated inverse CYNSN model, and use 3-ink and 1-ink printer data to test. It was found that for the 3-ink combination, the XYZ dimension reduction using one light source (D50), and the 4-ink combination, the reverse CYNSN model of XYZ dimension reduction under two light sources (D50, Lambda), not only in accuracy (more chromatic aberration under the light source), and have advantages in calculation speed.