CIE1931 correction table

STEP_NUM = 100
MAX_OUTPUT = 9990
INT_TYPE = 'unsigned int'
TABLE_NAME = 'cie_table'
TRIAC_DELAY = 10

def cie1931(L):
    L = L*100.0
    if L <= 8:
        return (L/902.3)
    else:
        return ((L+16.0)/116.0)**3

x = range(0,int(STEP_NUM+1))
y = [round(cie1931(float(L)/STEP_NUM)*MAX_OUTPUT) for L in x]


f = open('cie1931.h', 'w')
textOut = '// CIE1931 correction table\n'
textOut += '// Automatically generated\n'

textOut += '%s %s[%d] = {\n' % (INT_TYPE, TABLE_NAME, STEP_NUM+1)
textOut += '\t'
for i, L in enumerate(reversed(y)):
    val = MAX_OUTPUT - L
    textOut += '%d, ' % int(val)
    if i % 10 == 9:
        textOut += '\n\t'

textOut += '\n};\n\n'
print(textOut)