# Version 1.1 released by Yi-Chuan Lu on February 23, 2023.
#    Release 1.1 accommodates old Python 2 installations that
#    interpret some constants as integers instead of reals.
# Version 1.0 released by Yi-Chuan Lu on May 18, 2022.
# 
# When using this code, please cite:
# 
# @article{20heatindex,
#   Title   = {Extending the Heat Index},
#   Author  = {Yi-Chuan Lu and David M. Romps},
#   Journal = {Journal of Applied Meteorology and Climatology},
#   Year    = {2022},
#   Volume  = {61},
#   Number  = {10},
#   Pages   = {1367--1383},
#   Year    = {2022},
# }
#
# This headindex function returns the Heat Index in Kelvin. The inputs are:
# - T, the temperature in Kelvin
# - RH, the relative humidity, which is a value from 0 to 1
# - show_info is an optional logical flag. If true, the function returns the physiological state.

exec(open('heatindex.py').read())

error = 1e-8
if (all(list(map(lambda RH, HI: abs(heatindex(200,RH)/HI -1)<error, [0,.5,1], [199.9994020652,199.9997010342,200.0000000021]))) &
    all(list(map(lambda RH, HI: abs(heatindex(210,RH)/HI -1)<error, [0,.5,1], [209.9975943902,209.9987971085,209.9999998068]))) &
    all(list(map(lambda RH, HI: abs(heatindex(220,RH)/HI -1)<error, [0,.5,1], [219.9915822029,219.9957912306,219.9999999912]))) &
    all(list(map(lambda RH, HI: abs(heatindex(230,RH)/HI -1)<error, [0,.5,1], [229.9739691979,229.9869861009,230.0000001850]))) &
    all(list(map(lambda RH, HI: abs(heatindex(240,RH)/HI -1)<error, [0,.5,1], [239.9253828022,239.9626700074,240.0000000003]))) &
    all(list(map(lambda RH, HI: abs(heatindex(250,RH)/HI -1)<error, [0,.5,1], [249.7676757244,249.8837049107,250.0000000037]))) &
    all(list(map(lambda RH, HI: abs(heatindex(260,RH)/HI -1)<error, [0,.5,1], [259.3735990024,259.6864068902,259.9999999944]))) &
    all(list(map(lambda RH, HI: abs(heatindex(270,RH)/HI -1)<error, [0,.5,1], [268.5453870455,269.2745889562,270.0000002224]))) &
    all(list(map(lambda RH, HI: abs(heatindex(280,RH)/HI -1)<error, [0,.5,1], [277.2234200026,278.6369451963,280.0000000091]))) &
    all(list(map(lambda RH, HI: abs(heatindex(290,RH)/HI -1)<error, [0,.5,1], [285.7510545370,288.2813660100,290.7860610129]))) &
    all(list(map(lambda RH, HI: abs(heatindex(300,RH)/HI -1)<error, [0,.5,1], [297.5737503539,300.2922595865,305.3947127590]))) &
    all(list(map(lambda RH, HI: abs(heatindex(310,RH)/HI -1)<error, [0,.5,1], [305.5549530893,318.6225524695,359.9063248191]))) &
    all(list(map(lambda RH, HI: abs(heatindex(320,RH)/HI -1)<error, [0,.5,1], [313.0298872791,359.0538750602,407.5345212438]))) &
    all(list(map(lambda RH, HI: abs(heatindex(330,RH)/HI -1)<error, [0,.5,1], [320.5088548469,398.5759733823,464.9949352940]))) &
    all(list(map(lambda RH, HI: abs(heatindex(340,RH)/HI -1)<error, [0,.5,1], [328.0358006469,445.8599463105,530.5524786708]))) &
    all(list(map(lambda RH, HI: abs(heatindex(350,RH)/HI -1)<error, [0,.5,1], [333.2806160592,500.0421800191,601.9518435268]))) &
    all(list(map(lambda RH, HI: abs(heatindex(360,RH)/HI -1)<error, [0,.5,1], [343.6312984164,559.6640227151,677.2462089759]))) &
    all(list(map(lambda RH, HI: abs(heatindex(370,RH)/HI -1)<error, [0,.5,1], [354.1825692377,623.1960299857,755.0832658147])))):
   print('Success')
else:
   print('Failure')