/* See LICENSE file for copyright and license details. */
#include "libcolour.h"

#include <math.h>
#include <stdio.h>
#include <string.h>


#undef libcolour_unmarshal

#define libcolour_colour_t  libcolour_colour_lf_t
#define libcolour_rgb_t     libcolour_rgb_lf_t
#define libcolour_unmarshal libcolour_unmarshal_lf


#ifndef SKIP_2CONVERT
static int
test_2convert(libcolour_colour_t *c1, libcolour_colour_t *c2, libcolour_colour_t *c3)
{
	double ch1, ch2, ch3;
	if (libcolour_convert(c1, c2))
		return -1;
	if (libcolour_convert(c2, c3))
		return -1;
	if ((c1->model == LIBCOLOUR_CIELCHUV) &&
	    !isnan(c1->cielchuv.h) && !isnan(c3->cielchuv.h) &&
	    !isinf(c1->cielchuv.h) && !isinf(c3->cielchuv.h) &&
	    (c3->cielchuv.h / c3->cielchuv.one_revolution > 0.99999))
		c3->cielchuv.h = 0;
	ch1 = c1->srgb.R - c3->srgb.R;
	ch2 = c1->srgb.G - c3->srgb.G;
	ch3 = c1->srgb.B - c3->srgb.B;
	ch1 *= ch1;
	ch2 *= ch2;
	ch3 *= ch3;
	switch (c1->model) {
	case LIBCOLOUR_CIEXYY:
		if (c1->ciexyy.Y == 0)
			ch1 = ch2 = 0;
		if (c1->ciexyy.y == 0)
			ch1 = ch2 = ch3 = 0;
		break;
	case LIBCOLOUR_CIEXYZ:
		if (c1->ciexyz.Y == 0)
			ch1 = ch3 = 0;
		break;
	case LIBCOLOUR_CIE1960UCS:
		if (c1->cie1960ucs.v == 0)
			ch1 = ch2 = ch3 = 0;
		if (c1->cie1960ucs.Y == 0)
			ch1 = ch2 = 0;
		break;
	case LIBCOLOUR_CIEUVW:
		if (c1->cieuvw.W == 0)
			ch1 = ch2 = 0;
		break;
	case LIBCOLOUR_CIELCHUV:
		if (c1->cielchuv.C == 0)
			ch1 = ch3 = 0;
		break;
	case LIBCOLOUR_CIELAB:
		if (c1->cielchuv.L == 0)
			ch2 = ch3 = 0;
		break;
	default:
		break;
	}
	if (ch1 > 0.000001 ||
	    ch2 > 0.000001 ||
	    ch3 > 0.000001)
		return 0;
	return 1;
}
  

static int
test_2convert_11(libcolour_colour_t *c1, libcolour_model_t model)
{
	libcolour_colour_t c2, c3;
	int r1, r2;
	c3 = *c1;
	c2.model = model;
	switch (model) {
	case LIBCOLOUR_RGB:
		c2.rgb.with_transfer = 0;
		if (libcolour_get_rgb_colour_space(&c2.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB) < 0)
			return -1;
		return test_2convert(c1, &c2, &c3);
	case LIBCOLOUR_SRGB:
		c2.srgb.with_transfer = 0;
		if (r1 = test_2convert(c1, &c2, &c3), r1 < 0)
			return -1;
		c2.srgb.with_transfer = 1;
		if (r2 = test_2convert(c1, &c2, &c3), r2 < 0)
			return -1;
		return r1 & r2;
	case LIBCOLOUR_CIELCHUV:
		c2.cielchuv.one_revolution = 360.;
		/* fall though */
	case LIBCOLOUR_CIELUV:
		c2.cieluv.white.model = LIBCOLOUR_CIEXYZ;
		c2.cieluv.white.X = 1.0294;
		c2.cieluv.white.Y = 1;
		c2.cieluv.white.Z = 0.9118;
		return test_2convert(c1, &c2, &c3);
	case LIBCOLOUR_CIEUVW:
		c2.cieuvw.u0 = 0.37;
		c2.cieuvw.v0 = 0.30;
		return test_2convert(c1, &c2, &c3);
	default:
		return test_2convert(c1, &c2, &c3);
	}
}


static int
test_2convert_1n(libcolour_model_t model, const char *model_name, double ch1, double ch2, double ch3)
{
	libcolour_colour_t c1;
	int run, r, rc = 1;

	c1.model = model;
	switch (model) {
	case LIBCOLOUR_CIELAB:
		c1.cielab.L = ch1 * 100, c1.cielab.a = ch2 * 100, c1.cielab.b = ch3 * 100;
		break;
	case LIBCOLOUR_CIELCHUV:
		if (ch3 > 0.9999)
			return 1;
		/* fall though */
	default:
		c1.srgb.R = ch1, c1.srgb.G = ch2, c1.srgb.B = ch3;
		break;
	}

	for (run = 0;; run++) {
		switch (model) {
		case LIBCOLOUR_RGB:
			if (run == 2)
				return rc;
			c1.rgb.with_transfer = run;
			if (libcolour_get_rgb_colour_space(&c1.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB) < 0)
				return -1;
			break;
		case LIBCOLOUR_SRGB:
			if (run == 2)
				return rc;
			c1.srgb.with_transfer = run;
			break;
		case LIBCOLOUR_CIELCHUV:
			c1.cielchuv.one_revolution = run < 2 ? 360. : 2 * 3.14159265358979323846;
			c1.cielchuv.h *= c1.cielchuv.one_revolution;
			/* fall through */
		case LIBCOLOUR_CIELUV:
			c1.cieluv.white.model = LIBCOLOUR_CIEXYZ;
			if (run == 0) {
				c1.cieluv.white.X = 1.0294;
				c1.cieluv.white.Y = 1;
				c1.cieluv.white.Z = 0.9118;
			} else if (run == 1) {
				c1.cieluv.white.X = 1.03;
				c1.cieluv.white.Y = 0.8;
				c1.cieluv.white.Z = 0.92;
			} else if (run == 2 && model == LIBCOLOUR_CIELCHUV) {
				c1.cieluv.white.X = 1.0294;
				c1.cieluv.white.Y = 1;
				c1.cieluv.white.Z = 0.9118;
			} else if (run == 3 && model == LIBCOLOUR_CIELCHUV) {
				c1.cieluv.white.X = 1.03;
				c1.cieluv.white.Y = 0.8;
				c1.cieluv.white.Z = 0.92;
			} else {
				return rc;
			}
			break;
		case LIBCOLOUR_CIEUVW:
			if (run == 0) {
				c1.cieuvw.u0 = 0.37;
				c1.cieuvw.v0 = 0.30;
			} else if (run == 1) {
				c1.cieuvw.u0 = 0.47;
				c1.cieuvw.v0 = 0.31;
			} else {
				return rc;
			}
			break;
		default:
			if (run == 1)
				return rc;
			break;
		}

#define X(ENUM, T, N)\
		r = test_2convert_11(&c1, ENUM);\
		if (r < 0)\
			return -1;\
		if (!r)\
			printf("%s -> %s -> %s failed at run %i with (%lf, %lf, %lf)\n",\
			       model_name, #ENUM, model_name, run, ch1, ch2, ch3), rc = 0;
		LIBCOLOUR_LIST_MODELS(X,);
#undef X
	}
}


static int
test_2convert_nm(double ch1, double ch2, double ch3)
{
	int r, rc = 1;

#define X(ENUM, T, N)\
	r = test_2convert_1n(ENUM, #ENUM, ch1, ch2, ch3);\
	if (r < 0)\
		return -1;\
	if (!r)\
		rc = 0;
	LIBCOLOUR_LIST_MODELS(X,);
#undef X

	return rc;
}


static int
test_2convert_nm_all(void)
{
#define N 11
	int r, rc = 1;
	double ch1, ch2, ch3;
	int i, j, k;

	for (i = 0; i < N; i++) {
		for (j = 0; j < N; j++) {
			for (k = 0; k < N; k++) {
				ch1 = ((i + N / 2) % N) / (N - 1);
				ch2 = ((j + N / 2) % N) / (N - 1);
				ch3 = ((k + N / 2) % N) / (N - 1);

				if (ch1 > 0.999) ch1 = 1;
				if (ch2 > 0.999) ch2 = 1;
				if (ch3 > 0.999) ch3 = 1;

				r = test_2convert_nm(ch1, ch2, ch3);
				if (r < 0)
					return -1;
				if (!r)
					rc = 0;
			}
		}
	}

	return rc;
#undef N
}
#endif


static inline int
ftest(double a, double b, double l)
{
	return a - l <= b && b <= a + l;
}


#define test_rgb(COLOUR_SPACE, ...) test_rgb_(COLOUR_SPACE, #COLOUR_SPACE, __VA_ARGS__)
static int
test_rgb_(enum libcolour_rgb_colour_space colour_space, const char *name,
	  double r, double g, double b, double x, double y, double z,
	  double xx, double yy, double zz)
{
	libcolour_colour_t c1, c2;

	c1.model = LIBCOLOUR_RGB;
	c2.model = LIBCOLOUR_CIEXYZ;
	c1.rgb.with_transfer = 0;
	c1.rgb.R = r, c1.rgb.G = g, c1.rgb.B = b;
	if (libcolour_get_rgb_colour_space(&c1.rgb, colour_space) ||
	    libcolour_convert(&c1, &c2)) {
		printf("%s failed without gamma\n", name);
		return 0;
	} else if (!ftest(c2.ciexyz.X, x, 0.0001) ||
		   !ftest(c2.ciexyz.Y, y, 0.0001) ||
		   !ftest(c2.ciexyz.Z, z, 0.0001)) {
		printf("%s failed without gamma\n", name);
		return 0;
	}

	c1.model = LIBCOLOUR_RGB;
	c2.model = LIBCOLOUR_CIEXYZ;
	c1.rgb.with_transfer = 1;
	c1.rgb.R = r, c1.rgb.G = g, c1.rgb.B = b;
	if (libcolour_get_rgb_colour_space(&c1.rgb, colour_space) ||
	    libcolour_convert(&c1, &c2)) {
		printf("%s failed with gamma\n", name);
		return 0;
	} else if (!ftest(c2.ciexyz.X, xx, 0.0001) ||
		   !ftest(c2.ciexyz.Y, yy, 0.0001) ||
		   !ftest(c2.ciexyz.Z, zz, 0.0001)) {
		printf("%s failed with gamma\n", name);
		return 0;
	}

	if (libcolour_convert(&c2, &c1) ||
	    !ftest(c1.rgb.R, r, 0.0000001) ||
	    !ftest(c1.rgb.G, g, 0.0000001) ||
	    !ftest(c1.rgb.B, b, 0.0000001)) {
		printf("%s failed, encode\n", name);
		return 0;
	}

	return 1;
}


static int
test_1convert(libcolour_colour_t *c1, libcolour_model_t model, const char *name1,
	      const char *name2, double ch1, double ch2, double ch3, double d1,
	      double d2, double d3)
{
	libcolour_colour_t c2, c3;
	int rc = 1;
	c3 = *c1;
	c2.model = model;
	switch (model) {
	case LIBCOLOUR_CIELCHUV:
		c2.cielchuv.one_revolution = 360;
		/* fall-through */
	case LIBCOLOUR_CIELUV:
		c2.cieluv.white.X = 0.9504705586542831;
		c2.cieluv.white.Y = 1;
		c2.cieluv.white.Z = 1.0888287363958837;
		break;
	default:
		break;
	}
	if (libcolour_convert(c1, &c2))
		return -1;
	if (!ftest(c2.srgb.R, ch1, d1) ||
	    !ftest(c2.srgb.G, ch2, d2) ||
	    !ftest(c2.srgb.B, ch3, d3)) {
		printf("%s -> %s failed\n", name1, name2), rc = 0;
		rc = 0;
	}
	if (libcolour_convert(&c2, &c3))
		return -1;
	if (!ftest(c3.srgb.R, c1->srgb.R, 0.00000001) ||
	    !ftest(c3.srgb.G, c1->srgb.G, 0.00000001) ||
	    !ftest(c3.srgb.B, c1->srgb.B, 0.00000001)) {
		printf("%s -> %s failed\n", name2, name1), rc = 0;
		rc = 0;
	}
	return rc;
}


/**
 * Test libcolour
 * 
 * @return  0: All tests passed
 *          1: At least one test fail
 *          2: An error occurred
 */
int
main(int argc, char *argv[])
{
	int r, rc = 0;
	libcolour_colour_t c1, c2, c3, c4;
	double t1, t2;
	size_t n;
	char buf[sizeof(int) + sizeof(libcolour_rgb_t)];

#ifndef SKIP_2CONVERT
	r = test_2convert_nm_all();
	if (r < 0)
		goto fail;
	if (!r)
		rc = 1;
#else
	(void) r;
#endif

	c1.model = LIBCOLOUR_SRGB;
	c1.srgb.R = 0.02, c1.srgb.G = 0.5, c1.srgb.B = 0.9;
	c1.srgb.with_transfer = 0;
	c2.model = LIBCOLOUR_RGB;
	c2.rgb.with_transfer = 0;
	if (libcolour_get_rgb_colour_space(&c2.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB)) {
		printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
		goto colour_spaces_done;
	}
	if (c2.rgb.white_r != 1 || c2.rgb.white_g != 1 || c2.rgb.white_b != 1 ||
	    c2.rgb.red.model != LIBCOLOUR_CIEXYY || c2.rgb.green.model != LIBCOLOUR_CIEXYY ||
	    c2.rgb.blue.model != LIBCOLOUR_CIEXYY || c2.rgb.white.model != LIBCOLOUR_CIEXYY ||
	    c2.rgb.colour_space != LIBCOLOUR_RGB_COLOUR_SPACE_SRGB || c2.rgb.gamma != 2.4 ||
	    c2.rgb.encoding_type != LIBCOLOUR_ENCODING_TYPE_REGULAR || c2.rgb.slope != 12.92
	    || c2.rgb.offset != 0.055 || c2.rgb.transitioninv != c2.rgb.transition * c2.rgb.slope ||
	    !ftest(c2.rgb.transition, 0.00313067, 0.00000001)) {
		printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
		goto colour_spaces_done;
	}
	t1 = c2.rgb.transition * c2.rgb.slope;
	t2 = (1 + c2.rgb.offset) * pow(c2.rgb.transition, 1 / c2.rgb.gamma) - c2.rgb.offset;
	if (!ftest(t1, t2, 0.0000001) || c2.rgb.white.Y != 1 ||
	    c2.rgb.red.x != 0.64 || c2.rgb.red.y != 0.33 || !ftest(c2.rgb.red.Y, 0.21265, 0.00005) ||
	    c2.rgb.green.x != 0.30 || c2.rgb.green.y != 0.60 || !ftest(c2.rgb.green.Y, 0.71515, 0.00005) ||
	    c2.rgb.blue.x != 0.15 || c2.rgb.blue.y != 0.06 || !ftest(c2.rgb.blue.Y, 0.07215, 0.00005) ||
	    !ftest(c2.rgb.white.x, 0.312727, 0.000001) || !ftest(c2.rgb.white.y, 0.329023, 0.000001) ||
	    !ftest(c2.rgb.red.Y + c2.rgb.green.Y + c2.rgb.blue.Y, 1, 0.00000001)) {
		printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
		goto colour_spaces_done;
	}
	if (libcolour_convert(&c1, &c2) ||
	    !ftest(c1.srgb.R, c2.rgb.R, 0.0000001) ||
	    !ftest(c1.srgb.G, c2.rgb.G, 0.0000001) ||
	    !ftest(c1.srgb.B, c2.rgb.B, 0.0000001)) {
		printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
		goto colour_spaces_done;
	}
	if (libcolour_convert(&c2, &c1) ||
	    !ftest(c1.srgb.R, c2.rgb.R, 0.0000001) ||
	    !ftest(c1.srgb.G, c2.rgb.G, 0.0000001) ||
	    !ftest(c1.srgb.B, c2.rgb.B, 0.0000001)) {
		printf("LIBCOLOUR_RGB_COLOUR_SPACE_SRGB failed\n"), rc = 0;
		goto colour_spaces_done;
	}

	c1.model = c2.model = LIBCOLOUR_CIELAB;
	c1.cielab.L =  30, c2.cielab.L =  80;
	c1.cielab.a =  40, c2.cielab.a = -10;
	c1.cielab.b = -50, c2.cielab.b = -40;
	if (libcolour_delta_e(&c1, &c2, &t1) || !ftest(t1, 71.4142842854285, 0.0000001))
		printf("libcolour_delta_e failed\n"), rc = 0;

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_ADOBE_RGB, 0.3, 0.2, 0.9,
		       0.379497, 0.282430, 0.914245, 0.195429, 0.098925, 0.790020);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_APPLE_RGB, 0.3, 0.2, 0.9,
		       0.364212, 0.282789, 0.866008, 0.221570, 0.134028, 0.773782);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_BEST_RGB, 0.3, 0.2, 0.9,
		       0.345007, 0.246779, 0.736029, 0.151406, 0.064655, 0.647212);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_BETA_RGB, 0.3, 0.2, 0.9,
		       0.342837, 0.253386, 0.714198, 0.146437, 0.066824, 0.623382);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_BRUCE_RGB, 0.3, 0.2, 0.9,
		       0.368857, 0.276910, 0.915272, 0.191184, 0.096699, 0.791486);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_CIE_RGB, 0.3, 0.2, 0.9,
		       0.389293, 0.225188, 0.892857, 0.202678, 0.044608, 0.785312);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_COLORMATCH_RGB, 0.3, 0.2, 0.9,
		       0.337557, 0.274378, 0.651989, 0.186858, 0.123210, 0.581381);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_DON_RGB_4, 0.3, 0.2, 0.9,
		       0.344990, 0.251411, 0.727872, 0.150504, 0.066348, 0.638058);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB, 0.3, 0.2, 0.9,
		       0.353021, 0.286400, 0.695202, 0.196733, 0.134157, 0.630279);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB_V2, 0.3, 0.2, 0.9,
		       0.353021, 0.286400, 0.695202, 0.149594, 0.097238, 0.579927);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_EKTA_SPACE_PS5, 0.3, 0.2, 0.9,
		       0.320377, 0.229160, 0.713291, 0.127134, 0.043253, 0.622392);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_NTSC_RGB, 0.3, 0.2, 0.9,
		       0.397081, 0.310031, 1.017821, 0.235650, 0.156675, 0.906707);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_PAL_SECAM_RGB, 0.3, 0.2, 0.9,
		       0.357972, 0.272130, 0.877151, 0.208580, 0.120322, 0.769056);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_PROPHOTO_RGB, 0.3, 0.2, 0.9,
		       0.294559, 0.228864, 0.742689, 0.124735, 0.072340, 0.682655);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_C_RGB, 0.3, 0.2, 0.9,
		       0.363595, 0.281822, 0.890350, 0.216773, 0.131311, 0.783348);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_WIDE_GAMUT_RGB, 0.3, 0.2, 0.9,
		       0.367485, 0.237631, 0.706442, 0.170318, 0.052665, 0.614915);

	rc &= test_rgb(LIBCOLOUR_RGB_COLOUR_SPACE_LIGHTROOM_RGB, 0.3, 0.2, 0.9,
		       0.294559, 0.228864, 0.742689, 0.294559, 0.228864, 0.742689);

	/*
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_D65
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_DCI_P3_THEATER
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_625_LINE
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_601_525_LINE
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_709
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2020
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_PQ
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_PQ
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_PQ
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_EOTF_HLG
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OOTF_HLG
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_ITU_R_BT_2100_OETF_HLG
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_SGI_RGB
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_SMPTE_240M_RGB
  
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MEASUREMENTS
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_MATRIX
	  TODO test LIBCOLOUR_RGB_COLOUR_SPACE_CUSTOM_FROM_INV_MATRIX
	*/

	/* TODO test transfer functions more rigorously */

	c1.model = c2.model = LIBCOLOUR_RGB;
	c3.model = c4.model = LIBCOLOUR_SRGB;
	c1.rgb.with_transfer = c3.srgb.with_transfer = 0;
	c2.rgb.with_transfer = c4.srgb.with_transfer = 1;
	if (libcolour_get_rgb_colour_space(&c1.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB))
		goto fail;
	if (libcolour_get_rgb_colour_space(&c2.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_SRGB))
		goto fail;
	for (t1 = 0; t1 <= 1; t1 += 0.001) {
		c1.rgb.R = c3.srgb.R = t1 - 1;
		c1.rgb.G = c3.srgb.G = t1;
		c1.rgb.B = c3.srgb.B = t1 + 1;
		if (libcolour_convert(&c1, &c2))
			goto fail;
		if (libcolour_convert(&c3, &c4))
			goto fail;
		if (!ftest(c2.rgb.R, c4.srgb.R, 0.00000000001) ||
		    !ftest(c2.rgb.G, c4.srgb.G, 0.00000000001) ||
		    !ftest(c2.rgb.B, c4.srgb.B, 0.00000000001) ||
		    !ftest(libcolour_srgb_encode(c1.rgb.R), c2.srgb.R, 0.00000000001) ||
		    !ftest(libcolour_srgb_encode(c1.rgb.G), c2.srgb.G, 0.00000000001) ||
		    !ftest(libcolour_srgb_encode(c1.rgb.B), c2.srgb.B, 0.00000000001)) {
			printf("libcolour_srgb_encode failed\n"), rc = 0;
			goto colour_spaces_done;
		}
		if (libcolour_convert(&c2, &c1))
			goto fail;
		if (libcolour_convert(&c4, &c3))
			goto fail;
		if (!ftest(c1.rgb.R, c3.srgb.R, 0.00000000001) ||
		    !ftest(c1.rgb.G, c3.srgb.G, 0.00000000001) ||
		    !ftest(c1.rgb.B, c3.srgb.B, 0.00000000001) ||
		    !ftest(libcolour_srgb_decode(c2.rgb.R), c1.srgb.R, 0.00000000001) ||
		    !ftest(libcolour_srgb_decode(c2.rgb.G), c1.srgb.G, 0.00000000001) ||
		    !ftest(libcolour_srgb_decode(c2.rgb.B), c1.srgb.B, 0.00000000001)) {
			printf("libcolour_srgb_decode failed\n"), rc = 0;
			goto colour_spaces_done;
		}
	}

	c1.rgb.encoding_type = LIBCOLOUR_ENCODING_TYPE_SIMPLE;
	c1.rgb.gamma = 2.2;
	c1.rgb.with_transfer = 0;
	c1.rgb.R = 0.25;
	c1.rgb.G = 0.5;
	c1.rgb.B = 0.75;
	c4 = c3 = c2 = c1;
	c2.rgb.with_transfer = c3.srgb.with_transfer = 1;
	c3.rgb.gamma = c4.rgb.gamma = 1.8;
	if (libcolour_convert(&c1, &c2))
		goto fail;
	if (libcolour_convert(&c2, &c3))
		goto fail;
	if (libcolour_convert(&c3, &c4))
		goto fail;
	if (!ftest(c1.rgb.R, c4.srgb.R, 0.0000001) ||
	    !ftest(c1.rgb.G, c4.srgb.G, 0.0000001) ||
	    !ftest(c1.rgb.B, c4.srgb.B, 0.0000001)) {
		printf("libcolour_convert failed to convert between two transfer functions\n"), rc = 0;
		goto colour_spaces_done;
	}

	if (libcolour_convert(&c1, &c4))
		goto fail;
	if (!ftest(c1.rgb.R, c4.srgb.R, 0.0000001) ||
	    !ftest(c1.rgb.G, c4.srgb.G, 0.0000001) ||
	    !ftest(c1.rgb.B, c4.srgb.B, 0.0000001)) {
		printf("libcolour_convert failed to convert when two different transfer functions are not applied\n"), rc = 0;
		goto colour_spaces_done;
	}

#define TEST(TO, CH1, CH2, CH3, D1, D2, D3)\
	r = test_1convert(&c1, TO, "LIBCOLOUR_SRGB", #TO, CH1, CH2, CH3, D1, D2, D3);\
	if (r < 0)\
		goto fail;\
	if (!r)\
		rc = 1;

	c1.model = LIBCOLOUR_SRGB;
	c1.srgb.with_transfer = 0;
	c1.srgb.R = 0.3, c1.srgb.G = 0.2, c1.srgb.B = 0.7;
	TEST(LIBCOLOUR_CIEXYZ, 0.321558, 0.257355, 0.694851, 0.000001, 0.000001, 0.000001);
	TEST(LIBCOLOUR_CIEXYY, 0.252447, 0.202043, 0.257355, 0.000001, 0.000001, 0.000001);
	TEST(LIBCOLOUR_CIELAB, 57.7851, 30.3599, -44.9740, 0.0001, 0.0001, 0.0001);
	TEST(LIBCOLOUR_CIELUV, 57.7851, 5.5723, -74.1571, 0.0001, 0.0001, 0.0001);
	TEST(LIBCOLOUR_CIELCHUV, 57.7851, 74.3662, 274.2972, 0.0001, 0.0001, 0.0001);
	/*
	 * TODO convert test from sRGB to:
	 *     LChab  57.7851, 54.2622, 304.0215
	 *     YDbDr
	 *     YPbPr
	 *     YCgCo
	 *     CIE 1960 UCS
	 *     CIEUVW
	 */
	c1.srgb.with_transfer = 0;
	c1.srgb.R = 32 / 255., c1.srgb.G = 65 / 255., c1.srgb.B = 32 / 255.;
	TEST(LIBCOLOUR_YIQ, 0.201, -0.035, -0.068, 0.001, 0.001, 0.001);
	TEST(LIBCOLOUR_YUV, 0.201, -0.037, -0.067, 0.001, 0.001, 0.001);
	/* XXX YCbCr(60, 118, 116) */
	/* XXX HSV(120, 0.508, 0.255) */
	/* XXX HSI(120, 0.256, 0.169) */
	/* XXX HSL(120, 0.340, 0.190) */

	if (rc)
		goto colour_spaces_done;

#undef TEST
colour_spaces_done:

	memset(&c1, 0, sizeof(c1));
	memset(&c2, 0, sizeof(c2));
	c1.model = LIBCOLOUR_RGB;
	if (libcolour_get_rgb_colour_space(&c1.rgb, LIBCOLOUR_RGB_COLOUR_SPACE_ECI_RGB_V2))
		goto fail;
	c3 = c1;
	c1.rgb.to_encoded_red = NULL;
	c1.rgb.to_encoded_green = NULL;
	c1.rgb.to_encoded_blue = NULL;
	c1.rgb.to_decoded_red = NULL;
	c1.rgb.to_decoded_green = NULL;
	c1.rgb.to_decoded_blue = NULL;
	if (libcolour_marshal(&c1, NULL) > sizeof(buf)) {
		printf("libcolour_marshal failed\n"), rc = 0;
		goto marshal_done;
	}
	n = libcolour_marshal(&c1, buf);
	if (n > sizeof(buf)) {
		printf("libcolour_marshal failed\n"), rc = 0;
		goto marshal_done;
	}
	if (libcolour_unmarshal(NULL, buf) != n) {
		printf("libcolour_unmarshal failed\n"), rc = 0;
		goto marshal_done;
	}
	if (libcolour_unmarshal(&c2, buf) != n) {
		printf("libcolour_unmarshal failed\n"), rc = 0;
		goto marshal_done;
	}
	if (memcmp(&c2, &c3, sizeof(c2))) {
		printf("libcolour_(un)marshal failed\n"), rc = 0;
		goto marshal_done;
	}
 marshal_done:

	return rc;
 fail:
	perror(*argv);
	return 2;
	(void) argc;
}