- for (i = 0;i < 256;i++)
- {
- adjusted = (int) (i * 256.0 * prescale * scale + base * 65535.0);
- out[i] = bound(0, adjusted, 65535);
- }
- }
- else
- {
- invgamma = 1.0 / gamma;
- prescale /= 255.0f;
- for (i = 0;i < 256;i++)
- {
- d = pow((double) i * prescale, invgamma) * scale + base;
- adjusted = (int) (65535.0 * d);
- out[i] = bound(0, adjusted, 65535);
- }
+ adjusted = (int) (65535.0 * (pow((double) i * prescale, invgamma) * scale + base) + 0.5);
+ out[i] = bound(0, adjusted, 65535);