patterson-data.cpp

Go to the documentation of this file.

/* 
 * This file is part of PLaSK (https://plask.app) by Photonics Group at TUL
 * Copyright (c) 2022 Lodz University of Technology
 * 
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, version 3.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU General Public License for more details.
 */
#include "patterson-data.hpp"
 
namespace plask { namespace optical { namespace effective {
 
const double patterson_points[] = {
    0.000000000000000,
    0.774596669241483,
    0.434243749346803,
    0.960491268708020,
    0.223386686428967,
    0.621102946737226,
    0.888459232872257,
    0.993831963212755,
    0.112488943133187,
    0.331135393257977,
    0.531319743644376,
    0.702496206491527,
    0.836725938168869,
    0.929654857429740,
    0.981531149553740,
    0.999098124967668,
    0.056344313046593,
    0.168235251552207,
    0.277749822021824,
    0.383359324198730,
    0.483618026945841,
    0.577195710052046,
    0.662909660024781,
    0.739756044352695,
    0.806940531950218,
    0.863907938193690,
    0.910371156957004,
    0.946342858373403,
    0.972182874748582,
    0.988684757547429,
    0.997206259372222,
    0.999872888120358,
    0.028184648949746,
    0.084454040083711,
    0.140424233152560,
    0.195897502711100,
    0.250678730303483,
    0.304576441556714,
    0.357403837831532,
    0.408979821229889,
    0.459130011989832,
    0.507687757533717,
    0.554495132631933,
    0.599403930242243,
    0.642276642509759,
    0.682987431091079,
    0.721423085370099,
    0.757483966380514,
    0.791084933799848,
    0.822156254364980,
    0.850644494768350,
    0.876513414484705,
    0.899744899776940,
    0.920340025470012,
    0.938320397779593,
    0.953730006425761,
    0.966637851558417,
    0.977141514639706,
    0.985371499598520,
    0.991495721178106,
    0.995724104698407,
    0.998316635318407,
    0.999598799671911,
    0.999982430354892,
    0.014093886410782,
    0.042269164765364,
    0.070406976042855,
    0.098482396598119,
    0.126470584372302,
    0.154346811481378,
    0.182086496759252,
    0.209665238243181,
    0.237058845589830,
    0.264243372410927,
    0.291195148518247,
    0.317890812068477,
    0.344307341599438,
    0.370422087950078,
    0.396212806057616,
    0.421657686626163,
    0.446735387662028,
    0.471425065871659,
    0.495706407918761,
    0.519559661537457,
    0.542965666498311,
    0.565905885423654,
    0.588362434447663,
    0.610318113715186,
    0.631756437711194,
    0.652661665410017,
    0.673018830230419,
    0.692813769779115,
    0.712033155362252,
    0.730664521242181,
    0.748696293616937,
    0.766117819303760,
    0.782919394118283,
    0.799092290960841,
    0.814628787655137,
    0.829522194637401,
    0.843766882672709,
    0.857358310886232,
    0.870293055548114,
    0.882568840247342,
    0.894184568335559,
    0.905140358813262,
    0.915437587155765,
    0.925078932907076,
    0.934068436157726,
    0.942411565191083,
    0.950115297521295,
    0.957188216109861,
    0.963640621569812,
    0.969484659502459,
    0.974734459752403,
    0.979406281670863,
    0.983518657578633,
    0.987092527954034,
    0.990151370400770,
    0.992721344282789,
    0.994831502800621,
    0.996514145914890,
    0.997805354495957,
    0.998745614468095,
    0.999380338025024,
    0.999760490924432,
    0.999943996207054,
    0.999997596379748,
    0.007047138459337,
    0.021139853378331,
    0.035227882808441,
    0.049308104790869,
    0.063377399917322,
    0.077432652349857,
    0.091470750840355,
    0.105488589749542,
    0.119483070065440,
    0.133451100421162,
    0.147389598111940,
    0.161295490111305,
    0.175165714086311,
    0.188997219411722,
    0.202786968183065,
    0.216531936228473,
    0.230229114119222,
    0.243875508178893,
    0.257468141491070,
    0.271004054905513,
    0.284480308042726,
    0.297893980296858,
    0.311242171836872,
    0.324522004605922,
    0.337730623318886,
    0.350865196458001,
    0.363922917266550,
    0.376901004740559,
    0.389796704618471,
    0.402607290368737,
    0.415330064175322,
    0.427962357921063,
    0.440501534168876,
    0.452944987140767,
    0.465290143694635,
    0.477534464298829,
    0.489675444004456,
    0.501710613415392,
    0.513637539655989,
    0.525453827336443,
    0.537157119515795,
    0.548745098662529,
    0.560215487612728,
    0.571566050525743,
    0.582794593837319,
    0.593898967210122,
    0.604877064481584,
    0.615726824608993,
    0.626446232611720,
    0.637033320510492,
    0.647486168263572,
    0.657802904699714,
    0.667981708447750,
    0.678020808862644,
    0.687918486947839,
    0.697673076273711,
    0.707282963891961,
    0.716746591245747,
    0.726062455075390,
    0.735229108319492,
    0.744245161011347,
    0.753109281170558,
    0.761820195689839,
    0.770376691217077,
    0.778777615032823,
    0.787021875923539,
    0.795108445051101,
    0.803036356819269,
    0.810804709738147,
    0.818412667287926,
    0.825859458783650,
    0.833144380243173,
    0.840266795261030,
    0.847226135891581,
    0.854021903545469,
    0.860653669904300,
    0.867121077859315,
    0.873423842480859,
    0.879561752026556,
    0.885534668997285,
    0.891342531251320,
    0.896985353188317,
    0.902463227016166,
    0.907776324115059,
    0.912924896514371,
    0.917909278499078,
    0.922729888363349,
    0.927387230329537,
    0.931881896650954,
    0.936214569916451,
    0.940386025573670,
    0.944397134685867,
    0.948248866934137,
    0.951942293872574,
    0.955478592438184,
    0.958859048710200,
    0.962085061904651,
    0.965158148579916,
    0.968079947017760,
    0.970852221732792,
    0.973476868052507,
    0.975955916702012,
    0.978291538324759,
    0.980486047876721,
    0.982541908851081,
    0.984461737328815,
    0.986248305913008,
    0.987904547695124,
    0.989433560520241,
    0.990838611958294,
    0.992123145530863,
    0.993290788851685,
    0.994345364356723,
    0.995290903148810,
    0.996131662079315,
    0.996872143485260,
    0.997517116063472,
    0.998071634524930,
    0.998541055697168,
    0.998931050830811,
    0.999247618943342,
    0.999497112467187,
    0.999686286448318,
    0.999822363679788,
    0.999913081144678,
    0.999966730098486,
    0.999992298136258,
    0.999999672956734
};
 
const double patterson_weights_data[] = {
    // n = 1
    2.,
    // n = 3
    0.888888888888889,
    0.555555555555556,
    // n = 7
    0.450916538658474,
    0.268488089868333,
    0.401397414775962,
    0.104656226026467,
    // n = 15
    0.225510499798207,
    0.134415255243784,
    0.200628529376989,
    0.0516032829970797,
    0.219156858401587,
    0.171511909136391,
    0.0929271953151245,
    0.0170017196299403,
    // n = 31
    0.112755256720769,
    0.0672077542959907,
    0.100314278611796,
    0.0258075980961767,
    0.109578421055925,
    0.0857559200499903,
    0.046462893261758,
    0.00843456573932111,
    0.111956873020953,
    0.105669893580235,
    0.0936271099812645,
    0.0768796204990035,
    0.0569795094941234,
    0.0359571033071293,
    0.0164460498543878,
    0.00254478079156187,
    // n = 63
    0.0563776283603847,
    0.0336038771482077,
    0.0501571393058995,
    0.0129038001003513,
    0.0547892105279629,
    0.0428779600250077,
    0.0232314466399103,
    0.00421763044155885,
    0.0559784365104763,
    0.0528349467901165,
    0.046813554990628,
    0.0384398102494555,
    0.0284897547458336,
    0.0179785515681283,
    0.00822300795723593,
    0.00126515655623007,
    0.0562776998312543,
    0.0554814043565594,
    0.0539054993352661,
    0.0515832539520485,
    0.0485643304066732,
    0.0449145316536322,
    0.0407155101169443,
    0.0360644327807826,
    0.031073551111688,
    0.0258696793272147,
    0.0205942339159127,
    0.0154067504665595,
    0.0104982469096213,
    0.00611550682211725,
    0.00257904979468569,
    0.000363221481845531,
    // n = 127
    0.0281888141801924,
    0.0168019385741039,
    0.0250785696529498,
    0.00645190005017574,
    0.0273946052639814,
    0.0214389800125039,
    0.0116157233199551,
    0.00210881524572663,
    0.0279892182552382,
    0.0264174733950583,
    0.023406777495314,
    0.0192199051247278,
    0.0142448773729168,
    0.00898927578406414,
    0.00411150397865469,
    0.000632607319362634,
    0.0281388499156272,
    0.0277407021782797,
    0.026952749667633,
    0.0257916269760242,
    0.0242821652033366,
    0.0224572658268161,
    0.0203577550584722,
    0.0180322163903913,
    0.015536775555844,
    0.0129348396636074,
    0.0102971169579564,
    0.00770337523327974,
    0.00524912345480886,
    0.00305775341017553,
    0.00128952408261042,
    0.000180739564445388,
    0.0281763190330166,
    0.0280764557938172,
    0.0278772514766137,
    0.0275797495664819,
    0.0271855132296248,
    0.0266966229274504,
    0.0261156733767061,
    0.0254457699654648,
    0.0246905247444877,
    0.0238540521060385,
    0.0229409642293877,
    0.0219563663053178,
    0.020905851445812,
    0.0197954950480975,
    0.0186318482561388,
    0.0174219301594642,
    0.0161732187295777,
    0.0148936416648152,
    0.0135915710097655,
    0.0122758305600828,
    0.0109557333878379,
    0.00964117772970254,
    0.00834283875396816,
    0.00707248999543356,
    0.00584344987583564,
    0.00467105037211432,
    0.0035728927835173,
    0.00256876494379402,
    0.00168114286542147,
    0.000938369848542382,
    0.000377746646326985,
    5.05360952078625e-05,
    // n = 255
    0.0140944070900962,
    0.00840096928705193,
    0.0125392848264749,
    0.00322595002508787,
    0.0136973026319907,
    0.0107194900062519,
    0.00580786165997757,
    0.00105440762286332,
    0.0139946091276191,
    0.0132087366975291,
    0.011703388747657,
    0.00960995256236388,
    0.00712243868645839,
    0.00449463789203207,
    0.00205575198932735,
    0.000316303660822264,
    0.0140694249578136,
    0.0138703510891398,
    0.0134763748338165,
    0.0128958134880121,
    0.0121410826016683,
    0.011228632913408,
    0.0101788775292361,
    0.00901610819519564,
    0.00776838777792199,
    0.00646741983180369,
    0.00514855847897818,
    0.00385168761663987,
    0.00262456172740443,
    0.00152887670508777,
    0.000644762041305725,
    9.03727346587511e-05,
    0.0140881595165083,
    0.0140382278969086,
    0.0139386257383069,
    0.0137898747832409,
    0.0135927566148124,
    0.0133483114637252,
    0.013057836688353,
    0.0127228849827324,
    0.0123452623722438,
    0.0119270260530193,
    0.0114704821146939,
    0.0109781831526589,
    0.010452925722906,
    0.00989774752404875,
    0.0093159241280694,
    0.00871096507973209,
    0.00808660936478886,
    0.00744682083240759,
    0.00679578550488277,
    0.00613791528004138,
    0.00547786669391895,
    0.00482058886485127,
    0.00417141937698408,
    0.00353624499771678,
    0.00292172493791782,
    0.00233552518605716,
    0.00178644639175865,
    0.00128438247189701,
    0.000840571432710722,
    0.00046918492424785,
    0.000188873264506505,
    2.51578703842807e-05,
    0.0140928450691604,
    0.0140803519625537,
    0.01405538207265,
    0.0140179680394566,
    0.0139681588065169,
    0.0139060196013255,
    0.0138316319095064,
    0.0137450934430019,
    0.0136465181025713,
    0.0135360359349562,
    0.0134137930851101,
    0.0132799517439305,
    0.0131346900919602,
    0.0129782022395374,
    0.0128106981638774,
    0.0126324036435421,
    0.012443560190714,
    0.012244424981612,
    0.0120352707852796,
    0.0118163858908302,
    0.011588074033044,
    0.0113506543159806,
    0.0111044611340069,
    0.0108498440893373,
    0.0105871679048852,
    0.0103168123309476,
    0.0100391720440568,
    0.00975465653631741,
    0.00946368999383007,
    0.00916671116356079,
    0.0088641732094825,
    0.00855654356130769,
    0.00824430376303287,
    0.00792794933429485,
    0.00760798966571906,
    0.00728494798055381,
    0.00695936140939042,
    0.00663178124290189,
    0.00630277344908576,
    0.00597291956550817,
    0.00564281810138444,
    0.00531308660518706,
    0.00498436456476554,
    0.00465731729975685,
    0.00433264096809298,
    0.00401106872407502,
    0.00369337791702565,
    0.00338039799108692,
    0.00307301843470258,
    0.00277219576459345,
    0.00247895822665757,
    0.00219440692536384,
    0.00191971297101387,
    0.00165611272815445,
    0.00140490799565514,
    0.00116748411742996,
    0.000945361516858525,
    0.000740282804244503,
    0.000554295314930375,
    0.000389745284473282,
    0.000249212400482997,
    0.000135754910949229,
    5.32752936697806e-05,
    6.93793643241083e-06,
    // n = 511
    0.00704720354504809,
    0.00420048464352597,
    0.00626964241323744,
    0.00161297501254393,
    0.00684865131599536,
    0.00535974500312597,
    0.00290393082998878,
    0.000527203811431658,
    0.00699730456380954,
    0.00660436834876457,
    0.0058516943738285,
    0.00480497628118194,
    0.00356121934322919,
    0.00224731894601603,
    0.00102787599466367,
    0.000158151830411132,
    0.00703471247890679,
    0.00693517554456992,
    0.00673818741690826,
    0.00644790674400606,
    0.00607054130083415,
    0.00561431645670402,
    0.00508943876461804,
    0.00450805409759782,
    0.003884193888961,
    0.00323370991590184,
    0.00257427923948909,
    0.00192584380831994,
    0.00131228086370221,
    0.000764438352543883,
    0.000322381020652862,
    4.51863674126296e-05,
    0.00704407975825415,
    0.00701911394845431,
    0.00696931286915343,
    0.00689493739162047,
    0.0067963783074062,
    0.00667415573186259,
    0.00652891834417652,
    0.00636144249136619,
    0.00617263118612192,
    0.00596351302650963,
    0.00573524105734694,
    0.00548909157632946,
    0.00522646286145301,
    0.00494887376202437,
    0.0046579620640347,
    0.00435548253986604,
    0.00404330468239443,
    0.0037234104162038,
    0.00339789275244139,
    0.00306895764002069,
    0.00273893334695948,
    0.00241029443242563,
    0.00208570968849204,
    0.00176812249885839,
    0.00146086246895891,
    0.00116776259302858,
    0.000893223195879325,
    0.000642191235948505,
    0.000420285716355361,
    0.000234592462123925,
    9.44366322532705e-05,
    1.25792781889593e-05,
    0.0070464225345802,
    0.00704017598127683,
    0.00702769103632498,
    0.0070089840197283,
    0.00698407940325847,
    0.00695300980066273,
    0.00691581595475321,
    0.00687254672150095,
    0.00682325905128565,
    0.00676801796747811,
    0.00670689654255505,
    0.00663997587196526,
    0.00656734504598008,
    0.0064891011197687,
    0.00640534908193868,
    0.00631620182177104,
    0.00622178009535702,
    0.00612221249080599,
    0.00601763539263978,
    0.00590819294541512,
    0.00579403701652198,
    0.0056753271579903,
    0.00555223056700346,
    0.00542492204466866,
    0.0052935839524426,
    0.00515840616547381,
    0.00501958602202842,
    0.00487732826815871,
    0.00473184499691503,
    0.00458335558178039,
    0.00443208660474125,
    0.00427827178065384,
    0.00412215188151643,
    0.00396397466714742,
    0.00380399483285953,
    0.0036424739902769,
    0.00347968070469521,
    0.00331589062145094,
    0.00315138672454288,
    0.00298645978275408,
    0.00282140905069222,
    0.00265654330259353,
    0.00249218228238277,
    0.00232865864987843,
    0.00216632048404649,
    0.00200553436203751,
    0.00184668895851283,
    0.00169019899554346,
    0.00153650921735129,
    0.00138609788229673,
    0.00123947911332878,
    0.00109720346268192,
    0.000959856485506936,
    0.000828056364077226,
    0.000702453997827572,
    0.00058374205871498,
    0.000472680758429263,
    0.000370141402122252,
    0.000277147657465187,
    0.000194872642236641,
    0.000124606200241498,
    6.78774554733972e-05,
    2.66376412339001e-05,
    3.45456507169149e-06,
    0.00704700828844548,
    0.00704544633127951,
    0.00704232292709631,
    0.00703763909614153,
    0.00703139636865429,
    0.00702359678471226,
    0.00701424289402573,
    0.00700333775568107,
    0.00699088493783425,
    0.0069768885173552,
    0.00696135307942367,
    0.00694428371707783,
    0.00692568603071643,
    0.00690556612755588,
    0.00688393062104341,
    0.00686078663022781,
    0.00683614177908911,
    0.00681000419582895,
    0.00678238251212301,
    0.00675328586233753,
    0.00672272388271144,
    0.00669070671050613,
    0.00665724498312455,
    0.00662234983720169,
    0.00658603290766825,
    0.00654830632678944,
    0.00650918272318071,
    0.00646867522080231,
    0.00642679743793437,
    0.00638356348613414,
    0.0063389879691769,
    0.00629308598198199,
    0.00624587310952491,
    0.00619736542573666,
    0.00614757949239084,
    0.00609653235797889,
    0.00604424155657355,
    0.0059907251066801,
    0.0059360015100746,
    0.00588008975062789,
    0.00582300929311348,
    0.00576478008199711,
    0.00570542254020497,
    0.00564495756786715,
    0.00558340654103216,
    0.00552079131034779,
    0.0054571341997031,
    0.00539245800482556,
    0.00532678599182712,
    0.00526014189569259,
    0.00519254991870342,
    0.00512403472879005,
    0.0050546214578065,
    0.00498433569972103,
    0.00491320350871842,
    0.00484125139721057,
    0.00476850633375475,
    0.00469499574088179,
    0.00462074749284081,
    0.00454578991327213,
    0.00447015177282693,
    0.00439386228676004,
    0.0043169511125328,
    0.00423944834747438,
    0.0041613845265651,
    0.00408279062042158,
    0.00400369803358422,
    0.00392413860322996,
    0.00384414459846013,
    0.00376374872034296,
    0.00368298410292404,
    0.00360188431545532,
    0.00352048336613418,
    0.00343881570768791,
    0.00335691624518617,
    0.00327482034651234,
    0.00319256385597435,
    0.00311018311158428,
    0.00302771496658199,
    0.00294519681581858,
    0.00286266662764758,
    0.00278016298199139,
    0.00269772511525295,
    0.00261539297272236,
    0.00253320726907925,
    0.00245120955750556,
    0.00236944230779381,
    0.00228794899365196,
    0.00220677418916003,
    0.00212596367401473,
    0.00204556454679958,
    0.00196562534503151,
    0.00188619617015808,
    0.00180732881501809,
    0.00172907689054462,
    0.00165149594771915,
    0.00157464359003212,
    0.00149857957106457,
    0.00142336587141721,
    0.00134906674928353,
    0.00127574875977347,
    0.00120348074001266,
    0.00113233376051598,
    0.0010623810488534,
    0.000993697899638761,
    0.000926361595613111,
    0.000860451377808528,
    0.000796048517297551,
    0.000733236554224768,
    0.000672101776960108,
    0.000612734008012225,
    0.000555227733977308,
    0.000499683553312801,
    0.000446209810101403,
    0.000394924138246874,
    0.000345954492129904,
    0.000299439176850912,
    0.000255525589595237,
    0.000214368090034217,
    0.000176126765545083,
    0.000140970302204105,
    0.000109085545645742,
    8.06899228014035e-05,
    5.60319507856164e-05,
    3.5375137205519e-05,
    1.90213681905876e-05,
    7.36624069102322e-06,
    9.45715933950007e-07
};
 
const double* patterson_weights[] = {
    patterson_weights_data,
    patterson_weights_data + 1,
    patterson_weights_data + 3,
    patterson_weights_data + 7,
    patterson_weights_data + 15,
    patterson_weights_data + 31,
    patterson_weights_data + 63,
    patterson_weights_data + 127,
    patterson_weights_data + 255
};
 
}}} // namespace plask::optical::effective