src/gamma_lut.cc

   1 /*
   2     Copyright (C) 2012-2014 Carl Hetherington <cth@carlh.net>
   3
   4     This program is free software; you can redistribute it and/or modify
   5     it under the terms of the GNU General Public License as published by
   6     the Free Software Foundation; either version 2 of the License, or
   7     (at your option) any later version.
   8
   9     This program is distributed in the hope that it will be useful,
  10     but WITHOUT ANY WARRANTY; without even the implied warranty of
  11     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12     GNU General Public License for more details.
  13
  14     You should have received a copy of the GNU General Public License
  15     along with this program; if not, write to the Free Software
  16     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  17
  18 */
  19
  20 #include "gamma_lut.h"
  21 #include "lut_cache.h"
  22 #include <cmath>
  23
  24 using namespace dcp;
  25
  26 LUTCache<GammaLUT> GammaLUT::cache;
  27
  28 GammaLUT::GammaLUT (int bit_depth, float gamma, bool linearised)
  29         : _bit_depth (bit_depth)
  30         , _gamma (gamma)
  31         , _linearised (linearised)
  32 {
  33         _lut = new float[int(std::pow(2.0f, _bit_depth))];
  34         int const bit_length = pow (2, _bit_depth);
  35
  36         if (_linearised) {
  37                 for (int i = 0; i < bit_length; ++i) {
  38                         float const p = static_cast<float> (i) / (bit_length - 1);
  39                         if (p > 0.04045) {
  40                                 _lut[i] = pow ((p + 0.055) / 1.055, gamma);
  41                         } else {
  42                                 _lut[i] = p / 12.92;
  43                         }
  44                 }
  45         } else {
  46                 for (int i = 0; i < bit_length; ++i) {
  47                         _lut[i] = pow(float(i) / (bit_length - 1), gamma);
  48                 }
  49         }
  50 }