Add common.h

Signed-off-by: Mattias Andrée <maandree@kth.se>

1 files changed, 487 insertions, 0 deletions

diff --git a/common.h b/common.h
new file mode 100644
index 0000000..5c781b6
--- /dev/null
+++ b/common.h
@@ -0,0 +1,487 @@
+/* See LICENSE file for copyright and license details. */
+
+#if defined(__GNUC__) && !defined(__clang__)
+# pragma GCC diagnostic ignored "-Wfloat-equal"
+#elif defined(__clang__)
+# pragma clang diagnostic ignored "-Wcomma"
+# pragma clang diagnostic ignored "-Wfloat-equal"
+# pragma clang diagnostic ignored "-Wvla"
+# pragma clang diagnostic ignored "-Wtautological-compare"
+# pragma clang diagnostic ignored "-Wcovered-switch-default"
+# pragma clang diagnostic ignored "-Wfloat-conversion"
+# pragma clang diagnostic ignored "-Wconversion"
+# pragma clang diagnostic ignored "-Wdouble-promotion"
+# pragma clang diagnostic ignored "-Wswitch-enum"
+# pragma clang diagnostic ignored "-Wcast-align"
+#endif
+
+
+
+#define PI          D(3.14159265358979323846)
+#define PI2         (2 * D(3.14159265358979323846))
+#define WASDIV0(X)  (xisinf(X) || xisnan(X))
+
+#define REGULAR(S, T)    ((T) <= (S)->TRANSITION ? (S)->SLOPE * (T) : (1 + (S)->OFFSET) * xpow((T), 1 / (S)->GAMMA) - (S)->OFFSET)
+#define INVREGULAR(S, T) ((T) <= (S)->TRANSITIONINV ? (T) / (S)->SLOPE : xpow(((T) + (S)->OFFSET) / (1 + (S)->OFFSET), (S)->GAMMA))
+
+#define XTRANSFORM(X, Y, Z, A, B, C, R1C1, R1C2, R1C3, R2C1, R2C2, R2C3, R3C1, R3C2, R3C3)\
+	do {\
+		TYPE a__ = (A), b__ = (B), c__ = (C);\
+		(X) = (R1C1) * a__ + (R1C2) * b__ + (R1C3) * c__;\
+		(Y) = (R2C1) * a__ + (R2C2) * b__ + (R2C3) * c__;\
+		(Z) = (R3C1) * a__ + (R3C2) * b__ + (R3C3) * c__;\
+	} while (0)
+
+#define TRANSFORM(X, Y, Z, A, B, C, R1C1, R1C2, R1C3, R2C1, R2C2, R2C3, R3C1, R3C2, R3C3)\
+	XTRANSFORM((X), (Y), (Z), (A), (B), (C),\
+	           D(R1C1), D(R1C2), D(R1C3),\
+	           D(R2C1), D(R2C2), D(R2C3),\
+	           D(R3C1), D(R3C2), D(R3C3))
+
+
+
+#define SLOPE transfer.regular.slope
+#define TRANSITIONINV transfer.regular.transitioninv
+#define TRANSITION transfer.regular.transition
+#define GAMMA transfer.regular.gamma
+#define OFFSET transfer.regular.offset
+#define TO_ENCODED_RED transfer.custom.to_encoded_red
+#define TO_DECODED_RED transfer.custom.to_decoded_red
+#define TO_ENCODED_GREEN transfer.custom.to_encoded_green
+#define TO_DECODED_GREEN transfer.custom.to_decoded_green
+#define TO_ENCODED_BLUE transfer.custom.to_encoded_blue
+#define TO_DECODED_BLUE transfer.custom.to_decoded_blue
+
+
+
+#define CIEXYZ_TO_RGB(FROM_X, FROM_Y, FROM_Z, TO_R, TO_G, TO_B, TO_MINV)\
+	XTRANSFORM((TO_R), (TO_G), (TO_B), (FROM_X), (FROM_Y), (FROM_Z),\
+	           (TO_MINV)[0][0], (TO_MINV)[0][1], (TO_MINV)[0][2],\
+	           (TO_MINV)[1][0], (TO_MINV)[1][1], (TO_MINV)[1][2],\
+	           (TO_MINV)[2][0], (TO_MINV)[2][1], (TO_MINV)[2][2])
+
+#define CIEXYZ_TO_SRGB(FROM_X, FROM_Y, FROM_Z, TO_R, TO_G, TO_B)\
+	TRANSFORM((TO_R), (TO_G), (TO_B), (FROM_X), (FROM_Y), (FROM_Z),\
+	           3.240446254647737500675930277794, -1.537134761820080575134284117667, -0.498530193022728718155178739835,\
+	          -0.969266606244679751469561779231,  1.876011959788370209167851498933,  0.041556042214430065351304932619,\
+	           0.055643503564352832235773149705, -0.204026179735960239147729566866,  1.057226567722703292062647051353)
+
+#define YIQ_TO_SRGB(FROM_Y, FROM_I, FROM_Q, TO_R, TO_G, TO_B)\
+	TRANSFORM((TO_R), (TO_G), (TO_B), (FROM_Y), (FROM_I), (FROM_Q),\
+	          1.,  0.95629483232089407263032398986979387700557708740234,  0.62102512544472865396727456754888407886028289794922,\
+	          1., -0.27212147408397735492968649850809015333652496337891, -0.64738095351761570928061928498209454119205474853516,\
+	          1., -1.10698990856712820018969978264067322015762329101562,  1.70461497549882934343656870623817667365074157714844)
+	/* ⎛1    1000 cos 33° / 877                                  1000 sin 33° / 877                               ⎞
+	   ⎜1    9500 sin 33° / 24067 - 299000 * cos 33° / 514799    -9500 cos 33° / 24067 - 299000 * sin 33° / 514799⎟
+	   ⎝1    -250 sin 33° / 123                                  250 cos 33° / 123                                ⎠
+	*/
+
+#define YDBDR_TO_SRGB(FROM_Y, FROM_DB, FROM_DR, TO_R, TO_G, TO_B)\
+	TRANSFORM((TO_R), (TO_G), (TO_B), (FROM_Y), (FROM_DB), (FROM_DR),\
+	          1.,  0.000092303716148, -0.525912630661865,\
+	          1., -0.129132898890509,  0.267899328207599,\
+	          1.,  0.664679059978955, -0.000079202543533)
+
+#define YPBPR_TO_SRGB(FROM_Y, FROM_PB, FROM_PR, TO_R, TO_G, TO_B)\
+	do {\
+		TYPE y__ = (FROM_Y);\
+		TYPE r__ = (FROM_PR) + y__;\
+		TYPE b__ = (FROM_PB) + y__;\
+		(TO_R) = r__;\
+		(TO_B) = b__;\
+		(TO_G) = (y__ - r__ * D(0.2126) - b__ * D(0.0722)) / D(0.7152);\
+	} while (0)
+
+#define YUV_TO_SRGB(FROM_Y, FROM_U, FROM_V, TO_R, TO_G, TO_B)\
+	TRANSFORM((TO_R), (TO_G), (TO_B), (FROM_Y), (FROM_U), (FROM_V),\
+	          1.,  0.00028328010485821202317155420580263580632163211703,  1.14070449590558520291949662350816652178764343261719,\
+	          1., -0.39630886669497211727275498560629785060882568359375, -0.58107364288228224857846271333983168005943298339844,\
+	          1.,  2.03990003507541306504435851820744574069976806640625,  0.00017179031692307700847528739718228507626918144524)
+
+#define YCGCO_TO_SRGB(FROM_Y, FROM_CG, FROM_CO, TO_R, TO_G, TO_B)\
+	TRANSFORM((TO_R), (TO_G), (TO_B), (FROM_Y), (FROM_CG), (FROM_CO),\
+	           1., -1.,  1.,\
+	           1.,  1.,  0.,\
+	           1., -1., -1.)
+
+#define CIEXYZ_TO_CIEXYY(FROM_X, FROM_Y, FROM_Z, TO_X, TO_Y, TO_YY)\
+	do {\
+		TYPE x__ = (FROM_X), y__ = (FROM_Y), z__ = (FROM_Z);\
+		TYPE s__ = x__ + y__ + z__;\
+		x__ = x__ / s__;\
+		s__ = y__ / s__;\
+		if (WASDIV0(x__) || WASDIV0(s__))\
+			x__ = s__ = 0;\
+		(TO_X) = x__;\
+		(TO_Y) = s__;\
+		(TO_YY) = y__;\
+	} while (0)
+
+#define RGB_TO_CIEXYZ(FROM_R, FROM_G, FROM_B, TO_X, TO_Y, TO_Z, FROM_M)\
+	XTRANSFORM((TO_X), (TO_Y), (TO_Z), (FROM_R), (FROM_G), (FROM_B),\
+	           (FROM_M)[0][0], (FROM_M)[0][1], (FROM_M)[0][2],\
+	           (FROM_M)[1][0], (FROM_M)[1][1], (FROM_M)[1][2],\
+	           (FROM_M)[2][0], (FROM_M)[2][1], (FROM_M)[2][2])
+
+#define SRGB_TO_CIEXYZ(FROM_R, FROM_G, FROM_B, TO_X, TO_Y, TO_Z)\
+	TRANSFORM((TO_X), (TO_Y), (TO_Z), (FROM_R), (FROM_G), (FROM_B),\
+	          0.412457445582367576708548995157, 0.357575865245515878143578447634, 0.180437247826399665973085006954,\
+	          0.212673370378408277403536885686, 0.715151730491031756287156895269, 0.072174899130559869164791564344,\
+	          0.019333942761673460208893260415, 0.119191955081838593666354597644, 0.950302838552371742508739771438)
+
+#define CIEXYY_TO_CIEXYZ(FROM_X, FROM_Y, FROM_YY, TO_X, TO_Y, TO_Z)\
+	do {\
+		TYPE x__ = (FROM_X), y__ = (FROM_Y), Y__ = (FROM_YY);\
+		TYPE Yy__ = Y__ / y__;\
+		if (WASDIV0(Yy__)) {\
+			(TO_X) = (TO_Y) = (TO_Z) = Y__;\
+		} else {\
+			(TO_X) = x__ * Yy__;\
+			(TO_Y) = Y__;\
+			(TO_Z) = (1 - x__ - y__) * Yy__;\
+		}\
+	} while (0)
+
+#define CIELAB_TO_CIEXYZ(FROM_L, FROM_A, FROM_B, TO_X, TO_Y, TO_Z)\
+	do {\
+		TYPE Y__ = ((FROM_L) + 16) / 116;\
+		TYPE X__ = Y__ + (FROM_A) / 500;\
+		TYPE Z__ = Y__ - (FROM_B) / 200;\
+		(TO_X) = cielab_finv(X__) * D(0.95047);\
+		(TO_Y) = cielab_finv(Y__);\
+		(TO_Z) = cielab_finv(Z__) * D(1.08883);\
+	} while (0)
+
+#define CIELUV_TO_CIEXYZ(FROM_L, FROM_U, FROM_V, TO_X, TO_Y, TO_Z, FROM_WHITE_X, FROM_WHITE_Y, FROM_WHITE_Z)\
+	do {\
+		TYPE L__ = (FROM_L), X__ = (FROM_WHITE_X), Y__ = (FROM_WHITE_Y), Z__ = (FROM_WHITE_Z);\
+		TYPE L13__ = L__ * 13;\
+		TYPE t__ = X__ + 15 * Y__ + 3 * Z__;\
+		TYPE u__ = (FROM_U) / L13__ + 4 * X__ / t__;\
+		TYPE v__ = (FROM_V) / L13__ + 9 * Y__ / t__;\
+		if (L__ <= 8) {\
+			Y__ *= L__ * 27 / 24389;\
+		} else {\
+			L__ = (L__ + 16) / 116;\
+			Y__ *= L__ * L__ * L__;\
+		}\
+		(TO_X) = D(2.25) * Y__ * u__ / v__;\
+		(TO_Y) = Y__;\
+		(TO_Z) = Y__ * (3 / v__ - D(0.75) * u__ / v__ - 5);\
+	} while (0)
+
+#define YIQ_TO_CIEXYZ(FROM_Y, FROM_I, FROM_Q, TO_X, TO_Y, TO_Z)\
+	TRANSFORM((TO_X), (TO_Y), (TO_Z), (FROM_Y), (FROM_I), (FROM_Q),\
+	           0.95047055865428309306963683411595411598682403564453,\
+	           0.09738463974178063198294807989441324025392532348633,\
+	           0.33223466706854809515903070860076695680618286132812,\
+	           0.99999999999999988897769753748434595763683319091797,\
+	          -0.07112658301916767455974621725545148365199565887451,\
+	          -0.20786968876398304040264974901219829916954040527344,\
+	           1.08882873639588373393394249433185905218124389648438,\
+	          -1.06592139332461721679123911599162966012954711914062,\
+	           1.55474471255181900808395312196807935833930969238281)
+
+#define YDBDR_TO_CIEXYZ(FROM_Y, FROM_DB, FROM_DR, TO_X, TO_Y, TO_Z)\
+	TRANSFORM((TO_X), (TO_Y), (TO_Z), (FROM_Y), (FROM_DB), (FROM_DR),\
+	           0.95047055865428309306963683411595411598682403564453,\
+	           0.07379612357298576119646327242662664502859115600586,\
+	          -0.12113653724874726136384595065464964136481285095215,\
+	           0.99999999999999988897769753748434595763683319091797,\
+	          -0.04435684145428285540813106990754022262990474700928,\
+	           0.07973534004202506575431641522300196811556816101074,\
+	           1.08882873639588373393394249433185905218124389648438,\
+	           0.61625657933494271123464613992837257683277130126953,\
+	           0.02168821359337728266192257819966471288353204727173)
+
+#define YUV_TO_CIEXYZ(FROM_Y, FROM_U, FROM_V, TO_X, TO_Y, TO_Z)\
+	TRANSFORM((TO_X), (TO_Y), (TO_Z), (FROM_Y), (FROM_U), (FROM_V),\
+	           0.95047055865428309306963683411595411598682403564453,\
+	           0.22648030324549334180161963558930438011884689331055,\
+	           0.26274514929253273143316960158699657768011093139648,\
+	           0.99999999999999988897769753748434595763683319091797,\
+	          -0.13613114642319409930415474718756740912795066833496,\
+	          -0.17294595255115238763288232348713790997862815856934,\
+	           1.08882873639588373393394249433185905218124389648438,\
+	           1.89129144197893928058817891724174842238426208496094,\
+	          -0.04704173528403532422714761196402832865715026855469)
+
+#define YPBPR_TO_CIEXYZ(FROM_Y, FROM_PB, FROM_PR, TO_X, TO_Y, TO_Z)\
+	TRANSFORM((TO_X), (TO_Y), (TO_Z), (FROM_Y), (FROM_PB), (FROM_PR),\
+	           0.95047055865428298204733437160030007362365722656250,\
+	           0.14433968424876231217979238863335922360420227050781,\
+	           0.30616461986760712399302519770571961998939514160156,\
+	           0.99999999999999977795539507496869191527366638183594,\
+	          -0.00002022802471486617736928792510298080742359161377,\
+	           0.00008771894888734421691367515450110659003257751465,\
+	           1.08882873639588373393394249433185905218124389648438,\
+	           0.93827031735982591165168287261622026562690734863281,\
+	          -0.01609699914324668607035206946420657914131879806519)
+
+#define YCGCO_TO_CIEXYZ(FROM_Y, FROM_CG, FROM_CO, TO_X, TO_Y, TO_Z)\
+	TRANSFORM((TO_X), (TO_Y), (TO_Z), (FROM_Y), (FROM_CG), (FROM_CO),\
+	           0.95047055865428309306963683411595411598682403564453,\
+	          -0.23531882816325136453805555447615915909409523010254,\
+	           0.23202019775596791073546398820326430723071098327637,\
+	           0.99999999999999988897769753748434595763683319091797,\
+	           0.43030346098206362359661625305307097733020782470703,\
+	           0.14049847124784842211653312915586866438388824462891,\
+	           1.08882873639588373393394249433185905218124389648438,\
+	          -0.85044482623220662986796014592982828617095947265625,\
+	          -0.93096889579069830311652822274481877684593200683594)
+
+#define CIE1960UCS_TO_CIEXYZ(FROM_U, FROM_V, FROM_Y, TO_X, TO_Y, TO_Z)\
+	do {\
+		TYPE u__ = (FROM_U), v__ = (FROM_V), Y__ = (FROM_Y);\
+		(TO_X) = D(1.5) * Y__ * u__ / v__;\
+		(TO_Y) = Y__;\
+		(TO_Z) = (4 * Y__ - Y__ * u__ - 10 * Y__ * v__) / (2 * v__);\
+	} while (0)
+
+#define CIEXYZ_TO_CIELAB(FROM_X, FROM_Y, FROM_Z, TO_L, TO_A, TO_B)\
+	do {\
+		TYPE X__ = (FROM_X), Y__ = (FROM_Y), Z__ = (FROM_Z);\
+		X__ /= D(0.95047);\
+		Z__ /= D(1.08883);\
+		Y__ = cielab_f(Y__);\
+		(TO_L) = 116 * Y__ - 16;\
+		(TO_A) = 500 * (cielab_f(X__) - Y__);\
+		(TO_B) = 200 * (Y__ - cielab_f(Z__));\
+	} while (0)
+
+#define CIEXYZ_TO_CIELUV(FROM_X, FROM_Y, FROM_Z, TO_L, TO_U, TO_V, TO_WHITE_X, TO_WHITE_Y, TO_WHITE_Z)\
+	do {\
+		TYPE WX__ = (TO_WHITE_X), WY__ = (TO_WHITE_Y), WZ__ = (TO_WHITE_Z);\
+		TYPE X__ = (FROM_X), Y__ = (FROM_Y), Z__ = (FROM_Z);\
+		TYPE L2__, L__, u__, v__, t__;\
+		t__ = WX__ + 15 * WY__ + 3 * WZ__;\
+		u__ = 4 * WX__ / t__;\
+		v__ = 9 * WY__ / t__;\
+		t__ = X__ + 15 * Y__ + 3 * Z__;\
+		u__ = 4 * X__ / t__ - u__;\
+		v__ = 9 * Y__ / t__ - v__;\
+		L__ = Y__ / WY__;\
+		L2__ = L__ * 24389;\
+		L__ = L2__ <= 216 ? L2__ / 27 : xcbrt(L__) * 116 - 16;\
+		(TO_L) = L__;\
+		L__ *= 13;\
+		(TO_U) = u__ * L__;\
+		(TO_V) = v__ * L__;\
+	} while (0)
+
+#define CIELCHUV_TO_CIELUV(FROM_L, FROM_C, FROM_H, TO_L, TO_U, TO_V)\
+	do {\
+		TYPE C__ = (FROM_C), h__ = (FROM_H);\
+		(TO_L) = (FROM_L);\
+		(TO_U) = C__ * xcos(h__);\
+		(TO_V) = C__ * xsin(h__);\
+	} while (0)
+
+#define CIELUV_TO_CIELCHUV(FROM_L, FROM_U, FROM_V, TO_L, TO_C, TO_H, TO_ONE_REVOLUTION)\
+	do {\
+		TYPE u__ = (FROM_U), v__ = (FROM_V), rev__ = (TO_ONE_REVOLUTION), h__;\
+		(TO_L) = (FROM_L);\
+		(TO_C) = xsqrt(u__ * u__ + v__ * v__);\
+		h__ = xatan2(v__, u__) / PI2 * rev__;\
+		if (!WASDIV0(h__) && (h__ < 0))\
+			h__ += rev__;\
+		(TO_H) = h__;\
+	} while (0)
+
+#define SRGB_TO_YIQ(FROM_R, FROM_G, FROM_B, TO_Y, TO_I, TO_Q)\
+	TRANSFORM((TO_Y), (TO_I), (TO_Q), (FROM_R), (FROM_G), (FROM_B),\
+	           0.299, 0.587, 0.114,\
+	           0.59571613491277464191853141528554260730743408203125, /* (0.877 cos 33°)(1 - 0.299) - (0.492 sin 33°)(-0.299) */\
+	          -0.27445283783925644716106262421817518770694732666016, /* (0.877 cos 33°)(-0.587)    - (0.492 sin 33°)(-0.587) */\
+	          -0.32126329707351808373516632855171337723731994628906, /* (0.877 cos 33°)(-0.114)    - (0.492 sin 33°)(1 - 0.114) */\
+	           0.21145640212011795888713550084503367543220520019531, /* (0.877 sin 33°)(1 - 0.299) + (0.492 cos 33°)(-0.299) */\
+	          -0.52259104529161115593183239980135113000869750976562, /* (0.877 sin 33°)(-0.587)    + (0.492 cos 33°)(-0.587) */\
+	           0.31113464317149330806699936147197149693965911865234) /* (0.877 sin 33°)(-0.114)    + (0.492 cos 33°)(1 - 0.114) */
+
+#define CIEXYZ_TO_YIQ(FROM_X, FROM_Y, FROM_Z, TO_Y, TO_I, TO_Q)\
+	TRANSFORM((TO_Y), (TO_I), (TO_Q), (FROM_X), (FROM_Y), (FROM_Z),\
+	           0.40627729168038273499519164033699780702590942382812,\
+	           0.61835674212166968910509012857801280915737152099609,\
+	          -0.00414330221353725880462093300593551248311996459961,\
+	           2.17852787350219845308174626552499830722808837890625,\
+	          -1.36502666214454104753883711964590474963188171386719,\
+	          -0.64803574634025240541745915834326297044754028320312,\
+	           1.20905577682138853923277110880007967352867126464844,\
+	          -1.36890364998339797431015085749095305800437927246094,\
+	           0.20180559439597040016778350945969577878713607788086)
+
+#define SRGB_TO_YDBDR(FROM_R, FROM_G, FROM_B, TO_Y, TO_DB, TO_DR)\
+	TRANSFORM((TO_Y), (TO_DB), (TO_DR), (FROM_R), (FROM_G), (FROM_B),\
+	           0.299,  0.587, 0.114,\
+	          -0.450, -0.883, 1.333,\
+	          -1.333,  1.116, 0.217)
+
+#define YUV_TO_YDBDR(FROM_Y, FROM_U, FROM_V, TO_Y, TO_DB, TO_DR)\
+	do {\
+		(TO_Y)  = (FROM_Y);\
+		(TO_DB) = (FROM_U) *  D(3.069);\
+		(TO_DR) = (FROM_V) * D(-2.169);\
+	} while (0)
+
+#define CIEXYZ_TO_YDBDR(FROM_X, FROM_Y, FROM_Z, TO_Y, TO_DB, TO_DR)\
+	TRANSFORM((TO_Y), (TO_DB), (TO_DR), (FROM_X), (FROM_Y), (FROM_Z),\
+	           0.40627729168038273499519164033699780702590942382812,\
+	           0.61835674212166968910509012857801280915737152099609,\
+	          -0.00414330221353725880462093300593551248311996459961,\
+	          -0.52816561102614745237815441214479506015777587890625,\
+	          -1.23677481526212962315014465275453403592109680175781,\
+	           1.59692761635924940222253098909277468919754028320312,\
+	          -5.38914174974103143966885909321717917919158935546875,\
+	           4.09835630362728497999569299281574785709381103515625,\
+	           0.94033545560642795013706063400604762136936187744141)
+
+#define YDBDR_TO_YUV(FROM_Y, FROM_DB, FROM_DR, TO_Y, TO_U, TO_V)\
+	do {\
+		(TO_Y) = (FROM_Y);\
+		(TO_U) = (FROM_DB) /  D(3.069);\
+		(TO_V) = (FROM_DR) / D(-2.169);\
+	} while (0)
+
+#define SRGB_TO_YUV(FROM_R, FROM_G, FROM_B, TO_Y, TO_U, TO_V)\
+	TRANSFORM((TO_Y), (TO_U), (TO_V), (FROM_R), (FROM_G), (FROM_B),\
+	           0.29899999999999998800959133404830936342477798461914,\
+	           0.58699999999999996624922005139524117112159729003906,\
+	           0.11400000000000000410782519111307919956743717193604,\
+	          -0.14662756598240470062854967636667424812912940979004,\
+	          -0.28771586836102963635752871596196200698614120483398,\
+	           0.43434343434343436474165400795754976570606231689453,\
+	           0.61456892577224520035628074765554629266262054443359,\
+	          -0.51452282157676354490405401520547457039356231689453,\
+	          -0.10004610419548178035231700278018251992762088775635)
+
+#define CIEXYZ_TO_YUV(FROM_X, FROM_Y, FROM_Z, TO_Y, TO_U, TO_V)\
+	TRANSFORM((TO_Y), (TO_U), (TO_V), (FROM_X), (FROM_Y), (FROM_Z),\
+	           0.40627729168038273499519164033699780702590942382812,\
+	           0.61835674212166968910509012857801280915737152099609,\
+	          -0.00414330221353725880462093300593551248311996459961,\
+	          -0.17209697328971895746718701047939248383045196533203,\
+	          -0.40298951295605395239718404809536878019571304321289,\
+	           0.52034135430408912093014350830344483256340026855469,\
+	           2.48462044709130047692724474472925066947937011718750,\
+	          -1.88951420176453876997868519538315013051033020019531,\
+	          -0.43353409663735725798616726933687459677457809448242)
+
+#define SRGB_TO_YPBPR(FROM_R, FROM_G, FROM_B, TO_Y, TO_PB, TO_PR)\
+	do {\
+		TYPE R__ = (FROM_R), G__ = (FROM_G), B__ = (FROM_B);\
+		(TO_Y)  = R__ * D(0.2126) +\
+		          G__ * D(0.7152) +\
+		          B__ * D(0.0722);\
+		(TO_PB) = B__ - to->Y;\
+		(TO_PR) = R__ - to->Y;\
+	 } while (0)
+
+#define CIEXYZ_TO_YPBPR(FROM_X, FROM_Y, FROM_Z, TO_Y, TO_PB, TO_PR)\
+	TRANSFORM((TO_Y), (TO_PB), (TO_PR), (FROM_X), (FROM_Y), (FROM_Z),\
+	          -0.00028314209073960778378920011277841695118695497513,\
+	           1.00019821310075673892470149439759552478790283203125,\
+	           0.00006512054470741990286342115723527967929840087891,\
+	           0.05592664565509243568275365987574332393705844879150,\
+	          -1.20422439283671711685030913940863683819770812988281,\
+	           1.05716144717799576113748116767965257167816162109375,\
+	           3.24072939673847715269516811531502753496170043945312,\
+	          -2.53733297492083753610359053709544241428375244140625,\
+	          -0.49859531356743613805804216099204495549201965332031)
+
+#define SRGB_TO_YCGCO(FROM_R, FROM_G, FROM_B, TO_Y, TO_CG, TO_CO)\
+	TRANSFORM((TO_Y), (TO_CG), (TO_CO), (FROM_R), (FROM_G), (FROM_B),\
+		   0.25, 0.50,  0.25,\
+		  -0.25, 0.50, -0.25,\
+		   0.50, 0.,   -0.50)
+
+#define CIEXYZ_TO_YCGCO(FROM_X, FROM_Y, FROM_Z, TO_Y, TO_CG, TO_CO)\
+	TRANSFORM((TO_Y), (TO_CG), (TO_CO), (FROM_X), (FROM_Y), (FROM_Z),\
+	           0.33938913643068269188063368346774950623512268066406,\
+	           0.50271574450517486631895280879689380526542663574219,\
+	           0.16045211478220866574417868832824751734733581542969,\
+	          -1.30865574267536244335019546269904822111129760742188,\
+	           1.37329621528319534284889869013568386435508728027344,\
+	          -0.11889607256777862120955546743061859160661697387695,\
+	           1.59240137554169236544510113162687048316001892089844,\
+	          -0.66655429104206020962664069884340278804302215576172,\
+	          -0.77787838037271606062006412685150280594825744628906)
+
+#define CIEUVW_TO_CIE1960UCS(FROM_U, FROM_V, FROM_W, TO_U, TO_V, TO_Y, FROM_U0, FROM_V0)\
+	do {\
+		TYPE U__ = (FROM_U), V__ = (FROM_V), W__ = (FROM_W), Y__;\
+		Y__ = (W__ + 17) / 25;\
+		Y__ *= Y__ * Y__;\
+		W__ *= 13;\
+		(TO_U) = U__ / W__ + (FROM_U0);\
+		(TO_V) = V__ / W__ + (FROM_V0);\
+		(TO_Y) = Y__;\
+	} while (0)
+
+#define CIEXYZ_TO_CIE1960UCS(FROM_X, FROM_Y, FROM_Z, TO_U, TO_V, TO_Y)\
+	do {\
+		TYPE X__ = (FROM_X), Y__ = (FROM_Y), w__;\
+		w__ = X__ + 15 * Y__ + 3 * (FROM_Z);\
+		(TO_U) = 4 * X__ / w__;\
+		(TO_V) = 6 * Y__ / w__;\
+		(TO_Y) = Y__;\
+	} while (0)
+
+#define CIEUVW_TO_CIEUVW(FROM_U, FROM_V, FROM_W, TO_U, TO_V, TO_W, FROM_U0, FROM_V0, TO_U0, TO_V0)\
+	do {\
+		TYPE U__ = (FROM_U), V__ = (FROM_V), W__ = (FROM_W);\
+		TYPE w__ = W__ * 13;\
+		U__ += w__ * ((FROM_U0) - (TO_U0));\
+		V__ += w__ * ((FROM_V0) - (TO_V0));\
+		(TO_U) = U__;\
+		(TO_V) = V__;\
+		(TO_W) = W__;\
+	} while (0)
+
+#define CIE1960UCS_TO_CIEUVW(FROM_U, FROM_V, FROM_Y, TO_U, TO_V, TO_W, TO_U0, TO_V0)\
+	do {\
+		TYPE Y__ = 25 * xcbrt((FROM_Y)) - 17;\
+		TYPE w__ = Y__ * 13;\
+		to->U = w__ * ((FROM_U) - (TO_U0));\
+		to->V = w__ * ((FROM_V) - (TO_V0));\
+		to->W = Y__;\
+	} while (0)
+
+
+
+static inline TYPE
+srgb_encode(TYPE t)
+{
+	TYPE sign = 1;
+	if (t < 0) {
+		t = -t;
+		sign = -1;
+	}
+	t = t <= D(0.0031306684425217108) ? D(12.92) * t
+		: D(1.055) * xpow(t, 1 / D(2.4)) - D(0.055);
+	return t * sign;
+}
+
+static inline TYPE
+srgb_decode(TYPE t)
+{
+	TYPE sign = 1;
+	if (t < 0) {
+		t = -t;
+		sign = -1;
+	}
+	t = t <= D(0.0031306684425217108) * D(12.92) ? t / D(12.92)
+		: xpow((t + D(0.055)) / D(1.055), D(2.4));
+	return t * sign;
+}
+
+static inline TYPE
+cielab_finv(TYPE t)
+{
+	return (t > D(6.) / D(29.)) ? t * t * t : (t - D(4.) / D(29.)) * 108 / 841;
+}
+
+static inline TYPE
+cielab_f(TYPE t)
+{
+	return (t > D(216.) / D(24389.)) ? xcbrt(t) : t * D(841.) / D(108.) + D(4.) / D(29.);
+}