{
	Compute color rendering index based on CIE method.
	Used by scripts illumcal, reflcal.
}
sq(x) : x*x;

upi(i) = (10.872 + .404*c_r/c_k*ci(i) - 4*d_r/d_k*di(i)) /
	(16.518 + 1.481*c_r/c_k*ci(i) - d_r/d_k*di(i)) ;
vpi(i) = 5.520 / (16.518 + 1.481*c_r/c_k*ci(i) - d_r/d_k*di(i)) ;

Wi(i) = 25 * Yi(i)^(1/3) - 17;

dE(i) = sqrt( sq(U_ri(i) - 13*Wi(i)*(upi(i)-u_r)) +
		sq(V_ri(i) - 13*Wi(i)*(vpi(i)-v_r)) +
		sq(W_ri(i) - Wi(i)) );

Ri(i) = 100 - 4.6*dE(i);

cri = (Ri(1)+Ri(2)+Ri(3)+Ri(4)+Ri(5)+Ri(6)+Ri(7)+Ri(8))/8;

U_ri(i) = select(i, U_r1, U_r2, U_r3, U_r4, U_r5, U_r6, U_r7, U_r8);
V_ri(i) = select(i, V_r1, V_r2, V_r3, V_r4, V_r5, V_r6, V_r7, V_r8);
W_ri(i) = select(i, W_r1, W_r2, W_r3, W_r4, W_r5, W_r6, W_r7, W_r8);

{
	The following was generated by the command:

	lam ref.dat ref{1,2,3,4,5,6,7,8}.dat \
		| tabfunc {u,v,c,d}_r {U,V,W}_r1 {U,V,W}_r2 \
			{U,V,W}_r3 {U,V,W}_r4 {U,V,W}_r5 \
			{U,V,W}_r6 {U,V,W}_r7 {U,V,W}_r8

	then modified to eliminate the functional nature.
}
i_cct = fx`p17860(cct);		{ replaced in select calls below }

X`p17860(i)=select(i,2300,2400,2500,2600,2700,2800,2900,3000,
3100,3200,3300,3400,3500,3600,3700,3800,3900,4000,4100,4200,
4300,4400,4500,4600,4700,4800,4900,5000,5100,5200,5300,5400,
5500,5600,5700,5800,5900,6000,6100,6500,7000,7500,8000,8500,
9000,9500,10000,10500,11000,12000,13000,14000,15000,17000,
20000,25000);
fx`p17860(x)=if(x-2300,if(25000-x,fx2`p17860(x,56),56),1);
fx2`p17860(x,i)=if(x-X`p17860(i),
	i+(x-X`p17860(i))/(X`p17860(i+1)-X`p17860(i)),
	fx2`p17860(x,i-1));
u_r=select(i_cct,0.2836,0.2777,0.2722,0.2671,0.2625,
0.2582,0.2542,0.2506,0.2472,0.244,0.241,0.2383,0.2357,0.2333,
0.2311,0.2289,0.227,0.2251,0.2234,0.2217,0.2202,0.2187,0.2173,
0.216,0.2148,0.2136,0.2125,0.2092,0.2081,0.2071,0.2062,0.2053,
0.2044,0.2036,0.2029,0.2022,0.2021,0.2008,0.2001,0.1978,0.1955,
0.1935,0.1919,0.1906,0.1894,0.1884,0.1876,0.1868,0.1862,0.185,
0.1841,0.1834,0.1828,0.1819,0.1809,0.1798);
v_r=select(i_cct,0.3563,0.3552,0.3541,0.3528,0.3516,
0.3503,0.3489,0.3476,0.3462,0.3449,0.3435,0.3422,0.3408,0.3395,
0.3382,0.3369,0.3356,0.3344,0.3332,0.3319,0.3308,0.3296,0.3285,
0.3273,0.3262,0.3252,0.3241,0.3254,0.3244,0.3234,0.3224,0.3214,
0.3205,0.3196,0.3187,0.3178,0.3169,0.3161,0.3153,0.3122,0.3087,
0.3057,0.3029,0.3005,0.2983,0.2963,0.2945,0.2929,0.2914,0.2889,
0.2867,0.2849,0.2833,0.2808,0.278,0.2749);
c_r=select(i_cct,0.4292,0.4784,0.5286,0.5795,0.631,
0.6828,0.7349,0.787,0.839,0.8908,0.9423,0.9934,1.0441,1.0943,
1.1439,1.193,1.2413,1.289,1.336,1.3822,1.4278,1.4725,1.5166,
1.5598,1.6023,1.644,1.685,1.6497,1.6903,1.7295,1.7681,1.8059,
1.8431,1.8796,1.9155,1.9506,1.9851,2.019,2.0522,2.1785,2.3228,
2.4536,2.5717,2.6787,2.7758,2.8642,2.9449,3.0187,3.0863,3.2059,
3.308,3.3959,3.4722,3.5978,3.7381,3.8946);
d_r=select(i_cct,1.663,1.6877,1.7106,1.7317,1.7514,
1.7696,1.7867,1.8027,1.8176,1.8317,1.8449,1.8574,1.8691,1.8802,
1.8907,1.9007,1.9102,1.9192,1.9277,1.9359,1.9437,1.9511,1.9581,
1.9649,1.9714,1.9776,1.9836,1.9975,2.0033,2.0087,2.014,2.019,
2.0239,2.0285,2.033,2.0373,2.0415,2.0455,2.0494,2.0636,2.0789,
2.092,2.1034,2.1133,2.122,2.1297,2.1365,2.1427,2.1483,2.1579,
2.1658,2.1726,2.1783,2.1876,2.1978,2.2088);
U_r1=select(i_cct,39.69,39.75,39.76,39.71,39.62,39.5,
39.35,39.17,38.98,38.77,38.56,38.33,38.1,37.86,37.62,37.38,
37.14,36.9,36.66,36.43,36.19,35.96,35.53,35.53,35.29,35.07,
34.86,34.92,34.7,34.47,34.25,34.03,33.82,33.62,33.42,33.23,
33.04,32.86,32.68,32.03,31.31,30.69,30.15,29.68,29.26,28.89,
28.56,28.27,28,27.54,27.16,26.84,26.57,26.13,25.66,25.15);
V_r1=select(i_cct,0.37,0.74,1.11,1.47,1.83,2.18,2.52,
2.85,3.16,3.47,3.77,4.05,4.32,4.58,4.83,5.07,5.3,5.52,5.73,
5.93,6.12,6.3,6.64,6.64,6.8,6.95,7.1,6.83,7,7.13,7.26,7.39,
7.51,7.63,7.74,7.84,7.95,8.04,8.14,8.48,8.83,9.13,9.37,9.58,
9.75,9.9,10.02,10.13,10.23,10.38,10.51,10.6,10.68,10.8,10.91,
11.02);
W_r1=select(i_cct,63.95,63.74,63.55,63.37,63.2,63.05,
62.9,62.76,62.63,62.51,62.4,62.29,62.19,62.09,62,61.92,61.84,
61.76,61.68,61.61,61.55,61.48,61.37,61.37,61.31,61.26,61.21,
61.05,61,60.96,60.9,60.87,60.83,60.79,60.75,60.72,60.68,60.65,
60.62,60.5,60.37,60.27,60.17,60.09,60.02,59.96,59.91,59.86,
59.81,59.74,59.67,59.62,59.57,59.5,59.42,59.33);
U_r2=select(i_cct,15.82,16.13,16.39,16.61,16.8,16.95,
17.06,17.16,17.23,17.28,17.31,17.32,17.32,17.31,17.29,17.26,
17.21,17.16,17.11,17.05,16.98,16.91,16.84,16.77,16.69,16.61,
16.53,16.56,16.48,16.38,16.29,16.2,16.1,16.01,15.92,15.83,
15.74,15.66,15.57,15.24,14.87,14.52,14.21,13.92,13.67,13.43,
13.22,13.02,12.85,12.53,12.26,12.04,11.84,11.51,11.15,10.75);
V_r2=select(i_cct,4.7,5.36,6.04,6.71,7.39,8.07,8.74,
9.4,10.05,10.7,11.33,11.94,12.55,13.14,13.71,14.27,14.82,15.35,
15.86,16.37,16.85,17.33,17.78,18.23,18.66,19.08,19.49,18.9,
19.33,19.71,20.09,20.45,20.8,21.14,21.47,21.8,22.11,22.41,22.7,
23.8,24.99,26.02,26.92,27.7,28.39,28.99,29.53,59.21,30.44,
31.17,31.77,32.28,32.7,33.37,34.09,34.84);
W_r2=select(i_cct,61.6,61.51,61.42,61.33,61.25,61.18,
61.1,61.03,60.97,60.9,60.84,60.78,60.73,60.68,60.62,60.58,
60.53,60.48,60.44,60.4,60.36,60.32,60.28,60.25,60.21,60.18,
60.15,60.07,60.03,60,59.9,59.95,59.92,59.9,59.87,59.85,59.83,
59.8,59.78,59.7,59.61,59.53,59.46,59.4,59.34,59.29,59.25,60.77,
59.17,59.11,59.05,59.01,58.96,58.9,58.83,58.75);
U_r3=select(i_cct,-17.73,-17.15,-16.41,-15.8,-15.23,
-14.69,-14.18,-13.71,-13.28,-12.87,-12.5,-12.15,-11.83,-11.53,
-11.25,-11,-10.76,-10.55,-10.35,-10.17,-10,-9.84,-9.7,-9.57,
-9.45,-9.34,-9.24,-9.14,-9,-8.91,-8.83,-8.75,-8.67,-8.61,-8.55,
-8.4,-8.44,-8.39,-8.35,-8.21,-8.1,-8.04,-8.01,-8,-8.01,-8.02,
-8.05,-8.08,-8.11,-8.17,-8.24,-8.31,-8.37,-8.48,-8.62,-8.78);
V_r3=select(i_cct,9.35,10.2,11.07,11.96,12.84,13.74,
14.63,15.51,16.39,17.26,18.12,18.96,19.79,20.6,21.4,22.18,
22.95,23.7,24.43,25.14,25.83,26.51,27.17,27.81,28.43,29.04,
29.63,29.03,29.64,30.21,30.76,31.29,31.81,32.32,32.82,33.3,
33.76,34.22,34.66,36.31,38.14,39.74,41.14,42.38,43.47,44.44,
45.31,46.08,46.77,47.99,48.99,49.83,50.45,51.69,52.92,54.24);
W_r3=select(i_cct,60.84,60.9,60.96,61,61.04,61.07,
61.09,61.11,61.13,61.14,61.15,61.16,61.17,61.17,61.17,61.17,
61.17,61.17,61.17,61.16,61.16,61.16,61.15,61.15,61.14,61.13,
61.13,61.17,61.16,61.15,61.1,61.14,61.13,61.12,61.12,61.11,
61.1,61.09,61.08,61.05,61.01,60.97,60.93,60.89,60.86,60.83,
60.8,60.77,60.75,60.7,60.66,60.62,60.59,60.54,60.48,60.41);
U_r4=select(i_cct,-43.38,-43.01,-42.62,-42.22,-41.82,
-41.42,-41.03,-40.65,-40.27,-39.91,-39.56,-39.22,-38.9,-38.59,
-38.29,-38.01,-37.73,-37.47,-37.23,-36.99,-36.76,-36.55,-36.34,
-36.15,-36.96,-35.78,-35.61,-35.45,-35.25,-35.07,-34.89,-34.73,
-34.57,-34.42,-34.28,-34.14,-34.02,-33.89,-33.78,-33.39,-32.94,
-32.61,-32.33,-32.1,-31.92,-31.76,-31.62,-31.51,-31.41,-31.25,
-31.12,-31.02,-30.94,-30.82,-30.7,-30.59);
V_r4=select(i_cct,5.59,5.91,6.25,6.6,6.97,7.35,7.73,
8.13,8.52,8.93,9.33,9.73,10.13,10.53,10.92,11.31,11.7,12.08,
12.45,12.82,13.18,13.54,13.89,14.23,14.56,14.89,15.21,14.69,15,
15.3,15.6,15.89,16.17,16.45,16.71,16.98,17.23,17.48,17.73,
18.62,19.67,20.57,21.37,22.08,22.7,23.27,23.77,24.22,24.63,
25.34,25.94,26.44,26.86,27.55,28.3,29.1);
W_r4=select(i_cct,56.92,57.15,57.36,57.55,57.73,
57.89,58.04,58.18,58.31,58.43,58.54,58.65,58.75,58.84,58.92,59,
59.08,59.15,59.22,59.28,59.34,59.4,59.45,59.5,59.55,59.6,59.64,
59.8,59.83,59.87,59.9,59.94,59.97,60,60.03,60.06,60.09,60.11,
60.14,60.23,60.32,60.4,60.47,60.52,60.57,60.61,60.65,60.68,
60.71,60.76,60.8,60.84,60.87,60.91,60.96,61.01);
U_r5=select(i_cct,-37.38,-37.17,-36.92,-36.63,-36.33,
-36.01,-35.68,-35.34,-35.01,-34.67,-34.33,-34,-33.67,-33.34,
-33.03,-32.67,-32.41,-32.11,-31.83,-31.54,-31.27,-31.01,-30.75,
-30.5,-30.26,-30.02,-29.79,-29.76,-29.5,-29.26,-29.03,-28.81,
-28.6,-28.4,-28.2,-28.01,-27.82,-27.65,-27.47,-26.85,-26.17,
-25.64,-25.1,-24.68,-24.31,-23.98,-23.69,-23.44,-23.21,-22.82,
-22.49,-22.23,-22,-21.64,-21.25,-20.84);
V_r5=select(i_cct,-0.88,-1.22,-1.56,-1.88,-2.19,
-2.48,-2.76,-3.02,-3.26,-3.49,-3.71,-3.92,-4.11,-4.29,-4.46,
-4.5,-4.76,-4.9,-5.03,-5.15,-5.27,-5.37,-5.47,-5.57,-5.65,
-5.73,-5.81,-5.78,-5.85,-5.92,-5.99,-6.05,-6.11,-6.17,-6.22,
-6.27,-6.32,-6.36,-6.4,-6.54,-6.68,-6.76,-6.85,-6.91,-6.95,
-6.98,-7,-7.01,-7.02,-7.03,-7.03,-7.02,-7.01,-6.99,-6.95,-6.9);
W_r5=select(i_cct,58.06,58.27,58.47,58.65,58.81,
58.97,59.11,59.25,59.37,59.49,59.6,59.71,59.8,59.9,59.98,60.03,
60.15,60.22,60.29,60.36,60.42,60.48,60.54,60.6,60.65,60.7,
60.75,60.89,60.93,60.98,61,61.06,61.09,61.13,61.17,61.2,61.23,
61.26,61.29,61.41,61.52,61.61,61.71,61.79,61.86,61.92,61.97,
62.02,62.06,62.14,62.2,62.25,62.29,62.36,62.44,62.53);
U_r6=select(i_cct,-21.62,-22.11,-22.5,-22.8,-23.03,
-23.2,-23.31,-23.37,-23.4,-23.39,-23.35,-23.28,-23.19,-23.09,
-22.97,-22.84,-22.7,-22.55,-22.39,-22.23,-22.07,-21.9,-21.73,
-21.56,-21.39,-21.22,-21.06,-21.19,-20.94,-20.76,-20.58,-20.41,
-20.24,-20.07,-19.9,-19.74,-19.59,-19.43,-19.28,-18.72,-18.09,
-17.53,-17.04,-16.6,-16.21,-15.87,-15.55,-15.27,-15.02,-14.58,
-14.22,-13.91,-13.64,-13.22,-12.76,-12.26);
V_r6=select(i_cct,-8.16,-9.11,-10.05,-10.98,-11.88,
-12.77,-13.63,-14.46,-15.27,-16.05,-16.79,-17.51,-18.2,-18.85,
-19.49,-20.09,-20.67,-21.22,-21.75,-22.25,-22.73,-23.19,-23.63,
-24.05,-24.45,-24.83,-25.2,-24.69,-25.1,-25.44,-25.78,-26.1,
-26.4,-26.69,-26.97,-27.24,-27.5,-27.74,-27.98,-28.83,-29.72,
-30.45,-31.05,-31.55,-31.98,-32.34,-32.65,-32.91,-33.14,-33.52,
-33.82,-34.05,-34.24,-34.52,-34.8,-35.07);
W_r6=select(i_cct,57.43,57.58,57.72,57.86,57.98,
58.11,58.23,58.34,58.45,58.55,58.65,58.74,58.83,58.92,59,59.08,
59.16,59.23,59.3,59.37,59.43,59.5,59.56,59.61,59.67,59.72,
59.77,59.92,59.95,60,60.02,60.09,60.13,60.17,60.21,60.25,60.28,
60.32,60.35,60.48,60.63,60.75,60.86,60.96,61.05,61.13,61.2,
61.27,61.33,61.43,61.52,61.59,61.65,61.76,61.87,62);
U_r7=select(i_cct,21.04,20.23,19.46,18.74,18.07,
17.45,16.86,16.32,15.82,15.36,14.93,14.53,14.16,13.82,13.51,
13.23,12.95,12.7,12.47,12.26,12.07,11.89,11.72,11.56,11.42,
11.29,11.17,10.7,11.05,10.93,10.82,10.72,10.62,10.53,10.45,
10.37,10.3,10.23,10.17,9.96,9.78,9.65,9.56,9.5,9.45,9.42,9.41,
9.4,9.39,9.4,9.41,9.42,9.44,9.48,9.53,9.59);
V_r7=select(i_cct,-7.77,-8.45,-9.13,-9.82,-10.52,
-11.22,-11.91,-12.6,-13.28,-13.94,-14.59,-15.23,-15.85,-16.45,
-17.03,-17.59,-18.15,-18.68,-19.2,-19.7,-20.18,-20.64,-21.09,
-21.53,-21.94,-22.35,-22.73,-22.11,-22.44,-22.81,-23.16,-23.5,
-23.83,-24.15,-24.46,-24.76,-25.04,-25.32,-25.58,-26.56,-27.6,
-28.47,-29.22,-29.85,-30.4,-30.87,-31.28,-31.64,-31.96,-32.49,
-32.92,-33.24,-33.55,-34,-34.46,-34.93);
W_r7=select(i_cct,60.92,60.83,60.74,60.67,60.61,
60.55,60.5,60.46,60.42,60.39,60.36,60.33,60.31,60.29,60.27,
60.26,60.24,60.22,60.21,60.2,60.19,60.19,60.18,60.17,60.17,
60.16,60.16,60.1,60.15,60.15,60.14,60.14,60.14,60.14,60.14,
60.14,60.14,60.14,60.14,60.14,60.15,60.16,60.18,60.19,60.2,
60.2,60.23,60.25,60.26,60.28,60.3,60.32,60.34,60.37,60.4,60.44);
U_r8=select(i_cct,51.15,50.06,49,47.96,46.95,45.98,
45.04,44.14,43.28,42.46,41.67,40.92,40.21,39.53,38.88,38.26,
37.67,37.11,36.57,36.06,35.58,35.12,34.67,34.25,33.85,33.47,
33.1,33.17,32.81,32.44,32.08,31.75,31.43,31.12,30.82,30.54,
30.28,30.02,29.77,28.89,27.98,27.23,26.61,26.09,25.64,25.26,
24.94,24.65,24.4,23.98,23.65,23.38,23.16,22.81,22.45,22.09);
V_r8=select(i_cct,-6.3,-6.52,-6.77,-7.03,-7.31,-7.6,
-7.9,-8.21,-8.52,-8.84,-9.15,-9.47,-9.78,-10.1,-10.41,-10.71,
-11.01,-11.31,-11.6,-11.88,-12.16,-12.43,-12.7,-12.95,-13.21,
-13.45,-13.69,-13.35,-13.54,-13.77,-13.99,-14.21,-14.42,-14.63,
-14.82,-15.02,-15.2,-15.39,-15.56,-16.21,-16.92,-17.54,-18.06,
-18.52,-18.92,-19.27,-19.58,-19.85,-20.09,-20.5,-20.84,-21.11,
-21.34,-21.71,-22.09,-22.48);
W_r8=select(i_cct,65,64.75,64.53,64.32,64.14,63.97,
63.81,63.67,63.54,63.42,63.31,63.21,63.11,63.03,62.94,62.87,
62.8,62.73,62.67,62.61,62.55,62.5,62.45,62.41,62.37,62.33,
62.29,62.23,62.18,62.14,62.11,62.07,62.04,62.01,61.99,61.96,
61.94,61.91,61.89,61.81,61.73,61.67,61.62,61.58,61.54,61.51,
61.49,61.47,61.45,61.42,61.4,61.39,61.37,61.35,61.33,61.32);