SCALE_LOG
#define SCALE (1 << SCALE_LOG)
scale = ((y4 / SCALE_LOG) * (L + 1) -
(y2 / SCALE_LOG) * sumy2) * SCALE_LOG;