Home | History | Annotate | Download | only in src
      1 /* vim: set ts=8 sw=8 noexpandtab: */
      2 //  qcms
      3 //  Copyright (C) 2009 Mozilla Foundation
      4 //
      5 // Permission is hereby granted, free of charge, to any person obtaining
      6 // a copy of this software and associated documentation files (the "Software"),
      7 // to deal in the Software without restriction, including without limitation
      8 // the rights to use, copy, modify, merge, publish, distribute, sublicense,
      9 // and/or sell copies of the Software, and to permit persons to whom the Software
     10 // is furnished to do so, subject to the following conditions:
     11 //
     12 // The above copyright notice and this permission notice shall be included in
     13 // all copies or substantial portions of the Software.
     14 //
     15 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
     16 // EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO
     17 // THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
     18 // NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
     19 // LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
     20 // OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
     21 // WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
     22 
     23 #include "qcms.h"
     24 #include "qcmstypes.h"
     25 
     26 /* used as a lookup table for the output transformation.
     27  * we refcount them so we only need to have one around per output
     28  * profile, instead of duplicating them per transform */
     29 struct precache_output
     30 {
     31 	int ref_count;
     32 	/* We previously used a count of 65536 here but that seems like more
     33 	 * precision than we actually need.  By reducing the size we can
     34 	 * improve startup performance and reduce memory usage. ColorSync on
     35 	 * 10.5 uses 4097 which is perhaps because they use a fixed point
     36 	 * representation where 1. is represented by 0x1000. */
     37 #define PRECACHE_OUTPUT_SIZE 8192
     38 #define PRECACHE_OUTPUT_MAX (PRECACHE_OUTPUT_SIZE-1)
     39 	uint8_t data[PRECACHE_OUTPUT_SIZE];
     40 };
     41 
     42 #ifdef _MSC_VER
     43 #define ALIGN __declspec(align(16))
     44 #else
     45 #define ALIGN __attribute__(( aligned (16) ))
     46 #endif
     47 
     48 typedef struct _qcms_format_type {
     49 	int r;
     50 	int b;
     51 } qcms_format_type;
     52 
     53 struct _qcms_transform {
     54 	float ALIGN matrix[3][4];
     55 	float *input_gamma_table_r;
     56 	float *input_gamma_table_g;
     57 	float *input_gamma_table_b;
     58 
     59 	float *input_clut_table_r;
     60 	float *input_clut_table_g;
     61 	float *input_clut_table_b;
     62 	uint16_t input_clut_table_length;
     63 	float *r_clut;
     64 	float *g_clut;
     65 	float *b_clut;
     66 	uint16_t grid_size;
     67 	float *output_clut_table_r;
     68 	float *output_clut_table_g;
     69 	float *output_clut_table_b;
     70 	uint16_t output_clut_table_length;
     71 
     72 	float *input_gamma_table_gray;
     73 
     74 	float out_gamma_r;
     75 	float out_gamma_g;
     76 	float out_gamma_b;
     77 
     78 	float out_gamma_gray;
     79 
     80 	uint16_t *output_gamma_lut_r;
     81 	uint16_t *output_gamma_lut_g;
     82 	uint16_t *output_gamma_lut_b;
     83 
     84 	uint16_t *output_gamma_lut_gray;
     85 
     86 	size_t output_gamma_lut_r_length;
     87 	size_t output_gamma_lut_g_length;
     88 	size_t output_gamma_lut_b_length;
     89 
     90 	size_t output_gamma_lut_gray_length;
     91 
     92 	struct precache_output *output_table_r;
     93 	struct precache_output *output_table_g;
     94 	struct precache_output *output_table_b;
     95 
     96 	void (*transform_fn)(struct _qcms_transform *transform, unsigned char *src, unsigned char *dest, size_t length, struct _qcms_format_type output_format);
     97 };
     98 
     99 struct matrix {
    100 	float m[3][3];
    101 	bool invalid;
    102 };
    103 
    104 struct qcms_modular_transform;
    105 
    106 typedef void (*transform_module_fn_t)(struct qcms_modular_transform *transform, float *src, float *dest, size_t length);
    107 
    108 struct qcms_modular_transform {
    109 	struct matrix matrix;
    110 	float tx, ty, tz;
    111 
    112 	float *input_clut_table_r;
    113 	float *input_clut_table_g;
    114 	float *input_clut_table_b;
    115 	uint16_t input_clut_table_length;
    116 	float *r_clut;
    117 	float *g_clut;
    118 	float *b_clut;
    119 	uint16_t grid_size;
    120 	float *output_clut_table_r;
    121 	float *output_clut_table_g;
    122 	float *output_clut_table_b;
    123 	uint16_t output_clut_table_length;
    124 
    125 	uint16_t *output_gamma_lut_r;
    126 	uint16_t *output_gamma_lut_g;
    127 	uint16_t *output_gamma_lut_b;
    128 
    129 	size_t output_gamma_lut_r_length;
    130 	size_t output_gamma_lut_g_length;
    131 	size_t output_gamma_lut_b_length;
    132 
    133 	transform_module_fn_t transform_module_fn;
    134 	struct qcms_modular_transform *next_transform;
    135 };
    136 
    137 typedef int32_t s15Fixed16Number;
    138 typedef uint16_t uInt16Number;
    139 typedef uint8_t uInt8Number;
    140 
    141 struct XYZNumber {
    142 	s15Fixed16Number X;
    143 	s15Fixed16Number Y;
    144 	s15Fixed16Number Z;
    145 };
    146 
    147 struct curveType {
    148 	uint32_t type;
    149 	uint32_t count;
    150 	float parameter[7];
    151 	uInt16Number data[];
    152 };
    153 
    154 struct lutmABType {
    155 	uint8_t num_in_channels;
    156 	uint8_t num_out_channels;
    157 	// 16 is the upperbound, actual is 0..num_in_channels.
    158 	uint8_t num_grid_points[16];
    159 
    160 	s15Fixed16Number e00;
    161 	s15Fixed16Number e01;
    162 	s15Fixed16Number e02;
    163 	s15Fixed16Number e03;
    164 	s15Fixed16Number e10;
    165 	s15Fixed16Number e11;
    166 	s15Fixed16Number e12;
    167 	s15Fixed16Number e13;
    168 	s15Fixed16Number e20;
    169 	s15Fixed16Number e21;
    170 	s15Fixed16Number e22;
    171 	s15Fixed16Number e23;
    172 
    173 	// reversed elements (for mBA)
    174 	bool reversed;
    175 
    176 	float *clut_table;
    177 	struct curveType *a_curves[10];
    178 	struct curveType *b_curves[10];
    179 	struct curveType *m_curves[10];
    180 	float clut_table_data[];
    181 };
    182 
    183 /* should lut8Type and lut16Type be different types? */
    184 struct lutType { // used by lut8Type/lut16Type (mft2) only
    185 	uint8_t num_input_channels;
    186 	uint8_t num_output_channels;
    187 	uint8_t num_clut_grid_points;
    188 
    189 	s15Fixed16Number e00;
    190 	s15Fixed16Number e01;
    191 	s15Fixed16Number e02;
    192 	s15Fixed16Number e10;
    193 	s15Fixed16Number e11;
    194 	s15Fixed16Number e12;
    195 	s15Fixed16Number e20;
    196 	s15Fixed16Number e21;
    197 	s15Fixed16Number e22;
    198 
    199 	uint16_t num_input_table_entries;
    200 	uint16_t num_output_table_entries;
    201 
    202 	float *input_table;
    203 	float *clut_table;
    204 	float *output_table;
    205 
    206 	float table_data[];
    207 };
    208 #if 0
    209 /* this is from an intial idea of having the struct correspond to the data in
    210  * the file. I decided that it wasn't a good idea.
    211  */
    212 struct tag_value {
    213 	uint32_t type;
    214 	union {
    215 		struct {
    216 			uint32_t reserved;
    217 			struct {
    218 				s15Fixed16Number X;
    219 				s15Fixed16Number Y;
    220 				s15Fixed16Number Z;
    221 			} XYZNumber;
    222 		} XYZType;
    223 	};
    224 }; // I guess we need to pack this?
    225 #endif
    226 
    227 #define RGB_SIGNATURE  0x52474220
    228 #define GRAY_SIGNATURE 0x47524159
    229 #define XYZ_SIGNATURE  0x58595A20
    230 #define LAB_SIGNATURE  0x4C616220
    231 
    232 struct _qcms_profile {
    233 	uint32_t class;
    234 	uint32_t color_space;
    235 	uint32_t pcs;
    236 	qcms_intent rendering_intent;
    237 	struct XYZNumber redColorant;
    238 	struct XYZNumber blueColorant;
    239 	struct XYZNumber greenColorant;
    240 	struct curveType *redTRC;
    241 	struct curveType *blueTRC;
    242 	struct curveType *greenTRC;
    243 	struct curveType *grayTRC;
    244 	struct lutType *A2B0;
    245 	struct lutType *B2A0;
    246 	struct lutmABType *mAB;
    247 	struct lutmABType *mBA;
    248 	struct matrix chromaticAdaption;
    249 
    250 	struct precache_output *output_table_r;
    251 	struct precache_output *output_table_g;
    252 	struct precache_output *output_table_b;
    253 };
    254 
    255 #ifdef _MSC_VER
    256 #define inline _inline
    257 #endif
    258 
    259 /* produces the nearest float to 'a' with a maximum error
    260  * of 1/1024 which happens for large values like 0x40000040 */
    261 static inline float s15Fixed16Number_to_float(s15Fixed16Number a)
    262 {
    263 	return ((int32_t)a)/65536.f;
    264 }
    265 
    266 static inline s15Fixed16Number double_to_s15Fixed16Number(double v)
    267 {
    268 	return (int32_t)(v*65536);
    269 }
    270 
    271 static inline float uInt8Number_to_float(uInt8Number a)
    272 {
    273 	return ((int32_t)a)/255.f;
    274 }
    275 
    276 static inline float uInt16Number_to_float(uInt16Number a)
    277 {
    278 	return ((int32_t)a)/65535.f;
    279 }
    280 
    281 
    282 void precache_release(struct precache_output *p);
    283 qcms_bool set_rgb_colorants(qcms_profile *profile, qcms_CIE_xyY white_point, qcms_CIE_xyYTRIPLE primaries);
    284 
    285 void qcms_transform_data_rgb_out_lut_sse2(qcms_transform *transform,
    286                                           unsigned char *src,
    287                                           unsigned char *dest,
    288                                           size_t length,
    289                                           qcms_format_type output_format);
    290 void qcms_transform_data_rgba_out_lut_sse2(qcms_transform *transform,
    291                                           unsigned char *src,
    292                                           unsigned char *dest,
    293                                           size_t length,
    294                                           qcms_format_type output_format);
    295 void qcms_transform_data_rgb_out_lut_sse1(qcms_transform *transform,
    296                                           unsigned char *src,
    297                                           unsigned char *dest,
    298                                           size_t length,
    299                                           qcms_format_type output_format);
    300 void qcms_transform_data_rgba_out_lut_sse1(qcms_transform *transform,
    301                                           unsigned char *src,
    302                                           unsigned char *dest,
    303                                           size_t length,
    304                                           qcms_format_type output_format);
    305 
    306 extern qcms_bool qcms_supports_iccv4;
    307 
    308 
    309 #ifdef _MSC_VER
    310 
    311 long __cdecl _InterlockedIncrement(long volatile *);
    312 long __cdecl _InterlockedDecrement(long volatile *);
    313 #pragma intrinsic(_InterlockedIncrement)
    314 #pragma intrinsic(_InterlockedDecrement)
    315 
    316 #define qcms_atomic_increment(x) _InterlockedIncrement((long volatile *)&x)
    317 #define qcms_atomic_decrement(x) _InterlockedDecrement((long volatile*)&x)
    318 
    319 #else
    320 
    321 #define qcms_atomic_increment(x) __sync_add_and_fetch(&x, 1)
    322 #define qcms_atomic_decrement(x) __sync_sub_and_fetch(&x, 1)
    323 
    324 #endif
    325