The paper presents a comprehensive result for the 2+1 flavor QCD equation of state (EoS), extending previous studies by using finer lattices and including ensembles with different temporal extensions ($N_t = 6, 8, 10, 12, 16$). The authors employed a Symanzik improved gauge action and a stout-link improved staggered fermion action. They confirmed earlier findings and provided a parametrization of the EoS, making it accessible for further use. The main advancements include controlled continuum extrapolation, improved systematic error estimation, and a more accurate determination of the scale and line of constant physics (LCP). The results are consistent with the hadron resonance gas model at low temperatures and show good agreement with earlier studies, with minor differences noted in the high-temperature region. The paper also discusses the challenges and improvements in the simulation setup, emphasizing the computational efficiency and the well-behaved continuum extrapolation of the action.The paper presents a comprehensive result for the 2+1 flavor QCD equation of state (EoS), extending previous studies by using finer lattices and including ensembles with different temporal extensions ($N_t = 6, 8, 10, 12, 16$). The authors employed a Symanzik improved gauge action and a stout-link improved staggered fermion action. They confirmed earlier findings and provided a parametrization of the EoS, making it accessible for further use. The main advancements include controlled continuum extrapolation, improved systematic error estimation, and a more accurate determination of the scale and line of constant physics (LCP). The results are consistent with the hadron resonance gas model at low temperatures and show good agreement with earlier studies, with minor differences noted in the high-temperature region. The paper also discusses the challenges and improvements in the simulation setup, emphasizing the computational efficiency and the well-behaved continuum extrapolation of the action.