Analysis for the number of pixels along the equator N_eq of the planet (top panels) and moon (bottom panels). Shown are the maximum (solid line) and mean (dashed line) differences between results computed across the whole phase angle range and for consecutive values of N_eq values, for the reflected flux F(α) relative to F(α = 0°) (panels a and d), and degree of polarization P (panels b and e, note that for the moon, P = 0). Also shown is the computational time (in minutes) for the computation of a full phase curve (with 24 h temporal resolution) and the average disk integration (panels c and f). For (top panels), we used values of 10, 20, 30, 40, 50, 60, and 70, and for (bottom panels), 3, 5, 8, 11, 14, 16, and 19.