gauss_legendre.hh

Go to the documentation of this file.

#pragma once
 
#include <DiFfRG/common/cuda_prefix.hh>
 
namespace DiFfRG
{
  template <size_t N, typename ctype> struct GLQuadrature {
    static_assert(N, "Requested quadrature order has not been implemented.");
  };
 
  template <typename ctype> struct GLQuadrature<1, ctype> {
    __constant__ ctype x[1];
    __constant__ ctype w[1];
    consteval GLQuadrature()
    {
      x[0] = 0.50000000000000000000000000000000;
      w[0] = 1.0000000000000000000000000000000;
    }
  };
 
  template <typename ctype> struct GLQuadrature<2, ctype> {
    __constant__ ctype x[2];
    __constant__ ctype w[2];
    consteval GLQuadrature()
    {
      x[0] = 0.2113248654051871177454256097490;
      x[1] = 0.7886751345948128822545743902510;
      w[0] = 0.5000000000000000000000000000000;
      w[1] = 0.5000000000000000000000000000000;
    }
  };
 
  template <typename ctype> struct GLQuadrature<3, ctype> {
    __constant__ ctype x[3];
    __constant__ ctype w[3];
    consteval GLQuadrature()
    {
      x[0] = 0.1127016653792583114820734600218;
      x[1] = 0.500000000000000000000000000000;
      x[2] = 0.887298334620741688517926539978;
      w[0] = 0.2777777777777777777777777777778;
      w[1] = 0.4444444444444444444444444444444;
      w[2] = 0.2777777777777777777777777777778;
    }
  };
 
  template <typename ctype> struct GLQuadrature<4, ctype> {
    __constant__ ctype x[4];
    __constant__ ctype w[4];
    consteval GLQuadrature()
    {
      x[0] = 0.0694318442029737123880267555536;
      x[1] = 0.330009478207571867598667120448;
      x[2] = 0.669990521792428132401332879552;
      x[3] = 0.930568155797026287611973244446;
      w[0] = 0.1739274225687269286865319746110;
      w[1] = 0.3260725774312730713134680253890;
      w[2] = 0.3260725774312730713134680253890;
      w[3] = 0.1739274225687269286865319746110;
    }
  };
 
  template <typename ctype> struct GLQuadrature<5, ctype> {
    __constant__ ctype x[5];
    __constant__ ctype w[5];
    consteval GLQuadrature()
    {
      x[0] = 0.0469100770306680036011865608503;
      x[1] = 0.230765344947158454481842789650;
      x[2] = 0.500000000000000000000000000000;
      x[3] = 0.769234655052841545518157210350;
      x[4] = 0.953089922969331996398813439150;
      w[0] = 0.1184634425280945437571320203600;
      w[1] = 0.2393143352496832340206457574178;
      w[2] = 0.2844444444444444444444444444444;
      w[3] = 0.2393143352496832340206457574178;
      w[4] = 0.1184634425280945437571320203600;
    }
  };
 
  template <typename ctype> struct GLQuadrature<6, ctype> {
    __constant__ ctype x[6];
    __constant__ ctype w[6];
    consteval GLQuadrature()
    {
      x[0] = 0.0337652428984239860938492227530;
      x[1] = 0.169395306766867743169300202490;
      x[2] = 0.380690406958401545684749139160;
      x[3] = 0.619309593041598454315250860840;
      x[4] = 0.830604693233132256830699797510;
      x[5] = 0.966234757101576013906150777247;
      w[0] = 0.0856622461895851725201480710864;
      w[1] = 0.1803807865240693037849167569189;
      w[2] = 0.2339569672863455236949351719948;
      w[3] = 0.2339569672863455236949351719948;
      w[4] = 0.1803807865240693037849167569189;
      w[5] = 0.0856622461895851725201480710864;
    }
  };
 
  template <typename ctype> struct GLQuadrature<7, ctype> {
    __constant__ ctype x[7];
    __constant__ ctype w[7];
    consteval GLQuadrature()
    {
      x[0] = 0.0254460438286207377369051579761;
      x[1] = 0.129234407200302780068067613360;
      x[2] = 0.297077424311301416546696793962;
      x[3] = 0.500000000000000000000000000000;
      x[4] = 0.702922575688698583453303206038;
      x[5] = 0.87076559279969721993193238664;
      x[6] = 0.97455395617137926226309484202;
      w[0] = 0.06474248308443484663530571633954;
      w[1] = 0.1398526957446383339507338857119;
      w[2] = 0.1909150252525594724751848877445;
      w[3] = 0.2089795918367346938775510204082;
      w[4] = 0.1909150252525594724751848877445;
      w[5] = 0.1398526957446383339507338857119;
      w[6] = 0.06474248308443484663530571633954;
    }
  };
 
  template <typename ctype> struct GLQuadrature<8, ctype> {
    __constant__ ctype x[8];
    __constant__ ctype w[8];
    consteval GLQuadrature()
    {
      x[0] = 0.0198550717512318841582195657153;
      x[1] = 0.101666761293186630204223031762;
      x[2] = 0.237233795041835507091130475405;
      x[3] = 0.408282678752175097530261928820;
      x[4] = 0.591717321247824902469738071180;
      x[5] = 0.76276620495816449290886952459;
      x[6] = 0.89833323870681336979577696824;
      x[7] = 0.98014492824876811584178043428;
      w[0] = 0.05061426814518812957626567715498;
      w[1] = 0.1111905172266872352721779972131;
      w[2] = 0.1568533229389436436689811009933;
      w[3] = 0.1813418916891809914825752246386;
      w[4] = 0.1813418916891809914825752246386;
      w[5] = 0.1568533229389436436689811009933;
      w[6] = 0.1111905172266872352721779972131;
      w[7] = 0.05061426814518812957626567715498;
    }
  };
 
  template <typename ctype> struct GLQuadrature<9, ctype> {
    __constant__ ctype x[9];
    __constant__ ctype w[9];
    consteval GLQuadrature()
    {
      x[0] = 0.0159198802461869550822118985482;
      x[1] = 0.081984446336682102850285105965;
      x[2] = 0.193314283649704801345648980329;
      x[3] = 0.337873288298095535480730992678;
      x[4] = 0.500000000000000000000000000000;
      x[5] = 0.66212671170190446451926900732;
      x[6] = 0.80668571635029519865435101967;
      x[7] = 0.91801555366331789714971489403;
      x[8] = 0.98408011975381304491778810145;
      w[0] = 0.04063719418078720598594607905526;
      w[1] = 0.09032408034742870202923601562146;
      w[2] = 0.1303053482014677311593714347093;
      w[3] = 0.1561735385200014200343152032922;
      w[4] = 0.1651196775006298815822625346435;
      w[5] = 0.1561735385200014200343152032922;
      w[6] = 0.1303053482014677311593714347093;
      w[7] = 0.09032408034742870202923601562146;
      w[8] = 0.04063719418078720598594607905526;
    }
  };
 
  template <typename ctype> struct GLQuadrature<10, ctype> {
    __constant__ ctype x[10];
    __constant__ ctype w[10];
    consteval GLQuadrature()
    {
      x[0] = 0.0130467357414141399610179939578;
      x[1] = 0.067468316655507744633951655788;
      x[2] = 0.160295215850487796882836317443;
      x[3] = 0.283302302935376404600367028417;
      x[4] = 0.425562830509184394557586999435;
      x[5] = 0.57443716949081560544241300056;
      x[6] = 0.71669769706462359539963297158;
      x[7] = 0.83970478414951220311716368256;
      x[8] = 0.93253168334449225536604834421;
      x[9] = 0.98695326425858586003898200604;
      w[0] = 0.03333567215434406879678440494667;
      w[1] = 0.0747256745752902965728881698288;
      w[2] = 0.1095431812579910219977674671141;
      w[3] = 0.1346333596549981775456134607847;
      w[4] = 0.1477621123573764350869464973257;
      w[5] = 0.1477621123573764350869464973257;
      w[6] = 0.1346333596549981775456134607847;
      w[7] = 0.1095431812579910219977674671141;
      w[8] = 0.0747256745752902965728881698288;
      w[9] = 0.03333567215434406879678440494667;
    }
  };
 
  template <typename ctype> struct GLQuadrature<11, ctype> {
    __constant__ ctype x[11];
    __constant__ ctype w[11];
    consteval GLQuadrature()
    {
      x[0] = 0.0108856709269715035980309994386;
      x[1] = 0.056468700115952350462421115348;
      x[2] = 0.134923997212975337953291873984;
      x[3] = 0.240451935396594092037137165271;
      x[4] = 0.365228422023827513834234007300;
      x[5] = 0.50000000000000000000000000000;
      x[6] = 0.63477157797617248616576599270;
      x[7] = 0.75954806460340590796286283473;
      x[8] = 0.86507600278702466204670812602;
      x[9] = 0.94353129988404764953757888465;
      x[10] = 0.98911432907302849640196900056;
      w[0] = 0.02783428355808683324137686022127;
      w[1] = 0.06279018473245231231734714961197;
      w[2] = 0.09314510546386712571304882071583;
      w[3] = 0.1165968822959952399592618524216;
      w[4] = 0.1314022722551233310903444349453;
      w[5] = 0.1364625433889503153572417641682;
      w[6] = 0.1314022722551233310903444349453;
      w[7] = 0.1165968822959952399592618524216;
      w[8] = 0.09314510546386712571304882071583;
      w[9] = 0.06279018473245231231734714961197;
      w[10] = 0.02783428355808683324137686022127;
    }
  };
 
  template <typename ctype> struct GLQuadrature<12, ctype> {
    __constant__ ctype x[12];
    __constant__ ctype w[12];
    consteval GLQuadrature()
    {
      x[0] = 0.0092196828766403746547254549254;
      x[1] = 0.047941371814762571660767066940;
      x[2] = 0.115048662902847656481553083394;
      x[3] = 0.206341022856691276351648790530;
      x[4] = 0.31608425050090990312365423168;
      x[5] = 0.43738329574426554226377931527;
      x[6] = 0.56261670425573445773622068473;
      x[7] = 0.68391574949909009687634576832;
      x[8] = 0.79365897714330872364835120947;
      x[9] = 0.88495133709715234351844691661;
      x[10] = 0.95205862818523742833923293306;
      x[11] = 0.99078031712335962534527454507;
      w[0] = 0.02358766819325591359730798074251;
      w[1] = 0.05346966299765921548012735909700;
      w[2] = 0.0800391642716731131673262647717;
      w[3] = 0.1015837133615329608745322279049;
      w[4] = 0.1167462682691774043804249494624;
      w[5] = 0.1245735229067013925002812180215;
      w[6] = 0.1245735229067013925002812180215;
      w[7] = 0.1167462682691774043804249494624;
      w[8] = 0.1015837133615329608745322279049;
      w[9] = 0.0800391642716731131673262647717;
      w[10] = 0.05346966299765921548012735909700;
      w[11] = 0.02358766819325591359730798074251;
    }
  };
 
  template <typename ctype> struct GLQuadrature<13, ctype> {
    __constant__ ctype x[13];
    __constant__ ctype w[13];
    consteval GLQuadrature()
    {
      x[0] = 0.0079084726407059252635852755964;
      x[1] = 0.041200800388511017396726081750;
      x[2] = 0.099210954633345043602896755209;
      x[3] = 0.178825330279829889678007696502;
      x[4] = 0.27575362448177657356104357394;
      x[5] = 0.38477084202243260296723593945;
      x[6] = 0.50000000000000000000000000000;
      x[7] = 0.61522915797756739703276406055;
      x[8] = 0.72424637551822342643895642606;
      x[9] = 0.82117466972017011032199230350;
      x[10] = 0.90078904536665495639710324479;
      x[11] = 0.95879919961148898260327391825;
      x[12] = 0.99209152735929407473641472440;
      w[0] = 0.02024200238265793976001079610049;
      w[1] = 0.04606074991886422395721088797690;
      w[2] = 0.06943675510989361923180088843444;
      w[3] = 0.08907299038097286914002334599805;
      w[4] = 0.10390802376844425115626160965303;
      w[5] = 0.11314159013144861920604509301989;
      w[6] = 0.1162757766154369550972947576344;
      w[7] = 0.11314159013144861920604509301989;
      w[8] = 0.10390802376844425115626160965303;
      w[9] = 0.08907299038097286914002334599805;
      w[10] = 0.06943675510989361923180088843444;
      w[11] = 0.04606074991886422395721088797690;
      w[12] = 0.02024200238265793976001079610049;
    }
  };
 
  template <typename ctype> struct GLQuadrature<14, ctype> {
    __constant__ ctype x[14];
    __constant__ ctype w[14];
    consteval GLQuadrature()
    {
      x[0] = 0.0068580956515938305792013666480;
      x[1] = 0.035782558168213241331804430311;
      x[2] = 0.086399342465117503405102628675;
      x[3] = 0.156353547594157264925990098490;
      x[4] = 0.24237568182092295401735464072;
      x[5] = 0.34044381553605511978216408792;
      x[6] = 0.44597252564632816896687767489;
      x[7] = 0.55402747435367183103312232511;
      x[8] = 0.65955618446394488021783591208;
      x[9] = 0.75762431817907704598264535928;
      x[10] = 0.84364645240584273507400990151;
      x[11] = 0.91360065753488249659489737133;
      x[12] = 0.96421744183178675866819556969;
      x[13] = 0.99314190434840616942079863335;
      w[0] = 0.01755973016587593151591643806910;
      w[1] = 0.04007904357988010490281663853143;
      w[2] = 0.06075928534395159234470740453624;
      w[3] = 0.0786015835790967672848009693119;
      w[4] = 0.0927691987389689068708582950626;
      w[5] = 0.1025992318606478019829620328306;
      w[6] = 0.1076319267315788950979382216581;
      w[7] = 0.1076319267315788950979382216581;
      w[8] = 0.1025992318606478019829620328306;
      w[9] = 0.0927691987389689068708582950626;
      w[10] = 0.0786015835790967672848009693119;
      w[11] = 0.06075928534395159234470740453624;
      w[12] = 0.04007904357988010490281663853143;
      w[13] = 0.01755973016587593151591643806910;
    }
  };
 
  template <typename ctype> struct GLQuadrature<15, ctype> {
    __constant__ ctype x[15];
    __constant__ ctype w[15];
    consteval GLQuadrature()
    {
      x[0] = 0.0060037409897572857552171407067;
      x[1] = 0.031363303799647047846120526145;
      x[2] = 0.075896708294786391899675839613;
      x[3] = 0.137791134319914976291906972693;
      x[4] = 0.21451391369573057623138663137;
      x[5] = 0.30292432646121831505139631451;
      x[6] = 0.39940295300128273884968584830;
      x[7] = 0.50000000000000000000000000000;
      x[8] = 0.60059704699871726115031415170;
      x[9] = 0.69707567353878168494860368549;
      x[10] = 0.78548608630426942376861336863;
      x[11] = 0.86220886568008502370809302731;
      x[12] = 0.92410329170521360810032416039;
      x[13] = 0.96863669620035295215387947386;
      x[14] = 0.99399625901024271424478285929;
      w[0] = 0.01537662099805863417731419678860;
      w[1] = 0.03518302374405406235463370822533;
      w[2] = 0.05357961023358596750593477334293;
      w[3] = 0.06978533896307715722390239725551;
      w[4] = 0.08313460290849696677660043024060;
      w[5] = 0.09308050000778110551340028093321;
      w[6] = 0.09921574266355578822805916322192;
      w[7] = 0.1012891209627806364403100999838;
      w[8] = 0.09921574266355578822805916322192;
      w[9] = 0.09308050000778110551340028093321;
      w[10] = 0.08313460290849696677660043024060;
      w[11] = 0.06978533896307715722390239725551;
      w[12] = 0.05357961023358596750593477334293;
      w[13] = 0.03518302374405406235463370822533;
      w[14] = 0.01537662099805863417731419678860;
    }
  };
 
  template <typename ctype> struct GLQuadrature<16, ctype> {
    __constant__ ctype x[16];
    __constant__ ctype w[16];
    consteval GLQuadrature()
    {
      x[0] = 0.0052995325041750337019229132748;
      x[1] = 0.027712488463383711961005792233;
      x[2] = 0.067184398806084128059766051144;
      x[3] = 0.122297795822498483052449402576;
      x[4] = 0.19106187779867812577666411798;
      x[5] = 0.27099161117138630682879027851;
      x[6] = 0.35919822461037054338476974927;
      x[7] = 0.45249374508118127990734033229;
      x[8] = 0.54750625491881872009265966771;
      x[9] = 0.64080177538962945661523025073;
      x[10] = 0.72900838882861369317120972149;
      x[11] = 0.80893812220132187422333588202;
      x[12] = 0.87770220417750151694755059742;
      x[13] = 0.93281560119391587194023394886;
      x[14] = 0.97228751153661628803899420777;
      x[15] = 0.99470046749582496629807708673;
      w[0] = 0.01357622970587704742589028622801;
      w[1] = 0.03112676196932394643142191849719;
      w[2] = 0.04757925584124639240496255380112;
      w[3] = 0.06231448562776693602623814109601;
      w[4] = 0.0747979944082883660407508652737;
      w[5] = 0.0845782596975012690946560395152;
      w[6] = 0.0913017075224617944333818339846;
      w[7] = 0.0947253052275342481426983616041;
      w[8] = 0.0947253052275342481426983616041;
      w[9] = 0.0913017075224617944333818339846;
      w[10] = 0.0845782596975012690946560395152;
      w[11] = 0.0747979944082883660407508652737;
      w[12] = 0.06231448562776693602623814109601;
      w[13] = 0.04757925584124639240496255380112;
      w[14] = 0.03112676196932394643142191849719;
      w[15] = 0.01357622970587704742589028622801;
    }
  };
 
  template <typename ctype> struct GLQuadrature<20, ctype> {
    __constant__ ctype x[20];
    __constant__ ctype w[20];
    consteval GLQuadrature()
    {
      x[0] = 0.0034357004074525376069388057643;
      x[1] = 0.018014036361043104366166934401;
      x[2] = 0.043882785874337047066123779398;
      x[3] = 0.080441514088890588302735469149;
      x[4] = 0.12683404676992460369284746482;
      x[5] = 0.18197315963674248727358165189;
      x[6] = 0.24456649902458645099781797452;
      x[7] = 0.31314695564229021966372591149;
      x[8] = 0.38610707442917746095975190232;
      x[9] = 0.46173673943325133312267979530;
      x[10] = 0.53826326056674866687732020470;
      x[11] = 0.61389292557082253904024809768;
      x[12] = 0.68685304435770978033627408851;
      x[13] = 0.75543350097541354900218202548;
      x[14] = 0.81802684036325751272641834811;
      x[15] = 0.87316595323007539630715253518;
      x[16] = 0.91955848591110941169726453085;
      x[17] = 0.95611721412566295293387622060;
      x[18] = 0.98198596363895689563383306560;
      x[19] = 0.99656429959254746239306119424;
      w[0] = 0.00880700356957605915593098117593;
      w[1] = 0.02030071490019347066551997613747;
      w[2] = 0.03133602416705453178475326759352;
      w[3] = 0.04163837078835237436237907161102;
      w[4] = 0.05096505990862021751837506774017;
      w[5] = 0.05909726598075920865618868885569;
      w[6] = 0.0658443192245883134492472498741;
      w[7] = 0.0710480546591910256646491625336;
      w[8] = 0.0745864932363018733939143685010;
      w[9] = 0.0763766935653629253490421659775;
      w[10] = 0.0763766935653629253490421659775;
      w[11] = 0.0745864932363018733939143685010;
      w[12] = 0.0710480546591910256646491625336;
      w[13] = 0.0658443192245883134492472498741;
      w[14] = 0.05909726598075920865618868885569;
      w[15] = 0.05096505990862021751837506774017;
      w[16] = 0.04163837078835237436237907161102;
      w[17] = 0.03133602416705453178475326759352;
      w[18] = 0.02030071490019347066551997613747;
      w[19] = 0.00880700356957605915593098117593;
    }
  };
 
  template <typename ctype> struct GLQuadrature<24, ctype> {
    __constant__ ctype x[24];
    __constant__ ctype w[24];
    consteval GLQuadrature()
    {
      x[0] = 0.0024063900014893199100012951496;
      x[1] = 0.012635722014345250900804003496;
      x[2] = 0.030862723998633620738175499146;
      x[3] = 0.056792236497799482893422829009;
      x[4] = 0.08999900701304853902302506367;
      x[5] = 0.12993790421072281787808594845;
      x[6] = 0.17595317403151221537375210654;
      x[7] = 0.22728926430558023217081219139;
      x[8] = 0.28310324618697743075645788404;
      x[9] = 0.34247866015191831280660335434;
      x[10] = 0.40444056626319184542068008962;
      x[11] = 0.46797155356869718695747845869;
      x[12] = 0.53202844643130281304252154131;
      x[13] = 0.59555943373680815457931991038;
      x[14] = 0.65752133984808168719339664566;
      x[15] = 0.71689675381302256924354211596;
      x[16] = 0.77271073569441976782918780861;
      x[17] = 0.8240468259684877846262478935;
      x[18] = 0.8700620957892771821219140515;
      x[19] = 0.9100009929869514609769749363;
      x[20] = 0.9432077635022005171065771710;
      x[21] = 0.9691372760013663792618245009;
      x[22] = 0.9873642779856547490991959965;
      x[23] = 0.9975936099985106800899987049;
      w[0] = 0.006170614899993599773402833535019;
      w[1] = 0.01426569431446683159065390797594;
      w[2] = 0.02213871940870990308430137410567;
      w[3] = 0.02964929245771839037318387925005;
      w[4] = 0.03667324070554015286701680762656;
      w[5] = 0.04309508076597663795859260149187;
      w[6] = 0.04880932605205694413494033223212;
      w[7] = 0.05372213505798281739128867122330;
      w[8] = 0.05775283402686280067667224195339;
      w[9] = 0.0608352364639016956022315767381;
      w[10] = 0.0629187281734141480606876912556;
      w[11] = 0.0639690976733760784870280826123;
      w[12] = 0.0639690976733760784870280826123;
      w[13] = 0.0629187281734141480606876912556;
      w[14] = 0.0608352364639016956022315767381;
      w[15] = 0.05775283402686280067667224195339;
      w[16] = 0.05372213505798281739128867122330;
      w[17] = 0.04880932605205694413494033223212;
      w[18] = 0.04309508076597663795859260149187;
      w[19] = 0.03667324070554015286701680762656;
      w[20] = 0.02964929245771839037318387925005;
      w[21] = 0.02213871940870990308430137410567;
      w[22] = 0.01426569431446683159065390797594;
      w[23] = 0.006170614899993599773402833535019;
    }
  };
 
  template <typename ctype> struct GLQuadrature<32, ctype> {
    __constant__ ctype x[32];
    __constant__ ctype w[32];
    consteval GLQuadrature()
    {
      x[0] = 0.0013680690752592182275094356675;
      x[1] = 0.007194244227365832299912477685;
      x[2] = 0.017618872206246784613094035941;
      x[3] = 0.032546962031130155414540432582;
      x[4] = 0.05183942211697393801734637814;
      x[5] = 0.07531619313371501493315349752;
      x[6] = 0.10275810201602879651845135051;
      x[7] = 0.13390894062985515980628666745;
      x[8] = 0.16847786653489239951244241567;
      x[9] = 0.20614212137961883547962726180;
      x[10] = 0.24655004553388530498812626281;
      x[11] = 0.28932436193468232731794028191;
      x[12] = 0.33406569885893617511004159713;
      x[13] = 0.38035631887393146272769839542;
      x[14] = 0.42776401920860175325740681320;
      x[15] = 0.4758461671561308418825937148;
      x[16] = 0.5241538328438691581174062852;
      x[17] = 0.5722359807913982467425931868;
      x[18] = 0.6196436811260685372723016046;
      x[19] = 0.6659343011410638248899584029;
      x[20] = 0.7106756380653176726820597181;
      x[21] = 0.7534499544661146950118737372;
      x[22] = 0.7938578786203811645203727382;
      x[23] = 0.8315221334651076004875575843;
      x[24] = 0.8660910593701448401937133325;
      x[25] = 0.8972418979839712034815486495;
      x[26] = 0.9246838068662849850668465025;
      x[27] = 0.9481605778830260619826536219;
      x[28] = 0.9674530379688698445854595674;
      x[29] = 0.9823811277937532153869059641;
      x[30] = 0.9928057557726341677000875223;
      x[31] = 0.9986319309247407817724905643;
      w[0] = 0.003509305004735048300203531869427;
      w[1] = 0.00813719736545283530258528110319;
      w[2] = 0.01269603265463102972787629489461;
      w[3] = 0.01713693145651071655134386612619;
      w[4] = 0.02141794901111334032843932330306;
      w[5] = 0.02549902963118808809808162234476;
      w[6] = 0.02934204673926777357264181865009;
      w[7] = 0.03291111138818092341882503185347;
      w[8] = 0.03617289705442425311269967823924;
      w[9] = 0.03909694789353515323587045941415;
      w[10] = 0.04165596211347337761109953730217;
      w[11] = 0.04382604650220190557138573137590;
      w[12] = 0.04558693934788194235643428855582;
      w[13] = 0.04692219954040228281959011883406;
      w[14] = 0.0478193600396374297095410011021;
      w[15] = 0.0482700442573639002833824150318;
      w[16] = 0.0482700442573639002833824150318;
      w[17] = 0.0478193600396374297095410011021;
      w[18] = 0.04692219954040228281959011883406;
      w[19] = 0.04558693934788194235643428855582;
      w[20] = 0.04382604650220190557138573137590;
      w[21] = 0.04165596211347337761109953730217;
      w[22] = 0.03909694789353515323587045941415;
      w[23] = 0.03617289705442425311269967823924;
      w[24] = 0.03291111138818092341882503185347;
      w[25] = 0.02934204673926777357264181865009;
      w[26] = 0.02549902963118808809808162234476;
      w[27] = 0.02141794901111334032843932330306;
      w[28] = 0.01713693145651071655134386612619;
      w[29] = 0.01269603265463102972787629489461;
      w[30] = 0.00813719736545283530258528110319;
      w[31] = 0.003509305004735048300203531869427;
    }
  };
 
  template <typename ctype> struct GLQuadrature<48, ctype> {
    __constant__ ctype x[48];
    __constant__ ctype w[48];
    consteval GLQuadrature()
    {
      x[0] = 0.0006144963737869406997292542184;
      x[1] = 0.003234913866824621226035624575;
      x[2] = 0.007937708138586571127708199987;
      x[3] = 0.014704203726876374769294008100;
      x[4] = 0.02350614841978456963851966699;
      x[5] = 0.03430665464672283344291280995;
      x[6] = 0.04706043164221516358896258216;
      x[7] = 0.06171398986287605704715322260;
      x[8] = 0.07820586918780323464445507774;
      x[9] = 0.09646689798527868645872347849;
      x[10] = 0.11642048374212983037307228124;
      x[11] = 0.13798293453809267266275888325;
      x[12] = 0.16106381018366804739407435966;
      x[13] = 0.18556630161174318800241753347;
      x[14] = 0.2113876369580136480910953807;
      x[15] = 0.2384195126388834831608870654;
      x[16] = 0.2665485476245207977275355692;
      x[17] = 0.2956567590046416350418872521;
      x[18] = 0.3256220568539196309200910314;
      x[19] = 0.3563187563222722116320567693;
      x[20] = 0.3876181048026554693875672799;
      x[21] = 0.4193888219655541409717813046;
      x[22] = 0.4514976503952686505349730221;
      x[23] = 0.4838099145185653189833388784;
      x[24] = 0.5161900854814346810166611216;
      x[25] = 0.5485023496047313494650269779;
      x[26] = 0.5806111780344458590282186954;
      x[27] = 0.6123818951973445306124327201;
      x[28] = 0.6436812436777277883679432307;
      x[29] = 0.6743779431460803690799089686;
      x[30] = 0.7043432409953583649581127479;
      x[31] = 0.7334514523754792022724644308;
      x[32] = 0.7615804873611165168391129346;
      x[33] = 0.7886123630419863519089046193;
      x[34] = 0.8144336983882568119975824665;
      x[35] = 0.8389361898163319526059256403;
      x[36] = 0.8620170654619073273372411167;
      x[37] = 0.8835795162578701696269277188;
      x[38] = 0.9035331020147213135412765215;
      x[39] = 0.9217941308121967653555449223;
      x[40] = 0.9382860101371239429528467774;
      x[41] = 0.9529395683577848364110374178;
      x[42] = 0.9656933453532771665570871901;
      x[43] = 0.9764938515802154303614803330;
      x[44] = 0.9852957962731236252307059919;
      x[45] = 0.9920622918614134288722918000;
      x[46] = 0.9967650861331753787739643754;
      x[47] = 0.9993855036262130593002707458;
      w[0] = 0.001576673026152919316338655771946;
      w[1] = 0.003663776950638131051191989810893;
      w[2] = 0.005738617289617269744796333804546;
      w[3] = 0.00778965786147192436408847791723;
      w[4] = 0.00980808022867776390723035982611;
      w[5] = 0.01178538041966218957025965068922;
      w[6] = 0.01371325485417847410003691813125;
      w[7] = 0.01558361391639904445103287842318;
      w[8] = 0.01738861128238521944627429298190;
      w[9] = 0.01912067553291535315860862826186;
      w[10] = 0.02077254147173237460702941118053;
      w[11] = 0.02233728042834714020972429356293;
      w[12] = 0.02380832924624523741295331173946;
      w[13] = 0.02517951777692723747890380954393;
      w[14] = 0.02644509474259683354775252813235;
      w[15] = 0.02759975184999208143410174759582;
      w[16] = 0.02863864605020160785257511734235;
      w[17] = 0.02955741984919781787323740871676;
      w[18] = 0.03035221958294694002648461601391;
      w[19] = 0.03101971157994633195209889206880;
      w[20] = 0.03155709614312701282856301137512;
      w[21] = 0.03196211929232409331195310091276;
      w[22] = 0.03223308221797504110325209682885;
      w[23] = 0.0323688484063419612515124693683;
      w[24] = 0.0323688484063419612515124693683;
      w[25] = 0.03223308221797504110325209682885;
      w[26] = 0.03196211929232409331195310091276;
      w[27] = 0.03155709614312701282856301137512;
      w[28] = 0.03101971157994633195209889206880;
      w[29] = 0.03035221958294694002648461601391;
      w[30] = 0.02955741984919781787323740871676;
      w[31] = 0.02863864605020160785257511734235;
      w[32] = 0.02759975184999208143410174759582;
      w[33] = 0.02644509474259683354775252813235;
      w[34] = 0.02517951777692723747890380954393;
      w[35] = 0.02380832924624523741295331173946;
      w[36] = 0.02233728042834714020972429356293;
      w[37] = 0.02077254147173237460702941118053;
      w[38] = 0.01912067553291535315860862826186;
      w[39] = 0.01738861128238521944627429298190;
      w[40] = 0.01558361391639904445103287842318;
      w[41] = 0.01371325485417847410003691813125;
      w[42] = 0.01178538041966218957025965068922;
      w[43] = 0.00980808022867776390723035982611;
      w[44] = 0.00778965786147192436408847791723;
      w[45] = 0.005738617289617269744796333804546;
      w[46] = 0.003663776950638131051191989810893;
      w[47] = 0.001576673026152919316338655771946;
    }
  };
 
  template <typename ctype> struct GLQuadrature<64, ctype> {
    __constant__ ctype x[64];
    __constant__ ctype w[64];
    consteval GLQuadrature()
    {
      x[0] = 0.0003474791321139302715471878272;
      x[1] = 0.001829941614022360326537749662;
      x[2] = 0.004493314261627839630308808278;
      x[3] = 0.008331873057687021534350348922;
      x[4] = 0.01333658610504451812907324632;
      x[5] = 0.01949560017397314054069293905;
      x[6] = 0.02679431257079859196875925433;
      x[7] = 0.03521541393403021208925492272;
      x[8] = 0.04473893146074859712180966600;
      x[9] = 0.05534227700244294707329798086;
      x[10] = 0.06700030092295359011960830746;
      x[11] = 0.07968535187370981862415422765;
      x[12] = 0.09336734243860122012903833096;
      x[13] = 0.10801382052832929619488973739;
      x[14] = 0.1235900463697340516940681126;
      x[15] = 0.1400590749141945865755298911;
      x[16] = 0.1573818434728833787182208145;
      x[17] = 0.1755172643726713300711193840;
      x[18] = 0.1944223224138033748755735145;
      x[19] = 0.2140521768986829828580609417;
      x[20] = 0.2343602679900527271709930482;
      x[21] = 0.2552984271464735212607368465;
      x[22] = 0.2768169913732679560075261426;
      x[23] = 0.2988649210180041981521166144;
      x[24] = 0.3213899208311659420247786925;
      x[25] = 0.3443385640048945219212436507;
      x[26] = 0.3676564188956162918130179137;
      x[27] = 0.3912881781299964579251756255;
      x[28] = 0.4151777897880035909813431851;
      x[29] = 0.4392685903519397227648117683;
      x[30] = 0.4635034391061004802752285290;
      x[31] = 0.4878248536682877837455220786;
      x[32] = 0.5121751463317122162544779214;
      x[33] = 0.5364965608938995197247714710;
      x[34] = 0.5607314096480602772351882317;
      x[35] = 0.5848222102119964090186568149;
      x[36] = 0.6087118218700035420748243745;
      x[37] = 0.6323435811043837081869820863;
      x[38] = 0.6556614359951054780787563493;
      x[39] = 0.6786100791688340579752213075;
      x[40] = 0.7011350789819958018478833856;
      x[41] = 0.7231830086267320439924738574;
      x[42] = 0.7447015728535264787392631535;
      x[43] = 0.7656397320099472728290069518;
      x[44] = 0.7859478231013170171419390583;
      x[45] = 0.8055776775861966251244264855;
      x[46] = 0.8244827356273286699288806160;
      x[47] = 0.8426181565271166212817791855;
      x[48] = 0.8599409250858054134244701089;
      x[49] = 0.8764099536302659483059318874;
      x[50] = 0.8919861794716707038051102626;
      x[51] = 0.9066326575613987798709616690;
      x[52] = 0.9203146481262901813758457723;
      x[53] = 0.9329996990770464098803916925;
      x[54] = 0.9446577229975570529267020191;
      x[55] = 0.9552610685392514028781903340;
      x[56] = 0.9647845860659697879107450773;
      x[57] = 0.9732056874292014080312407457;
      x[58] = 0.9805043998260268594593070609;
      x[59] = 0.9866634138949554818709267537;
      x[60] = 0.9916681269423129784656496511;
      x[61] = 0.9955066857383721603696911917;
      x[62] = 0.9981700583859776396734622503;
      x[63] = 0.9996525208678860697284528122;
      w[0] = 0.000891640360848216473648039572486;
      w[1] = 0.002073516630281233817643767864276;
      w[2] = 0.003252228984489181428058680199991;
      w[3] = 0.004423379913181973861515457329865;
      w[4] = 0.005584069730065564409295246509604;
      w[5] = 0.00673152394835932129903038334298;
      w[6] = 0.00786301523801235966098299764877;
      w[7] = 0.00897585788784867154252265100056;
      w[8] = 0.01006741157676510468617015836427;
      w[9] = 0.01113508690419162707964916519208;
      w[10] = 0.01217635128435543666908877520453;
      w[11] = 0.01318873485752732933584589631261;
      w[12] = 0.01416983630712974161375565260012;
      w[13] = 0.01511732853620123943398702990977;
      w[14] = 0.01602896417742577679273375217395;
      w[15] = 0.01690258091857080469578274105536;
      w[16] = 0.01773610662844119190534657335762;
      w[17] = 0.01852756427012002302020755090479;
      w[18] = 0.01927507658930781456448124847340;
      w[19] = 0.01997687056636017069332846306417;
      w[20] = 0.02063128162131176430507814873682;
      w[21] = 0.02123675756182679450366988395441;
      w[22] = 0.02179186226466172668841393048687;
      w[23] = 0.02229527908187828153006735501547;
      w[24] = 0.02274581396370907223988549848563;
      w[25] = 0.02314239829065720864797662461613;
      w[26] = 0.02348409140810500866266314287729;
      w[27] = 0.02377008285741515433114110347211;
      w[28] = 0.02399969429822915386406308993567;
      w[29] = 0.02417238111740147858488476357901;
      w[30] = 0.02428773372075171346739953339199;
      w[31] = 0.0243454785045698601916826953674;
      w[32] = 0.0243454785045698601916826953674;
      w[33] = 0.02428773372075171346739953339199;
      w[34] = 0.02417238111740147858488476357901;
      w[35] = 0.02399969429822915386406308993567;
      w[36] = 0.02377008285741515433114110347211;
      w[37] = 0.02348409140810500866266314287729;
      w[38] = 0.02314239829065720864797662461613;
      w[39] = 0.02274581396370907223988549848563;
      w[40] = 0.02229527908187828153006735501547;
      w[41] = 0.02179186226466172668841393048687;
      w[42] = 0.02123675756182679450366988395441;
      w[43] = 0.02063128162131176430507814873682;
      w[44] = 0.01997687056636017069332846306417;
      w[45] = 0.01927507658930781456448124847340;
      w[46] = 0.01852756427012002302020755090479;
      w[47] = 0.01773610662844119190534657335762;
      w[48] = 0.01690258091857080469578274105536;
      w[49] = 0.01602896417742577679273375217395;
      w[50] = 0.01511732853620123943398702990977;
      w[51] = 0.01416983630712974161375565260012;
      w[52] = 0.01318873485752732933584589631261;
      w[53] = 0.01217635128435543666908877520453;
      w[54] = 0.01113508690419162707964916519208;
      w[55] = 0.01006741157676510468617015836427;
      w[56] = 0.00897585788784867154252265100056;
      w[57] = 0.00786301523801235966098299764877;
      w[58] = 0.00673152394835932129903038334298;
      w[59] = 0.005584069730065564409295246509604;
      w[60] = 0.004423379913181973861515457329865;
      w[61] = 0.003252228984489181428058680199991;
      w[62] = 0.002073516630281233817643767864276;
      w[63] = 0.000891640360848216473648039572486;
    }
  };
 
  template <typename ctype> struct GLQuadrature<96, ctype> {
    __constant__ ctype x[96];
    __constant__ ctype w[96];
    consteval GLQuadrature()
    {
      x[0] = 0.0001552480583846165861549471078;
      x[1] = 0.000817812068409161137925280237;
      x[2] = 0.002009078506395354674800457560;
      x[3] = 0.003728049838118687714053851149;
      x[4] = 0.00597293683518810025961861379;
      x[5] = 0.00874136821849266127647708406;
      x[6] = 0.01203041270743176677369948288;
      x[7] = 0.01583658576836789391317031678;
      x[8] = 0.02015585427562873034996598706;
      x[9] = 0.02498364110778118212195052578;
      x[10] = 0.03031483012362239153407131028;
      x[11] = 0.03614377163884565451765476365;
      x[12] = 0.04246428843955096289705776692;
      x[13] = 0.04926968234207382934038365135;
      x[14] = 0.0565527412987897919715612833;
      x[15] = 0.0643057470453517485631126178;
      x[16] = 0.0725204832826992722686065151;
      x[17] = 0.0811882443859064392528485916;
      x[18] = 0.0902998446310341622305001879;
      x[19] = 0.0998456279304295913856019280;
      x[20] = 0.1098154780662833911981977721;
      x[21] = 0.1201988294116762506485147794;
      x[22] = 0.1309846781277999335744171345;
      x[23] = 0.1421615938255161868874279260;
      x[24] = 0.1537177316789142193278771151;
      x[25] = 0.1656408449780419230237213962;
      x[26] = 0.1779182981075164466007938250;
      x[27] = 0.1905370799372657148068153565;
      x[28] = 0.2034838176112139596582221225;
      x[29] = 0.2167447907193014157978749025;
      x[30] = 0.2303059458378212818865987013;
      x[31] = 0.2441529114226661632072451273;
      x[32] = 0.2582710130397018201157971953;
      x[33] = 0.2726452889161284956821619260;
      x[34] = 0.2872605057963497273175903982;
      x[35] = 0.3021011750855456983574998784;
      x[36] = 0.3171515692638431824845522030;
      x[37] = 0.3323957385536872886918371876;
      x[38] = 0.3478175278227518234880351104;
      x[39] = 0.3634005937044754292563639221;
      x[40] = 0.3791284219180799938360340481;
      x[41] = 0.3949843447697163981985764071;
      x[42] = 0.4109515588161906986202986874;
      x[43] = 0.4270131426725515290054463333;
      x[44] = 0.4431520749446670395443959523;
      x[45] = 0.4593512522677872205027643685;
      x[46] = 0.4755935074319751344440208979;
      x[47] = 0.4918616275751985152104327182;
      x[48] = 0.5081383724248014847895672818;
      x[49] = 0.524406492568024865555979102;
      x[50] = 0.540648747732212779497235631;
      x[51] = 0.556847925055332960455604048;
      x[52] = 0.572986857327448470994553667;
      x[53] = 0.589048441183809301379701313;
      x[54] = 0.605015655230283601801423593;
      x[55] = 0.620871578081920006163965952;
      x[56] = 0.636599406295524570743636078;
      x[57] = 0.652182472177248176511964890;
      x[58] = 0.667604261446312711308162812;
      x[59] = 0.682848430736156817515447797;
      x[60] = 0.697898824914454301642500122;
      x[61] = 0.712739494203650272682409602;
      x[62] = 0.727354711083871504317838074;
      x[63] = 0.741728986960298179884202805;
      x[64] = 0.755847088577333836792754873;
      x[65] = 0.769694054162178718113401299;
      x[66] = 0.783255209280698584202125098;
      x[67] = 0.796516182388786040341777878;
      x[68] = 0.809462920062734285193184643;
      x[69] = 0.822081701892483553399206175;
      x[70] = 0.834359155021958076976278604;
      x[71] = 0.846282268321085780672122885;
      x[72] = 0.857838406174483813112572074;
      x[73] = 0.869015321872200066425582866;
      x[74] = 0.879801170588323749351485221;
      x[75] = 0.890184521933716608801802228;
      x[76] = 0.900154372069570408614398072;
      x[77] = 0.909700155368965837769499812;
      x[78] = 0.918811755614093560747151408;
      x[79] = 0.927479516717300727731393485;
      x[80] = 0.935694252954648251436887382;
      x[81] = 0.943447258701210208028438717;
      x[82] = 0.950730317657926170659616349;
      x[83] = 0.957535711560449037102942233;
      x[84] = 0.963856228361154345482345236;
      x[85] = 0.969685169876377608465928690;
      x[86] = 0.975016358892218817878049474;
      x[87] = 0.979844145724371269650034013;
      x[88] = 0.984163414231632106086829683;
      x[89] = 0.987969587292568233226300517;
      x[90] = 0.991258631781507338723522916;
      x[91] = 0.994027063164811899740381386;
      x[92] = 0.996271950161881312285946149;
      x[93] = 0.997990921493604645325199542;
      x[94] = 0.999182187931590838862074720;
      x[95] = 0.999844751941615383413845053;
      w[0] = 0.0003983960327760062147190717484718;
      w[1] = 0.000926980394473460866167962675447;
      w[2] = 0.001455365908967473204205308994700;
      w[3] = 0.001982277169222343336866707883710;
      w[4] = 0.002507101371463758846235097484515;
      w[5] = 0.003029272752117980841658371015865;
      w[6] = 0.003548235395576932634572080406072;
      w[7] = 0.004063438462849379608691213853928;
      w[8] = 0.004574335615391693316292301332604;
      w[9] = 0.005080385267504207878793818476913;
      w[10] = 0.005581051049919249295606631914281;
      w[11] = 0.00607580233554415981759067636833;
      w[12] = 0.00656411478348078631853183345130;
      w[13] = 0.00704547088615743045793081236232;
      w[14] = 0.00751936051349746900293813761105;
      w[15] = 0.00798528145128114569030822839575;
      w[16] = 0.00844273993212258622523877030343;
      w[17] = 0.00889125115802263041880711324304;
      w[18] = 0.00933033981370573369257837931107;
      w[19] = 0.00975954057007251120504261014060;
      w[20] = 0.01017839857716666229762260770864;
      w[21] = 0.01058646994609564949383693359550;
      w[22] = 0.01098332221937217459737819340078;
      w[23] = 0.01136853482916468700067392098875;
      w[24] = 0.01174169954296310992111796633381;
      w[25] = 0.01210242089618234564113366893634;
      w[26] = 0.01245031661124180514419109043417;
      w[27] = 0.01278501800267468074939858397180;
      w[28] = 0.01310617036783620695672898198223;
      w[29] = 0.01341343336279588109902836435708;
      w[30] = 0.01370648136301462141171054374455;
      w[31] = 0.01398500380842416721990928829451;
      w[32] = 0.01424870553254269282279975647290;
      w[33] = 0.01449730707527761827133939063984;
      w[34] = 0.01473054497908395298521816609143;
      w[35] = 0.01494817206816419299219403789720;
      w[36] = 0.01514995771041379689704438210325;
      w[37] = 0.01533568806183457450711441517810;
      w[38] = 0.01550516629315691871162488998532;
      w[39] = 0.01565821279843067790639213335753;
      w[40] = 0.01579466538536358427901037308498;
      w[41] = 0.01591437944720550326737686994277;
      w[42] = 0.01601722811599633160906948873511;
      w[43] = 0.01610310239701512533433355727862;
      w[44] = 0.01617191128428796421438741941447;
      w[45] = 0.01622358185703213468200639422442;
      w[46] = 0.01625805935693441799360274572389;
      w[47] = 0.0162753072461815831209807091486;
      w[48] = 0.0162753072461815831209807091486;
      w[49] = 0.01625805935693441799360274572389;
      w[50] = 0.01622358185703213468200639422442;
      w[51] = 0.01617191128428796421438741941447;
      w[52] = 0.01610310239701512533433355727862;
      w[53] = 0.01601722811599633160906948873511;
      w[54] = 0.01591437944720550326737686994277;
      w[55] = 0.01579466538536358427901037308498;
      w[56] = 0.01565821279843067790639213335753;
      w[57] = 0.01550516629315691871162488998532;
      w[58] = 0.01533568806183457450711441517810;
      w[59] = 0.01514995771041379689704438210325;
      w[60] = 0.01494817206816419299219403789720;
      w[61] = 0.01473054497908395298521816609143;
      w[62] = 0.01449730707527761827133939063984;
      w[63] = 0.01424870553254269282279975647290;
      w[64] = 0.01398500380842416721990928829451;
      w[65] = 0.01370648136301462141171054374455;
      w[66] = 0.01341343336279588109902836435708;
      w[67] = 0.01310617036783620695672898198223;
      w[68] = 0.01278501800267468074939858397180;
      w[69] = 0.01245031661124180514419109043417;
      w[70] = 0.01210242089618234564113366893634;
      w[71] = 0.01174169954296310992111796633381;
      w[72] = 0.01136853482916468700067392098875;
      w[73] = 0.01098332221937217459737819340078;
      w[74] = 0.01058646994609564949383693359550;
      w[75] = 0.01017839857716666229762260770864;
      w[76] = 0.00975954057007251120504261014060;
      w[77] = 0.00933033981370573369257837931107;
      w[78] = 0.00889125115802263041880711324304;
      w[79] = 0.00844273993212258622523877030343;
      w[80] = 0.00798528145128114569030822839575;
      w[81] = 0.00751936051349746900293813761105;
      w[82] = 0.00704547088615743045793081236232;
      w[83] = 0.00656411478348078631853183345130;
      w[84] = 0.00607580233554415981759067636833;
      w[85] = 0.005581051049919249295606631914281;
      w[86] = 0.005080385267504207878793818476913;
      w[87] = 0.004574335615391693316292301332604;
      w[88] = 0.004063438462849379608691213853928;
      w[89] = 0.003548235395576932634572080406072;
      w[90] = 0.003029272752117980841658371015865;
      w[91] = 0.002507101371463758846235097484515;
      w[92] = 0.001982277169222343336866707883710;
      w[93] = 0.001455365908967473204205308994700;
      w[94] = 0.000926980394473460866167962675447;
      w[95] = 0.0003983960327760062147190717484718;
    }
  };
 
  template <typename ctype> struct GLQuadrature<128, ctype> {
    __constant__ ctype x[128];
    __constant__ ctype w[128];
    consteval GLQuadrature()
    {
      x[0] = 0.0000875560264340427631959585091;
      x[1] = 0.000461270011312052494006121158;
      x[2] = 0.001133375687242990058921278326;
      x[3] = 0.002103620732509406567919395190;
      x[4] = 0.00337144354989353234828138946;
      x[5] = 0.00493609075413280833103484226;
      x[6] = 0.00679662863770689556438224038;
      x[7] = 0.00895194578214073198761718505;
      x[8] = 0.01140075426804630641731278674;
      x[9] = 0.01414159062643170954783081991;
      x[10] = 0.01717281678401736567708548671;
      x[11] = 0.02049262107315003595054072501;
      x[12] = 0.02409901932936780689110184920;
      x[13] = 0.02798985608488990893944429105;
      x[14] = 0.0321628058610418121084366226;
      x[15] = 0.0366153745605260783326877409;
      x[16] = 0.0413449009595197314817582099;
      x[17] = 0.0463485582991215930392570263;
      x[18] = 0.0516233559754209078067763012;
      x[19] = 0.0571661413273013912958537877;
      x[20] = 0.0629736015209841006522910226;
      x[21] = 0.0690422655302257697046838642;
      x[22] = 0.0753685062110155154181999355;
      x[23] = 0.0819485424695465764415623451;
      x[24] = 0.0887784415221780787677029071;
      x[25] = 0.0958541212460431699401788757;
      x[26] = 0.1031713526189033548783335267;
      x[27] = 0.1107257622467940165747529667;
      x[28] = 0.1185128349779526386101154283;
      x[29] = 0.1265279166014690094150587893;
      x[30] = 0.1347662166290455967641315285;
      x[31] = 0.1432228111582062933280700435;
      x[32] = 0.1518926458152428338074566988;
      x[33] = 0.1607705387761403703161221716;
      x[34] = 0.1698511838636769739470265826;
      x[35] = 0.1791291537188462214232375419;
      x[36] = 0.1885989030447075446192302007;
      x[37] = 0.1982547719207256878982134012;
      x[38] = 0.2080909891856184552249805714;
      x[39] = 0.2181016758866909580427846149;
      x[40] = 0.2282808487935948182779031864;
      x[41] = 0.2386224239744122607730260303;
      x[42] = 0.2491202204319277678551968533;
      x[43] = 0.2597679637979139870708621467;
      x[44] = 0.2705592900832239022754554348;
      x[45] = 0.2814877494814479185314785283;
      x[46] = 0.292546810223862492253863060;
      x[47] = 0.303729862483366278632175905;
      x[48] = 0.315030222325070486691704129;
      x[49] = 0.326441135701182245786918614;
      x[50] = 0.337955782487793312408373860;
      x[51] = 0.349567280561161398666422953;
      x[52] = 0.361268689911047798596884183;
      x[53] = 0.373053016788652839572191002;
      x[54] = 0.384913217886670006794506711;
      x[55] = 0.396842204548960391422970995;
      x[56] = 0.408832847007331408794808713;
      x[57] = 0.420877978642887533001262274;
      x[58] = 0.432970400269406107441212341;
      x[59] = 0.445102884436178126663512666;
      x[60] = 0.457268179747742250681750983;
      x[61] = 0.469459015197930215948106488;
      x[62] = 0.481668104515633253334892326;
      x[63] = 0.493888150519692117900973940;
      x[64] = 0.506111849480307882099026060;
      x[65] = 0.518331895484366746665107674;
      x[66] = 0.530540984802069784051893512;
      x[67] = 0.542731820252257749318249017;
      x[68] = 0.554897115563821873336487334;
      x[69] = 0.567029599730593892558787659;
      x[70] = 0.579122021357112466998737726;
      x[71] = 0.591167152992668591205191287;
      x[72] = 0.603157795451039608577029005;
      x[73] = 0.615086782113329993205493289;
      x[74] = 0.626946983211347160427808998;
      x[75] = 0.638731310088952201403115817;
      x[76] = 0.650432719438838601333577047;
      x[77] = 0.662044217512206687591626140;
      x[78] = 0.673558864298817754213081386;
      x[79] = 0.684969777674929513308295871;
      x[80] = 0.696270137516633721367824095;
      x[81] = 0.707453189776137507746136940;
      x[82] = 0.718512250518552081468521472;
      x[83] = 0.729440709916776097724544565;
      x[84] = 0.740232036202086012929137853;
      x[85] = 0.750879779568072232144803147;
      x[86] = 0.761377576025587739226973970;
      x[87] = 0.771719151206405181722096814;
      x[88] = 0.781898324113309041957215385;
      x[89] = 0.791909010814381544775019429;
      x[90] = 0.801745228079274312101786599;
      x[91] = 0.811401096955292455380769799;
      x[92] = 0.820870846281153778576762458;
      x[93] = 0.830148816136323026052973417;
      x[94] = 0.839229461223859629683877828;
      x[95] = 0.848107354184757166192543301;
      x[96] = 0.856777188841793706671929957;
      x[97] = 0.865233783370954403235868471;
      x[98] = 0.873472083398530990584941211;
      x[99] = 0.881487165022047361389884572;
      x[100] = 0.889274237753205983425247033;
      x[101] = 0.896828647381096645121666473;
      x[102] = 0.904145878753956830059821124;
      x[103] = 0.911221558477821921232297093;
      x[104] = 0.918051457530453423558437655;
      x[105] = 0.924631493788984484581800064;
      x[106] = 0.930957734469774230295316136;
      x[107] = 0.937026398479015899347708977;
      x[108] = 0.942833858672698608704146212;
      x[109] = 0.948376644024579092193223699;
      x[110] = 0.953651441700878406960742974;
      x[111] = 0.958655099040480268518241790;
      x[112] = 0.963384625439473921667312259;
      x[113] = 0.967837194138958187891563377;
      x[114] = 0.972010143915110091060555709;
      x[115] = 0.975900980670632193108898151;
      x[116] = 0.979507378926849964049459275;
      x[117] = 0.982827183215982634322914513;
      x[118] = 0.985858409373568290452169180;
      x[119] = 0.988599245731953693582687213;
      x[120] = 0.991048054217859268012382815;
      x[121] = 0.993203371362293104435617760;
      x[122] = 0.995063909245867191668965158;
      x[123] = 0.996628556450106467651718611;
      x[124] = 0.997896379267490593432080605;
      x[125] = 0.998866624312757009941078722;
      x[126] = 0.999538729988687947505993879;
      x[127] = 0.999912443973565957236804041;
      w[0] = 0.0002246904801460451881971461199944;
      w[1] = 0.0005229063396701743896564258000535;
      w[2] = 0.000821251509334514769395437797399;
      w[3] = 0.001119144215481309371811027136358;
      w[4] = 0.001416375735728995547642867323410;
      w[5] = 0.001712763020455107887168892330067;
      w[6] = 0.002008127491869321156597171743152;
      w[7] = 0.002302292128351477559145270990142;
      w[8] = 0.002595080916338165102535383567425;
      w[9] = 0.002886318771432849294667308813034;
      w[10] = 0.003175831580853594393607163913116;
      w[11] = 0.003463446283449406781713335180164;
      w[12] = 0.003748990962817364343835981344235;
      w[13] = 0.004032294945243028986464299348931;
      w[14] = 0.004313188899308374852489421890984;
      w[15] = 0.004591504935830437167239371844213;
      w[16] = 0.004867076707503402931774133046779;
      w[17] = 0.005139739507916078566607670163193;
      w[18] = 0.005409330369751538123829823138734;
      w[19] = 0.005675688162040208346640834226269;
      w[20] = 0.005938653686370139787945553462789;
      w[21] = 0.00619806977197546148441086409866;
      w[22] = 0.00645378136963367361022141700193;
      w[23] = 0.00670563564430816615724447580808;
      w[24] = 0.00695348206647599262214400369780;
      w[25] = 0.00719717250208342308841194600427;
      w[26] = 0.00743656130107365712619274926004;
      w[27] = 0.00767150538443257204299542687072;
      w[28] = 0.00790186432969967342948281584369;
      w[29] = 0.00812750045489259352582872823853;
      w[30] = 0.00834827890079460229454575397690;
      w[31] = 0.00856406771155568841534049380942;
      w[32] = 0.00877473791355885232435346281719;
      w[33] = 0.00898016359250434297009846376248;
      w[34] = 0.00918022196866567161064464549565;
      w[35] = 0.00937479347027235432545977423692;
      w[36] = 0.00956376180497547274325926583392;
      w[37] = 0.00974701402935330141151095934061;
      w[38] = 0.00992444061641543110997220663244;
      w[39] = 0.01009593552106502059033657920306;
      w[40] = 0.01026139624348003471614248389408;
      w[41] = 0.01042072389037557455679197421145;
      w[42] = 0.01057382323411067426850976759020;
      w[43] = 0.01072060276960423006855592693901;
      w[44] = 0.01086097476902603768763047888381;
      w[45] = 0.01099485533423024571706105329958;
      w[46] = 0.01112216444689988255231456680360;
      w[47] = 0.01124282601637248343591230197029;
      w[48] = 0.01135676792511823065485631796147;
      w[49] = 0.01146392207184342346020549360466;
      w[50] = 0.01156422441219351393964895120138;
      w[51] = 0.01165761499703138006120783563640;
      w[52] = 0.01174403800826795657651263664089;
      w[53] = 0.01182344179222380757182569615143;
      w[54] = 0.01189577889050170031939035494271;
      w[55] = 0.01196100606835172783622520440854;
      w[56] = 0.01201908434051202631879365840982;
      w[57] = 0.01206997899450964249885832694523;
      w[58] = 0.01211365961140762406004665422089;
      w[59] = 0.01215010008398593266172130318192;
      w[60] = 0.01217927863234531292663426012314;
      w[61] = 0.01220117781692479104664899484710;
      w[62] = 0.01221578454892502252742428071432;
      w[63] = 0.0122230900981312591056629263053;
      w[64] = 0.0122230900981312591056629263053;
      w[65] = 0.01221578454892502252742428071432;
      w[66] = 0.01220117781692479104664899484710;
      w[67] = 0.01217927863234531292663426012314;
      w[68] = 0.01215010008398593266172130318192;
      w[69] = 0.01211365961140762406004665422089;
      w[70] = 0.01206997899450964249885832694523;
      w[71] = 0.01201908434051202631879365840982;
      w[72] = 0.01196100606835172783622520440854;
      w[73] = 0.01189577889050170031939035494271;
      w[74] = 0.01182344179222380757182569615143;
      w[75] = 0.01174403800826795657651263664089;
      w[76] = 0.01165761499703138006120783563640;
      w[77] = 0.01156422441219351393964895120138;
      w[78] = 0.01146392207184342346020549360466;
      w[79] = 0.01135676792511823065485631796147;
      w[80] = 0.01124282601637248343591230197029;
      w[81] = 0.01112216444689988255231456680360;
      w[82] = 0.01099485533423024571706105329958;
      w[83] = 0.01086097476902603768763047888381;
      w[84] = 0.01072060276960423006855592693901;
      w[85] = 0.01057382323411067426850976759020;
      w[86] = 0.01042072389037557455679197421145;
      w[87] = 0.01026139624348003471614248389408;
      w[88] = 0.01009593552106502059033657920306;
      w[89] = 0.00992444061641543110997220663244;
      w[90] = 0.00974701402935330141151095934061;
      w[91] = 0.00956376180497547274325926583392;
      w[92] = 0.00937479347027235432545977423692;
      w[93] = 0.00918022196866567161064464549565;
      w[94] = 0.00898016359250434297009846376248;
      w[95] = 0.00877473791355885232435346281719;
      w[96] = 0.00856406771155568841534049380942;
      w[97] = 0.00834827890079460229454575397690;
      w[98] = 0.00812750045489259352582872823853;
      w[99] = 0.00790186432969967342948281584369;
      w[100] = 0.00767150538443257204299542687072;
      w[101] = 0.00743656130107365712619274926004;
      w[102] = 0.00719717250208342308841194600427;
      w[103] = 0.00695348206647599262214400369780;
      w[104] = 0.00670563564430816615724447580808;
      w[105] = 0.00645378136963367361022141700193;
      w[106] = 0.00619806977197546148441086409866;
      w[107] = 0.005938653686370139787945553462789;
      w[108] = 0.005675688162040208346640834226269;
      w[109] = 0.005409330369751538123829823138734;
      w[110] = 0.005139739507916078566607670163193;
      w[111] = 0.004867076707503402931774133046779;
      w[112] = 0.004591504935830437167239371844213;
      w[113] = 0.004313188899308374852489421890984;
      w[114] = 0.004032294945243028986464299348931;
      w[115] = 0.003748990962817364343835981344235;
      w[116] = 0.003463446283449406781713335180164;
      w[117] = 0.003175831580853594393607163913116;
      w[118] = 0.002886318771432849294667308813034;
      w[119] = 0.002595080916338165102535383567425;
      w[120] = 0.002302292128351477559145270990142;
      w[121] = 0.002008127491869321156597171743152;
      w[122] = 0.001712763020455107887168892330067;
      w[123] = 0.001416375735728995547642867323410;
      w[124] = 0.001119144215481309371811027136358;
      w[125] = 0.000821251509334514769395437797399;
      w[126] = 0.0005229063396701743896564258000535;
      w[127] = 0.0002246904801460451881971461199944;
    }
  };
 
} // namespace DiFfRG