062599d6569e877753e2321e9d6c27890c002bf0
[libdcp.git] / src / gamma_transfer_function.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_transfer_function.h"
21 #include <cmath>
22
23 using std::pow;
24 using boost::shared_ptr;
25 using boost::dynamic_pointer_cast;
26 using namespace dcp;
27
28 GammaTransferFunction::GammaTransferFunction (double gamma)
29         : _gamma (gamma)
30 {
31
32 }
33
34 double *
35 GammaTransferFunction::make_lut (int bit_depth, bool inverse) const
36 {
37         int const bit_length = int(std::pow(2.0f, bit_depth));
38         double* lut = new double[bit_length];
39         double const gamma = inverse ? (1 / _gamma) : _gamma;
40         for (int i = 0; i < bit_length; ++i) {
41                 lut[i] = pow(double(i) / (bit_length - 1), gamma);
42         }
43
44         return lut;
45 }
46
47 bool
48 GammaTransferFunction::about_equal (shared_ptr<const TransferFunction> other, double epsilon) const
49 {
50         shared_ptr<const GammaTransferFunction> o = dynamic_pointer_cast<const GammaTransferFunction> (other);
51         if (!o) {
52                 return false;
53         }
54
55         return fabs (_gamma - o->_gamma) < epsilon;
56 }