diff --git a/src/scgenerator/data/k_silica.txt b/src/scgenerator/data/k_silica.txt new file mode 100644 index 0000000..8c061ce --- /dev/null +++ b/src/scgenerator/data/k_silica.txt @@ -0,0 +1,1443 @@ +6.199e-10,1.503e-05 +6.351e-10,1.636e-05 +6.508e-10,1.781e-05 +6.665e-10,1.936e-05 +6.831e-10,6.3e-06 +6.997e-10,6.904e-06 +7.166e-10,7.564e-06 +7.345e-10,8.298e-06 +7.523e-10,9.083e-06 +7.705e-10,9.941e-06 +7.897e-10,1.09e-05 +8.087e-10,1.193e-05 +8.287e-10,1.308e-05 +8.492e-10,1.432e-05 +8.694e-10,1.562e-05 +8.906e-10,1.708e-05 +9.129e-10,1.872e-05 +9.35e-10,2.044e-05 +9.581e-10,2.238e-05 +9.816e-10,2.447e-05 +1.005e-09,2.673e-05 +1.03e-09,2.916e-05 +1.055e-09,3.189e-05 +1.081e-09,3.483e-05 +1.108e-09,3.812e-05 +1.134e-09,4.155e-05 +1.162e-09,4.537e-05 +1.191e-09,4.965e-05 +1.22e-09,5.423e-05 +1.25e-09,5.915e-05 +1.281e-09,6.468e-05 +1.312e-09,7.061e-05 +1.343e-09,7.686e-05 +1.376e-09,8.364e-05 +1.41e-09,9.121e-05 +1.445e-09,9.928e-05 +1.479e-09,0.000108 +1.516e-09,0.0001179 +1.554e-09,0.0001289 +1.592e-09,0.0001407 +1.631e-09,0.0001535 +1.671e-09,0.0001671 +1.71e-09,0.0001813 +1.754e-09,0.0001981 +1.797e-09,0.0002149 +1.839e-09,0.0002326 +1.884e-09,0.0002536 +1.931e-09,0.0002771 +1.977e-09,0.0003009 +2.026e-09,0.000327 +2.077e-09,0.000356 +2.127e-09,0.0003862 +2.179e-09,0.0004199 +2.234e-09,0.0004578 +2.287e-09,0.0004971 +2.344e-09,0.0001486 +2.403e-09,0.0001621 +2.46e-09,0.000176 +2.52e-09,0.000191 +2.583e-09,0.000208 +2.643e-09,0.000226 +2.707e-09,0.000245 +2.774e-09,0.000267 +2.844e-09,0.00029 +2.91e-09,0.000315 +2.987e-09,0.000344 +3.054e-09,0.000371 +3.131e-09,0.000404 +3.212e-09,0.00044 +3.289e-09,0.000477 +3.369e-09,0.000518 +3.453e-09,0.000563 +3.532e-09,0.000607 +3.625e-09,0.000662 +3.712e-09,0.000717 +3.803e-09,0.000777 +3.899e-09,0.000844 +3.986e-09,0.000909 +4.092e-09,0.000991 +4.188e-09,0.00107 +4.29e-09,0.00116 +4.396e-09,0.00125 +4.508e-09,0.00136 +4.609e-09,0.00146 +4.732e-09,0.00159 +4.843e-09,0.00171 +4.959e-09,0.00185 +5.081e-09,0.002 +5.209e-09,0.00216 +5.344e-09,0.00234 +5.462e-09,0.0025 +5.61e-09,0.00272 +5.74e-09,0.00292 +5.876e-09,0.00314 +6.018e-09,0.00337 +6.168e-09,0.00363 +6.325e-09,0.0039 +6.491e-09,0.00421 +6.63e-09,0.00447 +6.812e-09,0.00482 +6.965e-09,0.00511 +7.125e-09,0.00544 +7.293e-09,0.0058 +7.469e-09,0.00619 +7.653e-09,0.00660671 +7.68245e-09,0.00667199 +7.71201e-09,0.00673676 +7.74168e-09,0.00680055 +7.77147e-09,0.00686286 +7.80137e-09,0.00692316 +7.83138e-09,0.00698094 +7.86152e-09,0.00703564 +7.89177e-09,0.00708671 +7.92213e-09,0.00713359 +7.95261e-09,0.00717571 +7.98321e-09,0.00721249 +8.01393e-09,0.00724338 +8.04476e-09,0.0072678 +8.07571e-09,0.00728547 +8.10679e-09,0.00729822 +8.13798e-09,0.0073088 +8.16929e-09,0.00732001 +8.20072e-09,0.00733466 +8.23228e-09,0.00735512 +8.26395e-09,0.00738198 +8.29575e-09,0.00741543 +8.32767e-09,0.0074557 +8.35971e-09,0.00750303 +8.39188e-09,0.00755768 +8.42416e-09,0.00761994 +8.45658e-09,0.00769012 +8.48912e-09,0.00776856 +8.52178e-09,0.00785565 +8.55457e-09,0.00795179 +8.58748e-09,0.00805744 +8.62052e-09,0.00817301 +8.65369e-09,0.00829819 +8.68699e-09,0.00843201 +8.72041e-09,0.00857346 +8.75397e-09,0.00872146 +8.78765e-09,0.00887485 +8.82146e-09,0.00903238 +8.8554e-09,0.00919269 +8.88947e-09,0.00935436 +8.92368e-09,0.00951581 +8.95801e-09,0.00967541 +8.99248e-09,0.00983139 +9.02708e-09,0.00998224 +9.06181e-09,0.010129 +9.09668e-09,0.0102742 +9.13168e-09,0.0104204 +9.16682e-09,0.0105703 +9.20209e-09,0.0107268 +9.23749e-09,0.0108901 +9.27304e-09,0.0110501 +9.30871e-09,0.0111939 +9.34453e-09,0.0113082 +9.38049e-09,0.0113793 +9.41658e-09,0.0113939 +9.45281e-09,0.0113554 +9.48918e-09,0.0112851 +9.52569e-09,0.0112051 +9.56234e-09,0.0111369 +9.59914e-09,0.0111018 +9.63607e-09,0.011113 +9.67315e-09,0.0111511 +9.71037e-09,0.0111876 +9.74773e-09,0.0111941 +9.78523e-09,0.0111419 +9.82288e-09,0.0110056 +9.86068e-09,0.0107887 +9.89862e-09,0.0105163 +9.93671e-09,0.0102126 +9.97494e-09,0.00990006 +1.00133e-08,0.00959856 +1.00518e-08,0.00932001 +1.00905e-08,0.0090665 +1.01293e-08,0.00883923 +1.01683e-08,0.00863931 +1.02074e-08,0.00846784 +1.02467e-08,0.00832891 +1.02861e-08,0.00823321 +1.03257e-08,0.00819238 +1.03655e-08,0.00821871 +1.04053e-08,0.00832612 +1.04454e-08,0.00852018 +1.04856e-08,0.00877822 +1.05259e-08,0.00907073 +1.05664e-08,0.00936449 +1.06071e-08,0.00962189 +1.06479e-08,0.00981285 +1.06888e-08,0.00993346 +1.073e-08,0.00998539 +1.07713e-08,0.00997187 +1.08127e-08,0.00989745 +1.08543e-08,0.00976771 +1.08961e-08,0.00958899 +1.0938e-08,0.00936812 +1.09801e-08,0.00911211 +1.10223e-08,0.0088307 +1.10647e-08,0.0085869 +1.11073e-08,0.00848671 +1.115e-08,0.00863764 +1.11929e-08,0.00917269 +1.1236e-08,0.0102529 +1.12792e-08,0.011824 +1.13226e-08,0.0136312 +1.13662e-08,0.0152199 +1.14099e-08,0.0159479 +1.14538e-08,0.0154555 +1.14979e-08,0.0141406 +1.15422e-08,0.0125214 +1.15866e-08,0.0110014 +1.16311e-08,0.00981612 +1.16759e-08,0.00894896 +1.17208e-08,0.00830204 +1.17659e-08,0.00780469 +1.18112e-08,0.00740409 +1.18566e-08,0.00707182 +1.19023e-08,0.00680972 +1.1948e-08,0.00662284 +1.1994e-08,0.00651712 +1.20402e-08,0.00649628 +1.20865e-08,0.00654223 +1.2133e-08,0.00662961 +1.21797e-08,0.00673287 +1.22265e-08,0.00682613 +1.22736e-08,0.00689946 +1.23208e-08,0.00695735 +1.23682e-08,0.00700506 +1.24158e-08,0.00704808 +1.24636e-08,0.00709203 +1.25115e-08,0.00714272 +1.25597e-08,0.0072061 +1.2608e-08,0.00728837 +1.26565e-08,0.0073892 +1.27052e-08,0.00749106 +1.27541e-08,0.00757285 +1.28032e-08,0.00761267 +1.28524e-08,0.00760084 +1.29019e-08,0.00756848 +1.29515e-08,0.00755482 +1.30014e-08,0.00759898 +1.30514e-08,0.00772735 +1.31016e-08,0.00791519 +1.3152e-08,0.00812389 +1.32026e-08,0.00831101 +1.32534e-08,0.00844411 +1.33044e-08,0.008537 +1.33556e-08,0.00861526 +1.3407e-08,0.00870566 +1.34586e-08,0.00882831 +1.35104e-08,0.00898314 +1.35623e-08,0.0091667 +1.36145e-08,0.0093755 +1.36669e-08,0.00960703 +1.37195e-08,0.00986092 +1.37723e-08,0.010137 +1.38253e-08,0.0104351 +1.38785e-08,0.0107549 +1.39319e-08,0.0110962 +1.39855e-08,0.0114585 +1.40393e-08,0.0118412 +1.40933e-08,0.0122436 +1.41475e-08,0.0126649 +1.4202e-08,0.0131039 +1.42566e-08,0.0135593 +1.43115e-08,0.0140295 +1.43665e-08,0.0145126 +1.44218e-08,0.0150064 +1.44773e-08,0.0155085 +1.4533e-08,0.0160159 +1.45889e-08,0.0165255 +1.4645e-08,0.0170337 +1.47014e-08,0.0175365 +1.4758e-08,0.0180297 +1.48147e-08,0.0185088 +1.48717e-08,0.0189688 +1.4929e-08,0.0194046 +1.49864e-08,0.0198108 +1.50441e-08,0.0201825 +1.5102e-08,0.0205183 +1.51601e-08,0.0208194 +1.52184e-08,0.0210871 +1.52769e-08,0.0213231 +1.53357e-08,0.0215294 +1.53947e-08,0.0217082 +1.5454e-08,0.0218621 +1.55134e-08,0.0219935 +1.55731e-08,0.0221055 +1.5633e-08,0.0222008 +1.56932e-08,0.0222824 +1.57536e-08,0.0223535 +1.58142e-08,0.0224173 +1.5875e-08,0.0224769 +1.59361e-08,0.0225355 +1.59974e-08,0.0225965 +1.6059e-08,0.0226627 +1.61208e-08,0.0227347 +1.61828e-08,0.0228127 +1.62451e-08,0.0228969 +1.63076e-08,0.0229872 +1.63703e-08,0.0230838 +1.64333e-08,0.0231869 +1.64965e-08,0.0232966 +1.656e-08,0.0234131 +1.66237e-08,0.0235364 +1.66877e-08,0.0236667 +1.67519e-08,0.0238042 +1.68163e-08,0.0239491 +1.6881e-08,0.0241015 +1.6946e-08,0.0242616 +1.70112e-08,0.0244297 +1.70767e-08,0.0246059 +1.71424e-08,0.0247906 +1.72083e-08,0.0249844 +1.72745e-08,0.0251877 +1.7341e-08,0.0254009 +1.74077e-08,0.0256245 +1.74747e-08,0.0258591 +1.75419e-08,0.0261051 +1.76094e-08,0.0263633 +1.76772e-08,0.026634 +1.77452e-08,0.0269181 +1.78135e-08,0.027216 +1.7882e-08,0.0275286 +1.79508e-08,0.0278565 +1.80199e-08,0.0282005 +1.80892e-08,0.0285604 +1.81588e-08,0.0289342 +1.82287e-08,0.0293198 +1.82988e-08,0.0297147 +1.83692e-08,0.0301166 +1.84399e-08,0.0305229 +1.85109e-08,0.0309307 +1.85821e-08,0.0313371 +1.86536e-08,0.0317391 +1.87254e-08,0.0321336 +1.87974e-08,0.0325171 +1.88697e-08,0.0328862 +1.89423e-08,0.0332373 +1.90152e-08,0.0335668 +1.90884e-08,0.033873 +1.91618e-08,0.0341596 +1.92356e-08,0.0344306 +1.93096e-08,0.0346905 +1.93839e-08,0.034944 +1.94584e-08,0.0351958 +1.95333e-08,0.0354507 +1.96085e-08,0.0357138 +1.96839e-08,0.0359903 +1.97597e-08,0.0362854 +1.98357e-08,0.0366048 +1.9912e-08,0.0369541 +1.99886e-08,0.0373394 +2.00655e-08,0.0377658 +2.01427e-08,0.03823 +2.02202e-08,0.0387259 +2.0298e-08,0.0392466 +2.03761e-08,0.0397851 +2.04545e-08,0.040334 +2.05332e-08,0.0408854 +2.06122e-08,0.041431 +2.06916e-08,0.0419619 +2.07712e-08,0.0424692 +2.08511e-08,0.0429432 +2.09313e-08,0.043374 +2.10118e-08,0.0437518 +2.10927e-08,0.0440707 +2.11739e-08,0.0443375 +2.12553e-08,0.044561 +2.13371e-08,0.0447503 +2.14192e-08,0.044915 +2.15016e-08,0.0450648 +2.15843e-08,0.0452094 +2.16674e-08,0.0453589 +2.17508e-08,0.0455233 +2.18345e-08,0.0457129 +2.19185e-08,0.0459382 +2.20028e-08,0.0462098 +2.20875e-08,0.0465359 +2.21724e-08,0.0469134 +2.22578e-08,0.0473371 +2.23434e-08,0.0478015 +2.24294e-08,0.0483009 +2.25157e-08,0.0488296 +2.26023e-08,0.0493815 +2.26893e-08,0.0499503 +2.27766e-08,0.0505294 +2.28642e-08,0.0511119 +2.29522e-08,0.0516906 +2.30405e-08,0.052258 +2.31291e-08,0.0528109 +2.32181e-08,0.053351 +2.33075e-08,0.0538802 +2.33971e-08,0.0544006 +2.34872e-08,0.0549146 +2.35775e-08,0.0554244 +2.36683e-08,0.0559325 +2.37593e-08,0.0564413 +2.38507e-08,0.0569536 +2.39425e-08,0.0574719 +2.40346e-08,0.0579991 +2.41271e-08,0.0585387 +2.42199e-08,0.0590952 +2.43131e-08,0.0596734 +2.44067e-08,0.060278 +2.45006e-08,0.0609143 +2.45949e-08,0.0615875 +2.46895e-08,0.0623032 +2.47845e-08,0.0630673 +2.48798e-08,0.063886 +2.49756e-08,0.0647658 +2.50717e-08,0.0657114 +2.51681e-08,0.0667062 +2.5265e-08,0.067721 +2.53622e-08,0.068725 +2.54598e-08,0.0696855 +2.55577e-08,0.0705678 +2.56561e-08,0.0713362 +2.57548e-08,0.0719539 +2.58539e-08,0.0723838 +2.59534e-08,0.0725892 +2.60532e-08,0.0725374 +2.61535e-08,0.0722429 +2.62541e-08,0.0717632 +2.63551e-08,0.0711577 +2.64565e-08,0.0704846 +2.65583e-08,0.0698006 +2.66605e-08,0.0691601 +2.67631e-08,0.0686158 +2.6866e-08,0.0682188 +2.69694e-08,0.0680197 +2.70732e-08,0.0680648 +2.71774e-08,0.0683493 +2.72819e-08,0.0688396 +2.73869e-08,0.0695029 +2.74923e-08,0.0703062 +2.7598e-08,0.0712161 +2.77042e-08,0.0721976 +2.78108e-08,0.0732139 +2.79178e-08,0.0742264 +2.80253e-08,0.0751941 +2.81331e-08,0.0760912 +2.82413e-08,0.076926 +2.835e-08,0.0777112 +2.84591e-08,0.0784607 +2.85686e-08,0.0791888 +2.86785e-08,0.0799105 +2.87888e-08,0.0806413 +2.88996e-08,0.0813974 +2.90108e-08,0.0821955 +2.91224e-08,0.0830457 +2.92345e-08,0.0839365 +2.9347e-08,0.0848516 +2.94599e-08,0.085774 +2.95732e-08,0.0866862 +2.9687e-08,0.0875697 +2.98012e-08,0.0884054 +2.99159e-08,0.0891736 +3.0031e-08,0.089854 +3.01466e-08,0.0904356 +3.02626e-08,0.0909248 +3.0379e-08,0.0913303 +3.04959e-08,0.0916614 +3.06132e-08,0.0919274 +3.0731e-08,0.0921384 +3.08493e-08,0.0923041 +3.0968e-08,0.0924349 +3.10871e-08,0.0925409 +3.12067e-08,0.0926326 +3.13268e-08,0.0927202 +3.14473e-08,0.0928143 +3.15683e-08,0.0929254 +3.16898e-08,0.0930639 +3.18117e-08,0.0932405 +3.19341e-08,0.0934659 +3.2057e-08,0.0937509 +3.21803e-08,0.0941025 +3.23042e-08,0.0945236 +3.24284e-08,0.095017 +3.25532e-08,0.0955857 +3.26785e-08,0.096233 +3.28042e-08,0.0969624 +3.29304e-08,0.0977777 +3.30571e-08,0.0986825 +3.31843e-08,0.0996713 +3.3312e-08,0.100727 +3.34402e-08,0.101831 +3.35688e-08,0.102965 +3.3698e-08,0.104109 +3.38277e-08,0.105242 +3.39578e-08,0.106343 +3.40885e-08,0.10739 +3.42196e-08,0.108379 +3.43513e-08,0.109316 +3.44835e-08,0.11021 +3.46162e-08,0.111067 +3.47493e-08,0.111896 +3.4883e-08,0.112703 +3.50173e-08,0.113499 +3.5152e-08,0.114291 +3.52873e-08,0.115088 +3.5423e-08,0.1159 +3.55593e-08,0.116736 +3.56961e-08,0.117605 +3.58335e-08,0.118518 +3.59714e-08,0.119484 +3.61098e-08,0.120514 +3.62487e-08,0.121604 +3.63882e-08,0.122749 +3.65282e-08,0.123939 +3.66687e-08,0.125168 +3.68098e-08,0.126428 +3.69514e-08,0.127708 +3.70936e-08,0.129002 +3.72363e-08,0.1303 +3.73796e-08,0.131592 +3.75234e-08,0.132868 +3.76678e-08,0.134118 +3.78128e-08,0.135331 +3.79582e-08,0.136496 +3.81043e-08,0.137605 +3.82509e-08,0.13867 +3.83981e-08,0.13972 +3.85458e-08,0.140782 +3.86941e-08,0.141886 +3.8843e-08,0.143061 +3.89925e-08,0.144338 +3.91425e-08,0.145744 +3.92931e-08,0.147275 +3.94443e-08,0.14892 +3.95961e-08,0.150666 +3.97484e-08,0.152501 +3.99013e-08,0.15441 +4.00549e-08,0.156381 +4.0209e-08,0.158398 +4.03637e-08,0.160446 +4.0519e-08,0.162508 +4.06749e-08,0.164568 +4.08314e-08,0.166606 +4.09885e-08,0.168605 +4.11462e-08,0.170547 +4.13045e-08,0.172433 +4.14635e-08,0.17427 +4.1623e-08,0.176064 +4.17832e-08,0.177821 +4.19439e-08,0.17955 +4.21053e-08,0.181257 +4.22673e-08,0.182951 +4.24299e-08,0.184641 +4.25932e-08,0.186333 +4.27571e-08,0.188039 +4.29216e-08,0.189766 +4.30867e-08,0.191526 +4.32525e-08,0.193326 +4.34189e-08,0.195179 +4.3586e-08,0.197094 +4.37537e-08,0.199084 +4.39221e-08,0.201158 +4.40911e-08,0.203331 +4.42607e-08,0.205605 +4.4431e-08,0.207978 +4.4602e-08,0.210448 +4.47736e-08,0.213013 +4.49458e-08,0.215668 +4.51188e-08,0.218411 +4.52924e-08,0.221239 +4.54666e-08,0.224148 +4.56416e-08,0.227134 +4.58172e-08,0.230194 +4.59935e-08,0.233323 +4.61705e-08,0.236515 +4.63481e-08,0.239755 +4.65264e-08,0.243027 +4.67054e-08,0.246312 +4.68852e-08,0.249591 +4.70656e-08,0.252844 +4.72466e-08,0.25605 +4.74284e-08,0.259186 +4.76109e-08,0.26223 +4.77941e-08,0.265159 +4.7978e-08,0.267946 +4.81626e-08,0.270579 +4.83479e-08,0.273107 +4.85339e-08,0.275606 +4.87207e-08,0.278154 +4.89081e-08,0.280834 +4.90963e-08,0.283725 +4.92852e-08,0.286862 +4.94749e-08,0.290224 +4.96652e-08,0.29379 +4.98563e-08,0.297537 +5.00481e-08,0.301441 +5.02407e-08,0.305477 +5.0434e-08,0.309618 +5.06281e-08,0.313836 +5.08229e-08,0.318101 +5.10184e-08,0.322383 +5.12147e-08,0.326654 +5.14118e-08,0.330914 +5.16096e-08,0.335166 +5.18082e-08,0.339411 +5.20075e-08,0.343654 +5.22076e-08,0.347897 +5.24085e-08,0.352143 +5.26101e-08,0.356397 +5.28126e-08,0.360662 +5.30158e-08,0.364942 +5.32198e-08,0.369242 +5.34245e-08,0.373566 +5.36301e-08,0.377919 +5.38364e-08,0.382306 +5.40436e-08,0.386732 +5.42515e-08,0.391202 +5.44603e-08,0.395723 +5.46698e-08,0.4003 +5.48802e-08,0.404939 +5.50913e-08,0.409647 +5.53033e-08,0.414417 +5.55161e-08,0.419228 +5.57297e-08,0.424056 +5.59441e-08,0.428877 +5.61594e-08,0.433666 +5.63754e-08,0.438395 +5.65924e-08,0.443038 +5.68101e-08,0.447564 +5.70287e-08,0.451946 +5.72481e-08,0.456152 +5.74684e-08,0.460152 +5.76895e-08,0.463912 +5.79115e-08,0.467402 +5.81343e-08,0.470593 +5.8358e-08,0.473533 +5.85825e-08,0.476332 +5.88079e-08,0.4791 +5.90342e-08,0.481953 +5.92613e-08,0.48501 +5.94893e-08,0.488391 +5.97182e-08,0.492222 +5.9948e-08,0.496631 +6.01787e-08,0.501739 +6.04102e-08,0.507541 +6.06427e-08,0.513968 +6.0876e-08,0.520953 +6.11102e-08,0.528422 +6.13454e-08,0.5363 +6.15814e-08,0.544505 +6.18183e-08,0.552951 +6.20562e-08,0.561548 +6.2295e-08,0.570198 +6.25346e-08,0.578798 +6.27753e-08,0.587238 +6.30168e-08,0.595403 +6.32593e-08,0.603172 +6.35027e-08,0.610419 +6.3747e-08,0.617026 +6.39923e-08,0.622989 +6.42385e-08,0.628382 +6.44857e-08,0.633283 +6.47338e-08,0.637777 +6.49828e-08,0.641949 +6.52329e-08,0.645892 +6.54839e-08,0.649693 +6.57358e-08,0.653431 +6.59888e-08,0.657178 +6.62427e-08,0.661008 +6.64975e-08,0.664999 +6.67534e-08,0.669227 +6.70102e-08,0.673775 +6.72681e-08,0.6787 +6.75269e-08,0.683963 +6.77867e-08,0.689502 +6.80475e-08,0.69525 +6.83094e-08,0.70114 +6.85722e-08,0.707101 +6.8836e-08,0.71306 +6.91009e-08,0.718943 +6.93668e-08,0.724672 +6.96337e-08,0.730167 +6.99016e-08,0.735347 +7.01705e-08,0.740148 +7.04405e-08,0.744556 +7.07116e-08,0.748561 +7.09836e-08,0.752155 +7.12568e-08,0.755331 +7.15309e-08,0.758082 +7.18061e-08,0.760403 +7.20824e-08,0.762287 +7.23598e-08,0.763731 +7.26382e-08,0.76473 +7.29177e-08,0.765282 +7.31982e-08,0.765385 +7.34799e-08,0.765036 +7.37626e-08,0.764235 +7.40464e-08,0.762982 +7.43313e-08,0.761279 +7.46173e-08,0.759127 +7.49044e-08,0.756528 +7.51926e-08,0.753486 +7.54819e-08,0.750025 +7.57724e-08,0.746213 +7.60639e-08,0.742123 +7.63566e-08,0.737829 +7.66504e-08,0.733403 +7.69453e-08,0.728915 +7.72414e-08,0.724432 +7.75386e-08,0.720022 +7.78369e-08,0.71575 +7.81364e-08,0.71168 +7.8437e-08,0.707873 +7.87388e-08,0.704391 +7.90418e-08,0.701294 +7.93459e-08,0.698642 +7.96512e-08,0.696494 +7.99577e-08,0.69491 +8.02653e-08,0.69395 +8.05742e-08,0.693675 +8.08842e-08,0.69415 +8.11954e-08,0.695439 +8.15078e-08,0.697604 +8.18214e-08,0.700622 +8.21362e-08,0.704416 +8.24523e-08,0.708906 +8.27695e-08,0.714014 +8.3088e-08,0.719662 +8.34077e-08,0.725768 +8.37286e-08,0.732251 +8.40508e-08,0.739024 +8.43742e-08,0.745999 +8.46988e-08,0.753084 +8.50247e-08,0.760184 +8.53518e-08,0.7672 +8.56802e-08,0.774029 +8.60099e-08,0.780566 +8.63408e-08,0.786701 +8.6673e-08,0.792323 +8.70065e-08,0.797318 +8.73413e-08,0.801571 +8.76774e-08,0.804967 +8.80147e-08,0.807392 +8.83534e-08,0.808733 +8.86933e-08,0.808885 +8.90346e-08,0.807847 +8.93771e-08,0.805759 +8.9721e-08,0.802776 +9.00663e-08,0.799051 +9.04128e-08,0.794738 +9.07607e-08,0.789989 +9.11099e-08,0.784953 +9.14604e-08,0.779778 +9.18124e-08,0.774606 +9.21656e-08,0.769576 +9.25202e-08,0.764824 +9.28762e-08,0.760483 +9.32336e-08,0.756682 +9.35923e-08,0.753548 +9.39524e-08,0.751181 +9.43139e-08,0.749579 +9.46768e-08,0.748724 +9.50411e-08,0.7486 +9.54068e-08,0.749194 +9.57739e-08,0.750493 +9.61424e-08,0.752487 +9.65123e-08,0.755168 +9.68836e-08,0.758531 +9.72564e-08,0.762569 +9.76306e-08,0.767267 +9.80063e-08,0.77256 +9.83833e-08,0.778371 +9.87619e-08,0.784623 +9.91419e-08,0.791238 +9.95234e-08,0.798133 +9.99063e-08,0.805224 +1.00291e-07,0.812423 +1.00677e-07,0.819641 +1.01064e-07,0.826782 +1.01453e-07,0.83375 +1.01843e-07,0.840445 +1.02235e-07,0.846764 +1.02628e-07,0.852601 +1.03023e-07,0.857849 +1.0342e-07,0.8624 +1.03818e-07,0.866144 +1.04217e-07,0.868972 +1.04618e-07,0.870778 +1.05021e-07,0.871457 +1.05425e-07,0.870908 +1.0583e-07,0.869035 +1.06237e-07,0.865749 +1.06646e-07,0.86097 +1.07057e-07,0.854627 +1.07468e-07,0.846727 +1.07882e-07,0.837378 +1.08297e-07,0.826697 +1.08714e-07,0.814804 +1.09132e-07,0.801652 +1.09552e-07,0.785938 +1.09973e-07,0.765927 +1.10397e-07,0.740129 +1.10821e-07,0.707439 +1.11248e-07,0.669509 +1.11676e-07,0.631756 +1.12106e-07,0.598926 +1.12537e-07,0.574823 +1.1297e-07,0.562513 +1.13405e-07,0.561625 +1.13841e-07,0.569795 +1.14279e-07,0.585883 +1.14719e-07,0.610436 +1.1516e-07,0.644491 +1.15603e-07,0.689533 +1.16048e-07,0.747564 +1.16494e-07,0.819337 +1.16943e-07,0.902322 +1.17393e-07,0.991225 +1.17844e-07,1.07244 +1.18298e-07,1.12644 +1.18753e-07,1.14264 +1.1921e-07,1.1233 +1.19668e-07,1.06806 +1.20129e-07,0.9813 +1.20591e-07,0.888299 +1.21055e-07,0.800985 +1.21521e-07,0.709068 +1.21988e-07,0.610384 +1.22458e-07,0.517921 +1.22929e-07,0.441287 +1.23402e-07,0.381638 +1.23877e-07,0.333868 +1.24353e-07,0.294035 +1.24832e-07,0.259658 +1.25312e-07,0.229119 +1.25794e-07,0.201745 +1.26278e-07,0.177819 +1.26764e-07,0.157387 +1.27252e-07,0.139099 +1.27742e-07,0.121108 +1.28233e-07,0.102909 +1.28726e-07,0.0867397 +1.29222e-07,0.074374 +1.29719e-07,0.065736 +1.30218e-07,0.0592711 +1.30719e-07,0.0538154 +1.31222e-07,0.0488764 +1.31727e-07,0.044438 +1.32234e-07,0.040494 +1.32743e-07,0.0370279 +1.33253e-07,0.0340164 +1.33766e-07,0.0314246 +1.34281e-07,0.0291787 +1.34797e-07,0.0272061 +1.35316e-07,0.0254484 +1.35837e-07,0.0238583 +1.36359e-07,0.0223972 +1.36884e-07,0.0210335 +1.37411e-07,0.0197416 +1.37939e-07,0.0185013 +1.3847e-07,0.017315 +1.39003e-07,0.0162163 +1.39538e-07,0.0152324 +1.40075e-07,0.0143569 +1.40614e-07,0.0135555 +1.41155e-07,0.0127983 +1.41698e-07,0.0120577 +1.42243e-07,0.011308 +1.4279e-07,0.0105305 +1.4334e-07,0.00971361 +1.43891e-07,0.00885347 +1.44445e-07,0.00795408 +1.45001e-07,0.00703333 +1.45559e-07,0.00611857 +1.46119e-07,0.00523482 +1.46681e-07,0.00439974 +1.47245e-07,0.00362486 +1.47812e-07,0.00292099 +1.4838e-07,0.00229646 +1.48951e-07,0.00175642 +1.49524e-07,0.00130314 +1.501e-07,0.000937126 +1.50677e-07,0.00065561 +1.51257e-07,0.0004481 +1.51839e-07,0.000300489 +1.52423e-07,0.00019854 +1.5301e-07,0.0001298 +1.53598e-07,8.44066e-05 +1.54189e-07,5.51006e-05 +1.54783e-07,2.93033e-05 +1.55378e-07,8.87433e-06 +1.55976e-07,2.10754e-06 +1.56576e-07,5.35749e-07 +1.57179e-07,1.69996e-07 +1.57783e-07,8.31532e-08 +1.58391e-07,7.01033e-08 +1.59e-07,7e-08 +3.636e-06,7e-08 +3.64998e-06,7.00477e-08 +3.66402e-06,7.95381e-08 +3.67812e-06,1.49675e-07 +3.69226e-06,4.39428e-07 +3.70646e-06,1.54688e-06 +3.72072e-06,5.50399e-06 +3.73503e-06,1.66904e-05 +3.7494e-06,3.27504e-05 +3.76382e-06,3.83984e-05 +3.7783e-06,3.97792e-05 +3.79283e-06,4.16614e-05 +3.80742e-06,4.39037e-05 +3.82206e-06,4.62659e-05 +3.83676e-06,4.84523e-05 +3.85152e-06,5.01191e-05 +3.86633e-06,5.11057e-05 +3.8812e-06,5.15658e-05 +3.89613e-06,5.1703e-05 +3.91112e-06,5.17329e-05 +3.92616e-06,5.18735e-05 +3.94126e-06,5.22909e-05 +3.95642e-06,5.30591e-05 +3.97164e-06,5.42513e-05 +3.98691e-06,5.59554e-05 +4.00225e-06,5.82797e-05 +4.01764e-06,6.13075e-05 +4.03309e-06,6.49986e-05 +4.04861e-06,6.92787e-05 +4.06418e-06,7.4048e-05 +4.07981e-06,7.91691e-05 +4.0955e-06,8.44717e-05 +4.11125e-06,8.98251e-05 +4.12707e-06,9.51039e-05 +4.14294e-06,0.000100161 +4.15888e-06,0.000104828 +4.17487e-06,0.000108935 +4.19093e-06,0.000112584 +4.20705e-06,0.000116206 +4.22323e-06,0.000120306 +4.23947e-06,0.000125464 +4.25578e-06,0.000132373 +4.27215e-06,0.000141687 +4.28858e-06,0.000153473 +4.30507e-06,0.000167621 +4.32163e-06,0.000183924 +4.33825e-06,0.000202018 +4.35494e-06,0.000221302 +4.37169e-06,0.000240599 +4.38851e-06,0.000257936 +4.40539e-06,0.000270993 +4.42233e-06,0.000279198 +4.43934e-06,0.000283937 +4.45641e-06,0.000286751 +4.47355e-06,0.000287327 +4.49076e-06,0.000284324 +4.50803e-06,0.000276912 +4.52537e-06,0.000267002 +4.54278e-06,0.000257414 +4.56025e-06,0.000250347 +4.57779e-06,0.000246264 +4.5954e-06,0.000245121 +4.61307e-06,0.000246975 +4.63082e-06,0.000251994 +4.64863e-06,0.000260472 +4.66651e-06,0.000272931 +4.68445e-06,0.000290557 +4.70247e-06,0.000315147 +4.72056e-06,0.000349233 +4.73872e-06,0.000396513 +4.75694e-06,0.000462548 +4.77524e-06,0.000555512 +4.79361e-06,0.000683191 +4.81204e-06,0.00085231 +4.83055e-06,0.00106841 +4.84913e-06,0.00133303 +4.86778e-06,0.00163975 +4.8865e-06,0.00197029 +4.9053e-06,0.00230586 +4.92417e-06,0.00263969 +4.94311e-06,0.00296979 +4.96212e-06,0.00329899 +4.9812e-06,0.00363541 +5.00036e-06,0.00399277 +5.01959e-06,0.004379 +5.0389e-06,0.00475493 +5.05828e-06,0.00505588 +5.07774e-06,0.00521541 +5.09727e-06,0.00524639 +5.11687e-06,0.0052111 +5.13655e-06,0.00517433 +5.15631e-06,0.00516905 +5.17614e-06,0.00519168 +5.19605e-06,0.00523691 +5.21604e-06,0.00529961 +5.2361e-06,0.00537455 +5.25624e-06,0.00545634 +5.27646e-06,0.00553915 +5.29675e-06,0.00561563 +5.31712e-06,0.00567724 +5.33757e-06,0.00571524 +5.3581e-06,0.00572621 +5.37871e-06,0.00571822 +5.3994e-06,0.00570098 +5.42017e-06,0.0056835 +5.44101e-06,0.00566895 +5.46194e-06,0.00565729 +5.48295e-06,0.00564848 +5.50404e-06,0.00564247 +5.52521e-06,0.00563921 +5.54646e-06,0.00563868 +5.56779e-06,0.00564083 +5.58921e-06,0.00564564 +5.61071e-06,0.00565308 +5.63229e-06,0.00566313 +5.65395e-06,0.00567577 +5.6757e-06,0.00569098 +5.69753e-06,0.00570875 +5.71944e-06,0.00572907 +5.74144e-06,0.00575193 +5.76352e-06,0.00577734 +5.78569e-06,0.00580528 +5.80794e-06,0.00583577 +5.83028e-06,0.00586879 +5.85271e-06,0.00590438 +5.87522e-06,0.00594252 +5.89781e-06,0.0059832 +5.9205e-06,0.00602606 +5.94327e-06,0.00607053 +5.96613e-06,0.00611601 +5.98908e-06,0.0061619 +6.01211e-06,0.00620755 +6.03524e-06,0.00625234 +6.05845e-06,0.00629559 +6.08175e-06,0.00633662 +6.10514e-06,0.00637473 +6.12863e-06,0.00640921 +6.1522e-06,0.00643935 +6.17586e-06,0.00646461 +6.19961e-06,0.00648545 +6.22346e-06,0.00650261 +6.2474e-06,0.00651688 +6.27143e-06,0.00652905 +6.29555e-06,0.00653992 +6.31976e-06,0.00655029 +6.34407e-06,0.00656098 +6.36847e-06,0.0065728 +6.39296e-06,0.00658657 +6.41755e-06,0.00660312 +6.44224e-06,0.0066233 +6.46701e-06,0.00664796 +6.49189e-06,0.00667799 +6.51686e-06,0.00671428 +6.54192e-06,0.00675777 +6.56709e-06,0.00680945 +6.59234e-06,0.00687034 +6.6177e-06,0.00694156 +6.64315e-06,0.00702427 +6.6687e-06,0.00711976 +6.69435e-06,0.00722902 +6.7201e-06,0.00735184 +6.74595e-06,0.00748777 +6.7719e-06,0.00763633 +6.79794e-06,0.00779706 +6.82409e-06,0.00796942 +6.85034e-06,0.00815285 +6.87668e-06,0.00834672 +6.90313e-06,0.00855034 +6.92968e-06,0.0087629 +6.95634e-06,0.00898354 +6.98309e-06,0.00921124 +7.00995e-06,0.00944489 +7.03691e-06,0.00968324 +7.06398e-06,0.0099249 +7.09115e-06,0.0101684 +7.11842e-06,0.0104119 +7.1458e-06,0.0106538 +7.17329e-06,0.0108949 +7.20088e-06,0.0111456 +7.22857e-06,0.0114187 +7.25638e-06,0.0117282 +7.28429e-06,0.0120899 +7.3123e-06,0.0125218 +7.34043e-06,0.0130448 +7.36866e-06,0.0136838 +7.397e-06,0.0144693 +7.42545e-06,0.015438 +7.45402e-06,0.0166168 +7.48268e-06,0.0180248 +7.51147e-06,0.0196837 +7.54036e-06,0.0216176 +7.56936e-06,0.0238518 +7.59847e-06,0.0264114 +7.6277e-06,0.0293204 +7.65704e-06,0.0325988 +7.68649e-06,0.0362608 +7.71605e-06,0.0404369 +7.74573e-06,0.0463141 +7.77552e-06,0.0560878 +7.80543e-06,0.071371 +7.83545e-06,0.0932981 +7.86559e-06,0.122536 +7.89584e-06,0.158401 +7.92621e-06,0.200197 +7.95669e-06,0.247641 +7.9873e-06,0.300155 +8.01802e-06,0.356825 +8.04886e-06,0.415435 +8.07982e-06,0.472109 +8.11089e-06,0.523781 +8.14209e-06,0.569562 +8.17341e-06,0.609529 +8.20484e-06,0.644557 +8.2364e-06,0.676228 +8.26808e-06,0.706708 +8.29988e-06,0.737856 +8.3318e-06,0.769342 +8.36385e-06,0.800464 +8.39602e-06,0.831244 +8.42831e-06,0.861969 +8.46073e-06,0.892979 +8.49327e-06,0.924673 +8.52594e-06,0.95751 +8.55873e-06,0.992017 +8.59165e-06,1.0289 +8.6247e-06,1.06911 +8.65787e-06,1.11377 +8.69117e-06,1.16418 +8.7246e-06,1.22186 +8.75816e-06,1.28862 +8.79184e-06,1.36664 +8.82566e-06,1.45823 +8.8596e-06,1.56601 +8.89368e-06,1.69318 +8.92789e-06,1.84373 +8.96223e-06,2.0189 +8.9967e-06,2.20168 +9.0313e-06,2.36429 +9.06604e-06,2.4915 +9.10091e-06,2.58153 +9.13591e-06,2.63025 +9.17105e-06,2.62654 +9.20632e-06,2.565 +9.24173e-06,2.46299 +9.27728e-06,2.34381 +9.31296e-06,2.22772 +9.34878e-06,2.12384 +9.38474e-06,2.02754 +9.42084e-06,1.93384 +9.45707e-06,1.83864 +9.49345e-06,1.73866 +9.52996e-06,1.63154 +9.56661e-06,1.5174 +9.60341e-06,1.40043 +9.64035e-06,1.28467 +9.67743e-06,1.17325 +9.71465e-06,1.06847 +9.75201e-06,0.971888 +9.78952e-06,0.884144 +9.82717e-06,0.804163 +9.86497e-06,0.730624 +9.90292e-06,0.662502 +9.941e-06,0.599015 +9.97924e-06,0.539586 +1.00176e-05,0.48382 +1.00562e-05,0.431846 +1.00948e-05,0.384109 +1.01337e-05,0.340834 +1.01726e-05,0.302049 +1.02118e-05,0.267632 +1.0251e-05,0.237361 +1.02905e-05,0.210907 +1.033e-05,0.187644 +1.03698e-05,0.166946 +1.04097e-05,0.148338 +1.04497e-05,0.13146 +1.04899e-05,0.116048 +1.05302e-05,0.10191 +1.05707e-05,0.0890079 +1.06114e-05,0.0776218 +1.06522e-05,0.0679267 +1.06932e-05,0.0599456 +1.07343e-05,0.0536155 +1.07756e-05,0.0488427 +1.0817e-05,0.0455448 +1.08587e-05,0.0435994 +1.09004e-05,0.0427617 +1.09423e-05,0.0428521 +1.09844e-05,0.0437563 +1.10267e-05,0.0454015 +1.10691e-05,0.0477389 +1.11117e-05,0.050729 +1.11544e-05,0.0543309 +1.11973e-05,0.0584965 +1.12404e-05,0.063154 +1.12836e-05,0.0681959 +1.1327e-05,0.0734679 +1.13706e-05,0.0787631 +1.14143e-05,0.083962 +1.14582e-05,0.0892296 +1.15023e-05,0.0948216 +1.15465e-05,0.101061 +1.15909e-05,0.108354 +1.16355e-05,0.117218 +1.16803e-05,0.128112 +1.17252e-05,0.14096 +1.17703e-05,0.155465 +1.18156e-05,0.171302 +1.1861e-05,0.18836 +1.19066e-05,0.206507 +1.19524e-05,0.225535 +1.19984e-05,0.245159 +1.20445e-05,0.265005 +1.20909e-05,0.284542 +1.21374e-05,0.302712 +1.21841e-05,0.318125 +1.22309e-05,0.329376 +1.2278e-05,0.336472 +1.23252e-05,0.340559 +1.23726e-05,0.342938 +1.24202e-05,0.344028 +1.2468e-05,0.343246 +1.25159e-05,0.339997 +1.2564e-05,0.334064 +1.26124e-05,0.325787 +1.26609e-05,0.315576 +1.27096e-05,0.30385 +1.27585e-05,0.291016 +1.28075e-05,0.277457 +1.28568e-05,0.26353 +1.29062e-05,0.249624 +1.29559e-05,0.236104 +1.30057e-05,0.223268 +1.30557e-05,0.211469 +1.3106e-05,0.201171 +1.31564e-05,0.192762 +1.3207e-05,0.186434 +1.32578e-05,0.18171 +1.33088e-05,0.177986 +1.33599e-05,0.174739 +1.34113e-05,0.171743 +1.34629e-05,0.168993 +1.35147e-05,0.16649 +1.35667e-05,0.164233 +1.36189e-05,0.162225 +1.36712e-05,0.160467 +1.37238e-05,0.158961 +1.37766e-05,0.15771 +1.38296e-05,0.156718 +1.38828e-05,0.155974 +1.39362e-05,0.155459 +1.39898e-05,0.155161 +1.40436e-05,0.155062 +1.40976e-05,0.15515 +1.41518e-05,0.155413 +1.42063e-05,0.155836 +1.42609e-05,0.156409 +1.43158e-05,0.157118 +1.43708e-05,0.157954 +1.44261e-05,0.158903 +1.44816e-05,0.159955 +1.45373e-05,0.161096 +1.45932e-05,0.162316 +1.46493e-05,0.163602 +1.47057e-05,0.16494 +1.47622e-05,0.166319 +1.4819e-05,0.167723 +1.4876e-05,0.169144 +1.49332e-05,0.170583 +1.49907e-05,0.172042 +1.50483e-05,0.173524 +1.51062e-05,0.175031 +1.51643e-05,0.176568 +1.52226e-05,0.178137 +1.52812e-05,0.179741 +1.53399e-05,0.181383 +1.53989e-05,0.183067 +1.54582e-05,0.184797 +1.55176e-05,0.186576 +1.55773e-05,0.188408 +1.56372e-05,0.190297 +1.56974e-05,0.192247 +1.57578e-05,0.194262 +1.58184e-05,0.196347 +1.58792e-05,0.198506 +1.59403e-05,0.200744 +1.60016e-05,0.203067 +1.60631e-05,0.205479 +1.61249e-05,0.207978 +1.61869e-05,0.210565 +1.62492e-05,0.213237 +1.63117e-05,0.215994 +1.63744e-05,0.218834 +1.64374e-05,0.221755 +1.65006e-05,0.224755 +1.65641e-05,0.227834 +1.66278e-05,0.230988 +1.66918e-05,0.234215 +1.6756e-05,0.237514 +1.68204e-05,0.240882 +1.68851e-05,0.244316 +1.69501e-05,0.247812 +1.70152e-05,0.251369 +1.70807e-05,0.254982 +1.71464e-05,0.258649 +1.72123e-05,0.262364 +1.72785e-05,0.266124 +1.7345e-05,0.269925 +1.74117e-05,0.273762 +1.74787e-05,0.277631 +1.75459e-05,0.281525 +1.76134e-05,0.285439 +1.76811e-05,0.289369 +1.77491e-05,0.293307 +1.78174e-05,0.297247 +1.78859e-05,0.301183 +1.79547e-05,0.30512 +1.80238e-05,0.30908 +1.80931e-05,0.313089 +1.81627e-05,0.317172 +1.82326e-05,0.321357 +1.83027e-05,0.325672 +1.83731e-05,0.330147 +1.84438e-05,0.334815 +1.85147e-05,0.339709 +1.85859e-05,0.344866 +1.86574e-05,0.350321 +1.87292e-05,0.356117 +1.88012e-05,0.362296 +1.88735e-05,0.368905 +1.89461e-05,0.375992 +1.9019e-05,0.383632 +1.90921e-05,0.392136 +1.91656e-05,0.401996 +1.92393e-05,0.413762 +1.93133e-05,0.428058 +1.93876e-05,0.445614 +1.94621e-05,0.467304 +1.9537e-05,0.494201 +1.96121e-05,0.527658 +1.96876e-05,0.569035 +1.97633e-05,0.618589 +1.98393e-05,0.67614 +1.99156e-05,0.741208 +1.99922e-05,0.812845 +2.00691e-05,0.889735 +2.01463e-05,0.971396 +2.02238e-05,1.05765 +2.03016e-05,1.14822 +2.03796e-05,1.24271 +2.0458e-05,1.34066 +2.05367e-05,1.44218 +2.06157e-05,1.54798 +2.0695e-05,1.65904 +2.07746e-05,1.77662 +2.08545e-05,1.90223 +2.09347e-05,2.03335 +2.10152e-05,2.15884 +2.10961e-05,2.26473 +2.11772e-05,2.34386 +2.12587e-05,2.39692 +2.13404e-05,2.42589 +2.14225e-05,2.43235 +2.15049e-05,2.4178 +2.15876e-05,2.38428 +2.16706e-05,2.33419 +2.1754e-05,2.27017 +2.18377e-05,2.19497 +2.19217e-05,2.11128 +2.2006e-05,2.02168 +2.20906e-05,1.92856 +2.21756e-05,1.83402 +2.22609e-05,1.73993 +2.23465e-05,1.64784 +2.24324e-05,1.55904 +2.25187e-05,1.47455 +2.26053e-05,1.39516 +2.26923e-05,1.32146 +2.27796e-05,1.25381 +2.28672e-05,1.19191 +2.29551e-05,1.13495 +2.30434e-05,1.08225 +2.31321e-05,1.03321 +2.3221e-05,0.987302 +2.33103e-05,0.944064 +2.34e-05,0.9031 +2.349e-05,0.864061 +2.35803e-05,0.826644 +2.3671e-05,0.790586 +2.37621e-05,0.755664 +2.38535e-05,0.72169 +2.39452e-05,0.688628 +2.40373e-05,0.656572 +2.41298e-05,0.625599 +2.42226e-05,0.595774 +2.43158e-05,0.567142 +2.44093e-05,0.539737 +2.45032e-05,0.513577 +2.45974e-05,0.488671 +2.4692e-05,0.465018 +2.4787e-05,0.442608 +2.48823e-05,0.421424 +2.4978e-05,0.401443 +2.50741e-05,0.382634 +2.51705e-05,0.364946 +2.52674e-05,0.348317 +2.53645e-05,0.33269 +2.54621e-05,0.318011 +2.556e-05,0.304228 +2.56583e-05,0.291291 +2.5757e-05,0.279154 +2.58561e-05,0.267773 +2.59556e-05,0.257107 +2.60554e-05,0.247117 +2.61556e-05,0.237765 +2.62562e-05,0.229018 +2.63572e-05,0.220844 +2.64586e-05,0.213212 +2.65603e-05,0.206095 +2.66625e-05,0.199465 +2.6765e-05,0.193296 +2.6868e-05,0.187551 +2.69713e-05,0.18219 +2.70751e-05,0.177179 +2.71792e-05,0.172487 +2.72837e-05,0.168084 +2.73887e-05,0.163944 +2.7494e-05,0.160043 +2.75998e-05,0.156358 +2.77059e-05,0.152869 +2.78125e-05,0.149557 +2.79195e-05,0.146405 +2.80268e-05,0.143396 +2.81346e-05,0.140514 +2.82429e-05,0.137747 +2.83515e-05,0.13508 +2.84605e-05,0.132502 +2.857e-05,0.130001 +2.857e-05,0.13 +3.077e-05,0.092 +3.333e-05,0.067 +3.636e-05,0.056 +4e-05,0.046 +0.0001,0.0159 +0.000125,0.0119 +0.0001667,0.00862 +0.00025,0.00696 +0.0005,0.00796 diff --git a/src/scgenerator/data/materials.toml b/src/scgenerator/data/materials.toml index 08478db..a8a6bd0 100644 --- a/src/scgenerator/data/materials.toml +++ b/src/scgenerator/data/materials.toml @@ -1,271 +1,262 @@ [silica.sellmeier] -B = [0.6961663, 0.4079426, 0.8974794] -C = [4.67914826e-15, 1.35120631e-14, 9.79340025e-11] -kind = 1 -reference = [ - "I. H. Malitson. Interspecimen comparison of the refractive index of fused silica, J. Opt. Soc. Am. 55, 1205-1208 (1965)", - "C. Z. Tan. Determination of refractive index of silica glass for infrared wavelengths by IR spectroscopy, J. Non-Cryst. Solids 223, 158-163 (1998)", -] + B = [0.6961663, 0.4079426, 0.8974794] + C = [4.67914826e-15, 1.35120631e-14, 9.79340025e-11] + kind = 1 + reference = [ + "I. H. Malitson. Interspecimen comparison of the refractive index of fused silica, J. Opt. Soc. Am. 55, 1205-1208 (1965)", + "C. Z. Tan. Determination of refractive index of silica glass for infrared wavelengths by IR spectroscopy, J. Non-Cryst. Solids 223, 158-163 (1998)", + ] [nbk7.sellmeier] -B = [1.03961212, 0.231792344, 1.01046945] -C = [6.00069867e-15, 2.00179144e-14, 1.03560653e-10] -kind = 1 -reference = [ - "SCHOTT Zemax catalog 2017-01-20b (obtained from http://www.schott.com)", -] + B = [1.03961212, 0.231792344, 1.01046945] + C = [6.00069867e-15, 2.00179144e-14, 1.03560653e-10] + kind = 1 + reference = ["SCHOTT Zemax catalog 2017-01-20b (obtained from http://www.schott.com)"] [e7.sellmeier] -B = [0.0, 7.2e-15, 3e-28] -C = [] -kind = 3 -reference = [ - "J. Li, C. H. Wen, S. Gauza, R. Lu and S. T. Wu. Refractive indices of liquid crystals for display applications, J. Disp. Technol. 1, 51-61, 2005", -] + B = [0.0, 7.2e-15, 3e-28] + C = [] + kind = 3 + reference = [ + "J. Li, C. H. Wen, S. Gauza, R. Lu and S. T. Wu. Refractive indices of liquid crystals for display applications, J. Disp. Technol. 1, 51-61, 2005", + ] [d-zlaf52la.sellmeier] -B = [ - -8.48750283e9, - 3.1753525, - 3.71301651e-14, - -1.09609062e-27, - 2.4418582e-40, - -8.94583294e-54, -] -C = [] -kind = 3 -reference = ["CDGM Zemax catalog 2017-09 (obtained from http://www.cdgmgd.com)"] + B = [-8.48750283e9, 3.1753525, 3.71301651e-14, -1.09609062e-27, 2.4418582e-40, -8.94583294e-54] + C = [] + kind = 3 + reference = ["CDGM Zemax catalog 2017-09 (obtained from http://www.cdgmgd.com)"] -# ----- -# GASES -# ----- + # ----- + # GASES + # ----- [air] -a = 0.1358 -b = 3.64e-5 -atomic_mass = 4.809e-26 + a = 0.1358 + atomic_mass = 4.809e-26 + b = 3.64e-5 [air.sellmeier] -B = [57921050000.0, 1679170000.0] -C = [238018500000000.0, 57362000000000.0] -P0 = 101325 -T0 = 288.15 -kind = 2 + B = [57921050000.0, 1679170000.0] + C = [238018500000000.0, 57362000000000.0] + P0 = 101325 + T0 = 288.15 + kind = 2 [air.kerr] -P0 = 101325 -T0 = 293.15 -n2 = 4.0e-23 -source = "Pigeon, J. J., Tochitsky, S. Y., Welch, E. C., & Joshi, C. (2016). Measurements of the nonlinear refractive index of air, N 2, and O 2 at 10 μm using four-wave mixing. Optics letters, 41(17), 3924-3927." + P0 = 101325 + T0 = 293.15 + n2 = 4.0e-23 + source = "Pigeon, J. J., Tochitsky, S. Y., Welch, E. C., & Joshi, C. (2016). Measurements of the nonlinear refractive index of air, N 2, and O 2 at 10 μm using four-wave mixing. Optics letters, 41(17), 3924-3927." [nitrogen] -a = 0.137 -b = 1.709e-5 -atomic_mass = 2.3259e-26 + a = 0.137 + atomic_mass = 2.3259e-26 + b = 1.709e-5 [nitrogen.sellmeier] -B = [32431570000.0] -C = [144000000000000.0] -P0 = 101325 -T0 = 273.15 -const = 6.8552e-5 -kind = 2 + B = [32431570000.0] + C = [144000000000000.0] + P0 = 101325 + T0 = 273.15 + const = 6.8552e-5 + kind = 2 [nitrogen.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 2.2e-23 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 2.2e-23 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [helium] -a = 0.00346 -b = 2.38e-5 -ionization_energy = 3.9393356074281e-18 -atomic_number = 2 -atomic_mass = 6.646477e-27 + a = 0.00346 + atomic_mass = 6.646477e-27 + atomic_number = 2 + b = 2.38e-5 + ionization_energy = 3.9393356074281e-18 [helium.sellmeier] -B = [2.16463842e-5, 2.10561127e-7, 4.7509272e-5] -C = [-6.80769781e-16, 5.13251289e-15, 3.18621354e-15] -P0 = 101325 -T0 = 273.15 -kind = 1 -source = "A. Ermolov, K. F. Mak, M. H. Frosz, J. C. Travers, P. St. J. Russell, Supercontinuum generation in the vacuum ultraviolet through dispersive-wave and soliton-plasma interaction in a noble-gas-filled hollow-core photonic crystal fiber, Phys. Rev. A 92, 033821 (2015)" + B = [2.16463842e-5, 2.10561127e-7, 4.7509272e-5] + C = [-6.80769781e-16, 5.13251289e-15, 3.18621354e-15] + P0 = 101325 + T0 = 273.15 + kind = 1 + source = "A. Ermolov, K. F. Mak, M. H. Frosz, J. C. Travers, P. St. J. Russell, Supercontinuum generation in the vacuum ultraviolet through dispersive-wave and soliton-plasma interaction in a noble-gas-filled hollow-core photonic crystal fiber, Phys. Rev. A 92, 033821 (2015)" [helium.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 3.1e-25 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 3.1e-25 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [helium_alt] -a = 0.00346 -b = 2.38e-5 + a = 0.00346 + b = 2.38e-5 [helium_alt.sellmeier] -B = [14755297000.0] -C = [426297400000000.0] -P0 = 101325 -T0 = 273.15 -kind = 2 -source = " C. Cuthbertson and M. Cuthbertson. The refraction and dispersion of neon and helium. Proc. R. Soc. London A 135, 40-47 (1936)" + B = [14755297000.0] + C = [426297400000000.0] + P0 = 101325 + T0 = 273.15 + kind = 2 + source = " C. Cuthbertson and M. Cuthbertson. The refraction and dispersion of neon and helium. Proc. R. Soc. London A 135, 40-47 (1936)" [helium_alt.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 3.1e-25 + P0 = 30400.0 + T0 = 273.15 + n2 = 3.1e-25 [hydrogen] -a = 0.02453 -b = 2.651e-5 -atomic_mass = 1.674e-27 + a = 0.02453 + atomic_mass = 1.674e-27 + b = 2.651e-5 [hydrogen.sellmeier] -B = [0.0148956, 0.0049037] -C = [1.807e-10, 9.2e-11] -P0 = 101325 -T0 = 273.15 -kind = 2 -source = "E. R. Peck and S. Hung. Refractivity and dispersion of hydrogen in the visible and near infrared, J. Opt. Soc. Am. 67, 1550-1554 (1977)" + B = [0.0148956, 0.0049037] + C = [1.807e-10, 9.2e-11] + P0 = 101325 + T0 = 273.15 + kind = 2 + source = "E. R. Peck and S. Hung. Refractivity and dispersion of hydrogen in the visible and near infrared, J. Opt. Soc. Am. 67, 1550-1554 (1977)" [hydrogen.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 6.36e-24 -source = "Shelton, D. P., & Rice, J. E. (1994). Measurements and calculations of the hyperpolarizabilities of atoms and small molecules in the gas phase. Chemical Reviews, 94(1), 3-29" + P0 = 30400.0 + T0 = 273.15 + n2 = 6.36e-24 + source = "Shelton, D. P., & Rice, J. E. (1994). Measurements and calculations of the hyperpolarizabilities of atoms and small molecules in the gas phase. Chemical Reviews, 94(1), 3-29" [neon] -a = 0.02135 -b = 1.709e-5 -ionization_energy = 3.45501365359425e-18 -atomic_number = 10 -atomic_mass = 3.35092e-26 + a = 0.02135 + atomic_mass = 3.35092e-26 + atomic_number = 10 + b = 1.709e-5 + ionization_energy = 3.45501365359425e-18 [neon.sellmeier] -B = [1281450000.0, 22048600000.0] -C = [184661000000000.0, 376840000000000.0] -P0 = 101325 -T0 = 273.15 -kind = 2 + B = [1281450000.0, 22048600000.0] + C = [184661000000000.0, 376840000000000.0] + P0 = 101325 + T0 = 273.15 + kind = 2 [neon.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 8.7e-25 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 8.7e-25 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [argon] -a = 0.1355 -b = 3.201e-5 -ionization_energy = 2.5249661793774e-18 -atomic_number = 18 -atomic_mass = 6.6338e-26 + a = 0.1355 + atomic_mass = 6.6338e-26 + atomic_number = 18 + b = 3.201e-5 + ionization_energy = 2.5249661793774e-18 [argon.sellmeier] -B = [2501410000.0, 500283000.0, 52234300000.0] -C = [91012000000000.0, 87892000000000.0, 214020000000000.0] -P0 = 101325 -T0 = 273.15 -kind = 2 -source = "A. Bideau-Mehu, Y. Guern, R. Abjean, A. Johannin-Gilles. Measurement of refractive indices of neon, argon, krypton and xenon in the 253.7-140.4 nm wavelength range. Dispersion relations and estimated oscillator strengths of the resonance lines. J. Quant. Spectrosc. Rad. Transfer 25, 395-402 (1981)" + B = [2501410000.0, 500283000.0, 52234300000.0] + C = [91012000000000.0, 87892000000000.0, 214020000000000.0] + P0 = 101325 + T0 = 273.15 + kind = 2 + source = "A. Bideau-Mehu, Y. Guern, R. Abjean, A. Johannin-Gilles. Measurement of refractive indices of neon, argon, krypton and xenon in the 253.7-140.4 nm wavelength range. Dispersion relations and estimated oscillator strengths of the resonance lines. J. Quant. Spectrosc. Rad. Transfer 25, 395-402 (1981)" [argon.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 9.7e-24 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 9.7e-24 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [argon_alt] -a = 0.1355 -b = 3.201e-5 + a = 0.1355 + b = 3.201e-5 [argon_alt.sellmeier] -B = [0.0002033229, 0.0003445831] -C = [2.0612e-16, 8.066e-15] -P0 = 101325 -T0 = 273.15 -kind = 1 -source = "A. Börzsönyi, Z. Heiner, M. P. Kalashnikov, A. P. Kovács, and K. Osvay, Dispersion measurement of inert gases and gas mixtures at 800 nm, Appl. Opt. 47, 4856-4863 (2008)" + B = [0.0002033229, 0.0003445831] + C = [2.0612e-16, 8.066e-15] + P0 = 101325 + T0 = 273.15 + kind = 1 + source = "A. Börzsönyi, Z. Heiner, M. P. Kalashnikov, A. P. Kovács, and K. Osvay, Dispersion measurement of inert gases and gas mixtures at 800 nm, Appl. Opt. 47, 4856-4863 (2008)" [argon_alt.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 9.7e-24 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 9.7e-24 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [argon_alt2] -a = 0.1355 -b = 3.201e-5 + a = 0.1355 + b = 3.201e-5 [argon_alt2.sellmeier] -B = [0.030182943] -C = [144000000000000.0] -P0 = 101325 -T0 = 273.15 -const = 6.7867e-5 -kind = 2 -source = "E. R. Peck and D. J. Fisher. Dispersion of argon, J. Opt. Soc. Am. 54, 1362-1364 (1964)" + B = [0.030182943] + C = [144000000000000.0] + P0 = 101325 + T0 = 273.15 + const = 6.7867e-5 + kind = 2 + source = "E. R. Peck and D. J. Fisher. Dispersion of argon, J. Opt. Soc. Am. 54, 1362-1364 (1964)" [argon_alt2.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 9.7e-24 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 9.7e-24 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [krypton] -a = 0.2349 -b = 3.978e-5 -ionization_energy = 2.2429831039374e-18 -atomic_number = 36 -atomic_mass = 1.3915e-25 + a = 0.2349 + atomic_mass = 1.3915e-25 + atomic_number = 36 + b = 3.978e-5 + ionization_energy = 2.2429831039374e-18 [krypton.sellmeier] -B = [2536370000.0, 2736490000.0, 62080200000.0] -C = [65474200000000.0, 73698000000000.0, 181080000000000.0] -P0 = 101325 -T0 = 273.15 -kind = 2 + B = [2536370000.0, 2736490000.0, 62080200000.0] + C = [65474200000000.0, 73698000000000.0, 181080000000000.0] + P0 = 101325 + T0 = 273.15 + kind = 2 [krypton.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 2.2e-23 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 2.2e-23 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [xenon] -a = 0.425 -b = 5.105e-5 -ionization_energy = 1.94342415157935 -atomic_number = 54 -atomic_mass = 2.18017e-25 + a = 0.425 + atomic_mass = 2.18017e-25 + atomic_number = 54 + b = 5.105e-5 + ionization_energy = 1.94342415157935 [xenon.sellmeier] -B = [3228690000.0, 3553930000.0, 60676400000.0] -C = [46301000000000.0, 59578000000000.0, 112740000000000.0] -P0 = 101325 -T0 = 273.15 -kind = 2 + B = [3228690000.0, 3553930000.0, 60676400000.0] + C = [46301000000000.0, 59578000000000.0, 112740000000000.0] + P0 = 101325 + T0 = 273.15 + kind = 2 [xenon.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 5.8e-23 -source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." + P0 = 30400.0 + T0 = 273.15 + n2 = 5.8e-23 + source = "Wahlstrand, J. K., Cheng, Y. H., & Milchberg, H. M. (2012). High field optical nonlinearity and the Kramers-Kronig relations. Physical review letters, 109(11), 113904." [vacuum] -a = 0 -b = 0 -atomic_mass = 0 + a = 0 + atomic_mass = 0 + b = 0 [vacuum.sellmeier] -B = [] -C = [] -P0 = 101325 -T0 = 273.15 -kind = 1 + B = [] + C = [] + P0 = 101325 + T0 = 273.15 + kind = 1 [vacuum.kerr] -P0 = 30400.0 -T0 = 273.15 -n2 = 0 -source = "none" + P0 = 30400.0 + T0 = 273.15 + n2 = 0 + source = "none" diff --git a/src/scgenerator/data/n_silica.txt b/src/scgenerator/data/n_silica.txt new file mode 100644 index 0000000..89cf98b --- /dev/null +++ b/src/scgenerator/data/n_silica.txt @@ -0,0 +1,1198 @@ +6.199e-10,0.99993 +6.351e-10,0.99992 +6.508e-10,0.99992 +6.665e-10,0.99991 +6.831e-10,0.99991 +6.997e-10,0.99991 +7.166e-10,0.9999 +7.345e-10,0.9999 +7.523e-10,0.99989 +7.705e-10,0.99989 +7.897e-10,0.99988 +8.087e-10,0.99987 +8.287e-10,0.99987 +8.492e-10,0.99986 +8.694e-10,0.99985 +8.906e-10,0.99985 +9.129e-10,0.99984 +9.35e-10,0.99983 +9.581e-10,0.99982 +9.816e-10,0.99981 +1.005e-09,0.99981 +1.03e-09,0.9998 +1.055e-09,0.99979 +1.081e-09,0.99978 +1.108e-09,0.99976 +1.134e-09,0.99975 +1.162e-09,0.99974 +1.191e-09,0.99973 +1.22e-09,0.99971 +1.25e-09,0.9997 +1.281e-09,0.99968 +1.312e-09,0.99967 +1.343e-09,0.99965 +1.376e-09,0.99964 +1.41e-09,0.99962 +1.445e-09,0.9996 +1.479e-09,0.99958 +1.516e-09,0.99956 +1.554e-09,0.99954 +1.592e-09,0.99951 +1.631e-09,0.99949 +1.671e-09,0.99946 +1.71e-09,0.99944 +1.754e-09,0.99941 +1.797e-09,0.99938 +1.839e-09,0.99935 +1.884e-09,0.99932 +1.931e-09,0.99928 +1.977e-09,0.99925 +2.026e-09,0.99921 +2.077e-09,0.99917 +2.127e-09,0.99913 +2.179e-09,0.99909 +2.234e-09,0.99904 +2.287e-09,0.99899 +2.344e-09,0.99894 +2.403e-09,0.99889 +2.46e-09,0.99884 +2.52e-09,0.99878 +2.583e-09,0.99872 +2.643e-09,0.99866 +2.707e-09,0.99859 +2.774e-09,0.99852 +2.844e-09,0.99845 +2.91e-09,0.99837 +2.987e-09,0.99828 +3.054e-09,0.99821 +3.131e-09,0.99812 +3.212e-09,0.99802 +3.289e-09,0.99792 +3.369e-09,0.99782 +3.453e-09,0.99771 +3.532e-09,0.9976 +3.625e-09,0.99747 +3.712e-09,0.99735 +3.803e-09,0.99722 +3.899e-09,0.99708 +3.986e-09,0.99695 +4.092e-09,0.99678 +4.188e-09,0.99663 +4.29e-09,0.99646 +4.396e-09,0.99629 +4.508e-09,0.99609 +4.609e-09,0.99592 +4.732e-09,0.9957 +4.843e-09,0.99549 +4.959e-09,0.99527 +5.081e-09,0.99504 +5.209e-09,0.99478 +5.344e-09,0.99451 +5.462e-09,0.99427 +5.61e-09,0.99395 +5.74e-09,0.99367 +5.876e-09,0.99336 +6.018e-09,0.99304 +6.168e-09,0.99269 +6.325e-09,0.99231 +6.491e-09,0.9919 +6.63e-09,0.99155 +6.812e-09,0.99108 +6.965e-09,0.99068 +7.125e-09,0.99024 +7.293e-09,0.9898 +7.469e-09,0.9893 +7.653e-09,0.9887 +7.847e-09,0.9882 +8.051e-09,0.9875 +8.2e-09,0.9871 +8.4e-09,0.9865 +8.6e-09,0.986 +8.8e-09,0.9855 +9e-09,0.9853 +9.2e-09,0.9854 +9.4e-09,0.9858 +9.6e-09,0.9868 +9.8e-09,0.9872 +1e-08,0.9874 +1.02e-08,0.9865 +1.04e-08,0.9851 +1.06e-08,0.9841 +1.08e-08,0.9848 +1.1e-08,0.9844 +1.12e-08,0.9822 +1.14e-08,0.9828 +1.16e-08,0.9867 +1.18e-08,0.9858 +1.2e-08,0.9839 +1.22e-08,0.9823 +1.24e-08,0.9813 +1.26e-08,0.9803 +1.28e-08,0.9794 +1.3e-08,0.9789 +1.32e-08,0.9778 +1.34e-08,0.977 +1.36e-08,0.9761 +1.38e-08,0.9747 +1.5e-08,0.9634 +1.6e-08,0.9608 +1.7e-08,0.9562 +1.8e-08,0.9509 +1.9e-08,0.9458 +2e-08,0.9416 +2.1e-08,0.9386 +2.2e-08,0.9328 +2.3e-08,0.9271 +2.4e-08,0.9222 +2.5e-08,0.9164 +2.6e-08,0.9105 +2.7e-08,0.9207 +2.8e-08,0.9175 +2.9e-08,0.9137 +3e-08,0.913 +3.1e-08,0.907 +3.2e-08,0.901 +3.3e-08,0.895 +3.4e-08,0.888 +3.5e-08,0.882 +3.6e-08,0.877 +3.7e-08,0.87 +3.8e-08,0.866 +3.9e-08,0.858 +4e-08,0.851 +4.1e-08,0.845 +4.2e-08,0.839 +4.3e-08,0.833 +4.4e-08,0.827 +4.5e-08,0.822 +4.6e-08,0.817 +4.7e-08,0.813 +4.8e-08,0.808 +4.9e-08,0.804 +5e-08,0.803 +5.1e-08,0.804 +5.2e-08,0.806 +5.3e-08,0.810271 +5.32536e-08,0.81148 +5.35084e-08,0.812808 +5.37645e-08,0.814244 +5.40217e-08,0.81578 +5.42802e-08,0.817407 +5.45399e-08,0.819114 +5.48009e-08,0.820891 +5.50631e-08,0.822728 +5.53266e-08,0.824614 +5.55914e-08,0.82654 +5.58574e-08,0.828492 +5.61246e-08,0.830457 +5.63932e-08,0.832425 +5.6663e-08,0.834382 +5.69342e-08,0.836317 +5.72066e-08,0.838218 +5.74803e-08,0.840085 +5.77554e-08,0.841922 +5.80317e-08,0.843733 +5.83094e-08,0.845528 +5.85884e-08,0.847313 +5.88688e-08,0.849094 +5.91505e-08,0.850878 +5.94335e-08,0.852671 +5.97179e-08,0.854481 +6.00036e-08,0.856317 +6.02907e-08,0.858184 +6.05792e-08,0.86009 +6.08691e-08,0.862042 +6.11604e-08,0.864052 +6.1453e-08,0.86613 +6.17471e-08,0.86828 +6.20425e-08,0.87051 +6.23394e-08,0.872828 +6.26377e-08,0.875255 +6.29374e-08,0.877792 +6.32386e-08,0.880447 +6.35412e-08,0.883228 +6.38452e-08,0.886158 +6.41507e-08,0.889237 +6.44577e-08,0.892465 +6.47661e-08,0.895847 +6.5076e-08,0.899412 +6.53874e-08,0.903145 +6.57003e-08,0.90705 +6.60146e-08,0.911146 +6.63305e-08,0.915445 +6.66479e-08,0.919934 +6.69668e-08,0.924628 +6.72873e-08,0.929555 +6.76092e-08,0.934689 +6.79327e-08,0.940042 +6.82578e-08,0.94564 +6.85844e-08,0.951457 +6.89126e-08,0.957507 +6.92423e-08,0.96381 +6.95737e-08,0.970341 +6.99066e-08,0.977131 +7.02411e-08,0.984246 +7.05772e-08,0.991841 +7.09149e-08,1.00029 +7.12542e-08,1.00968 +7.15952e-08,1.02036 +7.19377e-08,1.0326 +7.2282e-08,1.04626 +7.26278e-08,1.06105 +7.29754e-08,1.0762 +7.33245e-08,1.09106 +7.36754e-08,1.10505 +7.40279e-08,1.11742 +7.43822e-08,1.1277 +7.47381e-08,1.13607 +7.50957e-08,1.14274 +7.5455e-08,1.14817 +7.58161e-08,1.15258 +7.61789e-08,1.15643 +7.65434e-08,1.15998 +7.69096e-08,1.16336 +7.72777e-08,1.16646 +7.76474e-08,1.16923 +7.8019e-08,1.17156 +7.83923e-08,1.17335 +7.87674e-08,1.17454 +7.91443e-08,1.17503 +7.9523e-08,1.17498 +7.99035e-08,1.17443 +8.02859e-08,1.17355 +8.067e-08,1.17245 +8.1056e-08,1.17126 +8.14439e-08,1.17011 +8.18336e-08,1.16914 +8.22252e-08,1.16847 +8.26186e-08,1.16829 +8.30139e-08,1.1687 +8.34112e-08,1.16983 +8.38103e-08,1.17196 +8.42113e-08,1.17511 +8.46143e-08,1.17937 +8.50191e-08,1.18481 +8.5426e-08,1.19127 +8.58347e-08,1.1987 +8.62454e-08,1.20715 +8.66581e-08,1.21648 +8.70728e-08,1.22663 +8.74894e-08,1.23759 +8.79081e-08,1.24928 +8.83287e-08,1.26162 +8.87514e-08,1.27455 +8.9176e-08,1.28791 +8.96027e-08,1.30135 +9.00315e-08,1.31456 +9.04623e-08,1.3272 +9.08951e-08,1.3389 +9.13301e-08,1.34927 +9.17671e-08,1.35797 +9.22062e-08,1.36472 +9.26474e-08,1.36956 +9.30907e-08,1.3727 +9.35362e-08,1.37437 +9.39837e-08,1.37487 +9.44334e-08,1.37447 +9.48853e-08,1.37343 +9.53393e-08,1.37203 +9.57955e-08,1.37054 +9.62539e-08,1.36918 +9.67145e-08,1.36815 +9.71773e-08,1.3677 +9.76423e-08,1.36801 +9.81095e-08,1.36933 +9.85789e-08,1.37188 +9.90506e-08,1.3759 +9.95246e-08,1.38161 +1.00001e-07,1.38906 +1.00479e-07,1.39817 +1.0096e-07,1.40888 +1.01443e-07,1.42113 +1.01929e-07,1.43482 +1.02416e-07,1.44988 +1.02906e-07,1.46623 +1.03399e-07,1.48379 +1.03894e-07,1.50246 +1.04391e-07,1.52216 +1.0489e-07,1.5428 +1.05392e-07,1.56428 +1.05896e-07,1.58649 +1.06403e-07,1.60935 +1.06912e-07,1.63275 +1.07424e-07,1.65645 +1.07938e-07,1.67981 +1.08454e-07,1.70207 +1.08973e-07,1.72245 +1.09495e-07,1.74016 +1.10019e-07,1.75437 +1.10545e-07,1.76418 +1.11074e-07,1.76864 +1.11606e-07,1.76697 +1.1214e-07,1.75827 +1.12676e-07,1.74162 +1.13215e-07,1.71586 +1.13757e-07,1.67996 +1.14301e-07,1.63501 +1.14848e-07,1.58664 +1.15398e-07,1.54184 +1.1595e-07,1.50807 +1.16505e-07,1.49116 +1.17062e-07,1.49425 +1.17622e-07,1.52134 +1.18185e-07,1.57393 +1.18751e-07,1.66289 +1.19319e-07,1.79851 +1.1989e-07,1.93716 +1.20464e-07,2.06021 +1.2104e-07,2.16289 +1.21619e-07,2.25357 +1.22201e-07,2.30964 +1.22786e-07,2.3342 +1.23373e-07,2.33671 +1.23964e-07,2.32575 +1.24557e-07,2.31009 +1.25153e-07,2.29183 +1.25752e-07,2.27153 +1.26354e-07,2.24954 +1.26958e-07,2.22635 +1.27566e-07,2.2024 +1.28176e-07,2.17817 +1.28789e-07,2.15415 +1.29406e-07,2.1308 +1.30025e-07,2.10842 +1.30647e-07,2.08693 +1.31272e-07,2.06637 +1.319e-07,2.04669 +1.32531e-07,2.02784 +1.33166e-07,2.00985 +1.33803e-07,1.99263 +1.34443e-07,1.9762 +1.35086e-07,1.9605 +1.35733e-07,1.94551 +1.36382e-07,1.93121 +1.37035e-07,1.91755 +1.3769e-07,1.9045 +1.38349e-07,1.89203 +1.39011e-07,1.8801 +1.39676e-07,1.86867 +1.40345e-07,1.85771 +1.41016e-07,1.84716 +1.41691e-07,1.83702 +1.42369e-07,1.82725 +1.4305e-07,1.81785 +1.43735e-07,1.8088 +1.44423e-07,1.8001 +1.45114e-07,1.79173 +1.45808e-07,1.78369 +1.46506e-07,1.77596 +1.47207e-07,1.76852 +1.47911e-07,1.76136 +1.48619e-07,1.75447 +1.4933e-07,1.74784 +1.50045e-07,1.74144 +1.50763e-07,1.73527 +1.51484e-07,1.72931 +1.52209e-07,1.72354 +1.52937e-07,1.71794 +1.53669e-07,1.7125 +1.54404e-07,1.7072 +1.55143e-07,1.70202 +1.55885e-07,1.69695 +1.56631e-07,1.69198 +1.57381e-07,1.68711 +1.58134e-07,1.68235 +1.58891e-07,1.67769 +1.59651e-07,1.67313 +1.60415e-07,1.66867 +1.61182e-07,1.66432 +1.61954e-07,1.66007 +1.62729e-07,1.65592 +1.63507e-07,1.65186 +1.6429e-07,1.64791 +1.65076e-07,1.64405 +1.65866e-07,1.6403 +1.66659e-07,1.63663 +1.67457e-07,1.63306 +1.68258e-07,1.62959 +1.69063e-07,1.62621 +1.69872e-07,1.62292 +1.70685e-07,1.61972 +1.71502e-07,1.61661 +1.72322e-07,1.61359 +1.73147e-07,1.61066 +1.73975e-07,1.60781 +1.74808e-07,1.60505 +1.75644e-07,1.60236 +1.76485e-07,1.59976 +1.77329e-07,1.59724 +1.78178e-07,1.59479 +1.7903e-07,1.59242 +1.79887e-07,1.59012 +1.80748e-07,1.58789 +1.81613e-07,1.58574 +1.82482e-07,1.58365 +1.83355e-07,1.58162 +1.84232e-07,1.57965 +1.85114e-07,1.57775 +1.85999e-07,1.5759 +1.86889e-07,1.57411 +1.87784e-07,1.57237 +1.88682e-07,1.57068 +1.89585e-07,1.56903 +1.90492e-07,1.56743 +1.91404e-07,1.56587 +1.9232e-07,1.56435 +1.9324e-07,1.56286 +1.94165e-07,1.5614 +1.95094e-07,1.55997 +1.96027e-07,1.55856 +1.96965e-07,1.55717 +1.97908e-07,1.55579 +1.98855e-07,1.55443 +1.99806e-07,1.55308 +2.00762e-07,1.55174 +2.01723e-07,1.55041 +2.02688e-07,1.54909 +2.03658e-07,1.54778 +2.04632e-07,1.54648 +2.05612e-07,1.5452 +2.06595e-07,1.54392 +2.07584e-07,1.54266 +2.08577e-07,1.54142 +2.09575e-07,1.54019 +2.10578e-07,1.53897 +2.11586e-07,1.53776 +2.12598e-07,1.53657 +2.13615e-07,1.5354 +2.14638e-07,1.53424 +2.15665e-07,1.5331 +2.16697e-07,1.53198 +2.17734e-07,1.53087 +2.18775e-07,1.52978 +2.19822e-07,1.52871 +2.20874e-07,1.52765 +2.21931e-07,1.52662 +2.22993e-07,1.52561 +2.2406e-07,1.52461 +2.25132e-07,1.52364 +2.26209e-07,1.52269 +2.27292e-07,1.52175 +2.28379e-07,1.52085 +2.29472e-07,1.51996 +2.3057e-07,1.5191 +2.31673e-07,1.51826 +2.32782e-07,1.51745 +2.33896e-07,1.51666 +2.35015e-07,1.5159 +2.3614e-07,1.51516 +2.37269e-07,1.51446 +2.38405e-07,1.51378 +2.39546e-07,1.51313 +2.40692e-07,1.5125 +2.41843e-07,1.51191 +2.43001e-07,1.51135 +2.44163e-07,1.51082 +2.45332e-07,1.51032 +2.46506e-07,1.50985 +2.47685e-07,1.50942 +2.4887e-07,1.50808 +2.50061e-07,1.50745 +2.51258e-07,1.50681 +2.5246e-07,1.50617 +2.53668e-07,1.50553 +2.54882e-07,1.50488 +2.56101e-07,1.5043 +2.57327e-07,1.5037 +2.58558e-07,1.50311 +2.59795e-07,1.50251 +2.61039e-07,1.50191 +2.62288e-07,1.50134 +2.63543e-07,1.50079 +2.64804e-07,1.50024 +2.66071e-07,1.49968 +2.67344e-07,1.49912 +2.68623e-07,1.49858 +2.69909e-07,1.49806 +2.712e-07,1.49754 +2.72498e-07,1.49702 +2.73802e-07,1.4965 +2.75112e-07,1.49598 +2.76428e-07,1.4955 +2.77751e-07,1.49501 +2.7908e-07,1.49452 +2.80415e-07,1.49403 +2.81757e-07,1.49355 +2.83105e-07,1.49309 +2.8446e-07,1.49263 +2.85821e-07,1.49217 +2.87189e-07,1.4917 +2.88563e-07,1.49125 +2.89944e-07,1.49082 +2.91331e-07,1.49039 +2.92725e-07,1.48995 +2.94126e-07,1.48951 +2.95533e-07,1.4891 +2.96947e-07,1.48869 +2.98368e-07,1.48828 +2.99796e-07,1.48786 +3.0123e-07,1.48745 +3.02672e-07,1.48706 +3.0412e-07,1.48667 +3.05575e-07,1.48628 +3.07037e-07,1.48589 +3.08507e-07,1.48551 +3.09983e-07,1.48514 +3.11466e-07,1.48477 +3.12956e-07,1.4844 +3.14454e-07,1.48403 +3.15959e-07,1.48368 +3.1747e-07,1.48333 +3.1899e-07,1.48298 +3.20516e-07,1.48263 +3.2205e-07,1.48229 +3.23591e-07,1.48195 +3.25139e-07,1.48162 +3.26695e-07,1.48128 +3.28258e-07,1.48095 +3.29829e-07,1.48064 +3.31407e-07,1.48032 +3.32993e-07,1.48 +3.34586e-07,1.47968 +3.36187e-07,1.47938 +3.37796e-07,1.47907 +3.39412e-07,1.47877 +3.41036e-07,1.47846 +3.42668e-07,1.47817 +3.44308e-07,1.47788 +3.45955e-07,1.47759 +3.47611e-07,1.4773 +3.49274e-07,1.47702 +3.50945e-07,1.47674 +3.52624e-07,1.47646 +3.54312e-07,1.47618 +3.56007e-07,1.47592 +3.57711e-07,1.47565 +3.59422e-07,1.47538 +3.61142e-07,1.47512 +3.6287e-07,1.47486 +3.64607e-07,1.47461 +3.66351e-07,1.47435 +3.68104e-07,1.4741 +3.69866e-07,1.47385 +3.71635e-07,1.4736 +3.73414e-07,1.47336 +3.752e-07,1.47312 +3.76996e-07,1.47288 +3.788e-07,1.47264 +3.80612e-07,1.47241 +3.82433e-07,1.47218 +3.84263e-07,1.47195 +3.86102e-07,1.47172 +3.8795e-07,1.4715 +3.89806e-07,1.47128 +3.91671e-07,1.47106 +3.93545e-07,1.47084 +3.95428e-07,1.47063 +3.9732e-07,1.47042 +3.99222e-07,1.4702 +4.01132e-07,1.47 +4.03051e-07,1.46979 +4.0498e-07,1.46959 +4.06918e-07,1.46938 +4.08865e-07,1.46919 +4.10821e-07,1.46899 +4.12787e-07,1.46879 +4.14762e-07,1.4686 +4.16747e-07,1.46841 +4.18741e-07,1.46821 +4.20745e-07,1.46803 +4.22758e-07,1.46784 +4.24781e-07,1.46766 +4.26813e-07,1.46747 +4.28856e-07,1.46729 +4.30908e-07,1.46711 +4.3297e-07,1.46693 +4.35041e-07,1.46676 +4.37123e-07,1.46659 +4.39215e-07,1.46641 +4.41316e-07,1.46625 +4.43428e-07,1.46608 +4.4555e-07,1.46591 +4.47682e-07,1.46574 +4.49824e-07,1.46558 +4.51976e-07,1.46542 +4.54139e-07,1.46526 +4.56312e-07,1.4651 +4.58496e-07,1.46494 +4.60689e-07,1.46479 +4.62894e-07,1.46463 +4.65109e-07,1.46448 +4.67334e-07,1.46433 +4.69571e-07,1.46418 +4.71817e-07,1.46403 +4.74075e-07,1.46388 +4.76344e-07,1.46374 +4.78623e-07,1.46359 +4.80913e-07,1.46345 +4.83214e-07,1.46331 +4.85526e-07,1.46316 +4.8785e-07,1.46303 +4.90184e-07,1.46289 +4.92529e-07,1.46275 +4.94886e-07,1.46262 +4.97254e-07,1.46248 +4.99634e-07,1.46235 +5.02024e-07,1.46222 +5.04427e-07,1.46209 +5.0684e-07,1.46196 +5.09265e-07,1.46183 +5.11702e-07,1.4617 +5.14151e-07,1.46158 +5.16611e-07,1.46145 +5.19083e-07,1.46133 +5.21567e-07,1.4612 +5.24062e-07,1.46108 +5.2657e-07,1.46096 +5.2909e-07,1.46084 +5.31621e-07,1.46072 +5.34165e-07,1.46061 +5.36721e-07,1.46049 +5.39289e-07,1.46038 +5.4187e-07,1.46026 +5.44463e-07,1.46015 +5.47068e-07,1.46004 +5.49686e-07,1.45993 +5.52316e-07,1.45981 +5.54959e-07,1.45971 +5.57614e-07,1.4596 +5.60282e-07,1.45949 +5.62963e-07,1.45938 +5.65657e-07,1.45928 +5.68364e-07,1.45917 +5.71083e-07,1.45907 +5.73816e-07,1.45896 +5.76562e-07,1.45886 +5.79321e-07,1.45876 +5.82093e-07,1.45866 +5.84878e-07,1.45856 +5.87677e-07,1.45846 +5.90489e-07,1.45836 +5.93314e-07,1.45826 +5.96153e-07,1.45817 +5.99006e-07,1.45807 +6.01872e-07,1.45798 +6.04752e-07,1.45788 +6.07646e-07,1.45779 +6.10553e-07,1.45769 +6.13475e-07,1.4576 +6.1641e-07,1.45751 +6.1936e-07,1.45742 +6.22323e-07,1.45733 +6.25301e-07,1.45724 +6.28293e-07,1.45715 +6.313e-07,1.45706 +6.3432e-07,1.45697 +6.37356e-07,1.45689 +6.40405e-07,1.4568 +6.4347e-07,1.45671 +6.46549e-07,1.45663 +6.49642e-07,1.45655 +6.52751e-07,1.45646 +6.55874e-07,1.45638 +6.59013e-07,1.45629 +6.62166e-07,1.45621 +6.65335e-07,1.45613 +6.68518e-07,1.45605 +6.71717e-07,1.45597 +6.74931e-07,1.45589 +6.78161e-07,1.45581 +6.81406e-07,1.45573 +6.84666e-07,1.45565 +6.87942e-07,1.45557 +6.91234e-07,1.4555 +6.94542e-07,1.45542 +6.97865e-07,1.45534 +7.01204e-07,1.45527 +7.0456e-07,1.45519 +7.07931e-07,1.45511 +7.11318e-07,1.45504 +7.14722e-07,1.45497 +7.18142e-07,1.45489 +7.21578e-07,1.45482 +7.25031e-07,1.45474 +7.285e-07,1.45467 +7.31986e-07,1.4546 +7.35489e-07,1.45453 +7.39008e-07,1.45446 +7.42544e-07,1.45438 +7.46097e-07,1.45431 +7.49667e-07,1.45424 +7.53254e-07,1.45417 +7.56859e-07,1.4541 +7.6048e-07,1.45403 +7.64119e-07,1.45397 +7.67775e-07,1.4539 +7.71449e-07,1.45383 +7.75141e-07,1.45376 +7.7885e-07,1.45369 +7.82576e-07,1.45362 +7.86321e-07,1.45356 +7.90084e-07,1.45349 +7.93864e-07,1.45342 +7.97663e-07,1.45336 +8.0148e-07,1.45329 +8.05315e-07,1.45323 +8.09168e-07,1.45316 +8.1304e-07,1.4531 +8.1693e-07,1.45303 +8.20839e-07,1.45297 +8.24767e-07,1.4529 +8.28714e-07,1.45284 +8.32679e-07,1.45277 +8.36663e-07,1.45271 +8.40667e-07,1.45264 +8.44689e-07,1.45258 +8.48731e-07,1.45252 +8.52792e-07,1.45246 +8.56873e-07,1.45239 +8.60973e-07,1.45233 +8.65093e-07,1.45227 +8.69232e-07,1.4522 +8.73392e-07,1.45214 +8.77571e-07,1.45208 +8.8177e-07,1.45202 +8.85989e-07,1.45196 +8.90229e-07,1.45189 +8.94488e-07,1.45183 +8.98769e-07,1.45177 +9.03069e-07,1.45171 +9.0739e-07,1.45165 +9.11732e-07,1.45159 +9.16095e-07,1.45153 +9.20478e-07,1.45147 +9.24883e-07,1.45141 +9.29308e-07,1.45134 +9.33755e-07,1.45128 +9.38223e-07,1.45122 +9.42712e-07,1.45116 +9.47223e-07,1.4511 +9.51756e-07,1.45104 +9.5631e-07,1.45098 +9.60886e-07,1.45092 +9.65484e-07,1.45086 +9.70104e-07,1.4508 +9.74746e-07,1.45074 +9.7941e-07,1.45068 +9.84096e-07,1.45062 +9.88805e-07,1.45056 +9.93536e-07,1.4505 +9.98291e-07,1.45044 +1.00307e-06,1.45038 +1.00787e-06,1.45032 +1.01269e-06,1.45026 +1.01754e-06,1.4502 +1.0224e-06,1.45014 +1.0273e-06,1.45008 +1.03221e-06,1.45002 +1.03715e-06,1.44996 +1.04211e-06,1.4499 +1.0471e-06,1.44983 +1.05211e-06,1.44977 +1.05715e-06,1.44971 +1.0622e-06,1.44965 +1.06729e-06,1.44959 +1.07239e-06,1.44953 +1.07752e-06,1.44947 +1.08268e-06,1.44941 +1.08786e-06,1.44935 +1.09307e-06,1.44929 +1.0983e-06,1.44922 +1.10355e-06,1.44916 +1.10883e-06,1.4491 +1.11414e-06,1.44904 +1.11947e-06,1.44898 +1.12483e-06,1.44891 +1.13021e-06,1.44885 +1.13562e-06,1.44879 +1.14105e-06,1.44873 +1.14651e-06,1.44866 +1.152e-06,1.4486 +1.15751e-06,1.44854 +1.16305e-06,1.44847 +1.16861e-06,1.44841 +1.1742e-06,1.44834 +1.17982e-06,1.44828 +1.18547e-06,1.44822 +1.19114e-06,1.44815 +1.19684e-06,1.44809 +1.20257e-06,1.44802 +1.20832e-06,1.44796 +1.2141e-06,1.44789 +1.21991e-06,1.44782 +1.22575e-06,1.44776 +1.23162e-06,1.44769 +1.23751e-06,1.44762 +1.24343e-06,1.44756 +1.24938e-06,1.44749 +1.25536e-06,1.44742 +1.26137e-06,1.44735 +1.2674e-06,1.44729 +1.27347e-06,1.44722 +1.27956e-06,1.44715 +1.28568e-06,1.44708 +1.29183e-06,1.44701 +1.29801e-06,1.44694 +1.30423e-06,1.44687 +1.31047e-06,1.4468 +1.31674e-06,1.44673 +1.32304e-06,1.44666 +1.32937e-06,1.44658 +1.33573e-06,1.44651 +1.34212e-06,1.44644 +1.34854e-06,1.44637 +1.355e-06,1.44629 +1.36148e-06,1.44622 +1.36799e-06,1.44615 +1.37454e-06,1.44607 +1.38112e-06,1.446 +1.38773e-06,1.44592 +1.39437e-06,1.44584 +1.40104e-06,1.44577 +1.40774e-06,1.44569 +1.41448e-06,1.44561 +1.42125e-06,1.44554 +1.42805e-06,1.44546 +1.43488e-06,1.44538 +1.44175e-06,1.4453 +1.44864e-06,1.44522 +1.45558e-06,1.44514 +1.46254e-06,1.44506 +1.46954e-06,1.44497 +1.47657e-06,1.44489 +1.48364e-06,1.44481 +1.49074e-06,1.44473 +1.49787e-06,1.44464 +1.50504e-06,1.44456 +1.51224e-06,1.44447 +1.51947e-06,1.44439 +1.52674e-06,1.4443 +1.53405e-06,1.44421 +1.54139e-06,1.44413 +1.54877e-06,1.44404 +1.55618e-06,1.44395 +1.56362e-06,1.44386 +1.57111e-06,1.44377 +1.57862e-06,1.44368 +1.58618e-06,1.44359 +1.59377e-06,1.44349 +1.60139e-06,1.4434 +1.60906e-06,1.44331 +1.61675e-06,1.44321 +1.62449e-06,1.44312 +1.63226e-06,1.44302 +1.64007e-06,1.44293 +1.64792e-06,1.44283 +1.65581e-06,1.44273 +1.66373e-06,1.44263 +1.67169e-06,1.44253 +1.67969e-06,1.44243 +1.68773e-06,1.44233 +1.6958e-06,1.44223 +1.70392e-06,1.44212 +1.71207e-06,1.44202 +1.72026e-06,1.44191 +1.72849e-06,1.44181 +1.73677e-06,1.4417 +1.74508e-06,1.44159 +1.75343e-06,1.44149 +1.76182e-06,1.44138 +1.77025e-06,1.44127 +1.77872e-06,1.44116 +1.78723e-06,1.44104 +1.79578e-06,1.44093 +1.80437e-06,1.44082 +1.81301e-06,1.4407 +1.82168e-06,1.44058 +1.8304e-06,1.44047 +1.83916e-06,1.44035 +1.84796e-06,1.44023 +1.8568e-06,1.44011 +1.86568e-06,1.43999 +1.87461e-06,1.43986 +1.88358e-06,1.43974 +1.8926e-06,1.43962 +1.90165e-06,1.43949 +1.91075e-06,1.43936 +1.91989e-06,1.43923 +1.92908e-06,1.4391 +1.93831e-06,1.43897 +1.94759e-06,1.43884 +1.9569e-06,1.43871 +1.96627e-06,1.43857 +1.97568e-06,1.43844 +1.98513e-06,1.4383 +1.99463e-06,1.43816 +2.00417e-06,1.43802 +2.01376e-06,1.43788 +2.0234e-06,1.43774 +2.03308e-06,1.4376 +2.04281e-06,1.43745 +2.05258e-06,1.43731 +2.06241e-06,1.43716 +2.07227e-06,1.43701 +2.08219e-06,1.43686 +2.09215e-06,1.43671 +2.10216e-06,1.43655 +2.11222e-06,1.4364 +2.12233e-06,1.43624 +2.13249e-06,1.43608 +2.14269e-06,1.43592 +2.15294e-06,1.43576 +2.16324e-06,1.4356 +2.1736e-06,1.43544 +2.184e-06,1.43527 +2.19445e-06,1.4351 +2.20495e-06,1.43493 +2.2155e-06,1.43476 +2.2261e-06,1.43459 +2.23675e-06,1.43442 +2.24745e-06,1.43424 +2.25821e-06,1.43406 +2.26901e-06,1.43388 +2.27987e-06,1.4337 +2.29078e-06,1.43352 +2.30174e-06,1.43333 +2.31275e-06,1.43314 +2.32382e-06,1.43296 +2.33494e-06,1.43276 +2.34611e-06,1.43257 +2.35734e-06,1.43238 +2.36862e-06,1.43218 +2.37995e-06,1.43198 +2.39134e-06,1.43178 +2.40278e-06,1.43158 +2.41428e-06,1.43137 +2.42583e-06,1.43116 +2.43744e-06,1.43095 +2.4491e-06,1.43074 +2.46082e-06,1.43053 +2.4726e-06,1.43031 +2.48443e-06,1.43009 +2.49632e-06,1.42987 +2.50826e-06,1.42965 +2.52026e-06,1.42942 +2.53232e-06,1.42919 +2.54444e-06,1.42896 +2.55662e-06,1.42873 +2.56885e-06,1.42849 +2.58114e-06,1.42826 +2.59349e-06,1.42802 +2.6059e-06,1.42777 +2.61837e-06,1.42753 +2.6309e-06,1.42728 +2.64349e-06,1.42703 +2.65614e-06,1.42677 +2.66885e-06,1.42652 +2.68162e-06,1.42626 +2.69445e-06,1.42599 +2.70734e-06,1.42573 +2.7203e-06,1.42546 +2.73331e-06,1.42519 +2.74639e-06,1.42492 +2.75953e-06,1.42464 +2.77274e-06,1.42436 +2.78601e-06,1.42407 +2.79934e-06,1.42379 +2.81273e-06,1.4235 +2.82619e-06,1.4232 +2.83971e-06,1.42291 +2.8533e-06,1.42261 +2.86695e-06,1.42231 +2.88067e-06,1.422 +2.89446e-06,1.42169 +2.90831e-06,1.42138 +2.92222e-06,1.42106 +2.93621e-06,1.42074 +2.95026e-06,1.42041 +2.96437e-06,1.42008 +2.97856e-06,1.41975 +2.99281e-06,1.41942 +3.00713e-06,1.41908 +3.02152e-06,1.41873 +3.03598e-06,1.41839 +3.0505e-06,1.41803 +3.0651e-06,1.41768 +3.07977e-06,1.41732 +3.0945e-06,1.41695 +3.10931e-06,1.41659 +3.12419e-06,1.41621 +3.13914e-06,1.41584 +3.15416e-06,1.41545 +3.16925e-06,1.41507 +3.18442e-06,1.41468 +3.19965e-06,1.41428 +3.21496e-06,1.41388 +3.23035e-06,1.41348 +3.24581e-06,1.41307 +3.26134e-06,1.41265 +3.27694e-06,1.41223 +3.29262e-06,1.41181 +3.30838e-06,1.41138 +3.32421e-06,1.41094 +3.34011e-06,1.4105 +3.3561e-06,1.41005 +3.37216e-06,1.4096 +3.38829e-06,1.40915 +3.4045e-06,1.40868 +3.42079e-06,1.40821 +3.43716e-06,1.40774 +3.45361e-06,1.40726 +3.47014e-06,1.40677 +3.48674e-06,1.40628 +3.50342e-06,1.40578 +3.52019e-06,1.40528 +3.53703e-06,1.40477 +3.55396e-06,1.40425 +3.57096e-06,1.40373 +3.58805e-06,1.40319 +3.60522e-06,1.40266 +3.62247e-06,1.40211 +3.6398e-06,1.40156 +3.65722e-06,1.401 +3.67472e-06,1.40044 +3.6923e-06,1.39986 +3.70997e-06,1.39928 +3.72772e-06,1.39869 +3.74556e-06,1.3981 +3.76348e-06,1.39749 +3.78149e-06,1.39688 +3.79958e-06,1.39626 +3.81777e-06,1.39563 +3.83603e-06,1.39499 +3.85439e-06,1.39435 +3.87283e-06,1.39369 +3.89136e-06,1.39303 +3.90998e-06,1.39236 +3.92869e-06,1.39168 +3.94749e-06,1.39099 +3.96638e-06,1.39029 +3.98536e-06,1.38958 +4.00443e-06,1.38886 +4.02359e-06,1.38813 +4.04284e-06,1.38739 +4.06219e-06,1.38664 +4.08163e-06,1.38588 +4.10116e-06,1.38511 +4.12078e-06,1.38433 +4.1405e-06,1.38354 +4.16031e-06,1.38273 +4.18022e-06,1.38192 +4.20022e-06,1.38109 +4.22032e-06,1.38025 +4.24051e-06,1.3794 +4.2608e-06,1.37854 +4.28119e-06,1.37766 +4.30168e-06,1.37677 +4.32226e-06,1.37587 +4.34294e-06,1.37495 +4.36372e-06,1.37403 +4.3846e-06,1.37308 +4.40558e-06,1.37213 +4.42666e-06,1.37115 +4.44785e-06,1.37017 +4.46913e-06,1.36917 +4.49051e-06,1.36815 +4.512e-06,1.36712 +4.53359e-06,1.36607 +4.55528e-06,1.36501 +4.57708e-06,1.36393 +4.59898e-06,1.36283 +4.62099e-06,1.36172 +4.6431e-06,1.36058 +4.66532e-06,1.35943 +4.68764e-06,1.35826 +4.71007e-06,1.35708 +4.73261e-06,1.35587 +4.75525e-06,1.35464 +4.77801e-06,1.3534 +4.80087e-06,1.35213 +4.82384e-06,1.35085 +4.84692e-06,1.34954 +4.87012e-06,1.34821 +4.89342e-06,1.34685 +4.91684e-06,1.34548 +4.94036e-06,1.34408 +4.964e-06,1.34266 +4.98775e-06,1.34121 +5.01162e-06,1.33974 +5.0356e-06,1.33824 +5.0597e-06,1.33672 +5.08391e-06,1.33517 +5.10823e-06,1.33359 +5.13268e-06,1.33198 +5.15724e-06,1.33035 +5.18191e-06,1.32868 +5.20671e-06,1.32699 +5.23162e-06,1.32526 +5.25666e-06,1.32351 +5.28181e-06,1.32172 +5.30708e-06,1.31989 +5.33248e-06,1.31803 +5.35799e-06,1.31614 +5.38363e-06,1.31421 +5.40939e-06,1.31225 +5.43528e-06,1.31024 +5.46128e-06,1.3082 +5.48742e-06,1.30612 +5.51367e-06,1.30399 +5.54006e-06,1.30182 +5.56657e-06,1.29961 +5.5932e-06,1.29736 +5.61996e-06,1.29505 +5.64686e-06,1.29271 +5.67388e-06,1.29031 +5.70103e-06,1.28786 +5.72831e-06,1.28536 +5.75572e-06,1.2828 +5.78326e-06,1.2802 +5.81093e-06,1.27753 +5.83873e-06,1.27481 +5.86667e-06,1.27202 +5.89474e-06,1.26918 +5.92295e-06,1.26627 +5.95129e-06,1.26329 +5.97977e-06,1.26024 +6.00838e-06,1.25713 +6.03713e-06,1.25394 +6.06602e-06,1.25067 +6.09505e-06,1.24733 +6.12421e-06,1.24391 +6.15351e-06,1.2404 +6.18296e-06,1.2368 +6.21254e-06,1.23312 +6.24227e-06,1.22934 +6.27214e-06,1.22547 +6.30215e-06,1.22149 +6.33231e-06,1.21741 +6.36261e-06,1.21323 +6.39305e-06,1.20893 +6.42364e-06,1.20451 +6.45438e-06,1.19997 +6.48527e-06,1.19531 +6.5163e-06,1.19051 +6.54748e-06,1.18558 +6.57881e-06,1.1805 +6.61029e-06,1.17528 +6.64192e-06,1.16989 +6.6737e-06,1.16435 +6.70563e-06,1.15863 +6.73772e-06,1.15273 +6.76996e-06,1.14665 +6.80235e-06,1.14037 +6.8349e-06,1.13389 +6.86761e-06,1.12719 +6.90047e-06,1.12026 +6.93349e-06,1.1131 +6.96666e-06,1.10568 +7e-06,1.098 diff --git a/src/scgenerator/helpers.py b/src/scgenerator/helpers.py index f167113..1289047 100644 --- a/src/scgenerator/helpers.py +++ b/src/scgenerator/helpers.py @@ -2,12 +2,14 @@ series of helper functions """ -from scgenerator.physics.materials import n_gas_2 -from scgenerator.physics.fiber import n_eff_marcatili, beta2, beta2_to_D -from scgenerator.physics.units import c import numpy as np -__all__ = ["capillary_dispersion"] +from scgenerator.math import all_zeros +from scgenerator.physics.fiber import beta2, n_eff_marcatili, n_eff_hasan +from scgenerator.physics.materials import n_gas_2 +from scgenerator.physics.units import c + +__all__ = ["capillary_dispersion", "capillary_zdw", "revolver_dispersion"] def capillary_dispersion( @@ -43,11 +45,100 @@ def capillary_dispersion( return beta2(w, n)[2:-2] -def extend_axis(wl): - dwl_left = wl[1] - wl[0] - dwl_right = wl[-1] - wl[-2] - wl = np.concatenate( - ([wl[0] - 2 * dwl_left, wl[0] - dwl_left], wl, [wl[-1] + dwl_right, wl[-1] + 2 * dwl_right]) +def capillary_zdw( + radius: float, + gas_name: str, + pressure=None, + temperature=None, + search_range: tuple[float, float] = (200e-9, 3000e-9), +) -> np.ndarray: + """find the zero dispersion wavelength of a capilally + + Parameters + ---------- + radius : float + in mnm + gas_name : str + gas name (case insensitive) + pressure : float, optional + pressure in Pa (multiply mbar by 100 to get Pa), by default atm pressure + temperature : float, optional + temperature in K, by default 20°C + search_range : (float, float), optional + range of wavelength (in m) in which to search for the ZDW, by default (200e-9, 3000e-9) + + Returns + ------- + np.ndarray + array of zero dispersion wavelength(s) + """ + wl = np.linspace(*search_range[:2], 1024) + disp = capillary_dispersion(wl, radius, gas_name, pressure, temperature) + return all_zeros(wl, disp) + + +def revolver_dispersion( + wl: np.ndarray, + core_radius: float, + gas_name: str, + capillary_num: int, + capillary_thickness: float, + capillary_spacing: float, + capillary_nested: int = 0, + capillary_resonance_strengths: list[float] = None, + pressure=None, + temperature=None, +) -> np.ndarray: + """computes the dispersion (beta2) of a capillary tube + + Parameters + ---------- + wl : np.ndarray + wavelength in m + radius : float + core radius in m + gas_name : str + gas name (case insensitive) + pressure : float, optional + pressure in Pa (multiply mbar by 100 to get Pa), by default atm pressure + temperature : float, optional + temperature in K, by default 20°C + + Returns + ------- + np.ndarray + D parameter + """ + capillary_resonance_strengths = capillary_resonance_strengths or [] + wl = extend_axis(wl) + if pressure is None: + pressure = 101325 + if temperature is None: + temperature = 293.15 + n = n_eff_hasan( + wl, + n_gas_2(wl, gas_name.lower(), pressure, temperature, False), + core_radius, + capillary_num, + capillary_nested, + capillary_thickness, + capillary_spacing, + capillary_resonance_strengths, + ) + w = 2 * np.pi * c / wl + return beta2(w, n)[2:-2] + + +def extend_axis(axis: np.ndarray) -> np.ndarray: + """add 4 values to an array, 2 on each 'side'""" + dwl_left = axis[1] - axis[0] + dwl_right = axis[-1] - axis[-2] + axis = np.concatenate( + ( + [axis[0] - 2 * dwl_left, axis[0] - dwl_left], + axis, + [axis[-1] + dwl_right, axis[-1] + 2 * dwl_right], + ) ) - return wl + return axis diff --git a/src/scgenerator/physics/fiber.py b/src/scgenerator/physics/fiber.py index 74fa9ad..6d66caa 100644 --- a/src/scgenerator/physics/fiber.py +++ b/src/scgenerator/physics/fiber.py @@ -4,15 +4,14 @@ import numpy as np from numpy import e from numpy.fft import fft from numpy.polynomial.chebyshev import Chebyshev, cheb2poly +from scgenerator import utils +from scgenerator.cache import np_cache +from scgenerator.math import argclosest, u_nm +from scgenerator.physics import materials as mat +from scgenerator.physics import units +from scgenerator.physics.units import c, pi from scipy.interpolate import interp1d -from .. import utils -from ..cache import np_cache -from ..math import argclosest, u_nm -from . import materials as mat -from . import units -from .units import c, pi - pipi = 2 * pi T = TypeVar("T") @@ -112,12 +111,12 @@ def gvd_from_n_eff(n_eff: np.ndarray, wl_for_disp: np.ndarray): def beta2_to_D(beta2, wl_for_disp): """returns the D parameter corresponding to beta2(wl_for_disp)""" - return -(pipi * c) / (wl_for_disp ** 2) * beta2 + return -(pipi * c) / (wl_for_disp**2) * beta2 def D_to_beta2(D, wl_for_disp): """returns the beta2 parameters corresponding to D(wl_for_disp)""" - return -(wl_for_disp ** 2) / (pipi * c) * D + return -(wl_for_disp**2) / (pipi * c) * D def A_to_C(A: np.ndarray, A_eff_arr: np.ndarray) -> np.ndarray: @@ -147,9 +146,9 @@ def plasma_dispersion(wl_for_disp, number_density, simple=False): w = units.m(wl_for_disp) if simple: w_pl = number_density * e2_me_e0 - return -(w_pl ** 2) / (c * w ** 2) + return -(w_pl**2) / (c * w**2) - beta2_arr = beta2(w, np.sqrt(1 - number_density * e2_me_e0 / w ** 2)) + beta2_arr = beta2(w, np.sqrt(1 - number_density * e2_me_e0 / w**2)) return beta2_arr @@ -224,7 +223,7 @@ def A_eff_marcatili(core_radius: float) -> float: float effective mode field area """ - return 1.5 * core_radius ** 2 + return 1.5 * core_radius**2 def capillary_spacing_hasan( @@ -261,7 +260,7 @@ def resonance_thickness( wl_for_disp / (4 * np.sqrt(n_si_2)) * (2 * order + 1) - * (1 - n_gas_2 / n_si_2 + wl_for_disp ** 2 / (4 * n_si_2 * core_radius ** 2)) ** -0.5 + * (1 - n_gas_2 / n_si_2 + wl_for_disp**2 / (4 * n_si_2 * core_radius**2)) ** -0.5 ) @@ -328,7 +327,9 @@ def n_eff_hasan( Reference ---------- - Hasan, Md Imran, Nail Akhmediev, and Wonkeun Chang. "Empirical formulae for dispersion and effective mode area in hollow-core antiresonant fibers." Journal of Lightwave Technology 36.18 (2018): 4060-4065. + Hasan, Md Imran, Nail Akhmediev, and Wonkeun Chang. "Empirical formulae for dispersion and + effective mode area in hollow-core antiresonant fibers." Journal of Lightwave Technology 36.18 + (2018): 4060-4065. """ u = u_nm(1, 1) alpha = 5e-12 @@ -340,7 +341,7 @@ def n_eff_hasan( if capillary_nested > 0: f2 += 0.0045 * np.exp(-4.1589 / (capillary_nested * Rg)) - R_eff = f1 * core_radius * (1 - f2 * wl_for_disp ** 2 / (core_radius * capillary_thickness)) + R_eff = f1 * core_radius * (1 - f2 * wl_for_disp**2 / (core_radius * capillary_thickness)) n_eff_2 = n_gas_2 - (u * wl_for_disp / (pipi * R_eff)) ** 2 @@ -350,8 +351,8 @@ def n_eff_hasan( for m, strength in enumerate(capillary_resonance_strengths): n_eff_2 += ( strength - * wl_for_disp ** 2 - / (alpha + wl_for_disp ** 2 - chi_sil * (2 * capillary_thickness / (m + 1)) ** 2) + * wl_for_disp**2 + / (alpha + wl_for_disp**2 - chi_sil * (2 * capillary_thickness / (m + 1)) ** 2) ) return np.sqrt(n_eff_2) @@ -374,7 +375,7 @@ def A_eff_hasan(core_radius, capillary_num, capillary_spacing): A_eff : float """ M_f = 1.5 / (1 - 0.5 * np.exp(-0.245 * capillary_num)) - return M_f * core_radius ** 2 * np.exp((capillary_spacing / 22e-6) ** 2.5) + return M_f * core_radius**2 * np.exp((capillary_spacing / 22e-6) ** 2.5) def V_eff_step_index( @@ -437,8 +438,8 @@ def V_parameter_koshiba(l: np.ndarray, pitch: float, pitch_ratio: float) -> floa V = A[0] + A[1] / (1 + A[2] * np.exp(A[3] * ratio_l)) - n_FSM2 = 1.45 ** 2 - (l * V / (pi2a)) ** 2 - V_eff = pi2a / l * np.sqrt(n_co ** 2 - n_FSM2) + n_FSM2 = 1.45**2 - (l * V / (pi2a)) ** 2 + V_eff = pi2a / l * np.sqrt(n_co**2 - n_FSM2) return V_eff @@ -627,7 +628,6 @@ def beta2(w_for_disp: np.ndarray, n_eff: np.ndarray) -> np.ndarray: return np.gradient(np.gradient(beta(w_for_disp, n_eff), w_for_disp), w_for_disp) - def frame_velocity(beta1_arr: np.ndarray, w0_ind: int) -> float: return 1.0 / beta1_arr[w0_ind] @@ -719,7 +719,7 @@ def dynamic_HCPCF_dispersion( in the fiber """ - A_eff = 1.5 * core_radius ** 2 + A_eff = 1.5 * core_radius**2 # defining function instead of storing every possilble value def pressure(r): @@ -811,7 +811,7 @@ def n_eff_pcf(wl_for_disp: np.ndarray, pitch: float, pitch_ratio: float) -> np.n V = A[0] + A[1] / (1 + A[2] * np.exp(A[3] * ratio_l)) W = B[0] + B[1] / (1 + B[2] * np.exp(B[3] * ratio_l)) - n_FSM2 = 1.45 ** 2 - (wl_for_disp * V / (pi2a)) ** 2 + n_FSM2 = 1.45**2 - (wl_for_disp * V / (pi2a)) ** 2 n_eff2 = (wl_for_disp * W / (pi2a)) ** 2 + n_FSM2 n_eff = np.sqrt(n_eff2) @@ -845,8 +845,8 @@ def saitoh_paramters(pitch_ratio: float) -> tuple[float, float]: di2 = np.array([9, 6.58, 10, 0.41]) di3 = np.array([10, 24.8, 15, 6]) - A = ai0 + ai1 * pitch_ratio ** bi1 + ai2 * pitch_ratio ** bi2 + ai3 * pitch_ratio ** bi3 - B = ci0 + ci1 * pitch_ratio ** di1 + ci2 * pitch_ratio ** di2 + ci3 * pitch_ratio ** di3 + A = ai0 + ai1 * pitch_ratio**bi1 + ai2 * pitch_ratio**bi2 + ai3 * pitch_ratio**bi3 + B = ci0 + ci1 * pitch_ratio**di1 + ci2 * pitch_ratio**di2 + ci3 * pitch_ratio**di3 return A, B @@ -960,7 +960,7 @@ def dispersion_from_coefficients( coef = np.array(beta2_coefficients) / np.cumprod([1] + list(range(1, len(beta2_coefficients)))) beta_arr = np.zeros_like(w_c) for k, b in reversed(list(enumerate(coef))): - beta_arr = beta_arr + b * w_c ** k + beta_arr = beta_arr + b * w_c**k return beta_arr @@ -985,12 +985,12 @@ def delayed_raman_t(t: np.ndarray, raman_type: str) -> np.ndarray: t_ = t - t[0] t = t_ if raman_type == "stolen": - hr_arr = (tau1 / tau2 ** 2 + 1 / tau1) * np.exp(-t_ / tau2) * np.sin(t_ / tau1) + hr_arr = (tau1 / tau2**2 + 1 / tau1) * np.exp(-t_ / tau2) * np.sin(t_ / tau1) elif raman_type == "agrawal": taub = 96e-15 - h_a = (tau1 / tau2 ** 2 + 1 / tau1) * np.exp(-t_ / tau2) * np.sin(t_ / tau1) - h_b = (2 * taub - t_) / taub ** 2 * np.exp(-t_ / taub) + h_a = (tau1 / tau2**2 + 1 / tau1) * np.exp(-t_ / tau2) * np.sin(t_ / tau1) + h_b = (2 * taub - t_) / taub**2 * np.exp(-t_ / taub) hr_arr = 0.79 * h_a + 0.21 * h_b elif raman_type == "measured": @@ -1096,12 +1096,12 @@ def effective_core_radius(wl_for_disp, core_radius, s=0.08, h=200e-9): ------- effective_core_radius : ndarray, shape (n, ) """ - return core_radius / (1 + s * wl_for_disp ** 2 / (core_radius * h)) + return core_radius / (1 + s * wl_for_disp**2 / (core_radius * h)) def effective_radius_HCARF(core_radius, t, f1, f2, wl_for_disp): """eq. 3 in Hasan 2018""" - return f1 * core_radius * (1 - f2 * wl_for_disp ** 2 / (core_radius * t)) + return f1 * core_radius * (1 - f2 * wl_for_disp**2 / (core_radius * t)) def capillary_loss(wl: np.ndarray, he_mode: tuple[int, int], core_radius: float) -> np.ndarray: @@ -1123,7 +1123,7 @@ def capillary_loss(wl: np.ndarray, he_mode: tuple[int, int], core_radius: float) """ chi_silica = abs(mat.sellmeier(wl, utils.load_material_dico("silica"))) nu_n = 0.5 * (chi_silica + 2) / np.sqrt(chi_silica) - return nu_n * (u_nm(*he_mode) * wl / pipi) ** 2 * core_radius ** -3 + return nu_n * (u_nm(*he_mode) * wl / pipi) ** 2 * core_radius**-3 def extinction_distance(loss: T, ratio=1 / e) -> T: diff --git a/src/scgenerator/physics/materials.py b/src/scgenerator/physics/materials.py index 41f630b..fec71f0 100644 --- a/src/scgenerator/physics/materials.py +++ b/src/scgenerator/physics/materials.py @@ -1,14 +1,15 @@ +from __future__ import annotations + import functools from dataclasses import dataclass, field from typing import Any, TypeVar import numpy as np - -from .. import utils -from ..cache import np_cache -from ..logger import get_logger -from . import units -from .units import NA, c, epsilon0, kB +from scgenerator import utils +from scgenerator.cache import np_cache +from scgenerator.logger import get_logger +from scgenerator.physics import units +from scgenerator.physics.units import NA, c, epsilon0, kB T = TypeVar("T", np.floating, np.ndarray) @@ -22,6 +23,21 @@ class Sellmeier: kind: int = 2 constant: float = 0 + @classmethod + def load(cls, name: str) -> Sellmeier: + mat_dico = utils.load_material_dico(name) + s = mat_dico.get("sellmeier", {}) + return cls( + **{ + newk: s.get(k, None) + for newk, k in zip( + ["B", "C", "pressure_ref", "temperature_ref", "kind", "constant"], + ["B", "C", "P0", "T0", "kind", "const"], + ) + if k in s + } + ) + def chi(self, wl: T) -> T: """n^2 - 1""" if isinstance(wl, np.ndarray): @@ -30,10 +46,10 @@ class Sellmeier: chi = 0 if self.kind == 1: for b, c_ in zip(self.B, self.C): - chi += wl ** 2 * b / (wl ** 2 - c_) + chi += wl**2 * b / (wl**2 - c_) elif self.kind == 2: # gives n-1 for b, c_ in zip(self.B, self.C): - chi += b / (c_ - 1 / wl ** 2) + chi += b / (c_ - 1 / wl**2) chi += self.constant chi = (chi + 1) ** 2 - 1 elif self.kind == 3: # Schott formula @@ -201,7 +217,7 @@ def pressure_from_gradient(ratio, p0, p1): ---------- the pressure (float) """ - return np.sqrt(p0 ** 2 - ratio * (p0 ** 2 - p1 ** 2)) + return np.sqrt(p0**2 - ratio * (p0**2 - p1**2)) def number_density_van_der_waals( @@ -241,7 +257,7 @@ def number_density_van_der_waals( pressure = 101325 if pressure is None else pressure temperature = 273.15 if temperature is None else temperature - ap = a / NA ** 2 + ap = a / NA**2 bp = b / NA # setup van der Waals equation for the number density @@ -306,11 +322,11 @@ def sellmeier( if kind == 1: logger.debug("materials : using Sellmeier 1st kind equation") for b, c_ in zip(B, C): - chi[ind] += temp_l ** 2 * b / (temp_l ** 2 - c_) + chi[ind] += temp_l**2 * b / (temp_l**2 - c_) elif kind == 2: # gives n-1 logger.debug("materials : using Sellmeier 2nd kind equation") for b, c_ in zip(B, C): - chi[ind] += b / (c_ - 1 / temp_l ** 2) + chi[ind] += b / (c_ - 1 / temp_l**2) chi += const chi = (chi + 1) ** 2 - 1 elif kind == 3: # Schott formula @@ -433,8 +449,8 @@ def gas_chi3(gas_name: str, wavelength: float, pressure: float, temperature: flo def n2_to_chi3(n2: float, n0: float) -> float: - return n2 * 4 * epsilon0 * n0 ** 2 * c / 3 + return n2 * 4 * epsilon0 * n0**2 * c / 3 def chi3_to_n2(chi3: float, n0: float) -> float: - return 3.0 * chi3 / (4.0 * epsilon0 * c * n0 ** 2) + return 3.0 * chi3 / (4.0 * epsilon0 * c * n0**2) diff --git a/src/scgenerator/utils.py b/src/scgenerator/utils.py index 0b3067e..f1627a1 100644 --- a/src/scgenerator/utils.py +++ b/src/scgenerator/utils.py @@ -14,16 +14,16 @@ from dataclasses import dataclass from functools import cache, lru_cache from pathlib import Path from string import printable as str_printable -from typing import Any, Callable, MutableMapping, Sequence, TypeVar, Set, Union +from typing import Any, Callable, MutableMapping, Sequence, TypeVar, Union import numpy as np import pkg_resources as pkg import tomli import tomli_w -from .const import PARAM_FN, PARAM_SEPARATOR, SPEC1_FN, Z_FN, ROOT_PARAMETERS -from .logger import get_logger -from .errors import DuplicateParameterError +from scgenerator.const import PARAM_FN, PARAM_SEPARATOR, ROOT_PARAMETERS, SPEC1_FN, Z_FN +from scgenerator.errors import DuplicateParameterError +from scgenerator.logger import get_logger T_ = TypeVar("T_") @@ -85,8 +85,8 @@ class Paths: class SubConfig: fixed: dict[str, Any] variable: list[dict[str, list]] - fixed_keys: Set[str] - variable_keys: Set[str] + fixed_keys: set[str] + variable_keys: set[str] def __init__(self, dico: dict[str, Any]): dico = dico.copy()