Why averaging for high harmonics in build_ellipse_model

[asterg1122]

photutils version: 1.12.0 In file isophote/model.py, Line 124, Col 64

There is a "/ 4.0" in the calculation of harm (contribution to isophote intensity from high harmonics) that takes the average of the term $A_3\sin(3\phi)+B_3\cos(3\phi)+A_4\sin(4\phi)+B_4\cos(4\phi)$. But according to the formula, that term is directly added into the intensity.

Even though the isophote intensity is weighted and added to the overlapping four pixels, the total intensity is intens + harm and seems not to need any averaging for harm. So I wonder why it gets averaged. Thank you very much!