Program Listing for File GKIntegratorData.hpp
↰ Return to documentation for file (src/utils/GKIntegratorData.hpp)
#ifndef GKINTEGRATORDATA_HPP
#define GKINTEGRATORDATA_HPP
#include <array>
#include <cstddef>
#include <cstdint>
namespace utils
{
enum class GKOrder : std::uint8_t
{
GK15, // NOLINT(readability-identifier-naming) -- named after the GSL's own gsl_integration_qk15, not a generic enumerator this codebase's usual camelBack convention fits
GK21, // NOLINT(readability-identifier-naming) -- see GK15
GK31, // NOLINT(readability-identifier-naming) -- see GK15
GK41, // NOLINT(readability-identifier-naming) -- see GK15
GK51, // NOLINT(readability-identifier-naming) -- see GK15
GK61 // NOLINT(readability-identifier-naming) -- see GK15
};
template <GKOrder Order>
struct GKRule;
template <>
struct GKRule<GKOrder::GK15>
{
static constexpr std::size_t ng = 4;
static constexpr std::size_t ngk = 8;
static constexpr std::array<double, ngk> xgk = {
0.991455371120812639206854697526329,
0.949107912342758524526189684047851,
0.864864423359769072789712788640926,
0.741531185599394439863864773280788,
0.586087235467691130294144838258730,
0.405845151377397166906606412076961,
0.207784955007898467600689403773245,
0.000000000000000000000000000000000
};
static constexpr std::array<double, ng> wg = {
0.129484966168869693270611432679082,
0.279705391489276667901467771423780,
0.381830050505118944950369775488975,
0.417959183673469387755102040816327
};
static constexpr std::array<double, ngk> wgk = {
0.022935322010529224963732008058970,
0.063092092629978553290700663189204,
0.104790010322250183839876322541518,
0.140653259715525918745189590510238,
0.169004726639267902826583426598550,
0.190350578064785409913256402421014,
0.204432940075298892414161999234649,
0.209482141084727828012999174891714
};
};
template <>
struct GKRule<GKOrder::GK21>
{
static constexpr std::size_t ng = 5;
static constexpr std::size_t ngk = 11;
static constexpr std::array<double, ngk> xgk = {
0.995657163025808080735527280689003,
0.973906528517171720077964012084452,
0.930157491355708226001207180059508,
0.865063366688984510732096688423493,
0.780817726586416897063717578345042,
0.679409568299024406234327365114874,
0.562757134668604683339000099272694,
0.433395394129247190799265943165784, // NOLINT(modernize-use-std-numbers) -- a Kronrod abscissa, not std::numbers::log10e, despite the numerical coincidence
0.294392862701460198131126603103866,
0.148874338981631210884826001129720,
0.000000000000000000000000000000000
};
static constexpr std::array<double, ng> wg = {
0.066671344308688137593568809893332,
0.149451349150580593145776339657697,
0.219086362515982043995534934228163,
0.269266719309996355091226921569469,
0.295524224714752870173892994651338
};
static constexpr std::array<double, ngk> wgk = {
0.011694638867371874278064396062192,
0.032558162307964727478818972459390,
0.054755896574351996031381300244580,
0.075039674810919952767043140916190,
0.093125454583697605535065465083366,
0.109387158802297641899210590325805,
0.123491976262065851077958109831074,
0.134709217311473325928054001771707,
0.142775938577060080797094273138717,
0.147739104901338491374841515972068,
0.149445554002916905664936468389821
};
};
template <>
struct GKRule<GKOrder::GK31>
{
static constexpr std::size_t ng = 8;
static constexpr std::size_t ngk = 16;
static constexpr std::array<double, ngk> xgk = {
0.998002298693397060285172840152271,
0.987992518020485428489565718586613,
0.967739075679139134257347978784337,
0.937273392400705904307758947710209,
0.897264532344081900882509656454496,
0.848206583410427216200648320774217,
0.790418501442465932967649294817947,
0.724417731360170047416186054613938,
0.650996741297416970533735895313275,
0.570972172608538847537226737253911,
0.485081863640239680693655740232351,
0.394151347077563369897207370981045,
0.299180007153168812166780024266389,
0.201194093997434522300628303394596,
0.101142066918717499027074231447392,
0.000000000000000000000000000000000
};
static constexpr std::array<double, ng> wg = {
0.030753241996117268354628393577204,
0.070366047488108124709267416450667,
0.107159220467171935011869546685869,
0.139570677926154314447804794511028,
0.166269205816993933553200860481209,
0.186161000015562211026800561866423,
0.198431485327111576456118326443839,
0.202578241925561272880620199967519
};
static constexpr std::array<double, ngk> wgk = {
0.005377479872923348987792051430128,
0.015007947329316122538374763075807,
0.025460847326715320186874001019653,
0.035346360791375846222037948478360,
0.044589751324764876608227299373280,
0.053481524690928087265343147239430,
0.062009567800670640285139230960803,
0.069854121318728258709520077099147,
0.076849680757720378894432777482659,
0.083080502823133021038289247286104,
0.088564443056211770647275443693774,
0.093126598170825321225486872747346,
0.096642726983623678505179907627589,
0.099173598721791959332393173484603,
0.100769845523875595044946662617570,
0.101330007014791549017374792767493
};
};
template <>
struct GKRule<GKOrder::GK41>
{
static constexpr std::size_t ng = 10;
static constexpr std::size_t ngk = 21;
static constexpr std::array<double, ngk> xgk = {
0.998859031588277663838315576545863,
0.993128599185094924786122388471320,
0.981507877450250259193342994720217,
0.963971927277913791267666131197277,
0.940822633831754753519982722212443,
0.912234428251325905867752441203298,
0.878276811252281976077442995113078,
0.839116971822218823394529061701521,
0.795041428837551198350638833272788,
0.746331906460150792614305070355642,
0.693237656334751384805490711845932, // NOLINT(modernize-use-std-numbers) -- a Kronrod abscissa, not std::numbers::ln2, despite the numerical coincidence
0.636053680726515025452836696226286,
0.575140446819710315342946036586425,
0.510867001950827098004364050955251,
0.443593175238725103199992213492640,
0.373706088715419560672548177024927,
0.301627868114913004320555356858592,
0.227785851141645078080496195368575,
0.152605465240922675505220241022678,
0.076526521133497333754640409398838,
0.000000000000000000000000000000000
};
static constexpr std::array<double, ng> wg = {
0.017614007139152118311861962351853,
0.040601429800386941331039952274932,
0.062672048334109063569506535187042,
0.083276741576704748724758143222046,
0.101930119817240435036750135480350,
0.118194531961518417312377377711382,
0.131688638449176626898494499748163,
0.142096109318382051329298325067165,
0.149172986472603746787828737001969,
0.152753387130725850698084331955098
};
static constexpr std::array<double, ngk> wgk = {
0.003073583718520531501218293246031,
0.008600269855642942198661787950102,
0.014626169256971252983787960308868,
0.020388373461266523598010231432755,
0.025882133604951158834505067096153,
0.031287306777032798958543119323801,
0.036600169758200798030557240707211,
0.041668873327973686263788305936895,
0.046434821867497674720231880926108,
0.050944573923728691932707670050345,
0.055195105348285994744832372419777,
0.059111400880639572374967220648594,
0.062653237554781168025870122174255,
0.065834597133618422111563556969398,
0.068648672928521619345623411885368,
0.071054423553444068305790361723210,
0.073030690332786667495189417658913,
0.074582875400499188986581418362488,
0.075704497684556674659542775376617,
0.076377867672080736705502835038061,
0.076600711917999656445049901530102
};
};
template <>
struct GKRule<GKOrder::GK51>
{
static constexpr std::size_t ng = 13;
static constexpr std::size_t ngk = 26;
static constexpr std::array<double, ngk> xgk = {
0.999262104992609834193457486540341,
0.995556969790498097908784946893902,
0.988035794534077247637331014577406,
0.976663921459517511498315386479594,
0.961614986425842512418130033660167,
0.942974571228974339414011169658471,
0.920747115281701561746346084546331,
0.894991997878275368851042006782805,
0.865847065293275595448996969588340,
0.833442628760834001421021108693570,
0.797873797998500059410410904994307,
0.759259263037357630577282865204361,
0.717766406813084388186654079773298,
0.673566368473468364485120633247622,
0.626810099010317412788122681624518,
0.577662930241222967723689841612654, // NOLINT(modernize-use-std-numbers) -- a Kronrod abscissa, not std::numbers::inv_sqrt3, despite the numerical coincidence
0.526325284334719182599623778158010,
0.473002731445714960522182115009192,
0.417885382193037748851814394594572,
0.361172305809387837735821730127641,
0.303089538931107830167478909980339,
0.243866883720988432045190362797452,
0.183718939421048892015969888759528,
0.122864692610710396387359818808037,
0.061544483005685078886546392366797,
0.000000000000000000000000000000000
};
static constexpr std::array<double, ng> wg = {
0.011393798501026287947902964113235,
0.026354986615032137261901815295299,
0.040939156701306312655623487711646,
0.054904695975835191925936891540473,
0.068038333812356917207187185656708,
0.080140700335001018013234959669111,
0.091028261982963649811497220702892,
0.100535949067050644202206890392686,
0.108519624474263653116093957050117,
0.114858259145711648339325545869556,
0.119455763535784772228178126512901,
0.122242442990310041688959518945852,
0.123176053726715451203902873079050
};
static constexpr std::array<double, ngk> wgk = {
0.001987383892330315926507851882843,
0.005561932135356713758040236901066,
0.009473973386174151607207710523655,
0.013236229195571674813656405846976,
0.016847817709128298231516667536336,
0.020435371145882835456568292235939,
0.024009945606953216220092489164881,
0.027475317587851737802948455517811,
0.030792300167387488891109020215229,
0.034002130274329337836748795229551,
0.037116271483415543560330625367620,
0.040083825504032382074839284467076,
0.042872845020170049476895792439495,
0.045502913049921788909870584752660,
0.047982537138836713906392255756915,
0.050277679080715671963325259433440,
0.052362885806407475864366712137873,
0.054251129888545490144543370459876,
0.055950811220412317308240686382747,
0.057437116361567832853582693939506,
0.058689680022394207961974175856788,
0.059720340324174059979099291932562,
0.060539455376045862945360267517565,
0.061128509717053048305859030416293,
0.061471189871425316661544131965264,
0.061580818067832935078759824240066
};
};
template <>
struct GKRule<GKOrder::GK61>
{
static constexpr std::size_t ng = 15;
static constexpr std::size_t ngk = 31;
static constexpr std::array<double, ngk> xgk = {
0.999484410050490637571325895705811,
0.996893484074649540271630050918695,
0.991630996870404594858628366109486,
0.983668123279747209970032581605663,
0.973116322501126268374693868423707,
0.960021864968307512216871025581798,
0.944374444748559979415831324037439,
0.926200047429274325879324277080474,
0.905573307699907798546522558925958,
0.882560535792052681543116462530226,
0.857205233546061098958658510658944,
0.829565762382768397442898119732502,
0.799727835821839083013668942322683,
0.767777432104826194917977340974503,
0.733790062453226804726171131369528,
0.697850494793315796932292388026640,
0.660061064126626961370053668149271,
0.620526182989242861140477556431189,
0.579345235826361691756024932172540,
0.536624148142019899264169793311073,
0.492480467861778574993693061207709,
0.447033769538089176780609900322854,
0.400401254830394392535476211542661,
0.352704725530878113471037207089374,
0.304073202273625077372677107199257,
0.254636926167889846439805129817805,
0.204525116682309891438957671002025,
0.153869913608583546963794672743256,
0.102806937966737030147096751318001,
0.051471842555317695833025213166723,
0.000000000000000000000000000000000
};
static constexpr std::array<double, ng> wg = {
0.007968192496166605615465883474674,
0.018466468311090959142302131912047,
0.028784707883323369349719179611292,
0.038799192569627049596801936446348,
0.048402672830594052902938140422808,
0.057493156217619066481721689402056,
0.065974229882180495128128515115962,
0.073755974737705206268243850022191,
0.080755895229420215354694938460530,
0.086899787201082979802387530715126,
0.092122522237786128717632707087619,
0.096368737174644259639468626351810,
0.099593420586795267062780282103569,
0.101762389748405504596428952168554,
0.102852652893558840341285636705415
};
static constexpr std::array<double, ngk> wgk = {
0.001389013698677007624551591226760,
0.003890461127099884051267201844516,
0.006630703915931292173319826369750,
0.009273279659517763428441146892024,
0.011823015253496341742232898853251,
0.014369729507045804812451432443580,
0.016920889189053272627572289420322,
0.019414141193942381173408951050128,
0.021828035821609192297167485738339,
0.024191162078080601365686370725232,
0.026509954882333101610601709335075,
0.028754048765041292843978785354334,
0.030907257562387762472884252943092,
0.032981447057483726031814191016854,
0.034979338028060024137499670731468,
0.036882364651821229223911065617136,
0.038678945624727592950348651532281,
0.040374538951535959111995279752468,
0.041969810215164246147147541285970,
0.043452539701356069316831728117073,
0.044814800133162663192355551616723,
0.046059238271006988116271735559374,
0.047185546569299153945261478181099,
0.048185861757087129140779492298305,
0.049055434555029778887528165367238,
0.049795683427074206357811569379942,
0.050405921402782346840893085653585,
0.050881795898749606492297473049805,
0.051221547849258772170656282604944,
0.051426128537459025933862879215781,
0.051494729429451567558340433647099
};
};
} // namespace utils
#endif // GKINTEGRATORDATA_HPP