From 67bc1de939bdfb2072074bc9554dfea9b32f64cd Mon Sep 17 00:00:00 2001 From: Mattias Andrée Date: Sun, 6 Oct 2019 20:36:00 +0200 Subject: Simplify MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: Mattias Andrée --- libclut.h | 2396 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 2396 insertions(+) create mode 100644 libclut.h (limited to 'libclut.h') diff --git a/libclut.h b/libclut.h new file mode 100644 index 0000000..f318df2 --- /dev/null +++ b/libclut.h @@ -0,0 +1,2396 @@ +/* See LICENSE file for copyright and license details. */ +#ifndef LIBCLUT_H +#define LIBCLUT_H + +#include +#include +#include + +/* Not documented, may be removed or modified in the future. */ +#define LIBCLUT_ILLUMINANT_A .white_x = 0.44757, .white_y = 0.40745, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_B .white_x = 0.34842, .white_y = 0.35161, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_C .white_x = 0.31006, .white_y = 0.31616, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_D50 .white_x = 0.34567, .white_y = 0.35850, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_D55 .white_x = 0.33242, .white_y = 0.34743, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_D65 .white_x = 0.31271, .white_y = 0.32902, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_D75 .white_x = 0.29902, .white_y = 0.31485, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_E .white_x = 1. / 3, .white_y = 1. / 3, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F1 .white_x = 0.31310, .white_y = 0.33727, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F2 .white_x = 0.37208, .white_y = 0.37529, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F3 .white_x = 0.40910, .white_y = 0.39430, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F4 .white_x = 0.44018, .white_y = 0.40329, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F5 .white_x = 0.31379, .white_y = 0.34531, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F6 .white_x = 0.37790, .white_y = 0.38835, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F7 .white_x = 0.31292, .white_y = 0.32933, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F8 .white_x = 0.34588, .white_y = 0.35875, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F9 .white_x = 0.37417, .white_y = 0.37281, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F10 .white_x = 0.34609, .white_y = 0.35986, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F11 .white_x = 0.38052, .white_y = 0.37713, .white_Y = 1 +#define LIBCLUT_ILLUMINANT_F12 .white_x = 0.43695, .white_y = 0.40441, .white_Y = 1 + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the sRGB colour space + * + * sRGB does not use a regular gamma function, but rather two different + * functions. `libclut_model_linear_to_standard1` and + * and `libclut_model_standard_to_linear1` can be used so to convert + * between sRGB and linear sRGB + */ +#define LIBCLUT_RGB_COLOUR_SPACE_SRGB_INITIALISER {\ + .red_x = 0.6400, .red_y = 0.3300, .red_Y = 0.212656,\ + .green_x = 0.3000, .green_y = 0.6000, .green_Y = 0.715158,\ + .blue_x = 0.1500, .blue_y = 0.0600, .blue_Y = 0.072186,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Adobe RGB (1998) colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_ADOBE_RGB_INITIALISER {\ + .red_x = 0.6400, .red_y = 0.3300, .red_Y = 0.297361,\ + .green_x = 0.2100, .green_y = 0.7100, .green_Y = 0.627355,\ + .blue_x = 0.1500, .blue_y = 0.0600, .blue_Y = 0.075285,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Apple RGB colour space + * + * This colour space's gamma is 1.8. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_APPLE_RGB_INITIALISER {\ + .red_x = 0.6250, .red_y = 0.3400, .red_Y = 0.244634,\ + .green_x = 0.2800, .green_y = 0.5950, .green_Y = 0.672034,\ + .blue_x = 0.1550, .blue_y = 0.0700, .blue_Y = 0.083332,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Best RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_BEST_RGB_INITIALISER {\ + .red_x = 0.7347, .red_y = 0.2653, .red_Y = 0.228457,\ + .green_x = 0.2150, .green_y = 0.7750, .green_Y = 0.737352,\ + .blue_x = 0.1300, .blue_y = 0.0350, .blue_Y = 0.034191,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Beta RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_BETA_RGB_INITIALISER {\ + .red_x = 0.6888, .red_y = 0.3112, .red_Y = 0.303273,\ + .green_x = 0.1986, .green_y = 0.7551, .green_Y = 0.663786,\ + .blue_x = 0.1265, .blue_y = 0.0352, .blue_Y = 0.032941,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Bruce RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_BRUCE_RGB_INITIALISER {\ + .red_x = 0.6400, .red_y = 0.3300, .red_Y = 0.240995,\ + .green_x = 0.2800, .green_y = 0.6500, .green_Y = 0.683554,\ + .blue_x = 0.1500, .blue_y = 0.0600, .blue_Y = 0.075452,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the CIE RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_CIE_RGB_INITIALISER {\ + .red_x = 0.7350, .red_y = 0.2650, .red_Y = 0.176204,\ + .green_x = 0.2740, .green_y = 0.7170, .green_Y = 0.812985,\ + .blue_x = 0.1670, .blue_y = 0.0090, .blue_Y = 0.010811,\ + LIBCLUT_ILLUMINANT_E} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ColorMatch RGB colour space + * + * This colour space's gamma is 1.8. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_COLORMATCH_RGB_INITIALISER {\ + .red_x = 0.6300, .red_y = 0.3400, .red_Y = 0.274884,\ + .green_x = 0.2950, .green_y = 0.6050, .green_Y = 0.658132,\ + .blue_x = 0.1500, .blue_y = 0.0750, .blue_Y = 0.066985,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the DCI-P3 D65 colour space + * + * This colour space's gamma is 2.6. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_DCI_P3_D65_INITIALISER {\ + .red_x = 0.680, .red_y = 0.320, .red_Y = 0.22897344,\ + .green_x = 0.265, .green_y = 0.690, .green_Y = 0.69175166,\ + .blue_x = 0.150, .blue_y = 0.060, .blue_Y = 0.07927490,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the DCI-P3 Theater colour space + * + * This colour space's gamma is 2.6. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_DCI_P3_THEATER_INITIALISER {\ + .red_x = 0.680, .red_y = 0.320, .red_Y = 0.20949168,\ + .green_x = 0.265, .green_y = 0.690, .green_Y = 0.72159525,\ + .blue_x = 0.150, .blue_y = 0.060, .blue_Y = 0.06891307,\ + .white_x = 0.314, .white_y = 0.351, .white_Y = 1} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Don RGB 4 colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_DON_RGB_4_INITIALISER {\ + .red_x = 0.6960, .red_y = 0.3000, .red_Y = 0.278350,\ + .green_x = 0.2150, .green_y = 0.7650, .green_Y = 0.687970,\ + .blue_x = 0.1300, .blue_y = 0.0350, .blue_Y = 0.033680,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ECI RGB v2 colour space + * + * This colour space's used the L* gamma function. It cannot + * be used with RGB colour space conversion unless the values are + * tranlated into using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_ECI_RGB_V2_INITIALISER {\ + .red_x = 0.6700, .red_y = 0.3300, .red_Y = 0.320250,\ + .green_x = 0.2100, .green_y = 0.7100, .green_Y = 0.602071,\ + .blue_x = 0.1400, .blue_y = 0.0800, .blue_Y = 0.077679,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Ekta Space PS5 colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_EKTA_SPACE_PS5_INITIALISER {\ + .red_x = 0.6950, .red_y = 0.3050, .red_Y = 0.260629,\ + .green_x = 0.2600, .green_y = 0.7000, .green_Y = 0.734946,\ + .blue_x = 0.1100, .blue_y = 0.0050, .blue_Y = 0.004425,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ITU-R Recommendation BT.601 (ITU-R BT.601), 625 line colour + * space + * + * This colour space uses a custom gamma function. It cannot be used + * with RGB colour space conversion unless the values are tranlated + * into using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_ITU_R_BT_601_625_LINE_INITIALISER {\ + .red_x = 0.640, .red_y = 0.330, .red_Y = 0.2220023,\ + .green_x = 0.290, .green_y = 0.600, .green_Y = 0.7066689,\ + .blue_x = 0.150, .blue_y = 0.060, .blue_Y = 0.0713288,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ITU-R Recommendation BT.601 (ITU-R BT.601), 525 line colour + * space + * + * This colour space uses a custom gamma function. It cannot be used + * with RGB colour space conversion unless the values are tranlated + * into using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_ITU_R_BT_601_525_LINE_INITIALISER {\ + .red_x = 0.630, .red_y = 0.340, .red_Y = 0.2220023,\ + .green_x = 0.310, .green_y = 0.595, .green_Y = 0.7066689,\ + .blue_x = 0.155, .blue_y = 0.070, .blue_Y = 0.0713288,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ITU-R Recommendation BT.709 (ITU-R BT.709) colour space + * + * This colour space's gamma is 2._. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_ITU_R_BT_709_INITIALISER {\ + .red_x = 0.6400, .red_y = 0.3300, .red_Y = 0.212656,\ + .green_x = 0.3000, .green_y = 0.6000, .green_Y = 0.715158,\ + .blue_x = 0.1500, .blue_y = 0.0600, .blue_Y = 0.072186,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ITU-R Recommendation BT.2020 (ITU-R BT.2020) colour space + * + * This colour space uses a custom gamma function. It cannot be used + * with RGB colour space conversion unless the values are tranlated + * into using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_ITU_R_BT_2020_INITIALISER {\ + .red_x = 0.7080, .red_y = 0.2920, .red_Y = 0.2627296,\ + .green_x = 0.1700, .green_y = 0.7970, .green_Y = 0.6767483,\ + .blue_x = 0.1310, .blue_y = 0.0460, .blue_Y = 0.0605221,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ITU-R Recommendation BT.2100 (ITU-R BT.2100) colour space + * + * This colour space uses a custom gamma function. It cannot be used + * with RGB colour space conversion unless the values are tranlated + * into using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_ITU_R_BT_2100_INITIALISER {\ + .red_x = 0.7080, .red_y = 0.2920, .red_Y = 0.2627296,\ + .green_x = 0.1700, .green_y = 0.7970, .green_Y = 0.6767483,\ + .blue_x = 0.1310, .blue_y = 0.0460, .blue_Y = 0.0605221,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Lightroom RGB colour space + * + * This colour space's gamma is 1 (linear). It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_LIGHTROOM_RGB_INITIALISER {\ + .red_x = 0.7347, .red_y = 0.2653, .red_Y = 0.288040,\ + .green_x = 0.1596, .green_y = 0.8404, .green_Y = 0.711874,\ + .blue_x = 0.0366, .blue_y = 0.0001, .blue_Y = 0.000086,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the NTSC RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_NTSC_RGB_INITIALISER {\ + .red_x = 0.6700, .red_y = 0.3300, .red_Y = 0.298839,\ + .green_x = 0.2100, .green_y = 0.7100, .green_Y = 0.586811,\ + .blue_x = 0.1400, .blue_y = 0.0800, .blue_Y = 0.114350,\ + LIBCLUT_ILLUMINANT_C} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the PAL/SECAM RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_PAL_SECAM_RGB_INITIALISER {\ + .red_x = 0.6400, .red_y = 0.3300, .red_Y = 0.222021,\ + .green_x = 0.2900, .green_y = 0.6000, .green_Y = 0.706645,\ + .blue_x = 0.1500, .blue_y = 0.0600, .blue_Y = 0.071334,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the ProPhoto RGB colour space + * + * This colour space's gamma is 1.8. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_PROPHOTO_RGB_INITIALISER {\ + .red_x = 0.7347, .red_y = 0.2653, .red_Y = 0.288040,\ + .green_x = 0.1596, .green_y = 0.8404, .green_Y = 0.711874,\ + .blue_x = 0.0366, .blue_y = 0.0001, .blue_Y = 0.000086,\ + LIBCLUT_ILLUMINANT_D50} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the SMPTE-C RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_SMPTE_C_RGB_INITIALISER {\ + .red_x = 0.6300, .red_y = 0.3400, .red_Y = 0.212395,\ + .green_x = 0.3100, .green_y = 0.5950, .green_Y = 0.701049,\ + .blue_x = 0.1550, .blue_y = 0.0700, .blue_Y = 0.086556,\ + LIBCLUT_ILLUMINANT_D65} + +/** + * Initialiser for `struct libclut_rgb_colour_space` with the values + * of the Wide Gamut RGB colour space + * + * This colour space's gamma is 2.2. It cannot be used with + * RGB colour space conversion unless the values are tranlated into + * using the sRGB gamma function. + */ +#define LIBCLUT_RGB_COLOUR_SPACE_WIDE_GAMUT_RGB_INITIALISER {\ + .red_x = 0.7350, .red_y = 0.2650, .red_Y = 0.258187,\ + .green_x = 0.1150, .green_y = 0.8260, .green_Y = 0.724938,\ + .blue_x = 0.1570, .blue_y = 0.0180, .blue_Y = 0.016875,\ + LIBCLUT_ILLUMINANT_D50} + +/* + * TODO gamma functions: + * https://en.wikipedia.org/wiki/Rec._2020#Transfer_characteristics + * http://www.itu.int/dms_pubrec/itu-r/rec/bt/R-REC-BT.2100-0-201607-I!!PDF-E.pdf + * https://en.wikipedia.org/wiki/Rec._601 + */ + +/** + * RGB colour space structure + */ +typedef struct libclut_rgb_colour_space { + /** + * The x-component of the red colour's xyY value + */ + double red_x; + + /** + * The y-component of the red colour's xyY value + */ + double red_y; + + /** + * The Y-component of the red colour's xyY value + */ + double red_Y; + + /** + * The x-component of the green colour's xyY value + */ + double green_x; + + /** + * The y-component of the green colour's xyY value + */ + double green_y; + + /** + * The Y-component of the green colour's xyY value + */ + double green_Y; + + /** + * The x-component of the blue colour's xyY value + */ + double blue_x; + + /** + * The y-component of the blue colour's xyY value + */ + double blue_y; + + /** + * The Y-component of the blue colour's xyY value + */ + double blue_Y; + + /** + * The x-component of the white point's xyY value + */ + double white_x; + + /** + * The y-component of the white point's xyY value + */ + double white_y; + + /** + * The Y-component of the white point's xyY value + */ + double white_Y; +} libclut_rgb_colour_space_t; + +/** + * Matrix date-type for colour space conversion + */ +typedef double libclut_colour_space_conversion_matrix_t[3][3]; + +/* This is to avoid warnings about comparing double, These are only + * used when it is safe, for example to test whether optimisations + * are possible. { */ +#if defined(__GNUC__) || defined(__clang__) +# pragma GCC diagnostic push +# pragma GCC diagnostic ignored "-Wfloat-equal" +#endif +static inline int libclut_eq__(double a, double b) { return a == b; } +static inline int libclut_1__(double x) { return libclut_eq__(x, 1); } +static inline int libclut_0__(double x) { return libclut_eq__(x, 0); } +#if defined(__GNUC__) || defined(__clang__) +# pragma GCC diagnostic pop +#endif +/* } */ + +#if defined(__clang__) +# pragma GCC diagnostic push +# pragma GCC diagnostic ignored "-Wdocumentation" +#endif + +/** + * Apply contrast correction on the colour curves using sRGB + * + * In this context, contrast is a measure of difference between + * the whitepoint and blackpoint, if the difference is 0 than + * they are both grey + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r The contrast parameter for the red curve + * @param g The contrast parameter for the green curve + * @param b The contrast parameter for the blue curve + */ +#define libclut_rgb_contrast(clut, max, type, r, g, b)\ + do {\ + const double h__ = (double)5 / 10;\ + if (!libclut_1__(r)) libclut__(clut, red, type, (LIBCLUT_VALUE - (max) * h__) * (r) + (max) * h__);\ + if (!libclut_1__(g)) libclut__(clut, green, type, (LIBCLUT_VALUE - (max) * h__) * (g) + (max) * h__);\ + if (!libclut_1__(b)) libclut__(clut, blue, type, (LIBCLUT_VALUE - (max) * h__) * (b) + (max) * h__);\ + } while (0) + +/** + * Apply contrast correction on the colour curves using CIE xyY + * + * In this context, contrast is a measure of difference between + * the whitepoint and blackpoint, if the difference is 0 than + * they are both grey + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r The contrast parameter for the red curve + * @param g The contrast parameter for the green curve + * @param b The contrast parameter for the blue curve + */ +#define libclut_cie_contrast(clut, max, type, r, g, b)\ + do {\ + const double h__ = (double)5 / 10;\ + libclut_cie__(clut, max, type, libclut_eq__((r), (g)) && libclut_eq__((g), (b)),\ + !libclut_1__(r), !libclut_1__(g), !libclut_1__(b),\ + (Y__ - h__) * (r) + h__, (Y__ - h__) * (g) + h__, (Y__ - h__) * (b) + h__);\ + } while (0) + +/** + * Apply brightness correction on the colour curves using sRGB + * + * In this context, brightness is a measure of the whiteness of the whitepoint + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r The brightness parameter for the red curve + * @param g The brightness parameter for the green curve + * @param b The brightness parameter for the blue curve + */ +#define libclut_rgb_brightness(clut, max, type, r, g, b)\ + do {\ + if (!libclut_1__(r)) libclut__(clut, red, type, LIBCLUT_VALUE * (r));\ + if (!libclut_1__(g)) libclut__(clut, green, type, LIBCLUT_VALUE * (g));\ + if (!libclut_1__(b)) libclut__(clut, blue, type, LIBCLUT_VALUE * (b));\ + } while (0) + +/** + * Apply brightness correction on the colour curves using CIE xyY + * + * In this context, brightness is a measure of the whiteness of the whitepoint + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r The brightness parameter for the red curve + * @param g The brightness parameter for the green curve + * @param b The brightness parameter for the blue curve + */ +#define libclut_cie_brightness(clut, max, type, r, g, b)\ + libclut_cie__(clut, max, type, libclut_eq__((r), (g)) && libclut_eq__((g), (b)),\ + !libclut_1__(r), !libclut_1__(g), !libclut_1__(b),\ + Y__ * (r), Y__ * (g), Y__ * (b)) + +/** + * Convert the curves from formatted in standard RGB to linear sRGB + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_standard_to_linear1` is not undefined + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Whether to convert the red colour curve + * @param g Whether to convert the green colour curve + * @param b Whether to convert the blue colour curve + */ +#define libclut_linearise(clut, max, type, r, g, b)\ + do {\ + double m__ = (double)(max);\ + if (r) libclut__(clut, red, type, m__ * libclut_model_standard_to_linear1(LIBCLUT_VALUE / m__));\ + if (g) libclut__(clut, green, type, m__ * libclut_model_standard_to_linear1(LIBCLUT_VALUE / m__));\ + if (b) libclut__(clut, blue, type, m__ * libclut_model_standard_to_linear1(LIBCLUT_VALUE / m__));\ + } while (0) + +/** + * Convert the curves from formatted in linear sRGB to standard RGB + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_linear_to_standard1` is not undefined + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Whether to convert the red colour curve + * @param g Whether to convert the green colour curve + * @param b Whether to convert the blue colour curve + */ +#define libclut_standardise(clut, max, type, r, g, b)\ + do {\ + double m__ = (double)(max);\ + if (r) libclut__(clut, red, type, m__ * libclut_model_linear_to_standard1(LIBCLUT_VALUE / m__));\ + if (g) libclut__(clut, green, type, m__ * libclut_model_linear_to_standard1(LIBCLUT_VALUE / m__));\ + if (b) libclut__(clut, blue, type, m__ * libclut_model_linear_to_standard1(LIBCLUT_VALUE / m__));\ + } while (0) + +/** + * Convert the curves between two RGB colour spaces + * + * Both RGB colour spaces must have same gamma functions as sRGB + * + * Requires that `clut->red_size`, `clut->green_size` + * and `clut->blue_size` are equal + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_linear_to_standard1`, + * `libclut_model_standard_to_linear1`, or + * `libclut_model_convert_rgb` is not undefined + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param m Conversion matrix. Can be created with + * `libclut_model_get_rgb_conversion_matrix` + * @param trunc Truncate values that are out of gamut + */ +#define libclut_convert_rgb_inplace(clut, max, type, m, trunc)\ + do {\ + double m__ = (double)(max), r__, g__, b__;\ + size_t i__, n__ = (clut)->red_size;\ + for (i__ = 0; i__ < n__; i__++) {\ + r__ = (clut)->red[i__] / m__;\ + g__ = (clut)->green[i__] / m__;\ + b__ = (clut)->blue[i__] / m__;\ + libclut_model_convert_rgb(r__, g__, b__, m, &r__, &g__, &b__);\ + r__ *= m__;\ + g__ *= m__;\ + b__ *= m__;\ + if (trunc) {\ + if (r__ < 0)\ + r__ = 0;\ + else if (r__ > m__)\ + r__ = m__;\ + if (g__ < 0)\ + g__ = 0;\ + else if (g__ > m__)\ + g__ = m__;\ + if (b__ < 0)\ + b__ = 0;\ + else if (b__ > m__)\ + b__ = m__;\ + }\ + (clut)->red[i__] = (type)r__;\ + (clut)->green[i__] = (type)g__;\ + (clut)->blue[i__] = (type)b__;\ + }\ + } while (0) + +/** + * Convert the curves between two RGB colour spaces + * + * Both RGB colour spaces must have same gamma functions as sRGB + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' if + * `libclut_model_linear_to_standard1`, + * `libclut_model_standard_to_linear1`, or + * `libclut_model_convert_rgb` is not undefined. + * Always requires linking with '-lm'. + * + * @param clut Pointer to the input gamma ramps, must have the + * arrays `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param m Conversion matrix. Can be created with + * `libclut_model_get_rgb_conversion_matrix` + * @param trunc Truncate values that are out of gamut + * @param out Pointer to the output gamma ramps, must have the + * arrays `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + */ +#define libclut_convert_rgb(clut, max, type, m, trunc, out)\ + do {\ + double m__ = (double)(max), r__, g__, b__, x__, y__;\ + size_t rn__ = (clut)->red_size;\ + size_t gn__ = (clut)->green_size;\ + size_t bn__ = (clut)->blue_size;\ + double w__;\ + size_t i__, j__, jj__;\ + for (i__ = 0; i__ < rn__; i__++) {\ + w__ = (double)i__ * gn__ / rn__;\ + j__ = (size_t)w__;\ + jj__ = j__ == gn__ ? j__ : (j__ + 1);\ + w__ = fmod(j__, (double)1);\ + x__ = (clut)->green[j__] / m__;\ + y__ = (clut)->green[jj__] / m__;\ + x__ = libclut_model_standard_to_linear1(x__);\ + y__ = libclut_model_standard_to_linear1(y__);\ + g__ = x__ * (1 - w__) + y__ * w__;\ + \ + w__ = (double)i__ * bn__ / rn__;\ + j__ = (size_t)w__;\ + jj__ = j__ == bn__ ? j__ : (j__ + 1);\ + w__ = fmod(j__, (double)1);\ + x__ = (clut)->blue[j__] / m__;\ + y__ = (clut)->blue[jj__] / m__;\ + x__ = libclut_model_standard_to_linear1(x__);\ + y__ = libclut_model_standard_to_linear1(y__);\ + b__ = x__ * (1 - w__) + y__ * w__;\ + \ + r__ = (M)[0][0] * r__ + (M)[0][1] * g__ + (M)[0][2] * b__;\ + r__ = libclut_model_linear_to_standard1(r__);\ + r__ *= m__;\ + if (trunc) {\ + if (r__ < 0)\ + r__ = 0;\ + else if (r__ > m__)\ + r__ = m__;\ + }\ + (out)->red[i__] = (type)r__;\ + }\ + for (i__ = 0; i__ < gn__; i__++) {\ + w__ = (double)i__ * rn__ / gn__;\ + j__ = (size_t)w__;\ + jj__ = j__ == rn__ ? j__ : (j__ + 1);\ + w__ = fmod(j__, (double)1);\ + x__ = (clut)->red[j__] / m__;\ + y__ = (clut)->red[jj__] / m__;\ + x__ = libclut_model_standard_to_linear1(x__);\ + y__ = libclut_model_standard_to_linear1(y__);\ + r__ = x__ * (1 - w__) + y__ * w__;\ + \ + w__ = (double)i__ * bn__ / gn__;\ + j__ = (size_t)w__;\ + jj__ = j__ == bn__ ? j__ : (j__ + 1);\ + w__ = fmod(j__, (double)1);\ + x__ = (clut)->blue[j__] / m__;\ + y__ = (clut)->blue[jj__] / m__;\ + x__ = libclut_model_standard_to_linear1(x__);\ + y__ = libclut_model_standard_to_linear1(y__);\ + b__ = x__ * (1 - w__) + y__ * w__;\ + \ + g__ = (M)[1][0] * r__ + (M)[0][1] * g__ + (M)[1][2] * b__;\ + g__ = libclut_model_linear_to_standard1(g__);\ + g__ *= m__;\ + if (trunc) {\ + if (g__ < 0)\ + g__ = 0;\ + else if (g__ > m__)\ + g__ = m__;\ + }\ + (out)->green[i__] = (type)g__;\ + }\ + for (i__ = 0; i__ < bn__; i__++) {\ + w__ = (double)i__ * rn__ / bn__;\ + j__ = (size_t)w__;\ + jj__ = j__ == rn__ ? j__ : (j__ + 1);\ + w__ = fmod(j__, (double)1);\ + x__ = (clut)->red[j__] / m__;\ + y__ = (clut)->red[jj__] / m__;\ + x__ = libclut_model_standard_to_linear1(x__);\ + y__ = libclut_model_standard_to_linear1(y__);\ + r__ = x__ * (1 - w__) + y__ * w__;\ + \ + w__ = (double)i__ * gn__ / bn__;\ + j__ = (size_t)w__;\ + jj__ = j__ == gn__ ? j__ : (j__ + 1);\ + w__ = fmod(j__, (double)1);\ + x__ = (clut)->green[j__] / m__;\ + y__ = (clut)->green[jj__] / m__;\ + x__ = libclut_model_standard_to_linear1(x__);\ + y__ = libclut_model_standard_to_linear1(y__);\ + g__ = x__ * (1 - w__) + y__ * w__;\ + \ + b__ = (M)[2][0] * r__ + (M)[2][1] * g__ + (M)[2][2] * b__;\ + b__ = libclut_model_linear_to_standard1(b__);\ + b__ *= m__;\ + if (trunc) {\ + if (b__ < 0)\ + b__ = 0;\ + else if (b__ > m__)\ + b__ = m__;\ + }\ + (out)->blue[i__] = (type)b__;\ + }\ + } while (0) + +/** + * Apply gamma correction on the colour curves + * + * None of the parameter may have side-effects + * + * Requires linking with '-lm' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r The gamma parameter the red colour curve + * @param g The gamma parameter the green colour curve + * @param b The gamma parameter the blue colour curve + */ +#define libclut_gamma(clut, max, type, r, g, b)\ + do {\ + double m__ = (double)(max);\ + if (!libclut_1__(r)) libclut__(clut, red, type, m__ * pow(LIBCLUT_VALUE / m__, 1 / (double)(r)));\ + if (!libclut_1__(g)) libclut__(clut, green, type, m__ * pow(LIBCLUT_VALUE / m__, 1 / (double)(g)));\ + if (!libclut_1__(b)) libclut__(clut, blue, type, m__ * pow(LIBCLUT_VALUE / m__, 1 / (double)(b)));\ + } while (0) + +/** + * Reverse the colour curves (negative image with gamma preservation) + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps. + * This parameter is not used, it is just a dummy, to unify + * the API with the other functions + * @param type The data type used for each stop in the ramps + * @param r Whether to invert the red colour curve + * @param g Whether to invert the green colour curve + * @param b Whether to invert the blue colour curve + */ +#define libclut_negative(clut, max, type, r, g, b)\ + do {\ + size_t i__, n__;\ + type t__;\ + if (r) {\ + for (i__ = 0, n__ = (clut)->red_size; i__ < (n__ >> 1); i__++) {\ + t__ = (clut)->red[i__];\ + (clut)->red[i__] = (clut)->red[n__ - i__ - 1];\ + (clut)->red[n__ - i__ - 1] = t__;\ + }\ + }\ + if (g) {\ + for (i__ = 0, n__ = (clut)->green_size; i__ < (n__ >> 1); i__++) {\ + t__ = (clut)->green[i__];\ + (clut)->green[i__] = (clut)->green[n__ - i__ - 1];\ + (clut)->green[n__ - i__ - 1] = t__;\ + }\ + }\ + if (b) {\ + for (i__ = 0, n__ = (clut)->blue_size; i__ < (n__ >> 1); i__++) {\ + t__ = (clut)->blue[i__];\ + (clut)->blue[i__] = (clut)->blue[n__ - i__ - 1];\ + (clut)->blue[n__ - i__ - 1] = t__;\ + }\ + }\ + } while (0) + +/** + * Invert the colour curves (negative image with gamma invertion), using sRGB + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Whether to invert the red colour curve + * @param g Whether to invert the green colour curve + * @param b Whether to invert the blue colour curve + */ +#define libclut_rgb_invert(clut, max, type, r, g, b)\ + do {\ + if (r) libclut__(clut, red, type, (max) - LIBCLUT_VALUE);\ + if (g) libclut__(clut, green, type, (max) - LIBCLUT_VALUE);\ + if (b) libclut__(clut, blue, type, (max) - LIBCLUT_VALUE);\ + } while (0) + +/** + * Invert the colour curves (negative image with gamma invertion), using CIE xyY + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Whether to invert the red colour curve + * @param g Whether to invert the green colour curve + * @param b Whether to invert the blue colour curve + */ +#define libclut_cie_invert(clut, max, type, r, g, b)\ + libclut_cie__(clut, max, type, (r) && (g) && (b), r, g, b, 1 - Y__, 1 - Y__, 1 - Y__) + +/** + * Apply S-curve correction on the colour curves. + * This is intended for fine tuning LCD monitors, + * 4.5 is good value start start testing at. + * You would probably like to use rgb_limits before + * this to adjust the blackpoint as that is the + * only way to adjust the blackpoint on many LCD + * monitors. + * + * None of the parameter may have side-effects + * + * Requires linking with '-lm' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param rp Pointer to the sigmoid parameter for the red curve. `NULL` for no adjustment + * @param gp Pointer to the sigmoid parameter for the green curve. `NULL` for no adjustment + * @param bp Pointer to the sigmoid parameter for the blue curve. `NULL` for no adjustment + */ +#define libclut_sigmoid(clut, max, type, rp, gp, bp)\ + do {\ + double *gcc_6_1_1_workaround, m__ = (double)(max);\ + const double h__ = (double)5 / 10;\ + gcc_6_1_1_workaround = rp;\ + if (gcc_6_1_1_workaround)\ + libclut_sigmoid__(clut, max, type, red);\ + gcc_6_1_1_workaround = gp;\ + if (gcc_6_1_1_workaround)\ + libclut_sigmoid__(clut, max, type, green);\ + gcc_6_1_1_workaround = bp;\ + if (gcc_6_1_1_workaround)\ + libclut_sigmoid__(clut, max, type, blue);\ + }\ + while (0) + +/** + * Apply S-curve correction on the colour curves. + * This is intended for fine tuning LCD monitors, + * 4.5 is good value start start testing at. + * You would probably like to use rgb_limits before + * this to adjust the blackpoint as that is the + * only way to adjust the blackpoint on many LCD + * monitors. + * + * None of the parameter may have side-effects + * + * Requires linking with '-lm' + * + * Intended for internal use + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param channel The channel, must be either "red", "green", or "blue" + */ +#define libclut_sigmoid__(clut, max, type, channel)\ + do {\ + double s__ = *gcc_6_1_1_workaround, l__;\ + size_t i__;\ + for (i__ = 0; i__ < (clut)->channel##_size; i__++) {\ + l__ = log(m__ / (clut)->channel[i__] - 1);\ + if (isnan(l__) || isinf(l__))\ + l__ = 37.024483 * (isinf(l__) > 0 ? +1 : -1);\ + (clut)->channel[i__] = (type)(m__ * (h__ - l__ / s__));\ + }\ + } while (0) + +/** + * Changes the blackpoint and the whitepoint, using sRGB + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param rmin The red component value of the blackpoint + * @param rmax The red component value of the whitepoint + * @param gmin The green component value of the blackpoint + * @param gmax The green component value of the whitepoint + * @param bmin The blue component value of the blackpoint + * @param bmax The blue component value of the whitepoint + */ +#define libclut_rgb_limits(clut, max, type, rmin, rmax, gmin, gmax, bmin, bmax)\ + do {\ + double diff__;\ + if (!libclut_0__(rmin) || !libclut_1__(rmax)) {\ + diff__ = (double)(rmax) - (double)(rmin);\ + libclut__(clut, red, type, LIBCLUT_VALUE / (double)(max) * diff__ + (rmin));\ + }\ + if (!libclut_0__(gmin) || !libclut_1__(gmax)) {\ + diff__ = (double)(gmax) - (double)(gmin);\ + libclut__(clut, green, type, LIBCLUT_VALUE / (double)(max) * diff__ + (gmin));\ + }\ + if (!libclut_0__(bmin) || !libclut_1__(bmax)) {\ + diff__ = (double)(bmax) - (double)(bmin);\ + libclut__(clut, blue, type, LIBCLUT_VALUE / (double)(max) * diff__ + (bmin));\ + }\ + } while (0) + +/** + * Changes the blackpoint and the whitepoint, using CIE xyY + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param rmin The red component value of the blackpoint + * @param rmax The red component value of the whitepoint + * @param gmin The green component value of the blackpoint + * @param gmax The green component value of the whitepoint + * @param bmin The blue component value of the blackpoint + * @param bmax The blue component value of the whitepoint + */ +#define libclut_cie_limits(clut, max, type, rmin, rmax, gmin, gmax, bmin, bmax)\ + do {\ + double rd__ = (rmax) - (rmin), gd__ = (gmax) - (gmin), bd__ = (bmax) - (bmin);\ + libclut_cie__(clut, max, type,\ + libclut_eq__((rmin), (gmin)) && libclut_eq__((gmin), (bmin)) &&\ + libclut_eq__((rmax), (gmax)) && libclut_eq__((gmax), (bmax)),\ + !libclut_0__(rmin) || !libclut_1__(rmax),\ + !libclut_0__(gmin) || !libclut_1__(gmax),\ + !libclut_0__(bmin) || !libclut_1__(bmax),\ + Y__ * rd__ + (rmin), Y__ * gd__ + (gmin), Y__ * bd__ + (bmin));\ + } while (0) + +/** + * Manipulate the colour curves using a function on the sRGB colour space + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Function to manipulate the red colour curve, should either + * be `NULL` or map a [0, 1] `double` to a [0, 1] `double` + * @param g Function to manipulate the green colour curve, should either + * be `NULL` or map a [0, 1] `double` to a [0, 1] `double` + * @param b Function to manipulate the blue colour curve, should either + * be `NULL` or map a [0, 1] `double` to a [0, 1] `double` + */ +#define libclut_manipulate(clut, max, type, r, g, b)\ + do {\ + double m__ = (double)(max);\ + double (*gcc_6_1_1_workaround__)(double);\ + gcc_6_1_1_workaround__ = r;\ + if (gcc_6_1_1_workaround__)\ + libclut__(clut, red, type, m__ * (gcc_6_1_1_workaround__)(LIBCLUT_VALUE / m__));\ + gcc_6_1_1_workaround__ = g;\ + if (gcc_6_1_1_workaround__)\ + libclut__(clut, green, type, m__ * (gcc_6_1_1_workaround__)(LIBCLUT_VALUE / m__));\ + gcc_6_1_1_workaround__ = b;\ + if (gcc_6_1_1_workaround__)\ + libclut__(clut, blue, type, m__ * (gcc_6_1_1_workaround__)(LIBCLUT_VALUE / m__));\ + } while (0) + +/** + * Manipulate the colour curves using a function on the CIE xyY colour space + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Function to manipulate the red colour curve, should either + * be `NULL` or map a [0, 1] `double` to a [0, 1] `double` + * @param g Function to manipulate the green colour curve, should either + * be `NULL` or map a [0, 1] `double` to a [0, 1] `double` + * @param b Function to manipulate the blue colour curve, should either + * be `NULL` or map a [0, 1] `double` to a [0, 1] `double` + */ +#define libclut_cie_manipulate(clut, max, type, r, g, b)\ + libclut_cie__(clut, max, type, (r) && (g) && (b), r, g, b, (r)(Y__), (g)(Y__), (b)(Y__)) + +/** + * Resets colour curvers to linear mappings + * (Identity mapping if imaginged to map from [0, 1] to [0, 1]) + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Whether to reset the red colour curve + * @param g Whether to reset the green colour curve + * @param b Whether to reset the blue colour curve + */ +#define libclut_start_over(clut, max, type, r, g, b)\ + do {\ + size_t i__;\ + double m__, max__ = (double)(max);\ + if (r) {\ + m__ = (double)((clut)->red_size - 1);\ + for (i__ = 0; i__ < (clut)->red_size; i__++)\ + (clut)->red[i__] = (type)(((double)i__ / m__) * max__);\ + }\ + if (g) {\ + m__ = (double)((clut)->green_size - 1);\ + for (i__ = 0; i__ < (clut)->green_size; i__++)\ + (clut)->green[i__] = (type)(((double)i__ / m__) * max__);\ + }\ + if (b) {\ + m__ = (double)((clut)->blue_size - 1);\ + for (i__ = 0; i__ < (clut)->blue_size; i__++)\ + (clut)->blue[i__] = (type)(((double)i__ / m__) * max__);\ + }\ + } while (0) + +/** + * Clip colour curves to only map to values between the minimum and maximum. + * This should be done, before apply the curves, and before applying changes + * with limited domain. + * + * Values below 0 are set to 0, and values above `max` are set to `max` + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param r Whether to clip the red colour curve + * @param g Whether to clip the green colour curve + * @param b Whether to clip the blue colour curve + */ +#define libclut_clip(clut, max, type, r, g, b)\ + do {\ + if (r) libclut__(clut, red, type, libclut_clip__(0, LIBCLUT_VALUE, max));\ + if (g) libclut__(clut, green, type, libclut_clip__(0, LIBCLUT_VALUE, max));\ + if (b) libclut__(clut, blue, type, libclut_clip__(0, LIBCLUT_VALUE, max));\ + } while (0) + +/** + * Truncates a value to fit a boundary + * + * None of the parameter may have side-effects + * + * Intended for internal use + * + * @param min The minimum allowed value + * @param val The current value + * @param max The maximum allowed value + * @return The value truncated into its boundary + */ +#define libclut_clip__(min, val, max)\ + (LIBCLUT_VALUE < (min) ? (min) : LIBCLUT_VALUE > (max) ? (max) : LIBCLUT_VALUE) + +/** + * Emulates low colour resolution + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param rx The desired emulated red encoding resolution, 0 for unchanged + * @param ry The desired emulated red output resolution, 0 for unchanged + * @param gx The desired emulated green encoding resolution, 0 for unchanged + * @param gy The desired emulated green output resolution, 0 for unchanged + * @param bx The desired emulated blue encoding resolution, 0 for unchanged + * @param by The desired emulated blue output resolution, 0 for unchanged + */ +#define libclut_lower_resolution(clut, max, type, rx, ry, gx, gy, bx, by)\ + do {\ + libclut_lower_resolution__(clut, red, max, type, rx, ry);\ + libclut_lower_resolution__(clut, green, max, type, gx, gy);\ + libclut_lower_resolution__(clut, blue, max, type, bx, by);\ + } while (0) + +/** + * Emulates low colour resolution of a channel + * + * None of the parameter may have side-effects + * + * Intended for internal use + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param channel The channel, must be either "red", "green", or "blue" + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param x The desired emulated encoding resolution, 0 for unchanged + * @param y The desired emulated output resolution, 0 for unchanged + */ +#define libclut_lower_resolution__(clut, channel, max, type, x, y)\ + do {\ + if ((x) || (y)) {\ + size_t x__, y__, i__, n__ = (clut)->channel##_size;\ + double xm__ = (double)((x) - 1), ym__ = (double)((y) - 1);\ + double m__ = (double)(max), nm__ = (double)(n__ - 1);\ + type c__[n__]; /* Do not use alloca! */\ + const double h__ = (double)5 / 10;\ + for (i__ = 0; i__ < n__; i__++) {\ + if ((x__ = i__), (x)) {\ + x__ = (size_t)((double)i__ * (double)(x) / (double)n__);\ + x__ = (size_t)((double)x__ * nm__ / xm__);\ + }\ + if (!(y)) {\ + c__[i__] = (clut)->channel[x__];\ + } else {\ + y__ = (size_t)((double)((clut)->channel[x__]) / (max) * ym__ + h__);\ + c__[i__] = (type)((double)y__ / ym__ * m__);\ + }\ + }\ + memcpy((clut)->channel, c__, n__ * sizeof(type));\ + }\ + } while (0) + +/** + * Translates a gamma ramp structure to another gamma ramp structure type + * + * None of the parameter may have side-effects + * + * @param dclut Pointer to the desired gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars `red_size`, + * `green_size`, and `blue_size`. Ramp structures from + * libgamma or libcoopgamma can be used. + * @param dmax The maximum value on each stop in the ramps in `dclut` + * @param dtype The data type used for each stop in the ramps in `dclut` + * @param sclut Pointer to the set gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars `red_size`, + * `green_size`, and `blue_size`. Ramp structures from + * libgamma or libcoopgamma can be used + * @param smax The maximum value on each stop in the ramps in `sclut` + * @param stype The data type used for each stop in the ramps in `sclut` + * (Not actually used) + */ +#define libclut_translate(dclut, dmax, dtype, sclut, smax, stype)\ + do {\ + libclut_translate__(dclut, dmax, dtype, sclut, smax, stype, red);\ + libclut_translate__(dclut, dmax, dtype, sclut, smax, stype, green);\ + libclut_translate__(dclut, dmax, dtype, sclut, smax, stype, blue);\ + } while (0) + +/** + * Translates a gamma ramp structure to another gamma ramp structure type + * + * None of the parameter may have side-effects + * + * This is intended for internal use + * + * @param dclut Pointer to the desired gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars `red_size`, + * `green_size`, and `blue_size`. Ramp structures from + * libgamma or libcoopgamma can be used. + * @param dmax The maximum value on each stop in the ramps in `dclut` + * @param dtype The data type used for each stop in the ramps in `dclut` + * @param sclut Pointer to the set gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars `red_size`, + * `green_size`, and `blue_size`. Ramp structures from + * libgamma or libcoopgamma can be used. + * @param smax The maximum value on each stop in the ramps in `sclut` + * @param stype The data type used for each stop in the ramps in `sclut`. + * (Not actually used) + * @param channel The channel, must be either "red", "green", or "blue" + */ +#define libclut_translate__(dclut, dmax, dtype, sclut, smax, stype, channel)\ + do {\ + size_t di__, si__, sj__;\ + size_t dn__ = (dclut)->channel##_size;\ + size_t sn__ = (sclut)->channel##_size;\ + double dm__ = (double)(dmax);\ + double sm__ = (double)(smax);\ + double dmsm__ = dm__ / sm__;\ + double x__, y__;\ + if (dn__ == sn__) {\ + for (di__ = 0; di__ < dn__; di__++) {\ + y__ = (double)((sclut)->channel[di__]) * dmsm__;\ + (dclut)->channel[di__] = (dtype)y__;\ + }\ + } else {\ + for (di__ = 0; di__ < dn__; di__++) {\ + x__ = di__ / (dn__ - 1) * (sn__ - 1);\ + si__ = (size_t)(x__);\ + sj__ = si__ + (si__ != sn__);\ + x__ -= (double)si__;\ + y__ = (double)((sclut)->channel[si__]) * (1 - x__);\ + y__ += (double)((sclut)->channel[sj__]) * (x__);\ + y__ *= dmsm__;\ + (dclut)->channel[di__] = (dtype)y__;\ + }\ + }\ + } while (0) + +/** + * Applies a filter or calibration + * + * None of the parameter may have side-effects + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param filter Same as `clut`, but for the filter to apply + * @param fmax Same as `max`, but for the filter to apply + * @param ftype Same as `type`, but for the filter to apply (Not actually used) + * @param r Whether to apply the filter for the red curve + * @param g Whether to apply the filter for the green curve + * @param b Whether to apply the filter for the blue curve + */ +#define libclut_apply(clut, max, type, filter, fmax, ftype, r, g, b)\ + do {\ + if (r) libclut_apply__(clut, max, type, filter, fmax, ftype, red);\ + if (g) libclut_apply__(clut, max, type, filter, fmax, ftype, green);\ + if (b) libclut_apply__(clut, max, type, filter, fmax, ftype, blue);\ + } while (0) + +/** + * Applies a filter or calibration for one channel + * + * None of the parameter may have side-effects + * + * Intended for internal use + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param filter Same as `clut`, but for the filter to apply + * @param fmax Same as `max`, but for the filter to apply + * @param ftype Same as `type`, but for the filter to apply (Not actually used) + * @param channel The channel, must be either "red", "green", or "blue" + */ +#define libclut_apply__(clut, max, type, filter, fmax, ftype, channel)\ + do {\ + size_t i__, rn__ = (clut)->channel##_size, fn__ = (filter)->channel##_size - 1;\ + double x__, rm__ = (double)(max), m__ = (double)(max) / (double)(fmax);\ + for (i__ = 0; i__ < rn__; i__++) {\ + x__ = (double)((clut)->channel[i__]) / rm__ * (double)fn__;\ + (clut)->channel[i__] = (type)((double)((filter)->channel[(size_t)x__]) * m__);\ + }\ + } while (0) + +/** + * Applies a filter or calibration, using CIE xyY + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param filter Same as `clut`, but for the filter to apply + * @param fmax Same as `max`, but for the filter to apply + * @param ftype Same as `type`, but for the filter to apply (Not actually used) + * @param r Whether to apply the filter for the red curve + * @param g Whether to apply the filter for the green curve + * @param b Whether to apply the filter for the blue curve + */ +#define libclut_cie_apply(clut, max, type, filter, fmax, ftype, r, g, b)\ + do {\ + size_t rfn__ = (filter)->red_size - 1, gfn__ = (filter)->green_size - 1;\ + size_t bfn__ = (filter)->blue_size - 1, x__;\ + size_t rm__ = (double)(max), fm__ = (double)(fmax);\ + libclut_cie__(clut, max, type, 0, r, g, b,\ + (x__ = (size_t)(Y__ / rm__ * rfn__), (double)((filter)->red[x__]) / fm__),\ + (x__ = (size_t)(Y__ / rm__ * gfn__), (double)((filter)->green[x__]) / fm__),\ + (x__ = (size_t)(Y__ / rm__ * bfn__), (double)((filter)->blue[x__]) / fm__));\ + } while (0) + +/** + * Modify a ramp + * + * None of the parameter may have side-effects + * + * This is intended for internal use + * + * @param clut Pointer to the gamma ramps, must have and array + * named `channel` and a scalar named `channel` followed + * by "_size" + * @param channel The channel, must be either "red", "green", or "blue" + * @param type The data type used for each stop in the ramps + * @param expr Expression that evalutes the value a stop should have. + * It can use the variable `LIBCLUT_VALUE` to get the + * current value of the stop. + */ +#define libclut__(clut, channel, type, expr)\ + do {\ + size_t i__, n__ = (clut)->channel##_size;\ + type LIBCLUT_VALUE;\ + for (i__ = 0; i__ < n__; i__++) {\ + LIBCLUT_VALUE = (clut)->channel[i__];\ + (clut)->channel[i__] = (type)(expr);\ + }\ + } while (0) + +/** + * Modify a ramp set in CIE xyY + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * This is intended for internal use + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param utest Whether all channels can be modified at the same time + * This test does not have to include the ramp size + * @param rtest Whether the red channel have to be modified + * @param gtest Whether the green channel have to be modified + * @param btest Whether the blue channel have to be modified + * @param rexpr Expression calculating the intensity of the red channel. + * The current value is stored in `Y__`. + * @param gexpr Expression calculating the intensity of the green channel. + * The current value is stored in `Y__`. + * @param bexpr Expression calculating the intensity of the blue channel. + * The current value is stored in `Y__`. + */ +#define libclut_cie__(clut, max, type, utest, rtest, gtest, btest, rexpr, gexpr, bexpr)\ + do {\ + size_t rn__ = (clut)->red_size;\ + size_t gn__ = (clut)->green_size;\ + size_t bn__ = (clut)->blue_size;\ + size_t i__;\ + double x__, y__, Y__, r__, g__, b__;\ + double m__ = (double)(max);\ + type* rs__ = (clut)->red;\ + type* gs__ = (clut)->green;\ + type* bs__ = (clut)->blue;\ + if ((rn__ == gn__) && (gn__ == bn__) && (utest)) {\ + if (!(rtest))\ + break;\ + for (i__ = 0; i__ < rn__; i__++) {\ + libclut_model_srgb_to_ciexyy(rs__[i__] / m__, gs__[i__] / m__,\ + bs__[i__] / m__, &x__, &y__, &Y__);\ + libclut_model_ciexyy_to_srgb(x__, y__, rexpr, &r__, &g__, &b__);\ + rs__[i__] = (type)(r__ * m__);\ + gs__[i__] = (type)(g__ * m__);\ + bs__[i__] = (type)(b__ * m__);\ + }\ + } else if ((rn__ == gn__) && (gn__ == bn__)) {\ + if (!(rtest) && !(gtest) && !(btest))\ + break;\ + for (i__ = 0; i__ < rn__; i__++) {\ + libclut_model_srgb_to_ciexyy(rs__[i__] / m__, gs__[i__] / m__,\ + bs__[i__] / m__, &x__, &y__, &Y__);\ + if (rtest) {\ + libclut_model_ciexyy_to_srgb(x__, y__, rexpr, &r__, &g__, &b__);\ + rs__[i__] = (type)(r__ * m__);\ + }\ + if (gtest) {\ + libclut_model_ciexyy_to_srgb(x__, y__, gexpr, &r__, &g__, &b__);\ + gs__[i__] = (type)(g__ * m__);\ + }\ + if (btest) {\ + libclut_model_ciexyy_to_srgb(x__, y__, bexpr, &r__, &g__, &b__);\ + bs__[i__] = (type)(b__ * m__);\ + }\ + }\ + } else {\ + if (rtest) {\ + for (i__ = 0; i__ < rn__; i__++) {\ + libclut_cie___(clut, max, type, r, rexpr, i__,\ + libclut_i__(i__, rn__, gn__),\ + libclut_i__(i__, rn__, bn__));\ + }\ + }\ + if (gtest) {\ + for (i__ = 0; i__ < rn__; i__++) {\ + libclut_cie___(clut, max, type, g, gexpr,\ + libclut_i__(i__, gn__, rn__), i__,\ + libclut_i__(i__, gn__, bn__));\ + }\ + }\ + if (btest) {\ + for (i__ = 0; i__ < rn__; i__++) {\ + libclut_cie___(clut, max, type, b, bexpr,\ + libclut_i__(i__, bn__, rn__),\ + libclut_i__(i__, bn__, gn__), i__);\ + }\ + }\ + }\ + } while (0) + +/** + * Modify a ramp stop in CIE xyY + * + * None of the parameter may have side-effects + * + * Requires linking with '-lclut' + * + * This is intended for internal use. + * Assumes the existence of variables defined in `libclut_cie__`. + * + * @param clut Pointer to the gamma ramps, must have the arrays + * `red`, `green`, and `blue`, and the scalars + * `red_size`, `green_size`, and `blue_size`. Ramp + * structures from libgamma or libcoopgamma can be used. + * @param max The maximum value on each stop in the ramps + * @param type The data type used for each stop in the ramps + * @param c Either "r" for red, "g" for green, or "b" for blue + * @param expr Expression calculating the intensity of the channel + * @param ri The index of the stop translated to the red channel + * @param gi The index of the stop translated to the green channel + * @param bi The index of the stop translated to the blue channel + */ +#define libclut_cie___(clut, max, type, c, expr, ri, gi, bi)\ + do {\ + for (i__ = 0; i__ < c##n__; i__++) {\ + libclut_model_srgb_to_ciexyy(rs__[(ri)] / m__, gs__[(gi)] / m__,\ + bs__[(bi)] / m__, &x__, &y__, &Y__);\ + libclut_model_ciexyy_to_srgb(x__, y__, expr, &r__, &g__, &b__);\ + c##s__[i__] = (type)(c##__ * m__);\ + }\ + } while (0) + +/** + * Translate an index from one channel to another + * + * @param i The index in the input channel + * @param in The size of the input channel + * @param out The size of the output channel + * @return The index in the output channel + */ +#define libclut_i__(i, in, out)\ + (size_t)((double)(i) * (double)(out) / (double)(in)) + +#if defined(__GNUC__) && !defined(__clang__) +# define LIBCLUT_GCC_ONLY__(x) x +#else +# define LIBCLUT_GCC_ONLY__(x) /* do nothing */ +#endif + +/** + * Convert one component from [0, 1] linear sRGB to [0, 1] sRGB + * + * If the macro variant is used, the argument must not have + * any side-effects. The macro variant requires linking with + * '-lm' + * + * @param c The linear sRGB value + * @return Corresponding sRGB value + */ +LIBCLUT_GCC_ONLY__(__attribute__((__const__, __leaf__))) +double (libclut_model_linear_to_standard1)(double); +#define libclut_model_linear_to_standard1(c)\ + (((double)(c) <= 0.0031308) ? (12.92 * (double)(c)) : ((1.055) * pow((double)(c), 1 / 2.4) - 0.055)) + +/** + * Convert [0, 1] linear sRGB to [0, 1] sRGB + * + * The macro variant requires linking with '-lm', + * if the 'libclut_model_linear_to_standard1' is defined, + * otherwise it requires linking with '-lclut' + * + * @param r Pointer to the linear red component, + * and output parameter for the red component + * @param g Pointer to the linear green component, + * and output parameter for the green component + * @param b Pointer to the linear blue component, + * and output parameter for the blue component + */ +void (libclut_model_linear_to_standard)(double *, double *, double *); +#define libclut_model_linear_to_standard(r, g, b)\ + do {\ + double *r__ = (r), *g__ = (g), *b__ = (b);\ + *r__ = libclut_model_linear_to_standard1(*r__);\ + *g__ = libclut_model_linear_to_standard1(*g__);\ + *b__ = libclut_model_linear_to_standard1(*b__);\ + } while (0) + +/** + * Convert one component from [0, 1] sRGB to [0, 1] linear sRGB + * + * If the macro variant is used, the argument must not have + * any side-effects. The macro variant requires linking with + * '-lm' + * + * @param c The sRGB value + * @return Corresponding linear sRGB value + */ +LIBCLUT_GCC_ONLY__(__attribute__((__const__, __leaf__))) +double (libclut_model_standard_to_linear1)(double); +#define libclut_model_standard_to_linear1(c) \ + (((double)(c) <= 0.04045) ? ((double)(c) / 12.92) : pow(((double)(c) + 0.055) / 1.055, 2.4)) + +/** + * Convert [0, 1] sRGB to [0, 1] linear sRGB + * + * The macro variant requires linking with '-lm', + * if the 'libclut_model_standard_to_linear1' is defined, + * otherwise it requires linking with '-lclut' + * + * @param r Pointer to the red component, and output + * parameter for the linear red component + * @param g Pointer to the green component, and output + * parameter for the linear green component + * @param b Pointer to the blue component, and output + * parameter for the linear blue component + */ +void (libclut_model_standard_to_linear)(double *, double *, double *); +#define libclut_model_standard_to_linear(r, g, b)\ + do {\ + double *r__ = (r), *g__ = (g), *b__ = (b);\ + *r__ = libclut_model_standard_to_linear1(*r__);\ + *g__ = libclut_model_standard_to_linear1(*g__);\ + *b__ = libclut_model_standard_to_linear1(*b__);\ + } while (0) + +/** + * Convert CIE xyY to CIE XYZ + * + * @param x The x parameter + * @param y The y parameter + * @param Y The Y parameter; this is also the Y (middle) parameter for the CIE XYZ colour + * @param X Output parameter for the X parameter + * @param Z Output parameter for the Z parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_ciexyy_to_ciexyz)(double, double, double, double *, double *); +#define libclut_model_ciexyy_to_ciexyz(x, y, Y, X, Z)\ + do {\ + double x__ = (x), y__ = (y), Y__ = (Y), *X__ = (X), *Z__ = (Z);\ + *X__ = libclut_0__(y__) ? Y__ : (Y__ * x__ / y__);\ + *Z__ = libclut_0__(y__) ? Y__ : (Y__ * (1 - x__ - y__) / y__);\ + } while (0) + +/** + * Convert CIE XYZ to CIE xyY + * + * @param X The X parameter + * @param Y The Y parameter; this is also the Y (last) parameter for the CIE xyY colour + * @param Z The Z parameter + * @param x Output parameter for the x parameter + * @param y Output parameter for the y parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_ciexyz_to_ciexyy)(double, double, double, double *, double *); +#define libclut_model_ciexyz_to_ciexyy(X, Y, Z, x, y)\ + do {\ + double X__ = (X), Y__ = (Y), Z__ = (Z), *x__ = (x), *y__ = (y);\ + double s__ = X__ + Y__ + Z__;\ + if (libclut_0__(s__))\ + *x__ = *y__ = 0;\ + else\ + *x__ = X__ / s__, *y__ = Y__ / s__;\ + } while (0) + +/** + * Convert CIE XYZ to [0, 1] linear sRGB + * + * @param X The X parameter + * @param Y The Y parameter + * @param Z The Z parameter + * @param r Output parameter for the red component + * @param g Output parameter for the green component + * @param b Output parameter for the blue component + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_ciexyz_to_linear)(double, double, double, double *, double *, double *); +#define libclut_model_ciexyz_to_linear(X, Y, Z, r, g, b)\ + do {\ + double X__ = (X), Y__ = (Y), Z__ = (Z);\ + *(r) = ( 3.2404500 * X__) + (-1.537140 * Y__) + (-0.4985320 * Z__);\ + *(g) = (-0.9692660 * X__) + ( 1.876010 * Y__) + ( 0.0415561 * Z__);\ + *(b) = ( 0.0556434 * X__) + (-0.204026 * Y__) + ( 1.0572300 * Z__);\ + } while (0) + +/** + * Convert [0, 1] linear sRGB to CIE XYZ + * + * @param r The red component + * @param g The green component + * @param b The blue component + * @param X Output parameter for the X parameter + * @param Y Output parameter for the Y parameter + * @param Z Output parameter for the Z parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_linear_to_ciexyz)(double, double, double, double *, double *, double *); +#define libclut_model_linear_to_ciexyz(r, g, b, X, Y, Z)\ + do {\ + double r__ = (r), g__ = (g), b__ = (b);\ + *(X) = (0.4124564 * r__) + (0.3575761 * g__) + (0.1804375 * b__);\ + *(Y) = (0.2126729 * r__) + (0.7151522 * g__) + (0.0721750 * b__);\ + *(Z) = (0.0193339 * r__) + (0.1191920 * g__) + (0.9503041 * b__);\ + } while (0) + +/** + * Convert [0, 1] sRGB to CIE xyY + * + * The macro variant requires linking with '-lclut' + * if any of `libclut_model_ciexyz_to_ciexyy`, + * `libclut_model_linear_to_ciexyz`, and + * `libclut_model_standard_to_linear` are undefined. + * The macro variant requires linking with '-lm' if + * neither `libclut_model_standard_to_linear` nor + * `libclut_model_standard_to_linear1` are undefined. + * + * @param r The red component + * @param g The green component + * @param b The blue component + * @param x Output parameter for the x parameter + * @param y Output parameter for the y parameter + * @param Y Output parameter for the Y parameter + */ +void (libclut_model_srgb_to_ciexyy)(double, double, double, double *, double *, double *); +#define libclut_model_srgb_to_ciexyy(r, g, b, x, y, Y)\ + do {\ + double r___ = (r), g___ = (g), b___ = (b);\ + double *x___ = (x), *y___ = (y), *Y___ = (Y);\ + double X___, Z___;\ + libclut_model_standard_to_linear(&r___, &g___, &b___);\ + libclut_model_linear_to_ciexyz(r___, g___, b___, &X___, Y___, &Z___);\ + libclut_model_ciexyz_to_ciexyy(X___, *Y___, Z___, x___, y___);\ + } while (0) + +/** + * Convert CIE xyY to [0, 1] sRGB + * + * The macro variant requires linking with '-lclut' + * if any of `libclut_model_ciexyy_to_ciexyz`, + * `libclut_model_ciexyz_to_linear`, and + * `libclut_model_linear_to_standard` are undefined. + * The macro variant requires linking with '-lm' if + * neither `libclut_model_linear_to_standard` nor + * `libclut_model_linear_to_standard1` are undefined. + * + * @param x The x parameter + * @param y The y parameter + * @param Y The Y parameter + * @param r Output parameter for the red component + * @param g Output parameter for the green component + * @param b Output parameter for the blue component + */ +void (libclut_model_ciexyy_to_srgb)(double, double, double, double *, double *, double *); +#define libclut_model_ciexyy_to_srgb(x, y, Y, r, g, b)\ + do {\ + double x___ = (x), y___ = (y), Y___ = (Y);\ + double *r___ = (r), *g___ = (g), *b___ = (b);\ + double X___, Z___;\ + libclut_model_ciexyy_to_ciexyz(x___, y___, Y___, &X___, &Z___);\ + libclut_model_ciexyz_to_linear(X___, Y___, Z___, r___, g___, b___);\ + libclut_model_linear_to_standard(r___, g___, b___);\ + } while(0) + +/** + * Convert from CIE XYZ to CIE L*a*b* + * + * The macro variant requires linking with '-lm' + * + * @param X The X parameter + * @param Y The Y parameter + * @param Z The Z parameter + * @param L Output parameter for the L* component + * @param a Output parameter for the a* component + * @param b Output parameter for the b* component + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_ciexyz_to_cielab)(double, double, double, double *, double *, double *); +#define libclut_model_ciexyz_to_cielab(X, Y, Z, L, a, b)\ + do {\ + double X__ = (X), Y__ = (Y), Z__ = (Z);\ + X__ /= 0.95047, Z__ /= 1.08883;\ + X__ = LIBCLUT_MODEL_CIEXYZ_TO_CIELAB__(X__);\ + Y__ = LIBCLUT_MODEL_CIEXYZ_TO_CIELAB__(Y__);\ + Z__ = LIBCLUT_MODEL_CIEXYZ_TO_CIELAB__(Z__);\ + *(L) = 116 * Y__ - 16;\ + *(a) = 500 * (X__ - Y__);\ + *(b) = 200 * (Y__ - Z__);\ + } while (0) +#define LIBCLUT_MODEL_CIEXYZ_TO_CIELAB__(C)\ + (((C) > 0.00885642) ? pow((C), 1.0 / 3) : ((7.78 + 703.0 / 99900) * (C) + 0.1379310)) + +/** + * Convert from CIE L*a*b* to CIE XYZ + * + * @param L The L* component + * @param a The a* component + * @param b The b* component + * @param X Output parameter for the X parameter + * @param Y Output parameter for the Y parameter + * @param Z Output parameter for the Z parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_cielab_to_ciexyz)(double, double, double, double *, double *, double *); +#define libclut_model_cielab_to_ciexyz(L, a, b, X, Y, Z)\ + do {\ + double L__ = (L), a__ = (a), b__ = (b);\ + double *X__ = (X), *Y__ = (Y), *Z__ = (Z);\ + *Y__ = (L__ + 16) / 116;\ + *X__ = a__ / 500 + *Y__;\ + *Z__ = *Y__ - b__ / 200;\ + *X__ = LIBCLUT_MODEL_CIELAB_TO_CIEXYZ__(*X__) * 0.95047;\ + *Y__ = LIBCLUT_MODEL_CIELAB_TO_CIEXYZ__(*Y__) * 1.08883;\ + *Z__ = LIBCLUT_MODEL_CIELAB_TO_CIEXYZ__(*Z__);\ + } while (0) +#define LIBCLUT_MODEL_CIELAB_TO_CIEXYZ__(C)\ + (((C)*(C)*(C) > 0.00885642) ? ((C)*(C)*(C)) : (((C) - 0.1379310) / (7.78 + 703.0 / 99900))) + +/** + * Convert from CIE XYZ to CIELUV + * + * Requires linking with `-lm` + * + * @param X The X component + * @param Y The Y component + * @param Z The Z component + * @param Xn The X component of the specified white object (white point) + * @param Yn The Y component of the specified white object (white point) + * @param Zn The Z component of the specified white object (white point) + * @param L Output parameter for the L* parameter + * @param u Output parameter for the u* parameter + * @param v Output parameter for the v* parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_ciexyz_to_cieluv)(double, double, double, double, double, double, double *, double *, double *); +#define libclut_model_ciexyz_to_cieluv(X, Y, Z, Xn, Yn, Zn, L, u, v)\ + do {\ + double xn__ = (Xn), yn__ = (Yn);\ + double t__ = xn__ + 15 * yn__ + 3 * (Zn);\ + double u__ = 4 * xn__ / t__;\ + double v__ = 9 * yn__ / t__;\ + double x__ = (X), y__ = (Y);\ + t__ = x__ + 15 * y__ + 3 * (Z);\ + u__ = 4 * x__ / t__ - u__;\ + v__ = 9 * y__ / t__ - v__;\ + y__ /= yn__;\ + if (y__ * 24389 <= (double)216)\ + y__ *= 24389, y__ /= 27;\ + else\ + y__ = cbrt(y__) * 116 - 16;\ + *(L) = y__;\ + y__ *= 13;\ + *(u) = y__ * u__;\ + *(v) = y__ * v__;\ + } while (0) + +/** + * Convert from CIELUV to CIE XYZ + * + * @param L The L* component + * @param u The u* component + * @param v The v* component + * @param Xn The X component of the specified white object (white point) + * @param Yn The Y component of the specified white object (white point) + * @param Zn The Z component of the specified white object (white point) + * @param X Output parameter for the X parameter + * @param Y Output parameter for the Y parameter + * @param Z Output parameter for the Z parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_cieluv_to_ciexyz)(double, double, double, double, double, double, double *, double *, double *); +#define libclut_model_cieluv_to_ciexyz(L, u, v, Xn, Yn, Zn, X, Y, Z)\ + do {\ + double x__ = (Xn), y__ = (Yn), l__ = (L), l13__ = l__ * 13;\ + double t__ = x__ + 15 * y__ + 3 * (Zn);\ + double u__ = 4 * x__ / t__;\ + double v__ = 9 * y__ / t__;\ + u__ = (u) / l13__ + u__;\ + v__ = (v) / l13__ + v__;\ + if (l__ <= (double)8) {\ + y__ *= l__ * 27 / 24389;\ + } else {\ + l__ += 16, l__ /= 116;\ + y__ *= l__ *= l__ * l__;\ + }\ + *(Y) = y__;\ + *(X) = y__ * 9 * u__ / (4 * v__);\ + *(Z) = y__ * (12 - 3 * u__ - 20 * v__) / (4 * v__);\ + } while (0) + +/** + * Convert from CIELCh to CIE L*u*v* + * + * Requires linking with `-lm` + * + * @param C The C*_uv component + * @param h The h_uv component + * @param u Output parameter for the u* parameter + * @param v Output parameter for the v* parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_cielch_to_cieluv)(double, double, double *, double *); +#define libclut_model_cielch_to_cieluv(C, h, u, v)\ + do {\ + double h__ = (h), C__ = (C);\ + *(v) = sin(h__) * C__;\ + *(u) = cos(h__) * C__;\ + } while (0) + +/** + * Convert from CIE L*u*v* to CIELCh + * + * Requires linking with `-lm` + * + * @param u The u* component + * @param v The v* component + * @param C Output parameter for the C*_uv parameter + * @param h Output parameter for the h_uv parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_cieluv_to_cielch)(double, double, double *, double *); +#define libclut_model_cieluv_to_cielch(u, v, C, h)\ + do {\ + double u__ = (u), v__ = (v);\ + *(C) = sqrt(u__ * u__ + v__ * v__);\ + *(h) = atan2(v__, u__); \ + } while (0) + +/** + * Convert from sRGB to YIQ + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_standard_to_linear1` is not undefined + * + * @param r The R component + * @param g The G component + * @param b The B component + * @param y Output parameter for the Y parameter + * @param i Output parameter for the I parameter + * @param q Output parameter for the Q parameter + */ +void (libclut_model_srgb_to_yiq)(double, double, double, double *, double *, double *); +#define libclut_model_srgb_to_yiq(r, g, b, y, i, q)\ + do {\ + double r__ = libclut_model_standard_to_linear1(r) / 100000000000000000ULL;\ + double g__ = libclut_model_standard_to_linear1(g) / 10000000000000000ULL;\ + double b__ = libclut_model_standard_to_linear1(b) / 100000000000000000ULL;\ + *(y) = r__ * 29893602129377540ULL + g__ * 5870430744511212ULL + b__ * 11402090425510336ULL;\ + *(i) = r__ * 59594574307079930ULL - g__ * 2743886357457892ULL - b__ * 32155710732501010ULL;\ + *(q) = r__ * 21149734030682846ULL - g__ * 5229106903029739ULL + b__ * 31141334999614540ULL;\ + } while (0) + +/** + * Convert from YIQ to sRGB + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_linear_to_standard1` is not undefined + * + * @param y The Y component + * @param i The I component + * @param q The Q component + * @param r Output parameter for the R parameter + * @param g Output parameter for the G parameter + * @param b Output parameter for the B parameter + */ +void (libclut_model_yiq_to_srgb)(double, double, double, double *, double *, double *); +#define libclut_model_yiq_to_srgb(y, i, q, r, g, b)\ + do {\ + double y__ = (y), i__ = (i), q__ = (q), r__, g__, b__;\ + r__ = y__ + i__ * 956 / 1000 + q__ * 621 / 1000;\ + g__ = y__ - i__ * 272 / 1000 - q__ * 647 / 1000;\ + b__ = y__ - i__ * 1106 / 1000 + q__ * 1703 / 1000;\ + *(r) = libclut_model_linear_to_standard1(r__);\ + *(g) = libclut_model_linear_to_standard1(g__);\ + *(b) = libclut_model_linear_to_standard1(b__);\ + } while (0) + +/** + * Convert from sRGB to YDbDr + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_standard_to_linear1` is not undefined + * + * @param r The R component + * @param g The G component + * @param b The B component + * @param y Output parameter for the Y parameter + * @param db Output parameter for the Db parameter + * @param dr Output parameter for the Dr parameter + */ +void (libclut_model_srgb_to_ydbdr)(double, double, double, double *, double *, double *); +#define libclut_model_srgb_to_ydbdr(r, g, b, y, db, dr)\ + do {\ + double r__ = libclut_model_standard_to_linear1(r);\ + double g__ = libclut_model_standard_to_linear1(g);\ + double b__ = libclut_model_standard_to_linear1(b);\ + *(y) = r__ * 299 / 1000 + g__ * 587 / 1000 + b__ * 114 / 1000;\ + *(db) = -r__ * 450 / 1000 - g__ * 883 / 1000 + b__ * 1333 / 1000;\ + *(dr) = -r__ * 1333 / 1000 + g__ * 1116 / 1000 + b__ * 217 / 1000;\ + } while (0) + +/** + * Convert from YDbDr to sRGB + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_linear_to_standard1` is not undefined + * + * @param y The Y component + * @param db The Db component + * @param dr The Dr component + * @param r Output parameter for the R parameter + * @param g Output parameter for the G parameter + * @param b Output parameter for the B parameter + */ +void (libclut_model_ydbdr_to_srgb)(double, double, double, double *, double *, double *); +#define libclut_model_ydbdr_to_srgb(y, db, dr, r, g, b)\ + do {\ + double y__ = (y), db__ = (db), dr__ = (dr), r__, g__, b__;\ + db__ /= 1000000000000000ULL;\ + dr__ /= 1000000000000000ULL;\ + r__ = y__ + db__ * 92303716148ULL - dr__ * 525912630661865ULL;\ + g__ = y__ - db__ * 129132898890509ULL + dr__ * 267899328207599ULL;\ + b__ = y__ + db__ * 664679059978955ULL - dr__ * 79202543533ULL;\ + *(r) = libclut_model_linear_to_standard1(r__);\ + *(g) = libclut_model_linear_to_standard1(g__);\ + *(b) = libclut_model_linear_to_standard1(b__);\ + } while (0) + +/** + * Convert from YUV to YDbDr + * + * @param u The U component + * @param v The V component + * @param db Output parameter for the Db parameter + * @param dr Output parameter for the Dr parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_yuv_to_ydbdr)(double, double, double *, double *); +#define libclut_model_yuv_to_ydbdr(u, v, db, dr)\ + do {\ + *(db) = 3069 * (u) / 1000;\ + *(dr) = -2169 * (v) / 1000;\ + } while (0) + +/** + * Convert from YDbDr to YUV + * + * @param db The Db component + * @param dr The Dr component + * @param u Output parameter for the U parameter + * @param v Output parameter for the V parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_ydbdr_to_yuv)(double, double, double *, double *); +#define libclut_model_ydbdr_to_yuv(db, dr, u, v)\ + do {\ + *(u) = (db) * 1000 / 3069;\ + *(v) = (dr) * 1000 / -2169;\ + } while (0) + +/** + * Convert from sRGB to YPbPr + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_standard_to_linear1` is not undefined + * + * @param r The R component + * @param g The G component + * @param b The B component + * @param y Output parameter for the Y parameter + * @param pb Output parameter for the Pb parameter + * @param pr Output parameter for the Pr parameter + */ +void (libclut_model_srgb_to_ypbpr)(double, double, double, double *, double *, double *); +#define libclut_model_srgb_to_ypbpr(r, g, b, y, pb, pr)\ + do {\ + double r__ = libclut_model_standard_to_linear1(r);\ + double g__ = libclut_model_standard_to_linear1(g);\ + double b__ = libclut_model_standard_to_linear1(b);\ + double y__;\ + y__ = r__ * 2126 / 10000;\ + y__ += g__ * 7152 / 10000;\ + y__ += b__ * 722 / 10000;\ + *(y) = y__;\ + *(pb) = b__ - y__;\ + *(pr) = r__ - y__;\ + } while (0) + +/** + * Convert from YPbPr to sRGB + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_linear_to_standard1` is not undefined + * + * @param y The Y component + * @param pb The Pb component + * @param pr The Pr component + * @param r Output parameter for the R parameter + * @param g Output parameter for the G parameter + * @param b Output parameter for the B parameter + */ +void (libclut_model_ypbpr_to_srgb)(double, double, double, double *, double *, double *); +#define libclut_model_ypbpr_to_srgb(y, pb, pr, r, g, b)\ + do {\ + double y__ = (y), r__ = (pr) + y__, g__, b__ = (pb) + y__;\ + y__ -= r__ * 2126 / 10000;\ + y__ -= b__ * 722 / 10000;\ + g__ = y__ * 10000 / 7152;\ + *(r) = libclut_model_linear_to_standard1(r__);\ + *(g) = libclut_model_linear_to_standard1(g__);\ + *(b) = libclut_model_linear_to_standard1(b__);\ + } while (0) + +/** + * Convert from sRGB to YCgCo + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_standard_to_linear1` is not undefined + * + * @param r The R component + * @param g The G component + * @param b The B component + * @param y Output parameter for the Y parameter + * @param cg Output parameter for the Cg parameter + * @param co Output parameter for the Co parameter + */ +void (libclut_model_srgb_to_ycgco)(double, double, double, double *, double *, double *); +#define libclut_model_srgb_to_ycgco(r, g, b, y, cg, co)\ + do {\ + double r__ = libclut_model_standard_to_linear1(r);\ + double g__ = libclut_model_standard_to_linear1(g);\ + double b__ = libclut_model_standard_to_linear1(b);\ + *(y) = r__ / 4 + g__ / 2 + b__ / 4;\ + *(cg) = r__ / -4 + g__ / 2 - b__ / 4;\ + *(co) = r__ / 2 - b__ / 2;\ + } while (0) + +/** + * Convert from YCgCo to sRGB + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_linear_to_standard1` is not undefined + * + * @param y The Y component + * @param cg The Cg component + * @param co The Co component + * @param r Output parameter for the R parameter + * @param g Output parameter for the G parameter + * @param b Output parameter for the B parameter + */ +void (libclut_model_ycgco_to_srgb)(double, double, double, double *, double *, double *); +#define libclut_model_ycgco_to_srgb(y, cg, co, r, g, b)\ + do {\ + double y__ = (y), cg__ = (cg), co__ = (co);\ + double r__ = y__ - cg__ + co__;\ + double g__ = y__ + cg__;\ + double b__ = y__ - cg__ - co__;\ + *(r) = libclut_model_linear_to_standard1(r__);\ + *(g) = libclut_model_linear_to_standard1(g__);\ + *(b) = libclut_model_linear_to_standard1(b__);\ + } while (0) + +/** + * Convert from CIE 1960 UCS to CIE XYZ + * + * @param u The u component + * @param v The v component + * @param Y The Y component + * @param x Output parameter for the X parameter + * @param y Output parameter for the Y parameter + * @param z Output parameter for the Z parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_cie_1960_ucs_to_ciexyz)(double, double, double, double *, double *, double *); +#define libclut_model_cie_1960_ucs_to_ciexyz(u, v, Y, x, y, z)\ + do {\ + double u__ = (u), v__ = (v), y__ = (Y);\ + *(y) = y__;\ + *(x) = 3 * y__ * u__ / (2 * v__);\ + *(z) = (4 * y__ - y__ * u__ - 10 * y__ * v__) / (2 * v__);\ + } while (0) + +/** + * Convert from CIE XYZ to CIE 1960 UCS + * + * @param x The X component + * @param y The Y component + * @param z The Z component + * @param u Output parameter for the u parameter + * @param v Output parameter for the v parameter + * @param Y Output parameter for the Y parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_ciexyz_to_cie_1960_ucs)(double, double, double, double *, double *, double *); +#define libclut_model_ciexyz_to_cie_1960_ucs(x, y, z, u, v, Y)\ + do {\ + double x__ = (x), y__ = (y);\ + double d__ = x__ + 15 * y__ + 3 * (z);\ + *(u) = 4 * x__ / d__;\ + *(v) = 6 * y__ / d__;\ + *(Y) = y__;\ + } while (0) + +/** + * Convert from CIEUVW to CIE 1960 UCS + * + * @param U:double The U* component + * @param V:double The V* component + * @param W:double The W* component + * @param u0:double The u parameter for the white point + * @param v0:double The v parameter for the white point + * @param u:double* Output parameter for the u parameter + * @param v:double* Output parameter for the v parameter + * @param Y:double* Output parameter for the Y parameter + */ +LIBCLUT_GCC_ONLY__(__attribute__((__leaf__))) +void (libclut_model_cieuvw_to_cie_1960_ucs)(double, double, double, double, double, double *, double *, double *); +#define libclut_model_cieuvw_to_cie_1960_ucs(U, V, W, u0, v0, u, v, Y)\ + do {\ + double w__ = (W), y__ = (w__ + 17) / 25;\ + *(Y) = y__ *= y__ * y__;\ + w__ *= 13;\ + *(u) = (U) / w__ + (u0);\ + *(v) = (V) / w__ + (v0);\ + } while (0) + +/** + * Convert from CIE 1960 UCS to CIEUVW + * + * Requires linking with `-lm` + * + * @param u:double The u component + * @param v:double The v component + * @param Y:double The Y component + * @param u0:double The u parameter for the white point + * @param v0:double The v parameter for the white point + * @param U:double* Output parameter for the U* parameter + * @param V:double* Output parameter for the V* parameter + * @param W:double* Output parameter for the W* parameter + */ +void (libclut_model_cie_1960_ucs_to_cieuvw)(double, double, double, double, double, double *, double *, double *); +#define libclut_model_cie_1960_ucs_to_cieuvw(u, v, Y, u0, v0, U, V, W)\ + do {\ + double w__ = 25 * cbrt(Y) - 17;\ + *(W) = w__;\ + w__ *= 13;\ + *(U) = w__ * ((u) - (u0));\ + *(V) = w__ * ((v) - (v0));\ + } while (0) + +/** + * Create a matrix for converting values between + * two RGB colour spaces + * + * @param from The input colour space, the Y-component is only necessary + * for the white point, `NULL` for CIE XYZ + * @param to The output colour space, the Y-component is only necessary + * for the white point, `NULL` for CIE XYZ + * @param M Output matrix for conversion from `from` to `to` + * @param Minv Output matrix for conversion from `to` to `from`, may be `NULL` + * @return Zero on success, -1 on error + * + * @throws EINVAL The colour space cannot be used + */ +int libclut_model_get_rgb_conversion_matrix(const libclut_rgb_colour_space_t *, const libclut_rgb_colour_space_t *, + libclut_colour_space_conversion_matrix_t, libclut_colour_space_conversion_matrix_t); + +/** + * Convert an RGB colour into another RGB colour space + * + * Both RGB colour spaces must have same gamma functions as sRGB + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_standard_to_linear1` or + * `libclut_model_linear_to_standard1` is not undefined + * + * @param r The red component of the colour to convert + * @param g The green component of the colour to convert + * @param b The blue component of the colour to convert + * @param M Conversion matrix, create with `libclut_model_get_rgb_conversion_matrix`, + * must not have side-effects + * @param out_r Output parameter for the new red component + * @param out_g Output parameter for the new green component + * @param out_b Output parameter for the new blue component + */ +void (libclut_model_convert_rgb)(double, double, double, libclut_colour_space_conversion_matrix_t, double *, double *, double *); +#define libclut_model_convert_rgb(r, g, b, M, out_r, out_g, out_b)\ + do {\ + double r__ = libclut_model_standard_to_linear1(r);\ + double g__ = libclut_model_standard_to_linear1(g);\ + double b__ = libclut_model_standard_to_linear1(b);\ + *(out_r) = libclut_model_linear_to_standard1((M)[0][0] * r__ + (M)[0][1] * g__ + (M)[0][2] * b__);\ + *(out_g) = libclut_model_linear_to_standard1((M)[1][0] * r__ + (M)[1][1] * g__ + (M)[1][2] * b__);\ + *(out_b) = libclut_model_linear_to_standard1((M)[2][0] * r__ + (M)[2][1] * g__ + (M)[2][2] * b__);\ + } while (0) + +/** + * Convert an RGB colour of a custom RGB colour space to CIE XYZ + * + * The RGB colour space must have same gamma functions as sRGB + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_standard_to_linear1` is not undefined + * + * @param r The red component + * @param g The green component + * @param b The blue component + * @param M Conversion matrix, create with `libclut_model_get_rgb_conversion_matrix`, + * must not have side-effects + * @param x Output parameter for the X component + * @param y Output parameter for the Y component + * @param z Output parameter for the Z component + */ +void (libclut_model_rgb_to_ciexyz)(double, double, double, libclut_colour_space_conversion_matrix_t, double *, double *, double *); +#define libclut_model_rgb_to_ciexyz(r, g, b, M, x, y, z)\ + do {\ + double r__ = libclut_model_standard_to_linear1(r);\ + double g__ = libclut_model_standard_to_linear1(g);\ + double b__ = libclut_model_standard_to_linear1(b);\ + *(x) = (M)[0][0] * r__ + (M)[0][1] * g__ + (M)[0][2] * b__;\ + *(y) = (M)[1][0] * r__ + (M)[1][1] * g__ + (M)[1][2] * b__;\ + *(z) = (M)[2][0] * r__ + (M)[2][1] * g__ + (M)[2][2] * b__;\ + } while (0) + +/** + * Convert a CIE XYZ colour to a custom RGB colour space + * + * The RGB colour space must have same gamma functions as sRGB + * + * Requires linking with '-lclut', or '-lm' if + * `libclut_model_linear_to_standard1` is not undefined + * + * @param x The X component + * @param y The Y component + * @param z The Z component + * @param M Conversion matrix, create with `libclut_model_get_rgb_conversion_matrix`, + * must not have side-effects. + * @param r Output parameter for the red component + * @param g Output parameter for the green component + * @param b Output parameter for the blue component + */ +void (libclut_model_ciexyz_to_rgb)(double, double, double, libclut_colour_space_conversion_matrix_t, double *, double *, double *); +#define libclut_model_ciexyz_to_rgb(x, y, z, M, r, g, b)\ + do {\ + double x__ = (x), y__ = (y), z__ = (z);\ + *(r) = libclut_model_linear_to_standard1((M)[0][0] * x__ + (M)[0][1] * y__ + (M)[0][2] * z__);\ + *(g) = libclut_model_linear_to_standard1((M)[1][0] * x__ + (M)[1][1] * y__ + (M)[1][2] * z__);\ + *(b) = libclut_model_linear_to_standard1((M)[2][0] * x__ + (M)[2][1] * y__ + (M)[2][2] * z__);\ + } while (0) + +#if defined(__clang__) +# pragma GCC diagnostic pop +#endif + +#endif -- cgit v1.2.3-70-g09d2