Repository navigation
Expand file tree
/
Copy pathCieDIlluminantData.cpp
More file actions
209 lines (196 loc) · 25.9 KB
/
Copy pathCieDIlluminantData.cpp
File metadata and controls
209 lines (196 loc) · 25.9 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
// SPDX-FileCopyrightText: 2026 Błażej Szczygieł <mumei6102@gmail.com>
// SPDX-License-Identifier: BSD-3-Clause
#include "CieDIlluminantData.hpp"
namespace CieDIlluminantData {
// Official CIE D-series basis functions are tabulated at 5 nm.
// These 1 nm values were generated from colour-science with Sprague interpolation.
static constexpr array<double, g_limitedVisSize> S0 = {
6.340000000000e+01, 6.338048000000e+01, 6.365568000000e+01, 6.408688000000e+01, 6.444288000000e+01, 6.460000000000e+01, 6.459200000000e+01, 6.437008000000e+01,
6.422336000000e+01, 6.460864000000e+01, 6.580000000000e+01, 6.772720000000e+01, 7.039616000000e+01, 7.359712000000e+01, 7.698048000000e+01, 8.030000000000e+01,
8.358304000000e+01, 8.683872000000e+01, 8.991200000000e+01, 9.260528000000e+01, 9.480000000000e+01, 9.647688000000e+01, 9.757000000000e+01, 9.828312000000e+01,
9.895584000000e+01, 9.980000000000e+01, 1.007780800000e+02, 1.019195200000e+02, 1.031013600000e+02, 1.041052000000e+02, 1.048000000000e+02, 1.052222400000e+02,
1.053565600000e+02, 1.053128800000e+02, 1.052744000000e+02, 1.053500000000e+02, 1.055172800000e+02, 1.057996800000e+02, 1.060768800000e+02, 1.061500800000e+02,
1.059000000000e+02, 1.053781600000e+02, 1.046172800000e+02, 1.036524000000e+02, 1.025419200000e+02, 1.013500000000e+02, 1.000812800000e+02, 9.870200000000e+01,
9.746832000000e+01, 9.675864000000e+01, 9.680000000000e+01, 9.756248000000e+01, 9.908752000000e+01, 1.011365600000e+02, 1.033120000000e+02, 1.053500000000e+02,
1.072881600000e+02, 1.091141600000e+02, 1.108137600000e+02, 1.124053600000e+02, 1.139000000000e+02, 1.152549600000e+02, 1.164217600000e+02, 1.174805600000e+02,
1.185649600000e+02, 1.197500000000e+02, 1.210176000000e+02, 1.223909600000e+02, 1.237387200000e+02, 1.248504800000e+02, 1.256000000000e+02, 1.260115200000e+02,
1.260639200000e+02, 1.258803200000e+02, 1.256659200000e+02, 1.255500000000e+02, 1.255104000000e+02, 1.255631200000e+02, 1.256466400000e+02, 1.256481600000e+02,
1.255000000000e+02, 1.252309600000e+02, 1.248506400000e+02, 1.243883200000e+02, 1.238928000000e+02, 1.234000000000e+02, 1.229067200000e+02, 1.224095200000e+02,
1.219455200000e+02, 1.215655200000e+02, 1.213000000000e+02, 1.211425600000e+02, 1.210907200000e+02, 1.211228800000e+02, 1.212030400000e+02, 1.213000000000e+02,
1.214142400000e+02, 1.215548800000e+02, 1.216475200000e+02, 1.215841600000e+02, 1.213000000000e+02, 1.208067200000e+02, 1.201003200000e+02, 1.192379200000e+02,
1.183147200000e+02, 1.174000000000e+02, 1.164872000000e+02, 1.155707200000e+02, 1.147150400000e+02, 1.140073600000e+02, 1.135000000000e+02, 1.131920800000e+02,
1.130986400000e+02, 1.131532000000e+02, 1.132449600000e+02, 1.133000000000e+02, 1.133286400000e+02, 1.133280000000e+02, 1.132901600000e+02, 1.132143200000e+02,
1.131000000000e+02, 1.129368000000e+02, 1.127125600000e+02, 1.124503200000e+02, 1.121884800000e+02, 1.119500000000e+02, 1.117302400000e+02, 1.115339200000e+02,
1.113372000000e+02, 1.111004800000e+02, 1.108000000000e+02, 1.104472800000e+02, 1.100512000000e+02, 1.096151200000e+02, 1.091446400000e+02, 1.086500000000e+02,
1.081332800000e+02, 1.075855200000e+02, 1.070721600000e+02, 1.066908000000e+02, 1.065000000000e+02, 1.064877600000e+02, 1.066569600000e+02, 1.069581600000e+02,
1.073089600000e+02, 1.076500000000e+02, 1.079872000000e+02, 1.083245600000e+02, 1.086123200000e+02, 1.087840800000e+02, 1.088000000000e+02, 1.086608800000e+02,
1.083563200000e+02, 1.079357600000e+02, 1.074816000000e+02, 1.070500000000e+02, 1.066337600000e+02, 1.062333600000e+02, 1.058665600000e+02, 1.055517600000e+02,
1.053000000000e+02, 1.051136000000e+02, 1.049957600000e+02, 1.049299200000e+02, 1.048884800000e+02, 1.048500000000e+02, 1.048158400000e+02, 1.047895200000e+02,
1.047388000000e+02, 1.046180800000e+02, 1.044000000000e+02, 1.040863200000e+02, 1.036700000000e+02, 1.031836800000e+02, 1.026817600000e+02, 1.022000000000e+02,
1.017331200000e+02, 1.012832800000e+02, 1.008490400000e+02, 1.004228000000e+02, 1.000000000000e+02, 9.958808000000e+01, 9.919544000000e+01, 9.881080000000e+01,
9.841536000000e+01, 9.800000000000e+01, 9.756784000000e+01, 9.711360000000e+01, 9.667016000000e+01, 9.628872000000e+01, 9.600000000000e+01, 9.579432000000e+01,
9.566816000000e+01, 9.560400000000e+01, 9.557264000000e+01, 9.555000000000e+01, 9.553696000000e+01, 9.553920000000e+01, 9.550864000000e+01, 9.537608000000e+01,
9.510000000000e+01, 9.468952000000e+01, 9.414584000000e+01, 9.350016000000e+01, 9.280448000000e+01, 9.210000000000e+01, 9.138448000000e+01, 9.065016000000e+01,
8.996584000000e+01, 8.942952000000e+01, 8.910000000000e+01, 8.897232000000e+01, 8.905984000000e+01, 8.929536000000e+01, 8.956688000000e+01, 8.980000000000e+01,
9.000528000000e+01, 9.017936000000e+01, 9.031744000000e+01, 9.042352000000e+01, 9.050000000000e+01, 9.053624000000e+01, 9.052104000000e+01, 9.047384000000e+01,
9.042704000000e+01, 9.040000000000e+01, 9.038880000000e+01, 9.039744000000e+01, 9.040568000000e+01, 9.037992000000e+01, 9.030000000000e+01, 9.016920000000e+01,
8.998208000000e+01, 8.976096000000e+01, 8.954304000000e+01, 8.935000000000e+01, 8.917712000000e+01, 8.902976000000e+01, 8.887920000000e+01, 8.867864000000e+01,
8.840000000000e+01, 8.805440000000e+01, 8.764800000000e+01, 8.719160000000e+01, 8.670320000000e+01, 8.620000000000e+01, 8.568240000000e+01, 8.514360000000e+01,
8.463680000000e+01, 8.424000000000e+01, 8.400000000000e+01, 8.390856000000e+01, 8.396960000000e+01, 8.414064000000e+01, 8.435088000000e+01, 8.455000000000e+01,
8.474336000000e+01, 8.493344000000e+01, 8.508432000000e+01, 8.514920000000e+01, 8.510000000000e+01, 8.493792000000e+01, 8.465792000000e+01, 8.429192000000e+01,
8.389312000000e+01, 8.350000000000e+01, 8.310880000000e+01, 8.271672000000e+01, 8.235744000000e+01, 8.207616000000e+01, 8.190000000000e+01, 8.183088000000e+01,
8.188104000000e+01, 8.200920000000e+01, 8.214656000000e+01, 8.225000000000e+01, 8.232704000000e+01, 8.237200000000e+01, 8.240776000000e+01, 8.247552000000e+01,
8.260000000000e+01, 8.276768000000e+01, 8.296856000000e+01, 8.320144000000e+01, 8.346432000000e+01, 8.375000000000e+01, 8.405632000000e+01, 8.439144000000e+01,
8.469656000000e+01, 8.488368000000e+01, 8.490000000000e+01, 8.474776000000e+01, 8.440696000000e+01, 8.394816000000e+01, 8.348896000000e+01, 8.310000000000e+01,
8.276720000000e+01, 8.250456000000e+01, 8.224232000000e+01, 8.186408000000e+01, 8.130000000000e+01, 8.057824000000e+01, 7.971312000000e+01, 7.873200000000e+01,
7.768048000000e+01, 7.660000000000e+01, 7.549072000000e+01, 7.433840000000e+01, 7.325648000000e+01, 7.241056000000e+01, 7.190000000000e+01, 7.171496000000e+01,
7.187848000000e+01, 7.227800000000e+01, 7.272592000000e+01, 7.310000000000e+01, 7.341888000000e+01, 7.367360000000e+01, 7.387992000000e+01, 7.408024000000e+01,
7.430000000000e+01, 7.450864000000e+01, 7.467192000000e+01, 7.482920000000e+01, 7.504608000000e+01, 7.535000000000e+01, 7.572832000000e+01, 7.620560000000e+01,
7.662328000000e+01, 7.673696000000e+01, 7.640000000000e+01, 7.565440000000e+01, 7.451048000000e+01, 7.306256000000e+01, 7.146784000000e+01, 6.985000000000e+01,
6.820192000000e+01, 6.650136000000e+01, 6.494760000000e+01, 6.382384000000e+01, 6.330000000000e+01, 6.336352000000e+01, 6.405600000000e+01, 6.517848000000e+01,
6.639936000000e+01, 6.750000000000e+01, 6.851232000000e+01, 6.942408000000e+01, 7.023744000000e+01, 7.098880000000e+01, 7.170000000000e+01, 7.232792000000e+01,
7.282304000000e+01, 7.325616000000e+01, 7.374528000000e+01, 7.435000000000e+01, 7.505168000000e+01, 7.588016000000e+01, 7.665264000000e+01, 7.708312000000e+01,
7.700000000000e+01, 7.643104000000e+01, 7.534224000000e+01, 7.391144000000e+01, 7.243504000000e+01, 7.110000000000e+01, 6.987440000000e+01, 6.878104000000e+01,
6.773328000000e+01, 6.657152000000e+01, 6.520000000000e+01, 6.368592000000e+01, 6.209528000000e+01, 6.039064000000e+01, 5.850960000000e+01, 5.645000000000e+01,
5.423344000000e+01, 5.180304000000e+01, 4.954904000000e+01, 4.806304000000e+01, 4.770000000000e+01, 4.838800000000e+01, 5.015384000000e+01, 5.268768000000e+01,
5.547312000000e+01, 5.815000000000e+01, 6.076016000000e+01, 6.331208000000e+01, 6.561240000000e+01, 6.742272000000e+01, 6.860000000000e+01, 6.912568000000e+01,
6.893720000000e+01, 6.825872000000e+01, 6.746384000000e+01, 6.680000000000e+01, 6.623168000000e+01, 6.577112000000e+01, 6.542696000000e+01, 6.517480000000e+01,
6.500000000000e+01
};
static constexpr array<double, g_limitedVisSize> S1 = {
3.850000000000e+01, 3.814512000000e+01, 3.783064000000e+01, 3.752416000000e+01, 3.717168000000e+01, 3.675000000000e+01, 3.626928000000e+01, 3.571016000000e+01,
3.519704000000e+01, 3.492192000000e+01, 3.500000000000e+01, 3.540720000000e+01, 3.615344000000e+01, 3.713768000000e+01, 3.819152000000e+01, 3.920000000000e+01,
4.017776000000e+01, 4.112408000000e+01, 4.200080000000e+01, 4.276752000000e+01, 4.340000000000e+01, 4.388376000000e+01, 4.419200000000e+01, 4.439024000000e+01,
4.458768000000e+01, 4.485000000000e+01, 4.516416000000e+01, 4.554304000000e+01, 4.592272000000e+01, 4.619640000000e+01, 4.630000000000e+01, 4.624728000000e+01,
4.602712000000e+01, 4.570096000000e+01, 4.537120000000e+01, 4.510000000000e+01, 4.487536000000e+01, 4.470856000000e+01, 4.454536000000e+01, 4.429416000000e+01,
4.390000000000e+01, 4.338352000000e+01, 4.275216000000e+01, 4.203280000000e+01, 4.127024000000e+01, 4.050000000000e+01, 3.972016000000e+01, 3.892400000000e+01,
3.817104000000e+01, 3.754608000000e+01, 3.710000000000e+01, 3.682928000000e+01, 3.674704000000e+01, 3.679280000000e+01, 3.686576000000e+01, 3.690000000000e+01,
3.690544000000e+01, 3.687760000000e+01, 3.682256000000e+01, 3.675952000000e+01, 3.670000000000e+01, 3.663240000000e+01, 3.654696000000e+01, 3.645352000000e+01,
3.636848000000e+01, 3.630000000000e+01, 3.624544000000e+01, 3.620912000000e+01, 3.616440000000e+01, 3.606968000000e+01, 3.590000000000e+01, 3.565888000000e+01,
3.534040000000e+01, 3.497192000000e+01, 3.459904000000e+01, 3.425000000000e+01, 3.391968000000e+01, 3.361232000000e+01, 3.330936000000e+01, 3.297840000000e+01,
3.260000000000e+01, 3.218128000000e+01, 3.172376000000e+01, 3.123824000000e+01, 3.074272000000e+01, 3.025000000000e+01, 2.975872000000e+01, 2.926824000000e+01,
2.878776000000e+01, 2.832928000000e+01, 2.790000000000e+01, 2.749992000000e+01, 2.713096000000e+01, 2.678400000000e+01, 2.644384000000e+01, 2.610000000000e+01,
2.575376000000e+01, 2.540520000000e+01, 2.504984000000e+01, 2.468248000000e+01, 2.430000000000e+01, 2.390184000000e+01, 2.348664000000e+01, 2.305944000000e+01,
2.262864000000e+01, 2.220000000000e+01, 2.177280000000e+01, 2.134704000000e+01, 2.092488000000e+01, 2.050872000000e+01, 2.010000000000e+01, 1.969992000000e+01,
1.931088000000e+01, 1.892784000000e+01, 1.854240000000e+01, 1.815000000000e+01, 1.775184000000e+01, 1.734624000000e+01, 1.694304000000e+01, 1.655784000000e+01,
1.620000000000e+01, 1.586632000000e+01, 1.555544000000e+01, 1.526256000000e+01, 1.497968000000e+01, 1.470000000000e+01, 1.442368000000e+01, 1.415256000000e+01,
1.387144000000e+01, 1.355832000000e+01, 1.320000000000e+01, 1.279928000000e+01, 1.235624000000e+01, 1.188120000000e+01, 1.139136000000e+01, 1.090000000000e+01,
1.040624000000e+01, 9.908000000000e+00, 9.424560000000e+00, 8.983120000000e+00, 8.600000000000e+00, 8.274880000000e+00, 8.013360000000e+00, 7.793840000000e+00,
7.580320000000e+00, 7.350000000000e+00, 7.106720000000e+00, 6.847840000000e+00, 6.582960000000e+00, 6.330080000000e+00, 6.100000000000e+00, 5.888080000000e+00,
5.692960000000e+00, 5.509840000000e+00, 5.330720000000e+00, 5.150000000000e+00, 4.968320000000e+00, 4.785840000000e+00, 4.599360000000e+00, 4.404880000000e+00,
4.200000000000e+00, 3.984720000000e+00, 3.758560000000e+00, 3.524400000000e+00, 3.287040000000e+00, 3.050000000000e+00, 2.812960000000e+00, 2.575600000000e+00,
2.341440000000e+00, 2.115280000000e+00, 1.900000000000e+00, 1.695840000000e+00, 1.504240000000e+00, 1.320640000000e+00, 1.137440000000e+00, 9.500000000000e-01,
7.592000000000e-01, 5.642400000000e-01, 3.688800000000e-01, 1.795200000000e-01, -2.775557561563e-17, -1.710400000000e-01, -3.339200000000e-01, -4.908000000000e-01,
-6.452800000000e-01, -8.000000000000e-01, -9.547200000000e-01, -1.109200000000e+00, -1.266080000000e+00, -1.428960000000e+00, -1.600000000000e+00, -1.777680000000e+00,
-1.959680000000e+00, -2.147680000000e+00, -2.344480000000e+00, -2.550000000000e+00, -2.763200000000e+00, -2.986880000000e+00, -3.201760000000e+00, -3.378640000000e+00,
-3.500000000000e+00, -3.570240000000e+00, -3.589760000000e+00, -3.571280000000e+00, -3.536000000000e+00, -3.500000000000e+00, -3.462080000000e+00, -3.420080000000e+00,
-3.394880000000e+00, -3.415680000000e+00, -3.500000000000e+00, -3.646480000000e+00, -3.858880000000e+00, -4.117280000000e+00, -4.388480000000e+00, -4.650000000000e+00,
-4.904800000000e+00, -5.152480000000e+00, -5.387360000000e+00, -5.604240000000e+00, -5.800000000000e+00, -5.972400000000e+00, -6.118160000000e+00, -6.245920000000e+00,
-6.370080000000e+00, -6.500000000000e+00, -6.634240000000e+00, -6.773520000000e+00, -6.916400000000e+00, -7.059280000000e+00, -7.200000000000e+00, -7.339600000000e+00,
-7.478000000000e+00, -7.616400000000e+00, -7.756800000000e+00, -7.900000000000e+00, -8.045600000000e+00, -8.194400000000e+00, -8.341200000000e+00, -8.478000000000e+00,
-8.600000000000e+00, -8.708400000000e+00, -8.803200000000e+00, -8.888000000000e+00, -8.968800000000e+00, -9.050000000000e+00, -9.131200000000e+00, -9.212000000000e+00,
-9.296800000000e+00, -9.391600000000e+00, -9.500000000000e+00, -9.620720000000e+00, -9.752400000000e+00, -9.894080000000e+00, -1.004416000000e+01, -1.020000000000e+01,
-1.036112000000e+01, -1.052968000000e+01, -1.068984000000e+01, -1.081800000000e+01, -1.090000000000e+01, -1.093880000000e+01, -1.093384000000e+01, -1.089688000000e+01,
-1.084752000000e+01, -1.080000000000e+01, -1.075296000000e+01, -1.070528000000e+01, -1.067000000000e+01, -1.066472000000e+01, -1.070000000000e+01, -1.077656000000e+01,
-1.089920000000e+01, -1.105184000000e+01, -1.120768000000e+01, -1.135000000000e+01, -1.148176000000e+01, -1.160064000000e+01, -1.171632000000e+01, -1.184600000000e+01,
-1.200000000000e+01, -1.217288000000e+01, -1.236072000000e+01, -1.256256000000e+01, -1.277680000000e+01, -1.300000000000e+01, -1.323136000000e+01, -1.347416000000e+01,
-1.370456000000e+01, -1.388696000000e+01, -1.400000000000e+01, -1.404544000000e+01, -1.401696000000e+01, -1.394048000000e+01, -1.385920000000e+01, -1.380000000000e+01,
-1.375808000000e+01, -1.373728000000e+01, -1.372128000000e+01, -1.368128000000e+01, -1.360000000000e+01, -1.348552000000e+01, -1.334248000000e+01, -1.317544000000e+01,
-1.299200000000e+01, -1.280000000000e+01, -1.259984000000e+01, -1.238784000000e+01, -1.219224000000e+01, -1.205464000000e+01, -1.200000000000e+01, -1.202424000000e+01,
-1.213024000000e+01, -1.229424000000e+01, -1.247664000000e+01, -1.265000000000e+01, -1.281760000000e+01, -1.297984000000e+01, -1.312368000000e+01, -1.323352000000e+01,
-1.330000000000e+01, -1.331968000000e+01, -1.328448000000e+01, -1.321528000000e+01, -1.314688000000e+01, -1.310000000000e+01, -1.307040000000e+01, -1.306248000000e+01,
-1.305376000000e+01, -1.300704000000e+01, -1.290000000000e+01, -1.274104000000e+01, -1.253376000000e+01, -1.228848000000e+01, -1.202240000000e+01, -1.175000000000e+01,
-1.147088000000e+01, -1.118128000000e+01, -1.091248000000e+01, -1.070968000000e+01, -1.060000000000e+01, -1.058088000000e+01, -1.065848000000e+01, -1.080208000000e+01,
-1.096048000000e+01, -1.110000000000e+01, -1.122560000000e+01, -1.133528000000e+01, -1.143016000000e+01, -1.151704000000e+01, -1.160000000000e+01, -1.167232000000e+01,
-1.172656000000e+01, -1.177280000000e+01, -1.182784000000e+01, -1.190000000000e+01, -1.198656000000e+01, -1.209200000000e+01, -1.218864000000e+01, -1.223328000000e+01,
-1.220000000000e+01, -1.209328000000e+01, -1.190864000000e+01, -1.167200000000e+01, -1.142656000000e+01, -1.120000000000e+01, -1.098784000000e+01, -1.079280000000e+01,
-1.060656000000e+01, -1.041232000000e+01, -1.020000000000e+01, -9.979040000000e+00, -9.760160000000e+00, -9.533280000000e+00, -9.281600000000e+00, -9.000000000000e+00,
-8.692480000000e+00, -8.350080000000e+00, -8.032480000000e+00, -7.830880000000e+00, -7.800000000000e+00, -7.927840000000e+00, -8.217280000000e+00, -8.622720000000e+00,
-9.068160000000e+00, -9.500000000000e+00, -9.924160000000e+00, -1.034272000000e+01, -1.072128000000e+01, -1.101584000000e+01, -1.120000000000e+01, -1.127200000000e+01,
-1.122208000000e+01, -1.108816000000e+01, -1.093344000000e+01, -1.080000000000e+01, -1.068192000000e+01, -1.058096000000e+01, -1.050080000000e+01, -1.044064000000e+01,
-1.040000000000e+01
};
static constexpr array<double, g_limitedVisSize> S2 = {
3.000000000000e+00, 2.772800000000e+00, 2.582880000000e+00, 2.418960000000e+00, 2.262240000000e+00, 2.100000000000e+00, 1.933920000000e+00, 1.763760000000e+00,
1.586400000000e+00, 1.399040000000e+00, 1.200000000000e+00, 9.903200000000e-01, 7.728000000000e-01, 5.452800000000e-01, 3.041600000000e-01, 5.000000000000e-02,
-2.156800000000e-01, -4.971200000000e-01, -7.649600000000e-01, -9.748000000000e-01, -1.100000000000e+00, -1.145040000000e+00, -1.106160000000e+00, -1.009280000000e+00,
-8.976000000000e-01, -8.000000000000e-01, -7.124800000000e-01, -6.360800000000e-01, -5.744800000000e-01, -5.288800000000e-01, -5.000000000000e-01, -4.914400000000e-01,
-5.072800000000e-01, -5.391200000000e-01, -5.729600000000e-01, -6.000000000000e-01, -6.217600000000e-01, -6.371200000000e-01, -6.504800000000e-01, -6.698400000000e-01,
-7.000000000000e-01, -7.399200000000e-01, -7.908800000000e-01, -8.478400000000e-01, -9.024000000000e-01, -9.500000000000e-01, -9.918400000000e-01, -1.026240000000e+00,
-1.063040000000e+00, -1.117840000000e+00, -1.200000000000e+00, -1.306720000000e+00, -1.437440000000e+00, -1.586160000000e+00, -1.742880000000e+00, -1.900000000000e+00,
-2.058080000000e+00, -2.218160000000e+00, -2.370240000000e+00, -2.500320000000e+00, -2.600000000000e+00, -2.669280000000e+00, -2.705040000000e+00, -2.718800000000e+00,
-2.729760000000e+00, -2.750000000000e+00, -2.777440000000e+00, -2.813600000000e+00, -2.852560000000e+00, -2.883520000000e+00, -2.900000000000e+00, -2.904320000000e+00,
-2.896560000000e+00, -2.880800000000e+00, -2.863840000000e+00, -2.850000000000e+00, -2.838560000000e+00, -2.830000000000e+00, -2.822640000000e+00, -2.813280000000e+00,
-2.800000000000e+00, -2.783760000000e+00, -2.765040000000e+00, -2.744320000000e+00, -2.722400000000e+00, -2.700000000000e+00, -2.677120000000e+00, -2.653520000000e+00,
-2.631120000000e+00, -2.612720000000e+00, -2.600000000000e+00, -2.592160000000e+00, -2.588320000000e+00, -2.588480000000e+00, -2.592640000000e+00, -2.600000000000e+00,
-2.610240000000e+00, -2.624480000000e+00, -2.634720000000e+00, -2.628960000000e+00, -2.600000000000e+00, -2.549200000000e+00, -2.475920000000e+00, -2.386640000000e+00,
-2.292160000000e+00, -2.200000000000e+00, -2.109280000000e+00, -2.019840000000e+00, -1.935600000000e+00, -1.861360000000e+00, -1.800000000000e+00, -1.752080000000e+00,
-1.719200000000e+00, -1.696320000000e+00, -1.675040000000e+00, -1.650000000000e+00, -1.622080000000e+00, -1.590720000000e+00, -1.557760000000e+00, -1.526800000000e+00,
-1.500000000000e+00, -1.476480000000e+00, -1.456160000000e+00, -1.437840000000e+00, -1.419520000000e+00, -1.400000000000e+00, -1.379520000000e+00, -1.357840000000e+00,
-1.336160000000e+00, -1.316480000000e+00, -1.300000000000e+00, -1.286320000000e+00, -1.275360000000e+00, -1.266400000000e+00, -1.258240000000e+00, -1.250000000000e+00,
-1.241760000000e+00, -1.233600000000e+00, -1.224640000000e+00, -1.213680000000e+00, -1.200000000000e+00, -1.183360000000e+00, -1.163040000000e+00, -1.140720000000e+00,
-1.119200000000e+00, -1.100000000000e+00, -1.082720000000e+00, -1.067920000000e+00, -1.052320000000e+00, -1.030720000000e+00, -1.000000000000e+00, -9.610400000000e-01,
-9.139200000000e-01, -8.608000000000e-01, -8.052800000000e-01, -7.500000000000e-01, -6.947200000000e-01, -6.392000000000e-01, -5.860800000000e-01, -5.389600000000e-01,
-5.000000000000e-01, -4.691200000000e-01, -4.468800000000e-01, -4.306400000000e-01, -4.160000000000e-01, -4.000000000000e-01, -3.830400000000e-01, -3.650400000000e-01,
-3.454400000000e-01, -3.238400000000e-01, -3.000000000000e-01, -2.736800000000e-01, -2.446400000000e-01, -2.136000000000e-01, -1.817600000000e-01, -1.500000000000e-01,
-1.182400000000e-01, -8.640000000000e-02, -5.536000000000e-02, -2.632000000000e-02, 1.283695372223e-16, 2.368000000000e-02, 4.464000000000e-02, 6.360000000000e-02,
8.176000000000e-02, 1.000000000000e-01, 1.182400000000e-01, 1.364000000000e-01, 1.553600000000e-01, 1.763200000000e-01, 2.000000000000e-01, 2.274400000000e-01,
2.609600000000e-01, 2.964800000000e-01, 3.272000000000e-01, 3.500000000000e-01, 3.660800000000e-01, 3.732800000000e-01, 3.852800000000e-01, 4.232800000000e-01,
5.000000000000e-01, 6.128000000000e-01, 7.628800000000e-01, 9.389600000000e-01, 1.122240000000e+00, 1.300000000000e+00, 1.473920000000e+00, 1.643760000000e+00,
1.806400000000e+00, 1.959040000000e+00, 2.100000000000e+00, 2.227840000000e+00, 2.340400000000e+00, 2.442960000000e+00, 2.544320000000e+00, 2.650000000000e+00,
2.759040000000e+00, 2.872160000000e+00, 2.986480000000e+00, 3.096800000000e+00, 3.200000000000e+00, 3.296960000000e+00, 3.387280000000e+00, 3.473600000000e+00,
3.560320000000e+00, 3.650000000000e+00, 3.742080000000e+00, 3.837200000000e+00, 3.931920000000e+00, 4.020640000000e+00, 4.100000000000e+00, 4.170640000000e+00,
4.231920000000e+00, 4.287200000000e+00, 4.342080000000e+00, 4.400000000000e+00, 4.460320000000e+00, 4.523600000000e+00, 4.587280000000e+00, 4.646960000000e+00,
4.700000000000e+00, 4.748160000000e+00, 4.793280000000e+00, 4.834400000000e+00, 4.869920000000e+00, 4.900000000000e+00, 4.925280000000e+00, 4.944000000000e+00,
4.968320000000e+00, 5.016640000000e+00, 5.100000000000e+00, 5.216000000000e+00, 5.365120000000e+00, 5.538240000000e+00, 5.720160000000e+00, 5.900000000000e+00,
6.078880000000e+00, 6.257440000000e+00, 6.427200000000e+00, 6.576960000000e+00, 6.700000000000e+00, 6.796560000000e+00, 6.865200000000e+00, 6.913840000000e+00,
6.955680000000e+00, 7.000000000000e+00, 7.045760000000e+00, 7.092640000000e+00, 7.146320000000e+00, 7.214000000000e+00, 7.300000000000e+00, 7.404720000000e+00,
7.529920000000e+00, 7.669120000000e+00, 7.811520000000e+00, 7.950000000000e+00, 8.085600000000e+00, 8.217920000000e+00, 8.347040000000e+00, 8.474160000000e+00,
8.600000000000e+00, 8.722960000000e+00, 8.841040000000e+00, 8.957120000000e+00, 9.076000000000e+00, 9.200000000000e+00, 9.328320000000e+00, 9.462320000000e+00,
9.593520000000e+00, 9.708720000000e+00, 9.800000000000e+00, 9.866960000000e+00, 9.904720000000e+00, 9.926480000000e+00, 9.954240000000e+00, 1.000000000000e+01,
1.006064000000e+01, 1.014048000000e+01, 1.021432000000e+01, 1.024216000000e+01, 1.020000000000e+01, 1.009536000000e+01, 9.930320000000e+00, 9.719280000000e+00,
9.486240000000e+00, 9.250000000000e+00, 9.009440000000e+00, 8.761280000000e+00, 8.535120000000e+00, 8.372960000000e+00, 8.300000000000e+00, 8.312880000000e+00,
8.414720000000e+00, 8.580560000000e+00, 8.768800000000e+00, 8.950000000000e+00, 9.127360000000e+00, 9.302160000000e+00, 9.454560000000e+00, 9.558960000000e+00,
9.600000000000e+00, 9.576080000000e+00, 9.480160000000e+00, 9.336240000000e+00, 9.184320000000e+00, 9.050000000000e+00, 8.929120000000e+00, 8.824240000000e+00,
8.727360000000e+00, 8.622480000000e+00, 8.500000000000e+00, 8.366080000000e+00, 8.225360000000e+00, 8.076640000000e+00, 7.917920000000e+00, 7.750000000000e+00,
7.573920000000e+00, 7.386640000000e+00, 7.209360000000e+00, 7.074080000000e+00, 7.000000000000e+00, 6.984240000000e+00, 7.030160000000e+00, 7.118080000000e+00,
7.215200000000e+00, 7.300000000000e+00, 7.375680000000e+00, 7.440880000000e+00, 7.496880000000e+00, 7.548880000000e+00, 7.600000000000e+00, 7.645840000000e+00,
7.681680000000e+00, 7.713520000000e+00, 7.751360000000e+00, 7.800000000000e+00, 7.857760000000e+00, 7.927520000000e+00, 7.991280000000e+00, 8.021040000000e+00,
8.000000000000e+00, 7.931040000000e+00, 7.811280000000e+00, 7.657520000000e+00, 7.497760000000e+00, 7.350000000000e+00, 7.211360000000e+00, 7.083520000000e+00,
6.961680000000e+00, 6.835840000000e+00, 6.700000000000e+00, 6.560160000000e+00, 6.423280000000e+00, 6.282400000000e+00, 6.125920000000e+00, 5.950000000000e+00,
5.757280000000e+00, 5.542000000000e+00, 5.342320000000e+00, 5.216640000000e+00, 5.200000000000e+00, 5.284560000000e+00, 5.471920000000e+00, 5.733280000000e+00,
6.020640000000e+00, 6.300000000000e+00, 6.575040000000e+00, 6.847280000000e+00, 7.093520000000e+00, 7.283760000000e+00, 7.400000000000e+00, 7.441440000000e+00,
7.401920000000e+00, 7.306400000000e+00, 7.196480000000e+00, 7.100000000000e+00, 7.013120000000e+00, 6.936800000000e+00, 6.874880000000e+00, 6.828960000000e+00,
6.800000000000e+00
};
static pair<double, double> CCT_to_xy_CIE_D(double CCT)
{
double xD;
if (CCT <= 7000.0)
{
xD = -4.607e9 / (CCT * CCT * CCT)
+ 2.9678e6 / (CCT * CCT)
+ 0.09911e3 / CCT
+ 0.244063;
}
else
{
xD = -2.0064e9 / (CCT * CCT * CCT)
+ 1.9018e6 / (CCT * CCT)
+ 0.24748e3 / CCT
+ 0.23704;
}
const double yD = -3.000 * xD * xD + 2.870 * xD - 0.275;
return {xD, yD};
}
VisData getSpd(double CCT)
{
const auto [xD, yD] = CCT_to_xy_CIE_D(CCT);
const double denom = 0.0241 + 0.2562 * xD - 0.7341 * yD;
const double M1 = (-1.3515 - 1.7703 * xD + 5.9114 * yD) / denom;
const double M2 = ( 0.0300 - 31.4424 * xD + 30.0717 * yD) / denom;
return views::iota(0, g_limitedVisSize) | views::transform([M1, M2](int i) {
return max(0.0, S0[i] + M1 * S1[i] + M2 * S2[i]);
}) | ranges::to<VisData>();
}
}