hi_res_mask
uround(pow(i / double(gamma_mask), m_gamma) * double(hi_res_mask));
uround(pow(i / double(hi_res_mask), inv_g) * double(gamma_mask));