#if !defined INCLUDE_ACES_MATRICES
#define INCLUDE_ACES_MATRICES

#include "/include/utility/color.glsl"

const mat3 ap0_to_xyz = mat3(
	 0.9525523959,  0.0000000000,  0.0000936786,
	 0.3439664498,  0.7281660966, -0.0721325464,
	 0.0000000000,  0.0000000000,  1.0088251844
);
const mat3 xyz_to_ap0 = mat3(
	 1.0498110175,  0.0000000000, -0.0000974845,
	-0.4959030231,  1.3733130458,  0.0982400361,
	 0.0000000000,  0.0000000000,  0.9912520182
);

const mat3 ap1_to_xyz = mat3(
	 0.6624541811,  0.1340042065,  0.1561876870,
	 0.2722287168,  0.6740817658,  0.0536895174,
	-0.0055746495,  0.0040607335,  1.0103391003
);
const mat3 xyz_to_ap1 = mat3(
	 1.6410233797, -0.3248032942, -0.2364246952,
	-0.6636628587,  1.6153315917,  0.0167563477,
	 0.0117218943, -0.0082844420,  0.9883948585
);

// Bradford chromatic adaptation matrix from d60 (ACES white point) to d65 (rec. 709 white point)
const mat3 d60_to_d65 = mat3(
	 0.9872240000, -0.0061132700,  0.0159533000,
	-0.0075983600,  1.0018600000,  0.0053300200,
	 0.0030725700, -0.0050959500,  1.0816800000
);
const mat3 d65_to_d60 = mat3(
	 1.0130349240,  0.0061053089, -0.0149709632,
	 0.0076982300,  0.9981648318, -0.0050320341,
	-0.0028413125,  0.0046851556,  0.9245066529
);

const mat3 ap0_to_ap1 = ap0_to_xyz * xyz_to_ap1;
const mat3 ap1_to_ap0 = ap1_to_xyz * xyz_to_ap0;

const mat3 rec709_to_ap0  = rec709_to_xyz * d65_to_d60 * xyz_to_ap0;
const mat3 ap0_to_rec709  = ap0_to_xyz * d60_to_d65 * xyz_to_rec709;
const mat3 rec709_to_ap1  = rec709_to_xyz * d65_to_d60 * xyz_to_ap1;
const mat3 ap1_to_rec709  = ap1_to_xyz * d60_to_d65 * xyz_to_rec709;
const mat3 rec2020_to_ap0 = rec2020_to_xyz * d65_to_d60 * xyz_to_ap0;
const mat3 ap0_to_rec2020 = ap0_to_xyz * d60_to_d65 * xyz_to_rec2020;
const mat3 rec2020_to_ap1 = rec2020_to_xyz * d65_to_d60 * xyz_to_ap1;
const mat3 ap1_to_rec2020 = ap1_to_xyz * d60_to_d65 * xyz_to_rec2020;

// Used for quantities that are multiplied by an illuminant, like the diffuse albedo. These take
// the white point of the light source
const mat3 rec709_to_ap1_unlit = rec709_to_xyz * xyz_to_ap1;
const mat3 rec2020_to_ap1_unlit = rec2020_to_xyz * xyz_to_ap1;

#endif // INCLUDE_ACES_MATRICES