libome-rs 0.0.1

/******************************************************************************
 * This file is part of libome                                                *
 * Copyright (C) 2025 Arnd Behring, Kay Schoenwald                            *
 * SPDX-License-Identifier: GPL-3.0-or-later                                  *
 ******************************************************************************/

#include <ome/AQqPS.h>

namespace ome
{

  /**************** AQqPS, regular part ****************/

  const ome_as<double> AQqPS_reg = 
    // AQqPS_reg
    ome_as<double>({
      // AQqPS_reg as^2
      ome_logm<double>({
        // AQqPS_reg as^2 LM^0
        ome_nf<double>({
          // AQqPS_reg as^2 LM^0 NF^0
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^2 LM^0 NF^0 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^2 LM^0 NF^0 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^0 NF^0 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    -11.0617283950617283950617283951,17.6438806007380114218610807789,-30.3561193992619885781389192211,
                    18.7160493827160493827160493827,3.95061728395061728395061728395,0.629135802469135802469135802469,
                    0.208592592592592592592592592593,0.0950728143109095490047871000252,0.051676204873483785048410898751,
                    0.0313865993193279663423276424032,0.0205728189766990472457844591883,0.0142547207956860033178776164186,
                    0.0103035950773936693740795363079,0.00769990042293496415132354950113,0.00591129590335512618357978126415,
                    0.00464032954794940062419835394114,0.00371145737027474337783076663244,0.00301630772134561885756481241734,
                    0.00248540255451718614453912340849,0.00207274401469949263344289259215,0.00174702411179932174566962224598,
                    0.00148641309377560593473567660343,0.00127536375790640981985714288104,0.00110259326023095637551427835171,
                    0.000959774920618354029164174360177,0.000840670359300234944986894497618,
                    0.000740541666545835030154743207916,0.000655745721899665455318412680772,
                    0.000583449429595313537968610728147,0.000521426717728167517633379340442,
                    0.000467911765072300876070976853717,0.000421491495538246849475173484843,
                    0.000381025887354731905312337548846,0.000345588243727238954595599830128,
                    0.000314419963416028313040757351598,0.000286895962959548049278985835598,
                    0.000262498005844367492028898348093,0.000240793958729487177455895650351,
                    0.000221421531341739717190038187014,0.000204075437314956578362390512602,
                    0.000188497186210704385828311958776,0.000174466914663459458294034522368,
                    0.000161796809139498515820602056938,0.000150325779419352845721516051238,
                    0.000139915121219315584350485332462,0.000130444965824015937732067133494,
                    0.000121811359512381878516477620456,0.000113923849723927034730736393782,
                    0.000106703481078251727053730125218,0.000100081124529595273300445125674,
                    0.000093996078579972056564157861636,0.0000883948936759754793021949996213
                  },-1),
                  // AQqPS_reg as^2 LM^0 NF^0 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    -5.33333333333333333333333333333,0,-20.1481481481481481481481481481,-2.37037037037037037037037037037,
                    -0.725925925925925925925925925926,-0.358518518518518518518518518519,-0.21671957671957671957671957672,
                    -0.146334089191232048374905517763,-0.105946082136558327034517510708,-0.0804820693709582598471487360376,
                    -0.0633295922184811073699962588851,-0.0511947760432608917457402305887,-0.0422782604600786418968237150055,
                    -0.0355257278334201411124488047565,-0.0302844083063863283643503423723,-0.0261317518460375603232746089889,
                    -0.0227840024898848428260192966075,-0.0200446409500619419881242257251,-0.0177739405971462467951672508624,
                    -0.0158702761566991284732490171268,-0.0142582577300952185532056261511,-0.0128809758884946854871666901742,
                    -0.0116948028862201645939647068953,-0.0106658281602430081550205282805,
                    -0.00976736489023823617396396224374,-0.00897817416078285643503034807383,
                    -0.00828117953588039058124528209998,-0.00766252283821134723318958836528,
                    -0.00711086121291071725907218213975,-0.00661683736632305041039918264688,
                    -0.00617267579639563836564335202669,-0.00577187183954559891981399749888,
                    -0.00540894988735779578651586976249,-0.00507927370420901557696451607095,
                    -0.00477889638002881161107018294813,-0.00450444071468952009403484111014,
                    -0.0042530031699032632739448799206,-0.00402207622059052436798604018455,
                    -0.00380948518005242319584362458261,-0.00361333649209208552911940936659,
                    -0.00343197516891602010590852301132,-0.00326394956971994449686299775532,
                    -0.00310798210640299572492678557926,-0.00296294476216573668252040401597,
                    -0.00282783853919766862061632371546,-0.00270177613029795113524214442395,
                    -0.00258396724864077065250607372547,-0.00247370615929324745837819039247,
                    -0.00237036104246471771756532855257,-0.0022733648870364954975090822743,
                    -0.00218220766764755788936889948723
                  },0),
                  // AQqPS_reg as^2 LM^0 NF^0 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    1.33333333333333333333333333333,6.66666666666666666666666666667,3.55555555555555555555555555556
                  },0),
                  // AQqPS_reg as^2 LM^0 NF^0 [below 0.25] y=x H[0, y]^3
                  ome_x<double>({
                    -0.888888888888888888888888888889,-0.888888888888888888888888888889
                  },0)
                },0)
              ),
              // AQqPS_reg as^2 LM^0 NF^0 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^2 LM^0 NF^0 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^0 NF^0 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -3.88727443033318664562764707759,-16.0542272743755863996413914931,-44.8525907724458029508207761453,
                    -118.780438761492441879239422676,-254.780521489156429841338011468,-529.046595299652055527863386375,
                    -1084.26548695276488456055478162,-2209.34471791947585488086284759,-4486.81464168941021365040770136,
                    -9091.20020925622669212914566518,-18389.193777612125199121887759,-37146.6767201752562191030719564,
                    -74955.2946779334628508556096307,-151109.000039289989984400158074,-304400.847068837562872660369669,
                    -612798.438614699524230331933406,-1.23294846482420663052503803088e6,-2.4794748085743247271223714092e6,
                    -4.9841195845613244918832728603e6,-1.00150542375765816496389818919e7,-2.01174488101128623575108624023e7,
                    -4.03983009886692604035618886484e7,-8.11031005948927601586804532456e7,-1.62782480494518513377325943617e8,
                    -3.26650885263190188416282014046e8,-6.55352276664231194229045723477e8,-1.31458388986278738785918515378e9,
                    -2.63652026956638550453966584195e9,-5.28699917210833249973652830013e9,
                    -1.06005406231012157466425171036e10,-2.12516264254404319496175656461e10,
                    -4.25996424863421558220657767403e10,-8.53833435071071131832959295439e10,
                    -1.71118591934098319563243618177e11,-3.4291079566667993814526845902e11,
                    -6.87112461111844191069013667293e11,-1.37670131207996443753153123148e12,
                    -2.75815812949240482025763361156e12,-5.52545727087918455486140084493e12,
                    -1.10685015582511475218821287662e13,-2.21708697525282901215442079089e13,
                    -4.44070099938503506453316875544e13,-8.89399156093602599832511820877e13,
                    -1.78122849200866756938148378747e14,-3.56715145025184348351481278739e14,
                    -7.14337778109029088813910982477e14,-1.43043101309443422930975467065e15,
                    -2.86426001813634030649138411386e15,-5.73510156103514323527351954356e15,
                    -1.14829584718141365573760091636e16,-2.29906518436578971770554198292e16
                  },0)
                },0)
              ),
              // AQqPS_reg as^2 LM^0 NF^0 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^2 LM^0 NF^0 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^0 NF^0 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    -6.10631502319054099451954755545,2.66666666666666666666666666667,-5.21742613430165210563065866656,
                    -6.18038909726461506859362162952,-7.06541441285244791655908802952,-7.51488417288472520896899887396,
                    -7.78522333227747082468543039141,-7.97421974691377874204045615206,-8.12082820615865699096589859012,
                    -8.2423644635158555179132744373,-8.34738626647327185841618033667,-8.44055366854941036969671159628,
                    -8.5246373340255429823012937049,-8.60141861964818214344106980186,-8.67212297139544026709499713337,
                    -8.73764245882000796051767904717,-8.79865694277361814838644908518,-8.85570373629650426618901941367,
                    -8.90921973915766361928560649022,-8.95956815886437501178103256523,-9.00705621529680910191492319638,
                    -9.05194734934821156892709621995,-9.09446994875159345150151754259,-9.13482378486569511872911400481,
                    -9.17318489307624747798368038913,-9.20970936140984781323326245209,-9.24453633125006043902052172891,
                    -9.27779041482295782993832163916,-9.30958367111356637293615220103,-9.34001724077920052363323338146,
                    -9.3691827131287953250595460685,-9.39716327939006403163567204809,-9.42403471327171572132792661301,
                    -9.44986621036498378039180506254,-9.47472111101869647639018329668,-9.49865752618513910873103043957,
                    -9.52172888185222958653435478677,-9.54398439470061365782725508821,-9.56546948931059011511587563879,
                    -9.58622616542380513662886217213,-9.6062933223175080539952033565,-9.6257070461872447024416085606,
                    -9.64450086549278327616562786926,-9.66270597845383123351656747893,-9.68035145625049224374106044169,
                    -9.69746442496077507820293295348,-9.71407022883245798635150264406,-9.73019257712259563118312147319,
                    -9.74585367643184487100548837509,-9.76107435020220101508102066423
                  },1)
                },0)
              )
            }
          )
        },0),
        // AQqPS_reg as^2 LM^1
        ome_nf<double>({
          // AQqPS_reg as^2 LM^1 NF^0
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^2 LM^1 NF^0 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^2 LM^1 NF^0 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^1 NF^0 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    -11.8518518518518518518518518519,10.6666666666666666666666666667,-32.,33.1851851851851851851851851852
                  },-1),
                  // AQqPS_reg as^2 LM^1 NF^0 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    -5.33333333333333333333333333333,-26.6666666666666666666666666667,-14.2222222222222222222222222222
                  },0),
                  // AQqPS_reg as^2 LM^1 NF^0 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    5.33333333333333333333333333333,5.33333333333333333333333333333
                  },0)
                },0)
              ),
              // AQqPS_reg as^2 LM^1 NF^0 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^2 LM^1 NF^0 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^1 NF^0 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -1.49384595028523357857654005328,-12.871872494022514052303647146,-49.7112998023598797063562968289,
                    -126.249897185091537102802751779,-271.651063209720463072572487533,-565.410823229367414998958816093,
                    -1166.3550620439949336490087641,-2391.91307129374776774694651932,-4883.645033310777969231058873,
                    -9937.33059160244098940807338078,-20167.2634116319207425099436693,-40842.883939329733175480983662,
                    -82576.8344289679088591986812865,-166728.221127689042798905985633,-336259.133543111024682198322062,
                    -677540.032107477952979190033015,-1.36413162571895448410958709039e6,-2.74466654052202399381759578949e6,
                    -5.51920491866683501573220554872e6,-1.10930517146641693414239447315e7,-2.22864618280652905915620562913e7,
                    -4.47579347371027323486005260608e7,-8.98581058604632844272920588502e7,-1.80351282759136852539819435465e8,
                    -3.61884469349241788303597893576e8,-7.25974467651174281697471791046e8,-1.45607495099396615781666985332e9,
                    -2.91988670531734095124783385815e9,-5.85431250043728826427132653917e9,-1.1736002684518208038431941737e10,
                    -2.35236570549175535734940739708e10,-4.71449366525686496458352787278e10,
                    -9.44746926658349619878982632157e10,-1.89299842056646771372102555429e11,
                    -3.79265221219868976195874623659e11,-7.59796127279807160483521502532e11,
                    -1.52200250184342223817546372225e12,-3.04860068768327593873471064737e12,
                    -6.10597468246234850522039989779e12,-1.22287171328852559898299881982e13,
                    -2.44895163017531518011464279984e13,-4.90404796489271816420014229823e13,
                    -9.81987664652439100497241502769e13,-1.96623609028571353825958426022e14,
                    -3.93681459918371389368424107199e14,-7.88197726431294493495880444457e14,
                    -1.57800165910734837556995816041e15,-3.15909619839191449975717847015e15,
                    -6.32415254718556331041921389465e15,-1.26597990761910337407871973731e16,
                    -2.53417796337157578782970792713e16
                  },0)
                },0)
              ),
              // AQqPS_reg as^2 LM^1 NF^0 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^2 LM^1 NF^0 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^1 NF^0 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    -4.44444444444444444444444444444,-7.11111111111111111111111111111,-7.64444444444444444444444444444,
                    -7.94074074074074074074074074074,-8.18201058201058201058201058201,-8.39576719576719576719576719577,
                    -8.58800705467372134038800705467,-8.76134038800705467372134038801,-8.91783870450537117203783870451,
                    -9.05941959275292608625941959275,-9.18784672118005451338784672118,-9.30470164755879041593327307613,
                    -9.41137606661416185225709035233,-9.50908079574746241412908079575,-9.59886311945135474547239253122,
                    -9.6816268214961025418541758411,-9.75815187118936688940128912609,-9.82911238362218754375617120715,
                    -9.89509237596809467663072262807,-9.95659929251662595717921197395,-10.0140754739440479598827106751,
                    -10.0679078204317132840487262348,-10.1184359090254903957959550906,-10.1659588077132627077230864882,
                    -10.2107407998403899199234642571,-10.2530162014228531758503319973,-10.2929934246609817124721842494,
                    -10.330858415084999419251281968,-10.3667775676595859866187484096,-10.4009002086636731313793561657,
                    -10.4333607148426754147652206502,-10.4642803287402466076769436631,-10.4937687188079236125316907939,
                    -10.5219253244639683944574989341,-10.5488405193851401581896111366,-10.5745966206812035297289114477,
                    -10.5992687669865439323501640648,-10.6229256847146224203528532591,-10.6456303586037560349398664566,
                    -10.6674406201112605634130895082,-10.6884096650859726668405660323,-10.7085865103847776082991918569,
                    -10.7280163976309519823796991746,-10.7467411510874771938055800983,-10.7647994955937266312531483301,
                    -10.7822273396539626217996176036,-10.799058028042330738232454968,-10.8153225676782104284356573214
                  },3)
                },0)
              )
            }
          )
        },0),
        // AQqPS_reg as^2 LM^2
        ome_nf<double>({
          // AQqPS_reg as^2 LM^2 NF^0
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^2 LM^2 NF^0 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^2 LM^2 NF^0 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^2 NF^0 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    -3.55555555555555555555555555556,-2.66666666666666666666666666667,2.66666666666666666666666666667,
                    3.55555555555555555555555555556
                  },-1),
                  // AQqPS_reg as^2 LM^2 NF^0 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    -5.33333333333333333333333333333,-5.33333333333333333333333333333
                  },0)
                },0)
              ),
              // AQqPS_reg as^2 LM^2 NF^0 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^2 LM^2 NF^0 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^2 NF^0 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -2.01037811107599308021769858389,-8.14122940743081942800301575889,-19.5555555555555555555555555556,
                    -46.2222222222222222222222222222,-96.,-197.688888888888888888888888889,-403.911111111111111111111111111,
                    -820.825396825396825396825396825,-1661.96825396825396825396825397,-3356.44444444444444444444444444,
                    -6765.98518518518518518518518519,-13620.2343434343434343434343434,-27389.4141414141414141414141414,
                    -55033.4358974358974358974358974,-110506.979242979242979242979243,-221782.146031746031746031746032,
                    -444916.622222222222222222222222,-892231.947712418300653594771242,-1.78874729411764705882352941176e6,
                    -3.58518967641325536062378167641e6,-7.18427860584795321637426900585e6,-1.43937861079365079365079365079e7,
                    -2.8833570354978354978354978355e7,-5.77513478296003513394817742644e7,-1.15657426241545893719806763285e8,
                    -2.31600146204444444444444444444e8,-4.6372798604034188034188034188e8,-9.28434882005698005698005698006e8,
                    -1.85869065300881834215167548501e9,-3.72077695005145046524356869184e9,-7.44790114241226053639846743295e9,
                    -1.49076929370838709677419354839e10,-2.98377073887885304659498207885e10,
                    -5.97173988176161616161616161616e10,-1.19513908738129530600118835413e11,
                    -2.39177144240068720821661998133e11,-4.78636609392965079365079365079e11,
                    -9.57807801873937937937937937938e11,-1.916629325547890627469574838e12,
                    -3.83518360853580926675663517769e12,-7.67402730975967179487179487179e12,
                    -1.53550225481172986449864498645e13,-3.07233260729231807975222609369e13,
                    -6.14719944985019771133259505353e13,-1.22992398650201783415550857411e14,
                    -2.46077366003103806060606060606e14,-4.92331915786013763607085346216e14,
                    -9.85003295003699678898139582691e14,-1.97065754971302067139479905437e15,
                    -3.94256446970912043537414965986e15,-7.88752881677069453641723356009e15
                  },0)
                },0)
              ),
              // AQqPS_reg as^2 LM^2 NF^0 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^2 LM^2 NF^0 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^2 LM^2 NF^0 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    -2.66666666666666666666666666667,0,-2.66666666666666666666666666667,-2.66666666666666666666666666667,
                    -2.75555555555555555555555555556,-2.84444444444444444444444444444,-2.92063492063492063492063492063,
                    -2.98412698412698412698412698413,-3.03703703703703703703703703704,-3.08148148148148148148148148148,
                    -3.11919191919191919191919191919,-3.15151515151515151515151515152,-3.17948717948717948717948717949,
                    -3.20390720390720390720390720391,-3.2253968253968253968253968254,-3.24444444444444444444444444444,
                    -3.26143790849673202614379084967,-3.27668845315904139433551198257,-3.29044834307992202729044834308,
                    -3.30292397660818713450292397661,-3.31428571428571428571428571429,-3.32467532467532467532467532468,
                    -3.33421168203776899429073342117,-3.34299516908212560386473429952,-3.35111111111111111111111111111,
                    -3.35863247863247863247863247863,-3.36562203228869895536562203229,-3.3721340388007054673721340388,
                    -3.37821565407772304324028461959,-3.38390804597701149425287356322,-3.38924731182795698924731182796,
                    -3.39426523297491039426523297491,-3.39898989898989898989898989899,-3.40344622697563874034462269756,
                    -3.4076563958916900093370681606,-3.41164021164021164021164021164,-3.41541541541541541541541541542,
                    -3.41899794531373478741899794531,-3.42240215924426450742240215924,-3.42564102564102564102564102564,
                    -3.42872628726287262872628726287,-3.4316686024003097173828881146,-3.43447766703580657069029162052,
                    -3.43716232088325111580925534414,-3.43973063973063973063973063973,-3.44219001610305958132045088567,
                    -3.4445472299311337239181827526,-3.44680851063829787234042553191,-3.44897959183673469387755102041,
                    -3.45106575963718820861678004535
                  },1)
                },0)
              )
            }
          )
        },0)
      },0),
      // AQqPS_reg as^3
      ome_logm<double>({
        // AQqPS_reg as^3 LM^0
        ome_nf<double>({
          // AQqPS_reg as^3 LM^0 NF^0
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^0 NF^0 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^0 NF^0 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^0 NF^0 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    1644.53378973794478252003472602,-2271.14031554448826915764721741,3830.47210124278586693893766206,
                    -3528.55630548144350900756942364,534.131050280579228724070255556,-272.930204280382789530771443699,
                    172.398006804625075158676254208,-149.489991117790996017266171761,155.306202012684467577935094603,
                    -126.357451915101666424687413467,182.322816574245012257319317206,-75.6855741264178537009004471991,
                    278.070069704375975351857325797,98.9437474115340568162086189036,609.729152239731789448748268997,
                    719.46541001842098625583937537,1786.47463175053838518971398277,2946.73733711194292521718947775,
                    6023.40389061982177960257674464,11016.7182900456640918054654811,21435.9070146398757210922863175,
                    40503.4368880245757655115644924,77959.1576939931887059805749901,149029.97143311761190849030483,
                    286665.261854106286029833401311,550967.398509668882633833372138,1.0617876221324772532673959342e6,
                    2.0476410829595642366251137843e6,3.95509878589456720322553479001e6,7.64703477579053537420903420035e6,
                    1.48027481145283125343377300746e7,2.86824913543887697748518946862e7,5.56314364597939925813318232101e7,
                    1.07997159643147735301799899906e8,2.09835394276963986500763554968e8,4.08032485121803042928546125573e8,
                    7.94042122616032929798576541022e8,1.54634567461485622511676782385e9,3.01347999610120314953504960931e9,
                    5.8764292252599662559387199551e9,1.14664368774695897196135522866e10,2.23872269172695857861921353765e10,
                    4.37337818031172191931906704918e10,8.5480589476614799825190403211e10,1.67163448882140526636470185826e11,
                    3.27060981933869178579781740472e11,6.40207320716451917251894046147e11,1.25374286243082603890488740419e12,
                    2.4563159350103347820185601051e12,4.81437980438997407600709337155e12,9.43995106266995875107975818965e12,
                    1.85167907285555520477292852358e13
                  },-1),
                  // AQqPS_reg as^3 LM^0 NF^0 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    295.166793174664762128675058891,1677.66605977860931081875625994,713.652908276707883682579613714,
                    1011.79621241575743813500835581,172.014729268223162024135366535,122.748155418777481583390053829,
                    -37.3004622159903985938959933835,58.3111803170903136671175754239,17.1469856355863775862998743548,
                    90.8386718388270820860165349193,117.47845431451624271687096121,249.441320927311437621352828658,
                    427.480596973950825268155256281,812.065416237019442090133836751,1485.45766648089739992658552532,
                    2785.93803882017209804918216461,5201.04060537263521605689215511,9789.83910781285332269948307817,
                    18456.5454376931274690634669277,34938.4611334761237536914556986,66299.2983450846711433359257104,
                    126160.368532642023364725216138,240608.82811241524244368892337,459882.50293338963119391541831,
                    880685.165062891564038449708813,1.68957633343469601962487624607e6,3.24676132852769957701795412387e6,
                    6.24866475212712465066072254624e6,1.20431028944604401009404871813e7,2.32413446889293316072731909802e7,
                    4.49072487221849976656920495351e7,8.68698372102966478144279950087e7,1.68223644693987740118263992189e8,
                    3.2609317657322822678163599176e8,6.3271295768593335565940505052e8,1.22873466704020740984261254311e9,
                    2.38821763404704959399881655202e9,4.64551595672451555195985742519e9,9.0431503734409326431413608821e9,
                    1.76162805019113020383279623395e10,3.43401018614847412961330182679e10,6.69834086145296668013519672156e10,
                    1.30736747851968012883609462194e11,2.55317337565124121687334093897e11,4.98888430230168743265705128618e11,
                    9.75340250390536304560043111009e11,1.90777990011372834988459449866e12,3.73344945916326706168029810439e12,
                    7.30959398246834459197141957081e12,1.43175544023003888668537559767e13,2.80562219849007360690490464369e13,
                    5.50005351401392862067413974679e13
                  },-1),
                  // AQqPS_reg as^3 LM^0 NF^0 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    -334.74793693608868690183406853,746.932659023622843578212367475,-614.045574208264839980762538622,
                    -203.341563786008230452674897119,-62.5003292181069958847736625514,8.95583539094650205761316872428,
                    -20.4694587497550460513423476386,42.1356952892856143045107519786,30.181582834649308152289591779,
                    120.641370816730899035425784397,172.551927707049293281990625799,374.064530209989108740936785694,
                    638.491207814867567012167641617,1230.87596541764970676989332598,2229.97527823493173297350940356,
                    4197.92844028191195662172656959,7790.38707017525539286095360377,14656.3981265825538843227699282,
                    27532.9971782305963007017984273,52038.9394225846950401489770863,98517.3269292150963457375503258,
                    187168.790060507957392254635033,356352.414308225397039598635603,680164.53674394121160884878139,
                    1.30080339405409323290828219342e6,2.49268235730708017060653429985e6,4.78491611517624099524906676152e6,
                    9.20010654296982709031069597159e6,1.7715722010733166138696860257e7,3.41609399477600815697032470363e7,
                    6.59569113005282297046078749054e7,1.27500950812332818831773913134e8,2.4674814523851518483712313388e8,
                    4.78025309629739252503105493344e8,9.26991831482256170341119659121e8,1.79929881642138974134551878692e9,
                    3.49550163870124477273985887064e9,6.79630861298553261887106312288e9,1.32243494198394039401676309536e10,
                    2.57510575767593847229784948244e10,5.01786129455986222099601368962e10,9.78429738101140491035549247583e10,
                    1.90903448589664653462523322729e11,3.72698641518703368885966805232e11,7.28030307795067396353926723131e11,
                    1.42290951363706085788826352931e12,2.78247081172673743234525903653e12,5.44376670515252494281720006623e12,
                    1.06555196806478267030607908279e13,2.08663884953039059101998152113e13,4.08798494176903889602060813121e13
                  },0),
                  // AQqPS_reg as^3 LM^0 NF^0 [below 0.25] y=x H[0, y]^3
                  ome_x<double>({
                    86.5701811790025987733317719502,-37.7754978333430802390139070621,60.3347050754458161865569272977,
                    27.4348422496570644718792866941,-1.15610425240054869684499314129,4.41569272976680384087791495199,
                    -0.345145992553399960807368214776,1.83336412754402172391590381008,-0.168728945624889187499416776313,
                    1.01525028850131730790167004159,-0.100858239459062504329994041928,0.647993137510533694596545326063,
                    -0.0673320444364215407986451757496,0.450525521750592975664200735426,-0.048230724339800448876557952667,
                    0.33178742088794998847908900819,-0.0362856335950017867229196205231,0.254700787762857906819630942831,
                    -0.028306646136195874525636732855,0.201768485108490513203610781481,-0.0227075956442257184365867379444,
                    0.163828286980612238200653683806,-0.0186250564484247065755734220924,0.135694591277515669513785519463,
                    -0.0155554329733423761289055694993,0.114249392809135160181859054645,-0.0131885451867724738886498937148,
                    0.0975246304246339419214098261617,-0.0113247048946355867459297715559,0.0842275406748567601699388904369,
                    -0.00983055774981527591565422730177,0.0734801214047635466363031705115,
                    -0.00861425352431056366000675554767,0.0646689017090533839299164914585,
                    -0.00761083497560144071392658765814,0.0573546416578705747096414910381,
                    -0.00677330134466075262146777172541,0.0512160899027297589602584656118,
                    -0.00606695787934274805263984655749,0.0460138580680601499630628909066,
                    -0.00546573823197736535385431442543,0.0415665825611404741992795840814,
                    -0.00494974928056773393229080666327,0.0377348156063265834055009664193,
                    -0.00450359471057406484024081184314,0.0344099064218961465024310458013,
                    -0.00411520709968715326139856185908,0.0315061797481395940579268103515,-0.0037750194379993652539003380652,
                    0.0289553390102478363247040863309,-0.00347536776699425886084676585004
                  },0),
                  // AQqPS_reg as^3 LM^0 NF^0 [below 0.25] y=x H[0, y]^4
                  ome_x<double>({
                    -13.6790123456790123456790123457,15.3580246913580246913580246914,9.21810699588477366255144032922
                  },0),
                  // AQqPS_reg as^3 LM^0 NF^0 [below 0.25] y=x H[0, y]^5
                  ome_x<double>({
                    1.18518518518518518518518518519,-2.54814814814814814814814814815
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^0 NF^0 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^0 NF^0 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^0 NF^0 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -9.00294907368899802725452883306,47.5548863266535362286990015128,832.534734486897607285260129405,
                    2470.29806337385411807977835362,4854.04727261027302129447323515,10011.0743560318416656774828131,
                    19058.7288332551684265046076982,37096.2899280883568466530961012,69733.9415961993547738752836726,
                    132649.887967804937363918104474,247235.241014886188037682386192,464199.210843145478287995156383,
                    858683.055298176957893674187081,1.59523693923357980414746793246e6,2.92575435094569437918971013275e6,
                    5.37509017414248931774901408086e6,9.7498719001921497617050588585e6,1.76663016719253870219021199494e7,
                    3.15544229527018123103207970327e7,5.60918192854854152581695504247e7,9.78942954147764875869126720357e7,
                    1.68982998294427065496717359021e8,2.83811465778742262183761488049e8,4.65235475645964564981888165549e8,
                    7.24761074168088516254958227109e8,1.0581438366886323169019897612e9,1.33306893117975851440948303542e9,
                    1.1683072925969101933687573554e9,-6.52102282996123664415801168786e8,-7.04445278892419358020303614623e9,
                    -2.55178898976247020455448517641e10,-7.30666706536744761427835798132e10,
                    -1.89925103950463511697614398867e11,-4.64527737269233137499992623171e11,
                    -1.0971599757460665421798282191e12,-2.52024970073227494000011843977e12,
                    -5.68691368441193906614082543325e12,-1.26300768552650285242392846234e13,
                    -2.77497759826569754988779412813e13,-6.03479257013015137363785997853e13,
                    -1.30298637715569343192805815363e14,-2.79332343180364460042188453497e14,
                    -5.95754692726814527461968144626e14,-1.26399036958804183760876216071e15,
                    -2.67142095762523152629964185165e15,-5.62356233094185394993751054809e15,
                    -1.1802519218380556749878472888e16,-2.46934114558390960462632646194e16,
                    -5.15396770811855235883126698416e16,-1.07303215426211549962534326736e17,
                    -2.22960126703502499335921171865e17
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^0 NF^0 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^0 NF^0 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^0 NF^0 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    -190.424439446550414781973926127,90.2743445019486368782422859915,-96.1487996303300003332095101537,
                    129.617592281080981843741723973,112.01820919887468555619174041,119.264346928043486077432749835,
                    124.495677679716005576153285159,122.122416317907918624040971254,109.835968757694332739744733821,
                    82.165915911192365009574942408,27.8427094396778461764221459965,-75.2446576280192604640508999074,
                    -269.703443833537069514216275244,-637.198763214277757964114864338,-1334.14562574116318432358722383,
                    -2660.38859035491506725499042369,-5191.50102650921535949047372203,-10033.9428345156488744379504689,
                    -19317.7546818594737010811097743,-37149.0059594428525654942179319,-71452.9975584386546657976346645,
                    -137544.903565683162657151996145,-265053.42649532850912018900235,-511358.255848942029407279259855,
                    -987691.823887962705498358501252,-1.90988196172300954239735616712e6,-3.69707279787603618745296624384e6,
                    -7.16392784354521595628829711325e6,-1.3895094912716944395886585873e7,-2.69752463112588685673587849023e7,
                    -5.24132218168522320638628292188e7,-1.01921862939878026962421868351e8,-1.98347150494653035020188108431e8,
                    -3.86277300016973738057373665383e8,-7.52784783207727175514314344142e8,-1.46799953369546459041217418665e9,
                    -2.86451141463973919700951422885e9,-5.59283304968218141015469447998e9,
                    -1.09259150529532389082983682236e10,-2.13558704169129753200200692682e10,
                    -4.17637802982464336587695682322e10,-8.17138462850932393785532515708e10,
                    -1.59954223133451305580047855036e11,-3.13250251586806440834228448151e11,
                    -6.13726370570461827294462582402e11,-1.20292413925710844819311832139e12,
                    -2.35871083163505738346003418277e12,-4.62676546563466565430231685785e12,
                    -9.07904749495739561411640642485e12,-1.78220294506559971553217025362e13
                  },1),
                  // AQqPS_reg as^3 LM^0 NF^0 [above 0.75] y=1-x H[0, y]^1
                  ome_x<double>({
                    -161.298283154550292495751046434,60.8221863861425907262344526424,-54.3961455872389894168235010346,
                    -113.300208026652122256627381561,-148.022876961205634784340408426,-188.145646888093504439391398063,
                    -229.787279440638877191747602956,-279.363467148685529017547934819,-348.750232145538389393468666684,
                    -458.54878370940356774706775874,-646.609698885678756902667703883,-983.866520264866802393202685708,
                    -1604.39168290543164765401975165,-2762.76784574006551479576484498,-4943.91015462785953537843100943,
                    -9073.71663893276736320642136897,-16923.8019395910705741571755825,-31890.3145532676078911322281049,
                    -60494.6840809315273212339891553,-115279.245392248003474997404395,-220400.974199409284845225787584,
                    -422449.642704440713766868397438,-811392.112909412407509057481547,-1.56116165402434498633271849643e6,
                    -3.00839262349603989534789558787e6,-5.80528457617956048954075328984e6,-1.12166380980776330529659689945e7,
                    -2.16974936986276100227271303375e7,-4.20172971666749304269062948637e7,-8.1449240737527373556033376817e7,
                    -1.58036949094616960588353422893e8,-3.06914943449763663932619545985e8,-5.96544329320851904757675869747e8,
                    -1.16041195061516648091592852034e9,-2.25895551609430316924306786065e9,-4.40060043235586938607633682949e9,
                    -8.5784550276725186087376879247e9,-1.67334089168606255204458780283e10,-3.2660616528678560069082405581e10,
                    -6.37846720115658684188884115619e10,-1.24637445214399362866966068973e11,
                    -2.43674502041214545942213772392e11,-4.76640320711817459057097846194e11,
                    -9.32783988410034614455269610095e11,-1.82629984530461372326978812255e12,
                    -3.57730044892523566663992605284e12,-7.0100842791452903310907707762e12,
                    -1.37425768151808352162301344816e13,-2.6951520833836730572360843706e13,
                    -5.28764158213760990107002831367e13
                  },1),
                  // AQqPS_reg as^3 LM^0 NF^0 [above 0.75] y=1-x H[0, y]^2
                  ome_x<double>({
                    -11.074606972731915198139803358,2.07407407407407407407407407407,-15.4202859850775942104854823704,
                    6.03238891204174735330052586007,6.69560846013861218277268330643,2.93825187654823421389093128778,
                    -6.97097497734159989153145083698,-27.5420233906885284974820334046,-67.1404931500439308352459144889,
                    -141.556651100552170655018637,-280.431630740913269293984257482,-539.513335124730198023554248852,
                    -1023.96934259192301036333357622,-1932.79998393197125613418667657,-3643.62040134519340491888843788,
                    -6874.82379686818820063920413709,-12996.3220532317676996812536806,-24626.0060670116044450785818081,
                    -46776.4427434589550938848890194,-89062.6270919181237994516798467,-169958.381024110587916681591303,
                    -325012.815935658941535108488355,-622730.022620001810699749812543,-1.19529330615290269953979433912e6,
                    -2.29807277759595396737235811833e6,-4.42499427864313327789255424368e6,-8.53241445578722032174916482034e6,
                    -1.64739316306359700670615856954e7,-3.18455711602871109654279793787e7,-6.16299348673497308547741892091e7,
                    -1.19396943003481168492809283641e8,-2.31539226023770705174435107395e8,-4.49427416329618642858214377516e8,
                    -8.73120693165888119933597148679e8,-1.69764526700182489672450019574e9,-3.30337306163409141340248637797e9,
                    -6.43262176532765919741428414481e9,-1.25349122333173367509734669811e10,
                    -2.4442293165296651411731894465e10,-4.76909118309758927165116492926e10,
                    -9.31084384575315068478385244597e10,-1.81882061720727561072038899434e11,
                    -3.55489383602475095540382008886e11,-6.95166221620875322574788077975e11,
                    -1.360084771811149718155931527e12,-2.66225264970102521076555055461e12,
                    -5.21350517412292757007217758894e12,-1.02140839687092833119438069332e13,
                    -2.00193713688546668507316732424e13,-3.92532507758470609000596965391e13
                  },1),
                  // AQqPS_reg as^3 LM^0 NF^0 [above 0.75] y=1-x H[0, y]^3
                  ome_x<double>({
                    -0.888888888888888888888888888889,-1.28395061728395061728395061728,-1.48148148148148148148148148148,
                    -2.35390946502057613168724279835,-2.92455418381344307270233196159,-3.38655692729766803840877914952,
                    -3.78155146832395509644186892864,-4.1273312617227961143305058649,-4.43457661070829795109630500577,
                    -4.71067563742460861802425588434,-4.96107437476685512600215181584,-5.18992997691089722515087273785,
                    -5.40049112717678719578054810723,-5.59534284434569334854235139135,-5.77657312965954941263583238892,
                    -5.94589050638433354482737198787,-6.1047090529853569710035172037,-6.25421116097351733218280272832,
                    -6.39539463550914032409595649768,-6.52910856544392856219320671933,-6.65608099911341675883750023347,
                    -6.77694055119317621979505446342,-6.89223345409584262458568638662,-7.00243714810741059660242849941,
                    -7.10797121215241394917638667109,-7.20920623017418857544890159818,-7.3064710396510249322019062608,
                    -7.4000587008988446487730431296,-7.49023144652898285854473968797,-7.57722481153337635678268829871,
                    -7.66125110027898946381347791565,-7.74250231322600687227569323363,-7.82115263061518218370929404372,
                    -7.89736053067317902284962506784,-7.97127060459366246125920492875,-8.04301511859351871007663721738,
                    -8.11271536392641277171876966065,-8.18048282827101051468747180867,-8.24642021595685607268809993817,
                    -8.31062233971323942162159353913,-8.37317690277123988564342239597,-8.43416518702165172179785049938,
                    -8.49366266038135140716275306138,-8.55173951443112567895523486666,-8.60846114166787061620171782203,
                    -8.663888560291772766375395781,-8.7180787932680224797020064699,-8.77108520741787058243865923798,
                    -8.82295781746958026353547424375,-8.87374355930730827315050193473
                  },1),
                  // AQqPS_reg as^3 LM^0 NF^0 [above 0.75] y=1-x H[0, y]^4
                  ome_x<double>({
                    -0.493827160493827160493827160494,0,-0.493827160493827160493827160494,-0.493827160493827160493827160494,
                    -0.51028806584362139917695473251,-0.526748971193415637860082304527,-0.540858318636096413874191651969,
                    -0.552616108171663727219282774838,-0.562414266117969821673525377229,-0.570644718792866941015089163237,
                    -0.57762813318368873924429479985,-0.583613916947250280583613916947,-0.588793922127255460588793922127,
                    -0.593316148871704427259982815538,-0.59729570840681951793062904174,-0.600823045267489711934156378601,
                    -0.60396998305495037521181312031,-0.606794157992415073025094811587,-0.609342285755541116164897841311,
                    -0.611652588260775395278319254927,-0.613756613756613756613756613757,-0.615680615680615680615680615681,
                    -0.617446607784772035979765448364,-0.619073179459652889604580425836,-0.620576131687242798353909465021,
                    -0.621968977524533080088635644191,-0.623263339312722028771411487461,-0.624469266444575086550395192371,
                    -0.625595491495874637637089744369,-0.62664963814389101745423584504,-0.627638391079251294305057745918,
                    -0.628567635736094517456524624983,-0.62944257388701833146277590722,-0.630267819810303470434189388438,
                    -0.63104748072068333506242002974,-0.631785224377816970409563002156,-0.632484336188039891743595447299,
                    -0.633147767650691627299814434317,-0.633778177637826760633778177638,-0.634377967711301044634377967711,
                    -0.634949312456087523838201344976,-0.635494185629686984700534836036,-0.636014382784408624201905855653,
                    -0.636511540904305762186899137804,-0.636987155505674024192542711061,-0.637442595574640663207490904753,
                    -0.637879116653913652577441250481,-0.638297872340425531914893617021,-0.638699924414210128495842781557,
                    -0.63908625178466448307718148988
                  },1)
                },0)
              )
            }
          ),
          // AQqPS_reg as^3 LM^0 NF^1
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^0 NF^1 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^0 NF^1 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^0 NF^1 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    -6.81002777386535765764427952013,83.4378680256578251381482321624,-99.1324370689804590584181875146,
                    22.6716962368658377154579315069,-3.55291071328077680434506921269,0.601329873207281993929953554987,
                    0.505861120642802855850810110571,0.355941410893315994995901280819,0.257345960564917347696881991439,
                    0.193517245205335624240967874115,0.150579610339223292708420876963,0.120493782504060944153248405712,
                    0.098646758966735417684842449363,0.0822964259959442844712149008988,0.0697447580362489867398027993335,
                    0.0598992752925464839262341935954,0.0520328014624871034653690011137,0.0456466844992287642638330621444,
                    0.0403899521280835018913552874792,0.0360098573252313742429758151704,0.0323206748532880305365191042663,
                    0.0291834600793266707800532013202,0.0264925889418015981672891318246,0.0241666056732327221817943019729,
                    0.0221418735309682782007405247968,0.0203680894663475369955573363517,0.01880506294973169142282187948,
                    0.0174203676980821324603099630328,0.0161876061108360728909846630147,0.0150851102933588238646383656576,
                    0.0140949585039657614274185951535,0.0132022224121087680344039699901,0.0123943852577757882531550545297,
                    0.0116608879457785050558031236684,0.0109927718912552815672468484564,0.0103823957318872406368088576791,
                    0.00982320893803406960656103516148,0.00930956961603120589940663093726,0.00883659690566364096164021054361,
                    0.00840005065724622856964549889355,0.0079962327695692312821918871499,0.00762190583983649040668415596135,
                    0.00727422573544535442913434123414,0.00695068542689497815889965727906,0.00664906798017884290945495440708,
                    0.0063674070384961434764189162734,0.00610395345832636586155384983493,0.00585714702696798699493497467734,
                    0.00562559239473542506533523386822,0.00540803851801301385986068787116,0.00520336103899023995500842956951,
                    0.00501054713151891432735709707651
                  },-1),
                  // AQqPS_reg as^3 LM^0 NF^1 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    50.5356811628617646416181293192,-11.5749860594735051709278246558,2.51058399169418226488705948991,
                    3.42386831275720164609053497942,0.519176954732510288065843621399,0.148345679012345679012345679012,
                    0.0439029142521206013269505332997,0.00404084032201946034145580630841,-0.0133247616114837077296134993288,
                    -0.0212047379537091471247849848673,-0.0246172857500953273493490478152,-0.0257836563780855179692033074341,
                    -0.0257814039085965812437993455629,-0.0251673885623981112538139556705,-0.0242422690582914567314375890009,
                    -0.023173119705243665939055205873,-0.0220542299056560970476863631511,-0.0209388841703807290195632537567,
                    -0.0198566873411289990140568223214,-0.0188233349001399883770294797773,-0.0178462749568190767311660208299,
                    -0.0169280611899892466983737303668,-0.0160683716701878061309257787709,-0.0152652388313202332271139851577,
                    -0.0145158039603767723381023458865,-0.013816780561257086033516584614,-0.0131647372550843552116568208694,
                    -0.0125562677994484185496819642252,-0.0119880901082932331690455050971,-0.0114571005453457076394233132579,
                    -0.01096040013474602484942689395,-0.0104953033105950922950232525648,-0.0100593360165915311685421970391,
                    -0.00965022753161044291489190324566,-0.00926589882792725778425876825854,
                    -0.00890444925134991353931093255066,-0.00856414264995450006609836460197,
                    -0.00824339364573932665803991725623,-0.00794075446149267277071900018817,
                    -0.00765490253191286438621220305718,-0.00738462900989398087002961034741,
                    -0.00712882820381452897732460701673,-0.00688648793505169075260671582254,
                    -0.00665668077710454469980088764119,-0.00643855612219106550825660095771,
                    -0.00623133301367376301145474040566,-0.00603429368028482187885384260033,
                    -0.00584677770893847031054098650287,-0.00566817679565486090279367669305,
                    -0.0054979300179288190383310865602,-0.0053355195762041031861752041347
                  },0),
                  // AQqPS_reg as^3 LM^0 NF^1 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    21.0986174691734849992437130534,-1.41990104934503351927480546512,-1.58024691358024691358024691358,
                    1.58024691358024691358024691358,0.483950617283950617283950617284,0.239012345679012345679012345679,
                    0.144479717813051146384479717813,0.0975560594608213655832703451751,0.0706307214243722180230116738053,
                    0.0536547129139721732314324906917,0.0422197281456540715799975059234,0.0341298506955072611638268203925,
                    0.028185506973385761264549143337,0.023683818555613427408299203171,0.0201896055375908855762335615815,
                    0.0174211678973583735488497393259,0.015189334993256561884012864405,0.0133630939667079613254161504834,
                    0.0118492937314308311967781672416,0.0105801841044660856488326780845,0.00950550515339681236880375076743,
                    0.00858731725899645699144446011614,0.00779653525748010972930980459684,0.00711055210682867210334701885369,
                    0.00651157659349215744930930816249,0.00598544944052190429002023204922,0.00552078635725359372083018806666,
                    0.00510834855880756482212639224352,0.00474057414194047817271478809317,0.00441122491088203360693278843126,
                    0.00411511719759709224376223468446,0.00384791455969706594654266499925,0.00360596659157186385767724650833,
                    0.00338618246947267705130967738063,0.00318593092001920774071345529875,0.00300296047645968006268989407343,
                    0.0028353354466021755159632532804,0.00268138414706034957865736012304,0.00253965678670161546389574972174,
                    0.00240889099472805701941293957772,0.00228798344594401340393901534088,0.00217596637981329633124199850355,
                    0.00207198807093533048328452371951,0.00197529650811049112168026934398,0.00188522569279844574707754914364,
                    0.00180118408686530075682809628263,0.00172264483242718043500404915031,0.00164913743952883163891879359498,
                    0.00158024069497647847837688570171,0.00151557659135766366500605484953,0.00145480511176503859291259965816
                  },0),
                  // AQqPS_reg as^3 LM^0 NF^1 [below 0.25] y=x H[0, y]^3
                  ome_x<double>({
                    3.22633744855967078189300411523,0.855967078189300411522633744856,-1.58024691358024691358024691358
                  },0),
                  // AQqPS_reg as^3 LM^0 NF^1 [below 0.25] y=x H[0, y]^4
                  ome_x<double>({
                    0.493827160493827160493827160494,0.493827160493827160493827160494
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^0 NF^1 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^0 NF^1 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^0 NF^1 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    2.56898329021512147545189358853,4.35687692540311668456103012534,-31.7843565937042448041794683949,
                    -75.5667956692759986998160215528,-115.746182569028318145972282983,-256.248259884302564003139298714,
                    -485.727476020352501360076832498,-1031.20869711815263541864797451,-2044.02029777572665129326537934,
                    -4255.1822816630500538455713715,-8519.090850054323241380451787,-17535.924305844045552889075441,
                    -35206.3583583608434699240579654,-71972.7799837872300966467141021,-144610.294261012455450864260238,
                    -294301.24641198619906507668135,-591378.640958886902347297964265,-1.1997851941118451620892290598e6,
                    -2.41052582542160608348610993984e6,-4.87938658731595388271964818712e6,-9.80099210787745113442352320768e6,
                    -1.98052897474423316638501739988e7,-3.97715297975301855145133790798e7,-8.02614025048745080845917054556e7,
                    -1.61133395504809502260647403944e8,-3.24833515292984762872424717897e8,-6.51978922725184762034797791323e8,
                    -1.31321474287216828179216617566e9,-2.63517984480836268930702842264e9,-5.30400741847195946656141053808e9,
                    -1.06411506496754846568798276235e10,-2.14054282421306847223395541347e10,
                    -4.29364560389972021902668326314e10,-8.63260826205896144182067535264e10,
                    -1.73128824303224606186849876574e11,-3.4793398636044750737597768423e11,
                    -6.97679363696511118281774853789e11,-1.40158732848506368238393965938e12,
                    -2.81006789956023957555446137824e12,-5.64337018899293526407551108151e12,
                    -1.13130094425143637962047720722e13,-2.27129953161457114597757216782e13,
                    -4.55262364234274831354377530094e13,-9.13790961947081373364239997237e13,
                    -1.83141256986328582583849288616e14,-3.67512602899146280399459669985e14,
                    -7.36490895739809492615437798579e14,-1.47762647450785934690717601923e15,
                    -2.96086800770665024580072767485e15,-5.93931556194900595203374356933e15,
                    -1.19001681469185770888394235342e16
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^0 NF^1 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^0 NF^1 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^0 NF^1 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    19.4155304255817812905077804577,-4.14814814814814814814814814815,-1.10269216708617636230063330346,
                    -3.98792474683281759915694339395,-3.68200690416932263798273268426,-3.51656587289003112607927337759,
                    -3.47506285544555335593180683316,-3.4937015789519623287290579104,-3.53943935594243677934199984348,
                    -3.59647224284919558677018100993,-3.6572044767268558161567240272,-3.71796245690720132184404355767,
                    -3.77699709826261603359702862094,-3.83352711746716758020267996521,-3.88726528679432024678118584103,
                    -3.93817668034416856903478366926,-3.98635234447924469961786256449,-4.03194223348657407836515080021,
                    -4.07511939356167103657287702383,-4.11606097310530711271098452427,-4.15493843083505039336578680221,
                    -4.19191281309321188335896769911,-4.22713282402123457144726328586,-4.26073441580457752477339049726,
                    -4.29284118085269294159067219153,-4.32356513945796887440329863593,-4.35300769391671337982941334067,
                    -4.38126062205148416837803784206,-4.40840704194907546232469705735,-4.43452231371520102968068030332,
                    -4.45967486352269617088057657706,-4.48392692619498302054597249261,-4.50733520852930373810895019555,
                    -4.52995147859515170841295182133,-4.55182308757192075247567927524,-4.57299343104878141777216123977,
                    -4.59350235653585836578769388981,-4.61338652348797633052004779371,-4.63267972157257992966656901873,
                    -4.65141315230931941228181169464,-4.66961567861869554230005692137,-4.68731404626620889674327203241,
                    -4.7045330806879892615165810461,-4.72129586223721435939413751022,-4.73762388249661819929874691732,
                    -4.75353718395755524186563864436,-4.76905448506587982598852445773,-4.78419329237445067489772162996,
                    -4.7989700013166021102831186268,-4.81339998691997532486941786493
                  },1),
                  // AQqPS_reg as^3 LM^0 NF^1 [above 0.75] y=1-x H[0, y]^1
                  ome_x<double>({
                    13.0141813477931493095651580575,-3.55555555555555555555555555556,9.29401673873965136717832678175,
                    9.76315254120878716964746258422,10.841885883803263888745450706,11.4755163457393043938516692805,
                    11.8791731078076508139142791186,12.1621926947416532791088072213,12.3759824459301996952457088603,
                    12.5465640732601483497037044733,12.688207798984992341522462239,12.8093167430137682748326960279,
                    12.9151516074410967621561629136,13.0091829602397369108487501697,13.0938053748940107498150517182,
                    13.1707356969620455001192533305,13.2412460406449555776358390619,13.3063059763552415282563727485,
                    13.3666726137815294507173666294,13.4229496132734830786933584071,13.4756270173937299633180862936,
                    13.5251088695352757780425315108,13.5717328341660143706030821826,13.6157844434451040683809353784,
                    13.6575076486966315457502142423,13.6971127764066393589931644197,13.7347826253761321853348000888,
                    13.7706772086684995054685365145,13.8049374912457137777979100379,13.8376883720664683887264103301,
                    13.8690410898877881201900716152,13.8990951838604116071354036128,13.9279401061327181052175159879,
                    13.9556565594922977889151023604,13.9823176155676396916871232162,14.0079896562763542055624335454,
                    14.0327331716808215343746653429,14.0566034402631688979494285407,14.0796511122090564327891162996,
                    14.1019227121358874976507883789,14.1234610744891834152683595294,14.1443057223255107741145127284,
                    14.1644931982299437974779600966,14.1840573545543044923416242422,14.2030296089155225789469098928,
                    14.2214391698910528596498195562,14.2393132370371916006771643647,14.2566771786958017125154727695,
                    14.2735546905142447346610756695,14.2899679371581264488226226541
                  },1),
                  // AQqPS_reg as^3 LM^0 NF^1 [above 0.75] y=1-x H[0, y]^2
                  ome_x<double>({
                    0.987654320987654320987654320988,0,-0.296296296296296296296296296296,-0.740740740740740740740740740741,
                    -0.802304526748971193415637860082,-0.824362139917695473251028806584,-0.841454606533971613336692701772,
                    -0.858167464516670865877215083564,-0.874773009340910575478476713045,-0.890979162584100855705793977399,
                    -0.906539177246247953318660389367,-0.921315771147420979070810720642,-0.935254763174991095219015446936,
                    -0.948354011113473872936632399392,-0.960641083392406143728895051646,-0.972159048084974010899936825863,
                    -0.982957877200476200860668373164,-0.993089415730280354997344189535,-1.0026045385639187251347453481,
                    -1.01155162811346614507768165907,-1.01997583955223151252440160986,-1.0279188289740517011010535374,
                    -1.03541874747218285428907394731,-1.04251038197063793320270905501,-1.04922537108779013626383543076,
                    -1.05559245328831508195584988656,-1.06163772230288392655150535985,-1.06738487563358092221481895608,
                    -1.07285544856984826307950776039,-1.07806903013523635140911308034,-1.08304345976519334217514225048,
                    -1.08779500489627425769538892898,-1.09233852041018776442864015805,-1.09668759126087278899561155337,
                    -1.10085465976873801553563648631,-1.10485113908682487363965336668,-1.10868751428735327128751742122,
                    -1.11237343242010135675690743697,-1.1159177827785696340124280289,-1.1193287684892431750870049816,
                    -1.12261397042126516405318980236,-1.12578040430273938674038339886,-1.12883457182778321713582446698,
                    -1.1317825064461243128303727588,-1.13462981444445479493780542489,-1.13738171185546193561453880899,
                    -1.14004305766574898954327651268,-1.14261838373694978374077380913,-1.14511192180441322361980131341,
                    -1.1475276278741037168813349622
                  },1),
                  // AQqPS_reg as^3 LM^0 NF^1 [above 0.75] y=1-x H[0, y]^3
                  ome_x<double>({
                    0.197530864197530864197530864198,0,0.197530864197530864197530864198,0.197530864197530864197530864198,
                    0.204115226337448559670781893004,0.210699588477366255144032921811,0.216343327454438565549676660788,
                    0.221046443268665490887713109935,0.224965706447187928669410150892,0.228257887517146776406035665295,
                    0.23105125327347549569771791994,0.233445566778900112233445566779,0.235517568850902184235517568851,
                    0.237326459548681770903993126215,0.238918283362727807172251616696,0.24032921810699588477366255144,
                    0.241587993221980150084725248124,0.242717663196966029210037924635,0.243736914302216446465959136524,
                    0.244661035304310158111327701971,0.245502645502645502645502645503,0.246272246272246272246272246272,
                    0.246978643113908814391906179346,0.247629271783861155841832170335,0.248230452674897119341563786008,
                    0.248787591009813232035454257676,0.249305335725088811508564594984,0.249787706577830034620158076948,
                    0.250238196598349855054835897748,0.250659855257556406981694338016,0.251055356431700517722023098367,
                    0.251427054294437806982609849993,0.251777029554807332585110362888,0.252107127924121388173675755375,
                    0.252418992288273334024968011896,0.252714089751126788163825200862,0.25299373447521595669743817892,
                    0.253259107060276650919925773727,0.253511271055130704253511271055,0.253751187084520417853751187085,
                    0.253979724982435009535280537991,0.254197674251874793880213934414,0.254405753113763449680762342261,
                    0.254604616361722304874759655121,0.254794862202269609677017084424,0.254977038229856265282996361901,
                    0.255151646661565461030976500192,0.255319148936170212765957446809,0.255479969765684051398337112623,
                    0.255634500713865793230872595952
                  },1)
                },0)
              )
            }
          )
        },0),
        // AQqPS_reg as^3 LM^1
        ome_nf<double>({
          // AQqPS_reg as^3 LM^1 NF^0
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^1 NF^0 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^1 NF^0 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^1 NF^0 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    1280.03548460689505488484005947,-1750.22859379067859853845854797,2382.00351626989170017765373851,
                    -1668.09792991314222892098356701,-187.654320987654320987654320988,-39.9297180023593364925568068134,
                    -0.395976685856216907386217480246,-9.76108938089916221976604766857,0.311206141801701578325825485206,
                    -3.7224097725312638014505048882,0.208225677035905484120005443057,-1.82111681981566042416325296273,
                    0.141257307501320048172918847247,-1.02833925644193207309445893501,0.103417792789054437945743901152,
                    -0.636507653884118869405380066935,0.0805011445707490286008675303063,-0.419716509404798244895284939889,
                    0.0655364111285738231324733482427,-0.289765952499791878488467280802,0.0551291150954218565912935292775,
                    -0.207072138747375291097938383055,0.0475164322571094137052099638009,-0.151958756300762003698156291779,
                    0.0417189355201355784709489578895,-0.113849488356673764955490216036,0.0371586762326039471495968199559,
                    -0.0866967625164512169576367799907,0.0334764346910276479662868796258,-0.0668647014395004482683709878642,
                    0.0304387735506696749442470317178,-0.0520747112398120291595711041154,0.0278880761924856919478211095887,
                    -0.0408480263979624779647034587281,0.0257143443268267652601399937835,-0.0321959072453307690139647385605,
                    0.0238385780066568727493958872525,-0.0254400148682423444762634279297,0.0222026024660341616375132398005,
                    -0.0201043962242411988142919031194,0.020762628361024897130460504377,-0.0158484383019883547064268622161,
                    0.0194850516827530735657491690987,-0.0124240772636414658416548726791,0.0183436402934982316304842132164,
                    -0.00964780652118426634491051300183,0.0173176035996650355842753866349,
                    -0.00738195886405131493256489586887,0.0163902392381131039337158239524,
                    -0.00552194084719308197848321458821,0.0155479655820088814654834831116,-0.00398737037825166595289711270734
                  },-1),
                  // AQqPS_reg as^3 LM^1 NF^0 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    265.481481481481481481481481481,1028.8814683484900969240883771,45.7343411023447991224624653356,
                    1309.89928331430218538051486893,74.6172839506172839506172839506,75.3466666666666666666666666667,
                    0.368691358024691358024691358025,27.2732422944486436549928613421,-1.07450071386579323087259595196,
                    14.5899570880674433206557167253,-0.960894898425762623293487491018,9.26215164525377768924345804256,
                    -0.77350200711327126998997216322,6.47020121680525250019728605116,-0.623110924597756853819879113674,
                    4.80622215090927493705830550101,-0.510033421809222727880789395994,3.72697164601646852063358535228,
                    -0.424726013951176847297745996839,2.98361552913919362980512823885,-0.359269093627128493143652494239,
                    2.44802602679995450995304717444,-0.308083642672727597388303608171,2.04835243211138343020872711632,
                    -0.26733471270121008974685611282,1.74158873080287523065841281475,-0.234365842381945966485328559264,
                    1.50063201198785938988170587595,-0.207305725292401448730400493994,1.30767146987425186948333123242,
                    -0.184810968170665415869870069548,1.15058917501353750471077790614,-0.165899033245078040998452744094,
                    1.02089544075776359771331815114,-0.149838632572761918207481002308,0.912491582964756540489844753225,
                    -0.136076645284643637679158928081,0.820900452593222263344142455536,-0.124188509485674400757922501003,
                    0.742772341742308033991898053016,-0.113843942507059949061965599805,0.67555875980198584456885603511,
                    -0.104782844657699677227676383311,0.617291711356178976166912693252,-0.0967980866255204389945459866843,
                    0.566431083138087825355161519387,-0.0897230283579883431116341022273,0.521757061151172902777948156655,
                    -0.0834223420206200821039993046505,0.482292960967520899831967149799,-0.0777851767985788437719628370692,
                    0.447248996887097829043385550767
                  },-1),
                  // AQqPS_reg as^3 LM^1 NF^0 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    -262.982948403583241858834366516,239.432719513337781886060696594,-226.666666666666666666666666667,
                    5.33333333333333333333333333333,-13.8444444444444444444444444444,1.32,-4.1853968253968253968253968254,
                    0.598639455782312925170068027211,-2.05498866213151927437641723356,0.342592592592592592592592592593,
                    -1.23079685746352413019079685746,0.222148760330578512396694214876,-0.822533022533022533022533022533,
                    0.15581854043392504930966469428,-0.589558426921064283701646339009,0.115374149659863945578231292517,
                    -0.443723155929038281979458450047,0.088884083044982698961937716263,-0.34624559604830350899676462719,
                    0.0705858212783420539653226633836,-0.277812074044761025370443652992,0.0574149659863945578231292517007,
                    -0.227897430381902431591872585661,0.0476183036760455560936743270478,-0.190357831624372267094384296653,
                    0.0401333333333333333333333333333,-0.161406793776024545255314486084,0.034285435994837704239413641123,
                    -0.138605507151392938658036267961,0.0296294498774539542332985512873,-0.120324705489831491902758827695,
                    0.0258619493583073187651751647589,-0.10544187320657143758080282534,0.0227703168044077134986225895317,
                    -0.0931629692875367615852044917789,0.0202019631381964550526092790057,-0.0829133415127812886916528373111,
                    0.0180450720991261531802072342613,-0.0742687274094027107428320275517,0.0162161866828391439514388587476,
                    -0.0669103327906527274553343478,0.0146519928613920285544318857823,-0.060594834431796689272760508637,
                    0.0133037530858746960114494689168,-0.0551338932085281463323648288494,0.0121334557698194061830425466789,
                    -0.0503798788178522455694692887845,0.0111111016027404645542680270863,-0.0462157417601900786163922008533,
                    0.0102127585728168818547827294183,-0.0425477096961233472981029146352
                  },0),
                  // AQqPS_reg as^3 LM^1 NF^0 [below 0.25] y=x H[0, y]^3
                  ome_x<double>({
                    69.9259259259259259259259259259,-74.0740740740740740740740740741,-37.9259259259259259259259259259
                  },0),
                  // AQqPS_reg as^3 LM^1 NF^0 [below 0.25] y=x H[0, y]^4
                  ome_x<double>({
                    -8.,16.
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^1 NF^0 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^1 NF^0 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^1 NF^0 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    119.613330334452320809793446158,324.538721060636300715239571333,484.710119509581528523861104091,
                    1234.60683543550711251338667645,2225.57226408307528186467929622,3872.63557891955776752137271281,
                    6392.99141527513968946336330403,9968.52160584063238202893336471,14423.9874403306691715742969342,
                    18069.825077265508124461181209,15508.79240593837162337737808,-9058.56645954098209472316851011,
                    -94618.726997940544359771682326,-337941.936805903742678983789711,-960021.746267311569767115894192,
                    -2.47336063827693949200052387251e6,-6.00598032680899508719398469603e6,-1.40787151289715463141032967721e7,
                    -3.21222865703474091601806258093e7,-7.19946149793632238360787298107e7,-1.58887577770845615511892958297e8,
                    -3.46944781561345736806555502852e8,-7.50073169856353045043299376068e8,-1.61031162103559142563819997378e9,
                    -3.43319544302504329668270782945e9,-7.28375628554233431085435046602e9,
                    -1.53741874870654072567439079598e10,-3.23340839345768128154299840685e10,
                    -6.77376609515789199367215552852e10,-1.41515560471340922295788024151e11,
                    -2.94739689319792330066829230855e11,-6.1254196360710353423630990109e11,
                    -1.26984647569539595576858259746e12,-2.62793911624688280942282313717e12,
                    -5.42735896769836948915652538112e12,-1.11930462531629105589166059083e13,
                    -2.30442512710228645585174397939e13,-4.73879690639430081777329140907e13,
                    -9.73063692581267903464807379738e13,-1.99612087828034989539774714755e14,
                    -4.08968415859454463526910206556e14,-8.37200001445574151509750014859e14,
                    -1.71197842523491571543620514925e15,-3.49828581835478404448711682451e15,
                    -7.14168896699226311257906707789e15,-1.4570561733021879270423802615e16,
                    -2.97022336119773606158665305745e16,-6.05154138468997528539205278901e16,
                    -1.23202964813332152922910428164e17,-2.50708438598134633588546332631e17,
                    -5.09834752860261033238979568471e17
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^1 NF^0 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^1 NF^0 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^1 NF^0 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    307.565772724117013179298876269,-151.44631017006396729314334874,66.9617077819614768102200581126,
                    23.5662460894144947815412735204,50.2201088901002062477337613964,59.9535228230086032171048121058,
                    61.5923739598077415677161632658,59.7206813493177174172625792576,56.23394636244321122984263351,
                    51.9581904587930625083644798277,47.2828503766668232795164481141,42.4060551483501867622969193227,
                    37.4357134807614375373990051175,32.4341444326970644252523114441,27.439116898523539557807647817,
                    22.4743896932612686622587190961,17.5552983559616800083148372397,12.6918761024426948509002617328,
                    7.89068428464600250881601559553,3.15593504675513405577836713747,-1.50979174773138511782580431114,
                    -6.10507399292340262434704631144,-10.6293233255897039666823413867,-15.08255792011103694827616498,
                    -19.4652395735166211022866883037,-23.7781547151472631547786569835,-28.0223261844227952383665328733,
                    -32.198947031853964479311152438,-36.30933037916475804456259155,-40.3548711761654071411661840257,
                    -44.337016889662784646164334589,-48.2572449745940395451648512707,-52.1170455439098065478945201654,
                    -55.9179080548971371270074061291,-59.661311118669055393458782838,-63.3487147509701336519529398582,
                    -66.9815545391834518533845224113,-70.5612373180047045333435497469,-74.089138035385586942258654599,
                    -77.5665975585517815886612897795,-80.9949212225176965391297723438,-84.3753779644096743615812822297,
                    -87.709199918892106734284490374,-90.9975823751515850155828996693,-94.2416840157899323430411088605,
                    -97.4426273737814698254126625677,-100.601499456255543302732831208,-103.719352493955864421469733843,
                    -106.797204783331089736080112102,-109.836041594735866534817748617
                  },1),
                  // AQqPS_reg as^3 LM^1 NF^0 [above 0.75] y=1-x H[0, y]^1
                  ome_x<double>({
                    51.6328270215766763740240053334,24.8888888888888888888888888889,15.8797406018235899542709189137,
                    -53.4289013734850520210377230617,-91.5506054933419611624797212378,-118.503914551470475242193324352,
                    -140.691471759485176113024398915,-160.144919848552604465323159085,-177.685673945650272072281655613,
                    -193.732407485193892325187737238,-208.542481189169099328496051804,-222.296137719859049160377839752,
                    -235.130395494775638098916889395,-247.154712986584261445854074701,-258.459240754520505152093974986,
                    -269.119706234829340011085578178,-279.20059487304675234942859153,-288.757377343522043763511152121,
                    -297.838152473745362098744493312,-306.484904647020025242204092455,-314.734491612361354205622667792,
                    -322.619435342955278476868158102,-330.168564280481833017392373086,-337.407540717088080253159700211,
                    -344.359297781814253589347225096,-351.044404284617646286303693701,-357.481371339355585715879035049,
                    -363.686911567017464488958156541,-369.676159374490918581073266666,-375.462859064741954961037502082,
                    -381.059526200790903059342007998,-386.477586610022935115509557228,-391.727496601957722803833035355,
                    -396.818847327975530791262143642,-401.760455696582857255423899564,-406.560443843634014700638947948,
                    -411.226308821719222638995580695,-415.764983900097504213601651084,-420.182892643334101325988841839,
                    -424.485996753294143308398811968,-428.679838507601027332154967395,-432.769578501988860911381416496,
                    -436.760029299328767477359151162,-440.655685500635728419080670909,-444.46075067997894907745969137,
                    -448.179161563439291002498559154,-451.81460978007224999467280225,-455.370561468611851827772396868,
                    -458.850274986058401901474026678,-462.256816932241349701294654913
                  },1),
                  // AQqPS_reg as^3 LM^1 NF^0 [above 0.75] y=1-x H[0, y]^2
                  ome_x<double>({
                    15.1111111111111111111111111111,7.11111111111111111111111111111,18.0740740740740740740740740741,
                    26.8148148148148148148148148148,32.1234567901234567901234567901,36.4641975308641975308641975309,
                    40.2162761400856638951877047115,43.5232552280171327790375409423,46.4728786988046247305506564766,
                    49.1286223789927493631197334901,51.5390530606692222853839015455,53.7422537058900695264331627968,
                    55.7685556608633531710454787378,57.6424228057927691627325326959,59.3838095590476542857495238448,
                    61.0091610858277524944191610858,62.5321606508746070453106208585,63.9642947471218648970795581074,
                    65.3152851468725432961789246378,66.5934226715448086729748055405,67.8058276144874994941337267742,
                    68.9586548709221701863734753446,70.0572569765321313530615431385,71.1063148013046667165562758821,
                    72.1099431627368989957045774338,73.0717768224174649444986412725,73.995041012846531980478747874,
                    74.8826096689400903648119046977,75.7370538144522417431606018965,76.560682009585224466520873938,
                    77.3555743541075684879975941341,78.1236112258477629457857780696,78.8664976925674447783534783301,
                    79.5857843478265897684499920629,80.282885175248989561873118517,80.9590929307630667010439505636,
                    81.6155924416224466653376551483,82.2534721488218035665090867199,82.8737341617770138647947147578,
                    83.4773030476872320441313318904,84.065033540431379602252774116,84.6377173233176249583898517889,
                    85.1960890150644084286576895932,85.7408314679253306163210955211,86.2725804700000025415685622968,
                    86.7919289298071898605163318991,87.2994306095876051486699387738,87.7956034641156468356608395561,
                    88.2809326336837084945975833914,88.7558731330993925998531390193
                  },1),
                  // AQqPS_reg as^3 LM^1 NF^0 [above 0.75] y=1-x H[0, y]^3
                  ome_x<double>({
                    2.96296296296296296296296296296,0,2.96296296296296296296296296296,2.96296296296296296296296296296,
                    3.06172839506172839506172839506,3.16049382716049382716049382716,3.24514991181657848324514991182,
                    3.31569664902998236331569664903,3.37448559670781893004115226337,3.42386831275720164609053497942,
                    3.4657687991021324354657687991,3.50168350168350168350168350168,3.53276353276353276353276353276,
                    3.55989689323022656355989689323,3.58377425044091710758377425044,3.6049382716049382716049382716,
                    3.62381989832970225127087872186,3.64076494795449043815056886952,3.65605371453324669698938704787,
                    3.66991552956465237166991552956,3.68253968253968253968253968254,3.69408369408369408369408369408,
                    3.70467964670863221587859269019,3.71443907675791733762748255502,3.72345679012345679012345679012,
                    3.73181386514719848053181386515,3.73958003587633217262846892477,3.74681559866745051930237115422,
                    3.75357294897524782582253846622,3.75989782886334610472541507024,3.76583034647550776583034647551,
                    3.7714058144165671047391477499,3.77665544332210998877665544332,3.78160691886182082260513633063,
                    3.78628488432410001037452017844,3.79071134626690182245737801293,3.79490601712823935046157268379,
                    3.7988866059041497637988866059,3.80266906582696056380266906583,3.80626780626780626780626780627,
                    3.80969587473652514302920806986,3.81296511377812190820320901622,3.81608629670645174521143513392,
                    3.81906924542583457312139482682,3.82192293303404414515525626637,3.82465557344784397924494542852,
                    3.82727469992348191546464750288,3.82978723404255319148936170213,3.83219954648526077097505668934,
                    3.83451751070798689846308893928
                  },1)
                },0)
              )
            }
          ),
          // AQqPS_reg as^3 LM^1 NF^1
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^1 NF^1 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^1 NF^1 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^1 NF^1 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    -17.9094650205761316872427983539,-42.8392885865989932573181607023,84.3020208888792201034760789022,
                    -24.0127895677881869775775028151,-0.395061728395061728395061728395,0.114567901234567901234567901235,
                    0.0725596707818930041152263374486,0.0494069035021415973796926177879,0.0377900154770903070222798113955,
                    0.0310636267048461555894158564856,0.0266150941400794428425233833817,0.0233738579387662814524019163002,
                    0.020853789055401135801094988933,0.0188089798991355347457158104423,0.017101806119341290722615950095,
                    0.0156480698995245905397211152913,0.0143922739878318691040360904888,0.0132955854855582066727962994967,
                    0.0123295364072352812229129281646,0.0114724965569316196624670268552,0.0107075729504780230673508011806,
                    0.0100212878929246308626932378055,0.00940270607849194727579087282845,0.00884283402806648721192142971199,
                    0.00833419239753301015885262772466,0.00787050253535107952139976927436,0.00744645124762318587135299929211,
                    0.00705751074193514223539826236477,0.00669979851899843875760438344472,0.00636996682578569681583092477863,
                    0.00606511439524318663946784925455,0.00578271525684551933937494463721,0.00552056080228751317054123972947,
                    0.0052767122655703905907953397459,0.00504946147022554672296684010668,0.00483729819898233046143510465347,
                    0.00463888291141742082245195887612,0.00445302381183968053102335975799,0.00427865747914485616929630982578,
                    0.00411483243076306218948647859983,0.00396069511687059716414069371107,0.00381547793785370331456848664355,
                    0.00367848895419464334044885897803,0.00354910301834449874953131593774,0.00342675410635090259699059209647,
                    0.00331092866572403508098937906124,0.00320115982729842459939979140124,0.00309702235424735442741634637906,
                    0.00299812822213810574450123476447,0.00290412274091728744610868871681,0.00281468114371904191930707822083,
                    0.00272950557897164143633424925778
                  },-1),
                  // AQqPS_reg as^3 LM^1 NF^1 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    -59.7402968519648994421755836046,4.25970314803510055782441639536,18.1728395061728395061728395062,
                    -3.16049382716049382716049382716,-0.967901234567901234567901234568,-0.478024691358024691358024691358,
                    -0.288959435626102292768959435626,-0.19511211892164273116654069035,-0.141261442848744436046023347611,
                    -0.107309425827944346462864981383,-0.0844394562913081431599950118469,-0.068259701391014522327653640785,
                    -0.056371013946771522529098286674,-0.047367637111226854816598406342,-0.0403792110751817711524671231631,
                    -0.0348423357947167470976994786519,-0.03037866998651312376802572881,-0.0267261879334159226508323009669,
                    -0.0236985874628616623935563344832,-0.021160368208932171297665356169,-0.0190110103067936247376075015349,
                    -0.0171746345179929139828889202323,-0.0155930705149602194586196091937,-0.0142211042136573442066940377074,
                    -0.013023153186984314898618616325,-0.0119708988810438085800404640984,-0.0110415727145071874416603761333,
                    -0.010216697117615129644252784487,-0.00948114828388095634542957618633,
                    -0.00882244982176406721386557686251,-0.00823023439519418448752446936892,
                    -0.00769582911939413189308532999851,-0.00721193318314372771535449301665,
                    -0.00677236493894535410261935476127,-0.00637186184003841548142691059751,
                    -0.00600592095291936012537978814686,-0.0056706708932043510319265065608,
                    -0.00536276829412069915731472024607,-0.00507931357340323092779149944348,
                    -0.00481778198945611403882587915545,-0.00457596689188802680787803068176,
                    -0.0043519327596265926624839970071,-0.00414397614187066096656904743901,
                    -0.00395059301622098224336053868796,-0.00377045138559689149415509828728,
                    -0.00360236817373060151365619256527,-0.00344528966485436087000809830062,
                    -0.00329827487905766327783758718995,-0.00316048138995295695675377140342,
                    -0.00303115318271532733001210969906,-0.00290961022353007718582519931631
                  },0),
                  // AQqPS_reg as^3 LM^1 NF^1 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    -15.4074074074074074074074074074,-8.2962962962962962962962962963,4.74074074074074074074074074074
                  },0),
                  // AQqPS_reg as^3 LM^1 NF^1 [below 0.25] y=x H[0, y]^3
                  ome_x<double>({
                    -2.37037037037037037037037037037,-2.37037037037037037037037037037
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^1 NF^1 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^1 NF^1 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^1 NF^1 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -13.0767536834002478168863376564,-50.9996287356508266042964342153,-108.570185802365407334704723402,
                    -235.668931522097998404442756519,-490.747594213294716503751106843,-1004.7263685295208509930410672,
                    -2044.11386346225445962264406453,-4136.07104790318340033380263091,-8325.9283394467167567102097638,
                    -16759.9841802412213497297099622,-33629.6148676876935586890113843,-67550.794642793504277306594348,
                    -135376.412457446844473112214264,-271630.029132821158972777037681,-544054.813382624372936330566246,
                    -1.09092938740704423791376490752e6,-2.18438199964070884693399370309e6,-4.37821313704022418268233388007e6,
                    -8.76492465931171304103080383498e6,-1.75623999477501852915063693514e7,-3.51545116208413751117745980854e7,
                    -7.04233432210055339610775559468e7,-1.40953092432458620187804865388e8,-2.82313433137012107583201946132e8,
                    -5.65014088889835230061250837085e8,-1.13149484347015851848262598956e9,-2.26441455193147120485669200163e9,
                    -4.53415708312672820001121378234e9,-9.07359941965909540628329689131e9,
                    -1.81666912547278802193191826029e10,-3.63531055714718691028625799524e10,
                    -7.27779472658180700467775609465e10,-1.45629842515579131957586235029e11,
                    -2.91525023544592660423523595815e11,-5.83328357242969582056193668337e11,
                    -1.1676424474358771754002941423e12,-2.33633565750174224282877272634e12,
                    -4.67634737061780781226138910356e12,-9.35667260890585294011540893727e12,
                    -1.87271008438641194276584412328e13,-3.74692960199080610811747098074e13,
                    -7.4990221680627192603395645423e13,-1.50037838847542522625386004518e14,
                    -3.00269993779777602061592007369e14,-6.00758637170470875367710134903e14,
                    -1.2022504358163729018106348169e15,-2.40533522753210652001226916378e15,
                    -4.81344299539071208481607805007e15,-9.63007476093509181096862450541e15,
                    -1.92706561238975100621262016092e16,-3.85535165167413975496599401925e16
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^1 NF^1 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^1 NF^1 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^1 NF^1 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    -22.7590373148713386099766806912,3.55555555555555555555555555556,-13.9689138580812151531865572344,
                    -12.8084200309207213260260634073,-13.6845525504699305554121173727,-14.3868579095253126243043441776,
                    -14.8769237706830302425816896454,-15.2219920924548831918852052585,-15.4724054547008471098280636739,
                    -15.6600536504048303346757149206,-15.8049266863595854021242391525,-15.9198071129095438386401904936,
                    -16.0130885982551756669983721943,-16.0904325616573916833700643707,-16.1557548676275248705251627989,
                    -16.211827200699448605643949199,-16.2606529412052183372103196253,-16.3037079656083798749317331424,
                    -16.3420985369241520340571492477,-16.3766670828487825379973217141,-16.4080644780563656181811303572,
                    -16.4368003436772631503220046127,-16.4632786391604605744718664343,-16.487823237360125221865168789,
                    -16.5106965645532346959608283558,-16.5321133654603239083574600534,-16.5522509927949020519496756561,
                    -16.5712571864776109699034185059,-16.5892560174342976722306575777,-16.6063524741215371031572753064,
                    -16.6226360346532434937447586629,-16.6381834732077307370444187101,-16.653061083006912392111410329,
                    -16.6673264508369366190338269145,-16.6810298839863951380208808575,-16.6942155656686916600712194639,
                    -16.7069224967698523723657165974,-16.7191852682545369396422720536,-16.7310346984655956580583471955,
                    -16.7424983619441821286087326778,-16.7536010306208621018178696375,-16.764365043810546449165915357,
                    -16.7748106200424093226588726385,-16.7849561211193635209806927828,-16.7948182767452550868645685006,
                    -16.8044123764444250370749707925,-16.8137524342250201717608193767,-16.8228513304271501210580943801,
                    -16.831720934391147881301924675,-16.8403722109351766395851461911
                  },1),
                  // AQqPS_reg as^3 LM^1 NF^1 [above 0.75] y=1-x H[0, y]^1
                  ome_x<double>({
                    -5.92592592592592592592592592593,0,1.77777777777777777777777777778,4.44444444444444444444444444444,
                    4.81382716049382716049382716049,4.94617283950617283950617283951,5.04872763920382968002015621063,
                    5.14900478710002519526329050139,5.24863805604546345287086027827,5.34587497550460513423476386439,
                    5.4392350634774877199119623362,5.52789462688452587442486432385,5.61152857904994657131409268161,
                    5.69012406668084323761979439635,5.76384650035443686237337030988,5.83295428850984406539962095518,
                    5.89774726320285720516401023898,5.95853649438168212998406513721,6.01562723138351235080847208862,
                    6.06930976868079687046608995444,6.11985503731338907514640965916,6.16751297384431020660632122442,
                    6.21251248483309712573444368385,6.25506229182382759921625433008,6.29535222652674081758301258458,
                    6.33355471972989049173509931938,6.36982633381730355930903215909,6.4043092538014855332889137365,
                    6.43713269141908957847704656236,6.46841418081141810845467848206,6.49826075859116005305085350288,
                    6.52677002937764554617233357389,6.55403112246112658657184094829,6.58012554756523673397366932025,
                    6.60512795861242809321381891784,6.62910683452094924183792020007,6.6521250857241196277251045273,
                    6.67424059452060814054144462181,6.69550669667141780407456817342,6.7159726109354590505220298896,
                    6.73568382252759098431913881415,6.75468242581643632044230039315,6.7730074309666993028149468019,
                    6.79069503867674587698223655282,6.80777888666672876962683254933,6.82429027113277161368723285396,
                    6.84025834599449393725965907611,6.85571030242169870244464285478,6.87067153082647934171880788046,
                    6.8851657672446223012880097732
                  },1),
                  // AQqPS_reg as^3 LM^1 NF^1 [above 0.75] y=1-x H[0, y]^2
                  ome_x<double>({
                    -1.77777777777777777777777777778,0,-1.77777777777777777777777777778,-1.77777777777777777777777777778,
                    -1.83703703703703703703703703704,-1.8962962962962962962962962963,-1.94708994708994708994708994709,
                    -1.98941798941798941798941798942,-2.02469135802469135802469135802,-2.05432098765432098765432098765,
                    -2.07946127946127946127946127946,-2.10101010101010101010101010101,-2.11965811965811965811965811966,
                    -2.13593813593813593813593813594,-2.15026455026455026455026455026,-2.16296296296296296296296296296,
                    -2.17429193899782135076252723312,-2.18445896877269426289034132171,-2.19363222871994801819363222872,
                    -2.20194931773879142300194931774,-2.20952380952380952380952380952,-2.21645021645021645021645021645,
                    -2.22280778802517932952715561411,-2.22866344605475040257648953301,-2.23407407407407407407407407407,
                    -2.23908831908831908831908831909,-2.24374802152579930357708135486,-2.24808935920047031158142269253,
                    -2.25214376938514869549352307973,-2.25593869731800766283524904215,-2.2594982078853046594982078853,
                    -2.26284348864994026284348864994,-2.26599326599326599326599326599,-2.26896415131709249356308179838,
                    -2.27177093059446000622471210707,-2.27442680776014109347442680776,-2.27694361027694361027694361028,
                    -2.27933196354248985827933196354,-2.2816014394961763382816014395,-2.28376068376068376068376068376,
                    -2.28581752484191508581752484192,-2.28777906826687314492192540973,-2.28965177802387104712686108035,
                    -2.29144154725550074387283689609,-2.29315375982042648709315375982,-2.29479334406870638754696725711,
                    -2.29636481995408914927878850173,-2.29787234042553191489361702128,-2.29931972789115646258503401361,
                    -2.30071050642479213907785336357
                  },1)
                },0)
              )
            }
          )
        },0),
        // AQqPS_reg as^3 LM^2
        ome_nf<double>({
          // AQqPS_reg as^3 LM^2 NF^0
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^2 NF^0 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^2 NF^0 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^2 NF^0 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    311.050714382376906611744854124,-628.490028971784280776710706882,905.952094134401981073733596135,
                    -556.570684794577825898543211719,4.21399176954732510288065843621,-18.6240329218106995884773662551,
                    -2.53461728395061728395061728395,-4.04319610313261106911900562694,-1.27834335144992741364623450791,
                    -1.72033840971230741953846010283,-0.756632155377011344089533389945,-0.950260021338210638622161255906,
                    -0.499554142932573704943871253432,-0.60254829324830031194373922353,-0.354588329784074210587867870756,
                    -0.416162464378863364493119353643,-0.264837903160596340432377411272,-0.304672986151613877213921153065,
                    -0.205400248561364635373375098921,-0.232690148990084918662943946483,-0.163990678350774596927542145226,
                    -0.18351938541935026940502566574,-0.133978768361777605095514816274,-0.148442457314835628698733763076,
                    -0.111528017270438909062644266548,-0.122542984315021610803461529012,-0.0942922520690634047556478202916,
                    -0.102876425771417338044232306454,-0.0807705884074694212418563501572,-0.0875914803991654466394970152468,
                    -0.0699663093087490824353999600035,-0.0754766648596154381269194848206,-0.0611962990935699022852490508756,
                    -0.0657121161303046841744211195835,-0.0539795733818581154643615145564,-0.0577269036960644832405545422295,
                    -0.0479694726622573944322771310372,-0.0511135241077174002574824938354,-0.0429110011412378997536297616277,
                    -0.0455748056242684027030455520371,-0.0386131858619359090306352106732,-0.0408899040815684248190472280805,
                    -0.0349306992660631665378460203751,-0.0368919279871886658790179205436,-0.0317513579057910853664743520977,
                    -0.0334528592454836660466976944785,-0.0289874429862402108320090414007,-0.0304731718492382317052730908249,
                    -0.0265695628181507410966416099207,-0.0278745470781708823401131836822,-0.0244422402581756512769500509456,
                    -0.0255946727970811284809152863938
                  },-1),
                  // AQqPS_reg as^3 LM^2 NF^0 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    71.1111111111111111111111111111,156.213164105673080942042527998,148.328270215766763740240053334,
                    275.75308641975308641975308642,-7.9012345679012345679012345679,15.0024691358024691358024691358,
                    -1.1950617283950617283950617284,4.4788712522045855379188712522,-0.487780297304106827916351725876,
                    2.18955236415553875871336188797,-0.268273564569860866157162453459,1.30881157251527621897992268363,
                    -0.170649253477536305819134101962,0.873750716174958599201023443448,-0.118419092778067137041496015855,
                    0.625877771665317452863240409028,-0.0871058394867918677442486966296,0.470869384790953418404398796556,
                    -0.0668154698335398066270807524172,0.36732810567435576710012318449,-0.0529009205223304282441633904226,
                    0.29467065975530118343291627379,-0.0429365862949822849572223005807,0.241692592981883401608603942502,
                    -0.0355527605341433605167350942685,0.201858874398256880928588553037,-0.0299272472026095214501011602461,
                    0.171144377074861405345735830066,-0.0255417427940378241106319612176,0.146957798400154823354158430888,
                    -0.0220561245544101680346639421563,0.127568633125509859556629275218,-0.0192395727984853297327133249963,
                    0.111784964338727779587994641758,-0.0169309123473633852565483869032,0.0987638585205954399621171142614,
                    -0.0150148023822984003134494703671,0.0878953988446674409948608516925,-0.0134069207353017478932868006152,
                    0.078729360387750079380768241374,-0.0120444549736402850970646978886,0.0709274868242644155221094755672,
                    -0.0108798318990664816562099925177,0.0642316301989952449818202353047,-0.00987648254055245560840134671989,
                    0.0584419964767518181594040234528,-0.00900592043432650378414048141317,0.0534019898052425934851255236597,
                    -0.00824568719764415819459396797488,0.048987461544270832829683456753,-0.00757788295678831832503027424765,
                    0.0450989584647161963802905894028
                  },-1),
                  // AQqPS_reg as^3 LM^2 NF^0 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    -94.2222222222222222222222222222,67.5555555555555555555555555556,9.48148148148148148148148148148
                  },0),
                  // AQqPS_reg as^3 LM^2 NF^0 [below 0.25] y=x H[0, y]^3
                  ome_x<double>({
                    16.5925925925925925925925925926,-26.0740740740740740740740740741
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^2 NF^0 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^2 NF^0 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^2 NF^0 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -28.4829604120560137045971232782,-98.9104335686272813806596789824,-10.9607968218768300018004710628,
                    90.3964479384806919403473222279,4.36934500210541410901613141079,-59.6019683199596112828216269481,
                    -664.247043887792381483116965809,-1855.41705878352929880668088369,-5466.96671553821897821507517864,
                    -13278.0995366930154036242630822,-32798.8958399170474768918828036,-75109.1996910456497786619280782,
                    -173052.787440911798413589576581,-383322.848985447350989033616325,-851096.711345955545651273008952,
                    -1.8454891005778309643648078642e6,-4.00528761230280838666977114915e6,-8.55831396468816762030369154998e6,
                    -1.82919500847102984152611076729e7,-3.86714371585967028570222161935e7,-8.17554395989601041994246490627e7,
                    -1.7146091287779262202692375568e8,-3.59546998816121284608275641674e8,-7.4938255273347861054843590007e8,
                    -1.56160466696648555262737143749e9,-3.2387211782962477030856834156e9,-6.71565555571539809521711280748e9,
                    -1.38724213304748729778199321052e10,-2.86502690079923198052728399058e10,
                    -5.89874875875174935794003447225e10,-1.21424792949854598167364080475e11,
                    -2.49310766200800019644204883781e11,-5.11795194824657976443006076786e11,
                    -1.04837448916434834071519473319e12,-2.14715764583165228755721200645e12,
                    -4.38952239219637770620179713229e12,-8.97229566897653967655612534268e12,
                    -1.83108404480344701128654820383e13,-3.73638314230941796905728005137e13,
                    -7.61384691786874746217723467869e13,-1.5513163476624921220002055829e14,
                    -3.15704698412722942784685676745e14,-6.42406396932217491126426200669e14,
                    -1.30582561355269181475435441553e15,-2.65407301310835431819507404353e15,
                    -5.38938292329829098520425585736e15,-1.0942621952643875359390661836e16,
                    -2.21996513966528986540025865964e16,-4.50329649702095448785365574417e16,
                    -9.12840142295039332669292662327e16,-1.85021270718384768755615943002e17
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^2 NF^0 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^2 NF^0 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^2 NF^0 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    -112.309085864084483273873799111,-10.6666666666666666666666666667,-48.0127895677881869775775028151,
                    -6.08686364186226105165157688918,-2.58012946696137345707699982253,-2.85611134144320191188513880525,
                    -3.58407564907413675775495206281,-4.38387457764287371989010130679,-5.23325490311444163569195610145,
                    -6.14213808861268348550301307917,-7.11302880100842685968221286779,-8.14116165578324020603923786967,
                    -9.2182611225272146105084147481,-10.3350888202924674542574664425,-11.4827598928025892474567366345,
                    -12.6533113857044560296314800875,-13.839876029109671639207081371,-15.0366654615722131278324994331,
                    -16.2388713093402675640759472453,-17.442538934345253767716012732,-18.644440125854696998884685068,
                    -19.8419562537604753269166767242,-21.0329759806093130605164516469,-22.2158080475781956238746787133,
                    -23.3891080474676958299946968371,-24.5518174961648878494086167047,-25.7031134006197875224928867167,
                    -26.8423666307109377157049253759,-27.9691075991882468053222608006,-29.0829979713944031689109673203,
                    -30.1838073336447063323826463237,-31.2713939333294536928153952495,-32.3456887615110228880849647074,
                    -33.4066823808754143959478508397,-34.4544140110490104167846707704,-35.4889624727498510372614389309,
                    -36.5104386652050564425331240823,-37.5189793106167414106224461412,-38.514741747683032144344940672,
                    -39.4978995953620944133808920399,-40.46863913992304709828044984,-41.4271563242593770787508982729,
                    -42.3736542395833597733229474427,-43.3083410368903899340005899058,-44.2314281897180347550203979875,
                    -45.1431290513198448050061983942,-46.0436576589062973333394775915,-46.9332277454594164464496913682,
                    -47.8120519261137416581624983963,-48.6803410314647786370874596838
                  },1),
                  // AQqPS_reg as^3 LM^2 NF^0 [above 0.75] y=1-x H[0, y]^1
                  ome_x<double>({
                    -58.6666666666666666666666666667,-14.2222222222222222222222222222,-65.7777777777777777777777777778,
                    -86.5185185185185185185185185185,-98.8444444444444444444444444444,-108.958024691358024691358024691,
                    -117.686067019400352733686067019,-125.346031746031746031746031746,-132.143209876543209876543209877,
                    -138.230970017636684303350970018,-143.728253968253968253968253968,-148.728726417615306504195393084,
                    -153.306977795866684755573644463,-157.523035694464265892837321409,-161.425704971101796498621895447,
                    -165.055059672837450615228393006,-168.444315379609497256556080085,-171.621247372009899243014711424,
                    -174.60927431531720030745371876,-177.428294013061814919400058719,-180.095333955389977798941384376,
                    -182.625062337584284161958293812,-185.030193028553263015260837065,-187.321809218834390459104417914,
                    -189.509624165537478230976682989,-191.602192867332186217635134044,-193.607085147734155875424078612,
                    -195.531028150116137205201155215,-197.380024407296945287127551148,-199.159450268543444005936357972,
                    -200.87413842383379215770858648,-202.528447470818746481224839569,-204.126320860273940884087697975,
                    -205.67133708458426441057564157,-207.166752606981057909968397551,-208.615538741808083461949771495,
                    -210.020413469401277023999686955,-211.383868989290854199280525627,-212.708195671865452868441883056,
                    -213.995502953403496080214136735,-215.24773762638701257221720788,-216.466699901587182572002072873,
                    -217.654057556929668666247530326,-218.811358437796514006520661101,-219.940041531998541094845348593,
                    -221.0414468084266617782797548,-222.116823979994032542768824507,-223.167340327826653659377909592,
                    -224.194087703880272832635395062,-225.198088812560654891622644457
                  },1)
                },0)
              )
            }
          ),
          // AQqPS_reg as^3 LM^2 NF^1
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^2 NF^1 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^2 NF^1 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^2 NF^1 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    24.7343459568466472519519893332,-26.2286170061163157110109736298,-1.58024691358024691358024691358,
                    1.58024691358024691358024691358,0.483950617283950617283950617284,0.239012345679012345679012345679,
                    0.144479717813051146384479717813,0.0975560594608213655832703451751,0.0706307214243722180230116738053,
                    0.0536547129139721732314324906917,0.0422197281456540715799975059234,0.0341298506955072611638268203925,
                    0.028185506973385761264549143337,0.023683818555613427408299203171,0.0201896055375908855762335615815,
                    0.0174211678973583735488497393259,0.015189334993256561884012864405,0.0133630939667079613254161504834,
                    0.0118492937314308311967781672416,0.0105801841044660856488326780845,0.00950550515339681236880375076743,
                    0.00858731725899645699144446011614,0.00779653525748010972930980459684,0.00711055210682867210334701885369,
                    0.00651157659349215744930930816249,0.00598544944052190429002023204922,0.00552078635725359372083018806666,
                    0.00510834855880756482212639224352,0.00474057414194047817271478809317,0.00441122491088203360693278843126,
                    0.00411511719759709224376223468446,0.00384791455969706594654266499925,0.00360596659157186385767724650833,
                    0.00338618246947267705130967738063,0.00318593092001920774071345529875,0.00300296047645968006268989407343,
                    0.0028353354466021755159632532804,0.00268138414706034957865736012304,0.00253965678670161546389574972174,
                    0.00240889099472805701941293957772,0.00228798344594401340393901534088,0.00217596637981329633124199850355,
                    0.00207198807093533048328452371951,0.00197529650811049112168026934398,0.00188522569279844574707754914364,
                    0.00180118408686530075682809628263,0.00172264483242718043500404915031,0.00164913743952883163891879359498,
                    0.00158024069497647847837688570171,0.00151557659135766366500605484953,0.00145480511176503859291259965816
                  },0),
                  // AQqPS_reg as^3 LM^2 NF^1 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    15.4074074074074074074074074074,8.2962962962962962962962962963,-4.74074074074074074074074074074
                  },0),
                  // AQqPS_reg as^3 LM^2 NF^1 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    3.55555555555555555555555555556,3.55555555555555555555555555556
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^2 NF^1 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^2 NF^1 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^2 NF^1 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    1.29273848016846738976737193701,5.06410027248028099941434463459,1.99375861722151856236232834667,
                    1.90512131666432113380035374816,14.6208965262911361239904078756,23.7832837257667068626487808024,
                    62.1343590168656212028550647373,113.95932105392000300643267393,249.058886575562193128743000821,
                    464.518295986722683996799711401,960.684743606294329746961335297,1808.48445212029717888767869252,
                    3656.09964994659490849474491446,6930.98006178378904898844287067,13860.2225630570221346975831217,
                    26422.8740216179537709754452176,52539.8911806159619248493154819,100621.01678777346109311308941,
                    199457.370717335837287316979024,383444.563481207082288406436541,758781.510366652305513550604693,
                    1.46339407969936787974403208246e6,2.89314230240124698833024045117e6,5.59501139433162716955742198063e6,
                    1.10562126169219386451824286164e7,2.14318696844006109973214473757e7,4.23434299025583751450404689445e7,
                    8.22487511637368070438213892341e7,1.62499772541552820782139640192e8,3.16210796270183805216059906081e8,
                    6.2480692061883533874448441503e8,1.21775552853674249622194610964e9,2.40660648073526828447531677293e9,
                    4.69714599795366161220800547615e9,9.28481734629161296587013382083e9,1.81448379162869871723426632088e10,
                    3.58754708932373884475132672624e10,7.01897295941348511461650767177e10,1.38812677513317699834613544794e11,
                    2.71865736684512640930001137274e11,5.37803412676503547383017975306e11,1.0542843303877438538820592568e12,
                    2.0861261591055806080767082762e12,4.09306194832584327598264766372e12,8.10109283448508798745673562288e12,
                    1.59072087452531711312037695136e13,3.14919462243769141436374899198e13,6.18821293577181519498688283485e13,
                    1.22540071752584363300771518214e14,2.40954550391423948722395485303e14,4.772567970846171959483969683e14
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^2 NF^1 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^2 NF^1 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^2 NF^1 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    5.92592592592592592592592592593,0,1.18518518518518518518518518519,0.296296296296296296296296296296,
                    0.282469135802469135802469135802,0.347654320987654320987654320988,0.405946082136558327034517510708,
                    0.448173343411438649533887629126,0.476699980403684107387811091515,0.495018616500097981579463060945,
                    0.505990739526093061446596800132,0.511718434950758183081415404648,0.513702568403423104277805132506,
                    0.513010365025017039669054321069,0.510404210721671039131356591674,0.506432908655130877353099575322,
                    0.501494816431379291817584781829,0.495881386615719564585385423525,0.48980734940939890879238732877,
                    0.483431820400661492038024183662,0.47687321333200737594073875955,0.470219887833440431513153424885,
                    0.46353783112960151418736343286,0.456876255130647923482896493094,0.450271712823586112947624142491,
                    0.443751152078951313413625006102,0.437334199409623053973277345651,0.431034880480416583944640928375,
                    0.424862925021564896504409603918,0.418824762578581504379509496615,0.412924286515230604694978770159,
                    0.407163443064803208080570536574,0.40154268743399035660497281851,0.396061338261594337810959788504,
                    0.390717853926187648473974944753,0.385510048455029687800412422635,0.380435260532640477734636230433,
                    0.37549048593352754594449634332,0.370672481319611484158874536468,0.365977845540955896379872283132,
                    0.361403083208246372307438823567,0.356944654257737388499759278986,0.35259901242394914174543055297,
                    0.34836263491310586188389135176,0.344232045087239218626300233722,0.340203829592213182183153878232,
                    0.336274651067990483575773143979,0.332441257347609712088435547624,0.328700487868407817102828764885,
                    0.325049277874184651002428441078
                  },1),
                  // AQqPS_reg as^3 LM^2 NF^1 [above 0.75] y=1-x H[0, y]^1
                  ome_x<double>({
                    3.55555555555555555555555555556,0,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.67407407407407407407407407407,3.79259259259259259259259259259,3.89417989417989417989417989418,
                    3.97883597883597883597883597884,4.04938271604938271604938271605,4.10864197530864197530864197531,
                    4.15892255892255892255892255892,4.20202020202020202020202020202,4.23931623931623931623931623932,
                    4.27187627187627187627187627188,4.30052910052910052910052910053,4.32592592592592592592592592593,
                    4.34858387799564270152505446623,4.36891793754538852578068264343,4.38726445743989603638726445744,
                    4.40389863547758284600389863548,4.41904761904761904761904761905,4.4329004329004329004329004329,
                    4.44561557605035865905431122822,4.45732689210950080515297906602,4.46814814814814814814814814815,
                    4.47817663817663817663817663818,4.48749604305159860715416270972,4.49617871840094062316284538507,
                    4.50428753877029739098704615946,4.51187739463601532567049808429,4.51899641577060931899641577061,
                    4.52568697729988052568697729988,4.53198653198653198653198653199,4.53792830263418498712616359675,
                    4.54354186118892001244942421413,4.54885361552028218694885361552,4.55388722055388722055388722055,
                    4.55866392708497971655866392708,4.56320287899235267656320287899,4.56752136752136752136752136752,
                    4.57163504968383017163504968383,4.57555813653374628984385081946,4.5793035560477420942537221607,
                    4.58288309451100148774567379219,4.58630751964085297418630751964,4.58958668813741277509393451422,
                    4.59272963990817829855757700346,4.59574468085106382978723404255,4.59863945578231292517006802721,
                    4.60142101284958427815570672714
                  },1)
                },0)
              )
            }
          )
        },0),
        // AQqPS_reg as^3 LM^3
        ome_nf<double>({
          // AQqPS_reg as^3 LM^3 NF^0
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^3 NF^0 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^3 NF^0 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^3 NF^0 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    71.9012345679012345679012345679,-105.421440792275276284117513086,96.0600406892062051973639683953,
                    -57.4156378600823045267489711934,-2.63374485596707818930041152263,-0.806584362139917695473251028807,
                    -0.398353909465020576131687242798,-0.240799529688418577307466196355,-0.162593432434702275972117241959,
                    -0.117717869040620363371686123009,-0.0894245215232869553857208178196,-0.0703662135760901192999958432057,
                    -0.0568830844925121019397113673208,-0.046975844955642935440915238895,-0.039473030926022379013832005285,
                    -0.0336493425626514759603892693026,-0.0290352798289306225814162322099,-0.0253155583220942698066881073417,
                    -0.0222718232778466022090269174724,-0.019748822885718051994630278736,-0.0176336401741101427480544634742,
                    -0.0158425085889946872813395846124,-0.0143121954316607616524074335269,-0.0129942254291335162155163409947,
                    -0.0118509201780477868389116980895,-0.0108526276558202624155155136042,
                    -0.00997574906753650715003372008203,-0.00920131059542265620138364677776,
                    -0.00851391426467927470354398707253,-0.00790095690323413028785798015528,
                    -0.00735204151813672267822131405209,-0.00685852866266182040627039114077,
                    -0.00641319093282844324423777499876,-0.00600994431928643976279541084721,
                    -0.00564363744912112841884946230106,-0.00530988486669867956785575883126,
                    -0.00500493412743280010448315678905,-0.004725559077670292526605422134,
                    -0.00446897357843391596442893353839,-0.0042327613111693591064929162029,
                    -0.00401481832454676169902156596287,-0.00381330574324002233989835890146,
                    -0.00362661063302216055206999750592,-0.00345331345155888413880753953251,
                    -0.00329216084685081853613378223996,-0.0031420428213307429117959152394,
                    -0.00300197347810883459471349380439,-0.00287107472071196739167341525052,
                    -0.00274856239921471939819798932496,-0.00263373449162746413062814283618,
                    -0.00252596098559610610834342474922,-0.00242467518627506432152099943026
                  },-1),
                  // AQqPS_reg as^3 LM^3 NF^0 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    14.2222222222222222222222222222,29.6296296296296296296296296296,23.7037037037037037037037037037,
                    6.32098765432098765432098765432
                  },-1),
                  // AQqPS_reg as^3 LM^3 NF^0 [below 0.25] y=x H[0, y]^2
                  ome_x<double>({
                    -13.037037037037037037037037037,18.962962962962962962962962963
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^3 NF^0 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^3 NF^0 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^3 NF^0 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    20.3866978962722508436820860639,70.9374555177035338623646445573,155.058524366895043723823425084,
                    344.73488103058584018331612463,633.853491466133890033699478846,1215.7430076985759971214131501,
                    2282.50375630291800218098690958,4331.78314135122992847165333662,8137.32915797510904153277279155,
                    15381.5499648294332073729393989,28889.6549661371400399522960701,54487.7278352901522200801927411,
                    102269.699506156573748894118392,192505.007179175654646266254678,360827.219599630201087268168166,
                    677605.059951778304576507584048,1.26731914112237429597495391433e6,2.37274836843090626174559577803e6,
                    4.42372791697308200356418001798e6,8.24958625501000483688057357801e6,1.53134344354900858219325507153e7,
                    2.84081513380704670704650229584e7,5.24213148325157528114850722167e7,9.65711264042318556981156487367e7,
                    1.76768143069777133812382496184e8,3.22569743037060895781430474203e8,5.83857620642649424093426945655e8,
                    1.05133975771859825374371695677e9,1.87239807010127849115301070407e9,3.30592723287534822389433053667e9,
                    5.74339717968268196779825614089e9,9.82655418938046604041645508799e9,1.6367336863451642355515207285e10,
                    2.64360864175766782014234398459e10,4.04045039676612773768593481903e10,5.68516330840805194538788492536e10,
                    6.62713233985885849601768733522e10,4.12028871905260061918992921823e10,
                    -9.84736334878332992131394515734e10,-5.45940741673170711702467612059e11,
                    -1.78315226217345832565341270107e12,-4.90238887196532790868915717325e12,
                    -1.24518622348014784655018711306e13,-3.00280965441440686302708199781e13,
                    -7.02111306184547856195022736746e13,-1.60109115899052618373209623991e14,
                    -3.5924064042406986888716758869e14,-7.94231903441717323675145891809e14,
                    -1.73864749983345018299799708638e15,-3.76918640213059425730324190495e15,
                    -8.11720666174177329207216517687e15
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^3 NF^0 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^3 NF^0 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^3 NF^0 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    24.2962962962962962962962962963,-2.96296296296296296296296296296,25.0864197530864197530864197531,
                    22.6172839506172839506172839506,21.8238683127572016460905349794,21.3267489711934156378600823045,
                    20.8571092634584698076761568825,20.3679180314100949021583942219,19.8585276634659350708733424783,
                    19.3367561565092429289960154158,18.8102705425264347823270382193,18.2849479240051630624021196412,
                    17.7649231848662048092247522447,17.2529695059609589524119438649,16.750868805895260921715948171,
                    16.2597022182207367392552577738,15.7800624205276262177454026743,15.312204678192118920044077035,
                    14.8561531853452746155779970106,14.4117757029760119700312459472,13.9788359429069046249251664751,
                    13.5570303625932895184455457688,13.1460140196644038931912074192,12.745418721204631363474856432,
                    12.354865723490063538821409628,11.9739745617836582674081365396,11.6023691222117200703013372472,
                    11.2396817430074702058087420972,10.8855559056213901939145101519,10.5396479168402851137200919792,
                    10.2016278703686417714945555728,9.87118009614031605795576861292,9.54800324822861244883570114283,
                    9.23181014090248969729226870199,8.92232741247776835989738354267,8.61929507487726596432068533958,
                    8.32246599094930590699944157518,8.03160530997470335773942507601,7.74648988325890123838737785493,
                    7.46690767542655575796336213726,7.19265718241013443610911861424,6.92354686371622980546726167104,
                    6.65939459404541761041842367361,6.40002713750057816753051652866,6.14527964627308462932177403993,
                    5.89499518471930521020346067939,5.64902427903735540924629571093,5.40722449225591120834502740874,
                    5.16946002390085259127382514889,4.93560133347240899148431085744
                  },1),
                  // AQqPS_reg as^3 LM^3 NF^0 [above 0.75] y=1-x H[0, y]^1
                  ome_x<double>({
                    -5.92592592592592592592592592593,0,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -6.12345679012345679012345679012,-6.32098765432098765432098765432,-6.49029982363315696649029982363,
                    -6.63139329805996472663139329806,-6.74897119341563786008230452675,-6.84773662551440329218106995885,
                    -6.9315375982042648709315375982,-7.00336700336700336700336700337,-7.06552706552706552706552706553,
                    -7.11979378646045312711979378646,-7.16754850088183421516754850088,-7.20987654320987654320987654321,
                    -7.24763979665940450254175744372,-7.28152989590898087630113773905,-7.31210742906649339397877409573,
                    -7.33983105912930474333983105913,-7.36507936507936507936507936508,-7.38816738816738816738816738817,
                    -7.40935929341726443175718538037,-7.42887815351583467525496511004,-7.44691358024691358024691358025,
                    -7.46362773029439696106362773029,-7.47916007175266434525693784953,-7.49363119733490103860474230845,
                    -7.50714589795049565164507693243,-7.51979565772669220945083014049,-7.53166069295101553166069295102,
                    -7.5428116288331342094782954998,-7.55331088664421997755331088664,-7.56321383772364164521027266125,
                    -7.57256976864820002074904035688,-7.58142269253380364491475602587,-7.58981203425647870092314536759,
                    -7.59777321180829952759777321181,-7.60533813165392112760533813165,-7.61253561253561253561253561254,
                    -7.61939174947305028605841613972,-7.62593022755624381640641803243,-7.63217259341290349042287026783,
                    -7.63813849085166914624278965364,-7.64384586606808829031051253273,-7.64931114689568795848989085704,
                    -7.65454939984696383092929500577,-7.65957446808510638297872340426,-7.66439909297052154195011337868,
                    -7.66903502141597379692617787856
                  },1)
                },0)
              )
            }
          ),
          // AQqPS_reg as^3 LM^3 NF^1
          ome_piecewisex<double>(
            {0.25,0.75},
            {
              // AQqPS_reg as^3 LM^3 NF^1 [below 0.25]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AQqPS_reg as^3 LM^3 NF^1 [below 0.25] y=x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^3 NF^1 [below 0.25] y=x H[0, y]^0
                  ome_x<double>({
                    -1.58024691358024691358024691358,-1.18518518518518518518518518519,1.18518518518518518518518518519,
                    1.58024691358024691358024691358
                  },-1),
                  // AQqPS_reg as^3 LM^3 NF^1 [below 0.25] y=x H[0, y]^1
                  ome_x<double>({
                    -2.37037037037037037037037037037,-2.37037037037037037037037037037
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^3 NF^1 [0.25,0.75]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AQqPS_reg as^3 LM^3 NF^1 [0.25,0.75] y=1/2-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^3 NF^1 [0.25,0.75] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -0.893501382700441368985643815062,-3.61832418108036419022356255951,-8.69135802469135802469135802469,
                    -20.5432098765432098765432098765,-42.6666666666666666666666666667,-87.8617283950617283950617283951,
                    -179.516049382716049382716049383,-364.811287477954144620811287478,-738.652557319223985890652557319,
                    -1491.75308641975308641975308642,-3007.10452674897119341563786008,-6053.43748597081930415263748597,
                    -12173.0729517396184062850729517,-24459.3048433048433048433048433,-49114.2129968796635463302129969,
                    -98569.8426807760141093474426808,-197740.720987654320987654320988,-396547.532316630355846042120552,
                    -794998.797385620915032679738562,-1.59341763396144682694390296729e6,-3.19301271371020142949967511371e6,
                    -6.39723827019400352733686067019e6,-1.28149201577681577681577681578e7,-2.56672657020446005953252330064e7,
                    -5.14033005517981749865807836822e7,-1.02933398313086419753086419753e8,-2.06101327129040835707502374169e8,
                    -4.12637725335865780310224754669e8,-8.26084734670585929845189104448e8,-1.65367864446731131788603052971e9,
                    -3.31017828551656023839931885909e9,-6.62564130537060931899641577061e9,
                    -1.32612032839060135404221425727e10,-2.65410661411627384960718294052e10,
                    -5.31172927725020136000528157391e10,-1.06300952995586098142960888059e11,
                    -2.12727381952428924162257495591e11,-4.25692356388416861305750194639e11,
                    -8.51835255799062501097588816887e11,-1.70452604823813745189183785675e12,
                    -3.41067880433763190883190883191e12,-6.82445446582991050888286660644e12,
                    -1.3654811587965858132232115972e13,-2.73208864437786564948115335712e13,
                    -5.46632882889785704069114921828e13,-1.09367718223601691582491582492e14,
                    -2.18814184793783894936482376096e14,-4.37779242223866523954728703418e14,
                    -8.75847799872453631731021801944e14,-1.75225087542627574905517762661e15,
                    -3.50556836300919757174099269337e15
                  },0)
                },0)
              ),
              // AQqPS_reg as^3 LM^3 NF^1 [above 0.75]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AQqPS_reg as^3 LM^3 NF^1 [above 0.75] y=1-x
                ome_logx<double>({
                  // AQqPS_reg as^3 LM^3 NF^1 [above 0.75] y=1-x H[0, y]^0
                  ome_x<double>({
                    -1.18518518518518518518518518519,0,-1.18518518518518518518518518519,-1.18518518518518518518518518519,
                    -1.22469135802469135802469135802,-1.26419753086419753086419753086,-1.29805996472663139329805996473,
                    -1.32627865961199294532627865961,-1.34979423868312757201646090535,-1.36954732510288065843621399177,
                    -1.38630751964085297418630751964,-1.40067340067340067340067340067,-1.41310541310541310541310541311,
                    -1.42395875729209062542395875729,-1.43350970017636684303350970018,-1.44197530864197530864197530864,
                    -1.44952795933188090050835148874,-1.45630597918179617526022754781,-1.46242148581329867879575481915,
                    -1.46796621182586094866796621183,-1.47301587301587301587301587302,-1.47763347763347763347763347763,
                    -1.48187185868345288635143707607,-1.48577563070316693505099302201,-1.48938271604938271604938271605,
                    -1.49272554605887939221272554606,-1.49583201435053286905138756991,-1.49872623946698020772094846169,
                    -1.50142917959009913032901538649,-1.5039591315453384418901660281,-1.5063321385902031063321385902,
                    -1.50856232576662684189565909996,-1.51066217732884399551066217733,-1.51264276754472832904205453225,
                    -1.51451395372964000414980807138,-1.51628453850676072898295120517,-1.51796240685129574018462907352,
                    -1.51955464236165990551955464236,-1.52106762633078422552106762633,-1.52250712250712250712250712251,
                    -1.52387834989461005721168322794,-1.52518604551124876328128360649,-1.52643451868258069808457405357,
                    -1.52762769817033382924855793073,-1.52876917321361765806210250655,-1.52986222937913759169797817141,
                    -1.53090987996939276618585900115,-1.53191489361702127659574468085,-1.53287981859410430839002267574,
                    -1.53380700428319475938523557571
                  },1)
                },0)
              )
            }
          )
        },0)
      },0)
    },2);


  const rpd_distribution<ome_as_view<double>, ome_as_plus_view<double>, ome_as_const_view<double>> AQqPS{
    AQqPS_reg.get_view(),
    std::nullopt,
    std::nullopt
  };
}


/* C interface */

double ome_AQqPS_reg(double as, double LM, double NF, double x)
{
  return(ome::AQqPS_reg(as, LM, NF, x));
}

double ome_AQqPS_reg_trunc_as(int trunc_order, double as, double LM, double NF, double x)
{
  auto trunc = ome::AQqPS_reg.truncate(trunc_order);
  return(trunc(as, LM, NF, x));
}

double ome_AQqPS_reg_coeff_as(int order_as, double LM, double NF, double x)
{
  auto& coeff = ome::AQqPS_reg[order_as];
  return(coeff(LM, NF, x));
}

double ome_AQqPS_reg_coeff_as_LM(int order_as, int order_LM, double NF, double x)
{
  auto& coeff = ome::AQqPS_reg[order_as][order_LM];
  return(coeff(NF, x));
}

double ome_AQqPS_reg_coeff_as_LM_NF(int order_as, int order_LM, int order_NF, double x)
{
  auto& coeff = ome::AQqPS_reg[order_as][order_LM][order_NF];
  return(coeff(x));
}

int ome_AQqPS_reg_min_power()
{
  return(ome::AQqPS_reg.min_power());
}

int ome_AQqPS_reg_max_power()
{
  return(ome::AQqPS_reg.max_power());
}

int ome_AQqPS_reg_coeff_as_min_power(int order_as)
{
  return(ome::AQqPS_reg[order_as].min_power());
}

int ome_AQqPS_reg_coeff_as_max_power(int order_as)
{
  return(ome::AQqPS_reg[order_as].max_power());
}

int ome_AQqPS_reg_coeff_as_LM_min_power(int order_as, int order_LM)
{
  return(ome::AQqPS_reg[order_as][order_LM].min_power());
}

int ome_AQqPS_reg_coeff_as_LM_max_power(int order_as, int order_LM)
{
  return(ome::AQqPS_reg[order_as][order_LM].max_power());
}