/* Convert the proposed ISO 12640-3 Reference Medium Gamut to Argyll Lab Gamut */
/* Copyright 2005 by Graeme W. Gill */
/* All rights reserved. */

#include <stdio.h>
#include <stdlib.h>
#include <stdarg.h>
#include <fcntl.h>
#include <string.h>
#include <math.h>
#include "aconfig.h"
#include "numlib.h"
#include "icc.h"
#include "gamut.h"
#if defined(__IBMC__) && defined(_M_IX86)
#include <float.h>
#endif

/* Gamut surface points in LCh space */
/* + black and white points */
#define NPOINTS (684+2+3)

double points[NPOINTS][3] = {
	{ 100,	0,	0 },			/* White */
	{ 3.1373,	0,	0 },		/* Black */
	{ 98.75, 30.75, 90 },		/* Three extra values along yellow ridge */
	{ 97.5, 61.5, 90 },
	{ 96.25, 92.25, 90 },
	{ 5,	10,	0 },
	{ 10,	25,	0 },
	{ 15,	38,	0 },
	{ 20,	50,	0 },
	{ 25,	62,	0 },
	{ 30,	73,	0 },
	{ 35,	80,	0 },
	{ 40,	86,	0 },
	{ 45,	89,	0 },
	{ 50,	86,	0 },
	{ 55,	85,	0 },
	{ 60,	81,	0 },
	{ 65,	73,	0 },
	{ 70,	65,	0 },
	{ 75,	55,	0 },
	{ 80,	45,	0 },
	{ 85,	34,	0 },
	{ 90,	23,	0 },
	{ 95,	12,	0 },
	{ 5,	10,	10 },
	{ 10,	24,	10 },
	{ 15,	37,	10 },
	{ 20,	50,	10 },
	{ 25,	61,	10 },
	{ 30,	71,	10 },
	{ 35,	80,	10 },
	{ 40,	88,	10 },
	{ 45,	90,	10 },
	{ 50,	89,	10 },
	{ 55,	85,	10 },
	{ 60,	81,	10 },
	{ 65,	73,	10 },
	{ 70,	65,	10 },
	{ 75,	55,	10 },
	{ 80,	45,	10 },
	{ 85,	34,	10 },
	{ 90,	22,	10 },
	{ 95,	11,	10 },
	{ 5,	10,	20 },
	{ 10,	22,	20 },
	{ 15,	36,	20 },
	{ 20,	50,	20 },
	{ 25,	62,	20 },
	{ 30,	73,	20 },
	{ 35,	80,	20 },
	{ 40,	90,	20 },
	{ 45,	93,	20 },
	{ 50,	91,	20 },
	{ 55,	86,	20 },
	{ 60,	83,	20 },
	{ 65,	73,	20 },
	{ 70,	70,	20 },
	{ 75,	55,	20 },
	{ 80,	50,	20 },
	{ 85,	35,	20 },
	{ 90,	25,	20 },
	{ 95,	11,	20 },
	{ 5,	9,	30 },
	{ 10,	22,	30 },
	{ 15,	34,	30 },
	{ 20,	48,	30 },
	{ 25,	61,	30 },
	{ 30,	74,	30 },
	{ 35,	86,	30 },
	{ 40,	98,	30 },
	{ 45,	99,	30 },
	{ 50,	101,	30 },
	{ 55,	95,	30 },
	{ 60,	90,	30 },
	{ 65,	81,	30 },
	{ 70,	73,	30 },
	{ 75,	65,	30 },
	{ 80,	52,	30 },
	{ 85,	35,	30 },
	{ 90,	30,	30 },
	{ 95,	12,	30 },
	{ 5,	8,	40 },
	{ 10,	20,	40 },
	{ 15,	33,	40 },
	{ 20,	46,	40 },
	{ 25,	60,	40 },
	{ 30,	73,	40 },
	{ 35,	83,	40 },
	{ 40,	93,	40 },
	{ 45,	96,	40 },
	{ 50,	101,	40 },
	{ 55,	99,	40 },
	{ 60,	97,	40 },
	{ 65,	90,	40 },
	{ 70,	83,	40 },
	{ 75,	71,	40 },
	{ 80,	60,	40 },
	{ 85,	44,	40 },
	{ 90,	30,	40 },
	{ 95,	15,	40 },
	{ 5,	8,	50 },
	{ 10,	20,	50 },
	{ 15,	30,	50 },
	{ 20,	42,	50 },
	{ 25,	54,	50 },
	{ 30,	66,	50 },
	{ 35,	77,	50 },
	{ 40,	88,	50 },
	{ 45,	95,	50 },
	{ 50,	99,	50 },
	{ 55,	101,	50 },
	{ 60,	100,	50 },
	{ 65,	97,	50 },
	{ 70,	92,	50 },
	{ 75,	83,	50 },
	{ 80,	72,	50 },
	{ 85,	55,	50 },
	{ 90,	37,	50 },
	{ 95,	19,	50 },
	{ 5,	7,	60 },
	{ 10,	17,	60 },
	{ 15,	26,	60 },
	{ 20,	36,	60 },
	{ 25,	46,	60 },
	{ 30,	55,	60 },
	{ 35,	64,	60 },
	{ 40,	75,	60 },
	{ 45,	80,	60 },
	{ 50,	90,	60 },
	{ 55,	95,	60 },
	{ 60,	100,	60 },
	{ 65,	101,	60 },
	{ 70,	102,	60 },
	{ 75,	98,	60 },
	{ 80,	90,	60 },
	{ 85,	72,	60 },
	{ 90,	50,	60 },
	{ 95,	25,	60 },
	{ 5,	6,	70 },
	{ 10,	15,	70 },
	{ 15,	25,	70 },
	{ 20,	32,	70 },
	{ 25,	42,	70 },
	{ 30,	50,	70 },
	{ 35,	58,	70 },
	{ 40,	68,	70 },
	{ 45,	75,	70 },
	{ 50,	82,	70 },
	{ 55,	90,	70 },
	{ 60,	95,	70 },
	{ 65,	100,	70 },
	{ 70,	104,	70 },
	{ 75,	107,	70 },
	{ 80,	109,	70 },
	{ 85,	100,	70 },
	{ 90,	74,	70 },
	{ 95,	34,	70 },
	{ 5,	5,	80 },
	{ 10,	13,	80 },
	{ 15,	23,	80 },
	{ 20,	30,	80 },
	{ 25,	37,	80 },
	{ 30,	47,	80 },
	{ 35,	55,	80 },
	{ 40,	62,	80 },
	{ 45,	70,	80 },
	{ 50,	77,	80 },
	{ 55,	85,	80 },
	{ 60,	90,	80 },
	{ 65,	95,	80 },
	{ 70,	103,	80 },
	{ 75,	106,	80 },
	{ 80,	110,	80 },
	{ 85,	113,	80 },
	{ 90,	110,	80 },
	{ 95,	70,	80 },
	{ 5,	5,	90 },
	{ 10,	13,	90 },
	{ 15,	21,	90 },
	{ 20,	30,	90 },
	{ 25,	36,	90 },
	{ 30,	45,	90 },
	{ 35,	55,	90 },
	{ 40,	60,	90 },
	{ 45,	64,	90 },
	{ 50,	73,	90 },
	{ 55,	80,	90 },
	{ 60,	88,	90 },
	{ 65,	95,	90 },
	{ 70,	100,	90 },
	{ 75,	106,	90 },
	{ 80,	111,	90 },
	{ 85,	117,	90 },
	{ 90,	120,	90 },
	{ 95,	123,	90 },
	{ 5,	5,	100 },
	{ 10,	14,	100 },
	{ 15,	21,	100 },
	{ 20,	30,	100 },
	{ 25,	35,	100 },
	{ 30,	45,	100 },
	{ 35,	53,	100 },
	{ 40,	60,	100 },
	{ 45,	65,	100 },
	{ 50,	73,	100 },
	{ 55,	80,	100 },
	{ 60,	87,	100 },
	{ 65,	94,	100 },
	{ 70,	100,	100 },
	{ 75,	106,	100 },
	{ 80,	109,	100 },
	{ 85,	113,	100 },
	{ 90,	108,	100 },
	{ 95,	90,	100 },
	{ 5,	6,	110 },
	{ 10,	14,	110 },
	{ 15,	22,	110 },
	{ 20,	30,	110 },
	{ 25,	37,	110 },
	{ 30,	46,	110 },
	{ 35,	52,	110 },
	{ 40,	61,	110 },
	{ 45,	67,	110 },
	{ 50,	75,	110 },
	{ 55,	82,	110 },
	{ 60,	89,	110 },
	{ 65,	96,	110 },
	{ 70,	100,	110 },
	{ 75,	103,	110 },
	{ 80,	106,	110 },
	{ 85,	107,	110 },
	{ 90,	102,	110 },
	{ 95,	75,	110 },
	{ 5,	6,	120 },
	{ 10,	15,	120 },
	{ 15,	22,	120 },
	{ 20,	31,	120 },
	{ 25,	40,	120 },
	{ 30,	49,	120 },
	{ 35,	57,	120 },
	{ 40,	65,	120 },
	{ 45,	73,	120 },
	{ 50,	80,	120 },
	{ 55,	87,	120 },
	{ 60,	93,	120 },
	{ 65,	98,	120 },
	{ 70,	100,	120 },
	{ 75,	101,	120 },
	{ 80,	96,	120 },
	{ 85,	88,	120 },
	{ 90,	73,	120 },
	{ 95,	50,	120 },
	{ 5,	6,	130 },
	{ 10,	15,	130 },
	{ 15,	24,	130 },
	{ 20,	34,	130 },
	{ 25,	44,	130 },
	{ 30,	53,	130 },
	{ 35,	63,	130 },
	{ 40,	72,	130 },
	{ 45,	80,	130 },
	{ 50,	87,	130 },
	{ 55,	93,	130 },
	{ 60,	97,	130 },
	{ 65,	101,	130 },
	{ 70,	99,	130 },
	{ 75,	94,	130 },
	{ 80,	84,	130 },
	{ 85,	72,	130 },
	{ 90,	55,	130 },
	{ 95,	30,	130 },
	{ 5,	7,	140 },
	{ 10,	18,	140 },
	{ 15,	27,	140 },
	{ 20,	37,	140 },
	{ 25,	47,	140 },
	{ 30,	57,	140 },
	{ 35,	67,	140 },
	{ 40,	77,	140 },
	{ 45,	85,	140 },
	{ 50,	95,	140 },
	{ 55,	98,	140 },
	{ 60,	101,	140 },
	{ 65,	97,	140 },
	{ 70,	93,	140 },
	{ 75,	85,	140 },
	{ 80,	75,	140 },
	{ 85,	61,	140 },
	{ 90,	42,	140 },
	{ 95,	21,	140 },
	{ 5,	7,	150 },
	{ 10,	18,	150 },
	{ 15,	29,	150 },
	{ 20,	40,	150 },
	{ 25,	51,	150 },
	{ 30,	61,	150 },
	{ 35,	71,	150 },
	{ 40,	81,	150 },
	{ 45,	92,	150 },
	{ 50,	97,	150 },
	{ 55,	99,	150 },
	{ 60,	96,	150 },
	{ 65,	91,	150 },
	{ 70,	85,	150 },
	{ 75,	76,	150 },
	{ 80,	66,	150 },
	{ 85,	52,	150 },
	{ 90,	37,	150 },
	{ 95,	19,	150 },
	{ 5,	7,	160 },
	{ 10,	20,	160 },
	{ 15,	30,	160 },
	{ 20,	42,	160 },
	{ 25,	53,	160 },
	{ 30,	66,	160 },
	{ 35,	79,	160 },
	{ 40,	92,	160 },
	{ 45,	96,	160 },
	{ 50,	99,	160 },
	{ 55,	97,	160 },
	{ 60,	90,	160 },
	{ 65,	87,	160 },
	{ 70,	79,	160 },
	{ 75,	69,	160 },
	{ 80,	59,	160 },
	{ 85,	46,	160 },
	{ 90,	32,	160 },
	{ 95,	19,	160 },
	{ 5,	8,	170 },
	{ 10,	20,	170 },
	{ 15,	30,	170 },
	{ 20,	41,	170 },
	{ 25,	52,	170 },
	{ 30,	63,	170 },
	{ 35,	73,	170 },
	{ 40,	83,	170 },
	{ 45,	91,	170 },
	{ 50,	96,	170 },
	{ 55,	93,	170 },
	{ 60,	89,	170 },
	{ 65,	82,	170 },
	{ 70,	75,	170 },
	{ 75,	65,	170 },
	{ 80,	55,	170 },
	{ 85,	42,	170 },
	{ 90,	29,	170 },
	{ 95,	15,	170 },
	{ 5,	8,	180 },
	{ 10,	20,	180 },
	{ 15,	29,	180 },
	{ 20,	40,	180 },
	{ 25,	51,	180 },
	{ 30,	62,	180 },
	{ 35,	72,	180 },
	{ 40,	81,	180 },
	{ 45,	86,	180 },
	{ 50,	92,	180 },
	{ 55,	90,	180 },
	{ 60,	86,	180 },
	{ 65,	79,	180 },
	{ 70,	71,	180 },
	{ 75,	61,	180 },
	{ 80,	51,	180 },
	{ 85,	38,	180 },
	{ 90,	25,	180 },
	{ 95,	13,	180 },
	{ 5,	8,	190 },
	{ 10,	20,	190 },
	{ 15,	29,	190 },
	{ 20,	40,	190 },
	{ 25,	49,	190 },
	{ 30,	60,	190 },
	{ 35,	68,	190 },
	{ 40,	76,	190 },
	{ 45,	81,	190 },
	{ 50,	87,	190 },
	{ 55,	85,	190 },
	{ 60,	82,	190 },
	{ 65,	75,	190 },
	{ 70,	69,	190 },
	{ 75,	60,	190 },
	{ 80,	50,	190 },
	{ 85,	38,	190 },
	{ 90,	26,	190 },
	{ 95,	13,	190 },
	{ 5,	8,	200 },
	{ 10,	20,	200 },
	{ 15,	28,	200 },
	{ 20,	38,	200 },
	{ 25,	47,	200 },
	{ 30,	56,	200 },
	{ 35,	63,	200 },
	{ 40,	70,	200 },
	{ 45,	75,	200 },
	{ 50,	80,	200 },
	{ 55,	78,	200 },
	{ 60,	78,	200 },
	{ 65,	72,	200 },
	{ 70,	66,	200 },
	{ 75,	58,	200 },
	{ 80,	49,	200 },
	{ 85,	38,	200 },
	{ 90,	27,	200 },
	{ 95,	14,	200 },
	{ 5,	8,	210 },
	{ 10,	20,	210 },
	{ 15,	28,	210 },
	{ 20,	37,	210 },
	{ 25,	45,	210 },
	{ 30,	53,	210 },
	{ 35,	60,	210 },
	{ 40,	66,	210 },
	{ 45,	72,	210 },
	{ 50,	79,	210 },
	{ 55,	80,	210 },
	{ 60,	76,	210 },
	{ 65,	69,	210 },
	{ 70,	64,	210 },
	{ 75,	57,	210 },
	{ 80,	49,	210 },
	{ 85,	38,	210 },
	{ 90,	27,	210 },
	{ 95,	14,	210 },
	{ 5,	8,	220 },
	{ 10,	20,	220 },
	{ 15,	27,	220 },
	{ 20,	37,	220 },
	{ 25,	44,	220 },
	{ 30,	51,	220 },
	{ 35,	57,	220 },
	{ 40,	64,	220 },
	{ 45,	70,	220 },
	{ 50,	76,	220 },
	{ 55,	73,	220 },
	{ 60,	71,	220 },
	{ 65,	68,	220 },
	{ 70,	63,	220 },
	{ 75,	56,	220 },
	{ 80,	48,	220 },
	{ 85,	38,	220 },
	{ 90,	27,	220 },
	{ 95,	14,	220 },
	{ 5,	8,	230 },
	{ 10,	20,	230 },
	{ 15,	28,	230 },
	{ 20,	38,	230 },
	{ 25,	45,	230 },
	{ 30,	52,	230 },
	{ 35,	57,	230 },
	{ 40,	65,	230 },
	{ 45,	69,	230 },
	{ 50,	75,	230 },
	{ 55,	72,	230 },
	{ 60,	71,	230 },
	{ 65,	66,	230 },
	{ 70,	61,	230 },
	{ 75,	54,	230 },
	{ 80,	46,	230 },
	{ 85,	36,	230 },
	{ 90,	26,	230 },
	{ 95,	13,	230 },
	{ 5,	9,	240 },
	{ 10,	20,	240 },
	{ 15,	29,	240 },
	{ 20,	40,	240 },
	{ 25,	48,	240 },
	{ 30,	55,	240 },
	{ 35,	61,	240 },
	{ 40,	66,	240 },
	{ 45,	71,	240 },
	{ 50,	74,	240 },
	{ 55,	70,	240 },
	{ 60,	66,	240 },
	{ 65,	61,	240 },
	{ 70,	56,	240 },
	{ 75,	49,	240 },
	{ 80,	41,	240 },
	{ 85,	32,	240 },
	{ 90,	23,	240 },
	{ 95,	12,	240 },
	{ 5,	11,	250 },
	{ 10,	22,	250 },
	{ 15,	32,	250 },
	{ 20,	42,	250 },
	{ 25,	50,	250 },
	{ 30,	59,	250 },
	{ 35,	65,	250 },
	{ 40,	70,	250 },
	{ 45,	73,	250 },
	{ 50,	70,	250 },
	{ 55,	68,	250 },
	{ 60,	62,	250 },
	{ 65,	58,	250 },
	{ 70,	52,	250 },
	{ 75,	45,	250 },
	{ 80,	38,	250 },
	{ 85,	29,	250 },
	{ 90,	21,	250 },
	{ 95,	11,	250 },
	{ 5,	14,	260 },
	{ 10,	27,	260 },
	{ 15,	38,	260 },
	{ 20,	48,	260 },
	{ 25,	56,	260 },
	{ 30,	63,	260 },
	{ 35,	67,	260 },
	{ 40,	72,	260 },
	{ 45,	73,	260 },
	{ 50,	69,	260 },
	{ 55,	65,	260 },
	{ 60,	60,	260 },
	{ 65,	56,	260 },
	{ 70,	50,	260 },
	{ 75,	43,	260 },
	{ 80,	35,	260 },
	{ 85,	27,	260 },
	{ 90,	20,	260 },
	{ 95,	10,	260 },
	{ 5,	16,	270 },
	{ 10,	32,	270 },
	{ 15,	44,	270 },
	{ 20,	55,	270 },
	{ 25,	62,	270 },
	{ 30,	70,	270 },
	{ 35,	75,	270 },
	{ 40,	74,	270 },
	{ 45,	72,	270 },
	{ 50,	68,	270 },
	{ 55,	66,	270 },
	{ 60,	59,	270 },
	{ 65,	54,	270 },
	{ 70,	48,	270 },
	{ 75,	41,	270 },
	{ 80,	33,	270 },
	{ 85,	26,	270 },
	{ 90,	18,	270 },
	{ 95,	9,	270 },
	{ 5,	21,	280 },
	{ 10,	42,	280 },
	{ 15,	55,	280 },
	{ 20,	68,	280 },
	{ 25,	73,	280 },
	{ 30,	81,	280 },
	{ 35,	80,	280 },
	{ 40,	78,	280 },
	{ 45,	76,	280 },
	{ 50,	70,	280 },
	{ 55,	67,	280 },
	{ 60,	60,	280 },
	{ 65,	55,	280 },
	{ 70,	49,	280 },
	{ 75,	41,	280 },
	{ 80,	33,	280 },
	{ 85,	26,	280 },
	{ 90,	18,	280 },
	{ 95,	9,	280 },
	{ 5,	26,	290 },
	{ 10,	52,	290 },
	{ 15,	66,	290 },
	{ 20,	83,	290 },
	{ 25,	84,	290 },
	{ 30,	89,	290 },
	{ 35,	87,	290 },
	{ 40,	84,	290 },
	{ 45,	80,	290 },
	{ 50,	75,	290 },
	{ 55,	69,	290 },
	{ 60,	63,	290 },
	{ 65,	57,	290 },
	{ 70,	50,	290 },
	{ 75,	42,	290 },
	{ 80,	34,	290 },
	{ 85,	26,	290 },
	{ 90,	18,	290 },
	{ 95,	9,	290 },
	{ 5,	25,	300 },
	{ 10,	65,	300 },
	{ 15,	79,	300 },
	{ 20,	95,	300 },
	{ 25,	93,	300 },
	{ 30,	92,	300 },
	{ 35,	90,	300 },
	{ 40,	87,	300 },
	{ 45,	85,	300 },
	{ 50,	77,	300 },
	{ 55,	73,	300 },
	{ 60,	66,	300 },
	{ 65,	59,	300 },
	{ 70,	52,	300 },
	{ 75,	44,	300 },
	{ 80,	36,	300 },
	{ 85,	28,	300 },
	{ 90,	19,	300 },
	{ 95,	10,	300 },
	{ 5,	20,	310 },
	{ 10,	50,	310 },
	{ 15,	67,	310 },
	{ 20,	91,	310 },
	{ 25,	97,	310 },
	{ 30,	100,	310 },
	{ 35,	98,	310 },
	{ 40,	95,	310 },
	{ 45,	90,	310 },
	{ 50,	84,	310 },
	{ 55,	77,	310 },
	{ 60,	70,	310 },
	{ 65,	63,	310 },
	{ 70,	55,	310 },
	{ 75,	47,	310 },
	{ 80,	39,	310 },
	{ 85,	30,	310 },
	{ 90,	20,	310 },
	{ 95,	10,	310 },
	{ 5,	18,	320 },
	{ 10,	41,	320 },
	{ 15,	60,	320 },
	{ 20,	79,	320 },
	{ 25,	90,	320 },
	{ 30,	102,	320 },
	{ 35,	101,	320 },
	{ 40,	98,	320 },
	{ 45,	94,	320 },
	{ 50,	89,	320 },
	{ 55,	83,	320 },
	{ 60,	76,	320 },
	{ 65,	68,	320 },
	{ 70,	60,	320 },
	{ 75,	51,	320 },
	{ 80,	42,	320 },
	{ 85,	32,	320 },
	{ 90,	21,	320 },
	{ 95,	11,	320 },
	{ 5,	16,	330 },
	{ 10,	35,	330 },
	{ 15,	53,	330 },
	{ 20,	71,	330 },
	{ 25,	82,	330 },
	{ 30,	91,	330 },
	{ 35,	99,	330 },
	{ 40,	104,	330 },
	{ 45,	102,	330 },
	{ 50,	98,	330 },
	{ 55,	90,	330 },
	{ 60,	84,	330 },
	{ 65,	76,	330 },
	{ 70,	67,	330 },
	{ 75,	57,	330 },
	{ 80,	47,	330 },
	{ 85,	36,	330 },
	{ 90,	24,	330 },
	{ 95,	12,	330 },
	{ 5,	14,	340 },
	{ 10,	31,	340 },
	{ 15,	45,	340 },
	{ 20,	61,	340 },
	{ 25,	71,	340 },
	{ 30,	82,	340 },
	{ 35,	91,	340 },
	{ 40,	101,	340 },
	{ 45,	103,	340 },
	{ 50,	99,	340 },
	{ 55,	95,	340 },
	{ 60,	89,	340 },
	{ 65,	80,	340 },
	{ 70,	71,	340 },
	{ 75,	61,	340 },
	{ 80,	50,	340 },
	{ 85,	38,	340 },
	{ 90,	26,	340 },
	{ 95,	13,	340 },
	{ 5,	12,	350 },
	{ 10,	28,	350 },
	{ 15,	41,	350 },
	{ 20,	54,	350 },
	{ 25,	65,	350 },
	{ 30,	76,	350 },
	{ 35,	83,	350 },
	{ 40,	93,	350 },
	{ 45,	95,	350 },
	{ 50,	92,	350 },
	{ 55,	90,	350 },
	{ 60,	85,	350 },
	{ 65,	77,	350 },
	{ 70,	68,	350 },
	{ 75,	58,	350 },
	{ 80,	48,	350 },
	{ 85,	36,	350 },
	{ 90,	24,	350 },
	{ 95,	12,	350 }
};

/* selected "cusp" points */
double cpoints[6][3] = {
	{ 50,	101,	36 },
	{ 95,	123,	90 },
	{ 60,	101,	140 },
	{ 55,	80,	210 },
	{ 20,	95,	300 },
	{ 40,	104,	330 },
};

/* Our idealised cusp point angles for Lab in degrees: */ 
/* 36.0,		Red */
/* 101.0,		Yellow */
/* 149.0,		Green */
/* 225.0,		Cyan */
/* 300.0,		Blue */
/* 337.0		Magenta */

#define GAMRES 1.0		/* Default surface resolution */

int
main(int argc, char *argv[]) {
	int i;
	char out_name[MAXNAMEL+1];			/* VRML output file */

	double gamres = GAMRES;				/* Surface resolution */
	gamut *gam;


#if defined(__IBMC__) && defined(_M_IX86)
	_control87(EM_UNDERFLOW, EM_UNDERFLOW);
#endif

	if (argc < 2) {
		strcpy(out_name, "RefMediumGamut.gam");
	} else {
		strncpy(out_name,argv[1],MAXNAMEL); out_name[MAXNAMEL] = '\000';
	}

	/* Creat a gamut surface */
	gam = new_gamut(gamres, 0, 0);

	/* For all the supplied surface points, */
	/* add them to the gamut */
	for (i = 0; i < NPOINTS; i++) {
		double lab[3];

		icmLCh2Lab(lab, points[i]);
//printf("~1 LCh %f %f %f -> Lab %f %f %f\n", points[i][0], points[i][1], points[i][2], lab[0], lab[1], lab[2]);
		gam->expand(gam, lab);
	}

	/* Add selected cusp points to set cusp values */
	gam->setcusps(gam, 0, NULL);
	for (i = 0; i < 6; i++) {
		double lab[3];

		icmLCh2Lab(lab, cpoints[i]);
//printf("~1 LCh %f %f %f -> Lab %f %f %f\n", points[i][0], points[i][1], points[i][2], lab[0], lab[1], lab[2]);
		gam->setcusps(gam, 1, lab);
	}
	gam->setcusps(gam, 2, NULL);

	if (gam->write_gam(gam, out_name))
		error ("write gamut failed on '%s'",out_name);

	gam->del(gam);

	return 0;
}