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/AgqQ.h>

namespace ome
{

  /**************** AgqQ, regular part ****************/

  const ome_as<double> AgqQ_reg = 
    // AgqQ_reg
    ome_as<double>({
      // AgqQ_reg as^2
      ome_logm<double>({
        // AgqQ_reg as^2 LM^0
        ome_nf<double>({
          // AgqQ_reg as^2 LM^0 NF^0
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^2 LM^0 NF^0 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^2 LM^0 NF^0 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^0 NF^0 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    11.0617283950617283950617283951,-16.987654320987654320987654321,13.2345679012345679012345679012,
                    -3.75308641975308641975308641975,-1.13580246913580246913580246914,-0.543209876543209876543209876543,
                    -0.301234567901234567901234567901,-0.174955908289241622574955908289,-0.099823633156966490299823633157,
                    -0.0513815402704291593180482069371,-0.0184244562022339800117577895356,0.00486558708780931003153225375448,
                    0.0217882528993640104751215862327,0.034341789897345452901008456564,0.0437991755452072912390372707833,
                    0.0510047037031164015290999417984,0.0565370490767316164141560966958,0.0608038125031589083876665576012,
                    0.0640987935650245018436739569637,0.0666379838951791709502981146141,0.0685829069408934103349881356991,
                    0.0700561669071178137962746543021,0.0711520542532564566258426577674,0.0719439308331907185469829583019,
                    0.0724894679187234570111665197671,0.0728344228389097595361336555755,0.073015402432016463652705546048,
                    0.0730619122738709527834612827244,0.0729978947953814952776715502107,0.0728428965969446797881172762139,
                    0.0726129633629285444736017154855,0.0723213323541607306335237874201,0.0719789728803945327184977930086,
                    0.0715950114840991667976591916733,0.0711770688967372219937896346492,0.0707315289064228961737815027031,
                    0.0702637542655826296103405604862,0.0697782611033434504289275782553,0.0692788606021867369509846245905,
                    0.0687687746830291083734606934359,0.0682507309287155258945650752981,0.0677270408292406061100305419861,
                    0.0671996645571522184331829611923,0.0666702648093942532400909350538,0.0661402517319204046307221083808,
                    0.0656108205387265327012572652536,0.0650829831200797962669925222594,0.0645575946852057276730036572958,
                    0.0640353762871551452472290773173,0.0635169339203822903102172967719,0.0630027747558720479705586189184,
                    0.0624933209776789359260392115673
                  },-1)
                },0)
              ),
              // AgqQ_reg as^2 LM^0 NF^0 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^2 LM^0 NF^0 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^0 NF^0 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    10.6257609527036173890333144336,36.0297318477967518516906667563,81.5984932429369003956740485743,
                    180.660921136027535395785839907,351.655764263892252842839651265,708.372731161425308227888121379,
                    1419.35511469503772497791180454,2820.07381441318329072396683246,5692.7484188562664576362174619,
                    11262.0708884645275214277198799,22792.920185013752781878095037,45022.3296408593794224273585843,
                    91200.6645502016478738150942537,180061.53669307333572216785011,364817.686329154112527579728712,
                    720271.224406785361574035315945,1.45914512353389683088772652022e6,2.88143939103311008962008174208e6,
                    5.835731463011793220802542156e6,1.15276418196987813086059961834e7,2.33389190541505058395169506151e7,
                    4.6118864385224818788697202808e7,9.33387944282811926635306225895e7,1.84509390736817045480418506483e8,
                    3.73287560320159592182559847383e8,7.38171449955975463725254532629e8,1.49288637976900322070901973986e9,
                    2.95320398123015575519976251299e9,5.97053060156614688371323994114e9,1.18147997499775288964288494881e10,
                    2.38782506063087985122395608337e10,4.72667468539848498486070952932e10,9.54983008842111581575442258659e10,
                    1.89095604891760237675700320831e11,3.81937522376337112445766634741e11,7.56490726671082647813570989992e11,
                    1.52753945577286706611854819853e12,3.02637250850646295470257399222e12,6.10936131245665983230545903716e12,
                    1.21070390056769134593224399285e13,2.44344326675621093812410711307e13,4.8434015989553725309967839147e13,
                    9.77263301126973314170383959517e13,1.93758248111635267066882293275e14,3.90862143262071024406653704289e14,
                    7.7511704823738789767424958222e14,1.56328489584347726904072033805e15,3.10078699439764637046745086435e15,
                    6.25251847543256773811364220232e15,1.24043583831241067308705665974e16,2.50077144379083408392719885287e16
                  },0)
                },0)
              ),
              // AgqQ_reg as^2 LM^0 NF^0 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^2 LM^0 NF^0 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    8.49382716049382716049382716049,2.56790123456790123456790123457,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951,
                    11.0617283950617283950617283951,11.0617283950617283950617283951,11.0617283950617283950617283951
                  },0),
                  // AgqQ_reg as^2 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    4.74074074074074074074074074074,1.18518518518518518518518518519,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593,
                    5.92592592592592592592592592593,5.92592592592592592592592592593,5.92592592592592592592592592593
                  },0),
                  // AgqQ_reg as^2 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^2
                  ome_x<double>({
                    0.888888888888888888888888888889,0.888888888888888888888888888889,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778,
                    1.77777777777777777777777777778,1.77777777777777777777777777778,1.77777777777777777777777777778
                  },0)
                },0)
              )
            }
          )
        },0),
        // AgqQ_reg as^2 LM^1
        ome_nf<double>({
          // AgqQ_reg as^2 LM^1 NF^0
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^2 LM^1 NF^0 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^2 LM^1 NF^0 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^1 NF^0 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    11.8518518518518518518518518519,-18.962962962962962962962962963,13.037037037037037037037037037,
                    -2.37037037037037037037037037037,-1.18518518518518518518518518519,-0.82962962962962962962962962963,
                    -0.651851851851851851851851851852,-0.541798941798941798941798941799,-0.465608465608465608465608465608,
                    -0.409171075837742504409171075838,-0.365432098765432098765432098765,-0.330415263748597081930415263749,
                    -0.301683501683501683501683501684,-0.277648277648277648277648277648,-0.257224257224257224257224257224,
                    -0.239641839641839641839641839642,-0.224338624338624338624338624339,-0.21089324618736383442265795207,
                    -0.19898329702251270878721859114,-0.188357604250277108894239957191,-0.178817413905133203378817413905,
                    -0.170203285992759676970203285993,-0.162385762385762385762385762386,-0.155258590041198736850910763954,
                    -0.148733713951105255453081540038,-0.142737520128824476650563607085,-0.137207977207977207977207977208,
                    -0.132092434314656536878759100981,-0.127345905123682901460679238457,-0.122929716799448600214883739788,
                    -0.118810436051815362160189746397,-0.114959007951221521855559675359,-0.111350059737156511350059737157,
                    -0.107961333767785380688606495058,-0.104773222420281243810655575361,-0.10176838412132529779588603118,
                    -0.0989314244216205000518726008922,-0.096248629581962915296248629582,-0.0937077428305498480937077428305,
                    -0.0912977755083018240912977755083,-0.0890088469035837456890088469036,-0.0868320478076575637551247307345,
                    -0.0847593237837140276164666408569,-0.0827833748990867538797204140369,-0.0808975692696622929181068715952,
                    -0.0790958682431550648604912170804,-0.0773727614307324452251988483873,-0.0757232101049554025193864848506,
                    -0.074142597731866927056566279508,-0.0726266866084495263826570817452,-0.0711715797430083144368858654573,
                    -0.0697736872526788493175047796897
                  },-1)
                },0)
              ),
              // AgqQ_reg as^2 LM^1 NF^0 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^2 LM^1 NF^0 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^1 NF^0 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    10.4312843209486342866616404019,38.4520405431006204470970375696,87.3824418762934816568567207941,
                    198.468587456290667017417145292,366.122360097766519220019475769,779.652127602940445847446358945,
                    1481.55610705773274354674456974,3095.31115591440728603528808128,5960.3577615642643075203116123,
                    12328.8707259503981035767784274,23924.506883999561639252322287,49180.328695047383660098359501,
                    95922.365134202511427940826746,196342.79812647288092373972135,384336.891328240837142554737776,
                    784250.919001618019421454611896,1.53930537662757466318153356242e6,3.13355248048826479477265318148e6,
                    6.16335152597289850748285175148e6,1.25232484291578904293294997037e7,2.46731247673476197150592593164e7,
                    5.00573306472544151471889725454e7,9.87573126984530342385818892948e7,2.00111011585918344445961486518e8,
                    3.95246084126894513514964987382e8,8.00045185286961329767450103588e8,1.58172051401232721793945861174e9,
                    3.19881774503697479042149398827e9,6.32941284897052625519744010877e9,1.27905592728938794478961472207e10,
                    2.53264452432517253252759653044e10,5.11457863337303117719625514633e10,1.01336622352304590091117456571e11,
                    2.04525207359667288143568680926e11,4.05455538151123841488167268765e11,8.1789521348459530485336817927e11,
                    1.6222105058895420154542205374e12,3.27084618756080190115964465933e12,6.49023392565418831665818962247e12,
                    1.30807438749898358415124295463e13,2.59659530685368262789162970892e13,5.23134307889020337405162158865e13,
                    1.03881991757030691613178515139e14,2.09219056929037974234349059256e14,4.15594145126307544249317154945e14,
                    8.36749759137065969191871592275e14,1.66261851126929215488693939789e15,3.34653578118334764497562386493e15,
                    6.65136192196947778900904743521e15,1.33844399031786069524250680616e16,2.66087177787927529548879290181e16
                  },0)
                },0)
              ),
              // AgqQ_reg as^2 LM^1 NF^0 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^2 LM^1 NF^0 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^1 NF^0 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    9.48148148148148148148148148148,2.37037037037037037037037037037,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519,
                    11.8518518518518518518518518519,11.8518518518518518518518518519,11.8518518518518518518518518519
                  },0),
                  // AgqQ_reg as^2 LM^1 NF^0 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    3.55555555555555555555555555556,3.55555555555555555555555555556,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111,
                    7.11111111111111111111111111111,7.11111111111111111111111111111,7.11111111111111111111111111111
                  },0)
                },0)
              )
            }
          )
        },0),
        // AgqQ_reg as^2 LM^2
        ome_nf<double>({
          // AgqQ_reg as^2 LM^2 NF^0
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^2 LM^2 NF^0 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^2 LM^2 NF^0 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^2 NF^0 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    3.55555555555555555555555555556,-3.55555555555555555555555555556,1.77777777777777777777777777778
                  },-1)
                },0)
              ),
              // AgqQ_reg as^2 LM^2 NF^0 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^2 LM^2 NF^0 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^2 NF^0 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    4.44444444444444444444444444444,12.4444444444444444444444444444,28.4444444444444444444444444444,
                    56.8888888888888888888888888889,113.777777777777777777777777778,227.555555555555555555555555556,
                    455.111111111111111111111111111,910.222222222222222222222222222,1820.44444444444444444444444444,
                    3640.88888888888888888888888889,7281.77777777777777777777777778,14563.5555555555555555555555556,
                    29127.1111111111111111111111111,58254.2222222222222222222222222,116508.444444444444444444444444,
                    233016.888888888888888888888889,466033.777777777777777777777778,932067.555555555555555555555556,
                    1.86413511111111111111111111111e6,3.72827022222222222222222222222e6,7.45654044444444444444444444444e6,
                    1.49130808888888888888888888889e7,2.98261617777777777777777777778e7,5.96523235555555555555555555556e7,
                    1.19304647111111111111111111111e8,2.38609294222222222222222222222e8,4.77218588444444444444444444444e8,
                    9.54437176888888888888888888889e8,1.90887435377777777777777777778e9,3.81774870755555555555555555556e9,
                    7.63549741511111111111111111111e9,1.52709948302222222222222222222e10,3.05419896604444444444444444444e10,
                    6.10839793208888888888888888889e10,1.22167958641777777777777777778e11,2.44335917283555555555555555556e11,
                    4.88671834567111111111111111111e11,9.77343669134222222222222222222e11,1.95468733826844444444444444444e12,
                    3.90937467653688888888888888889e12,7.81874935307377777777777777778e12,1.56374987061475555555555555556e13,
                    3.12749974122951111111111111111e13,6.25499948245902222222222222222e13,1.25099989649180444444444444444e14,
                    2.50199979298360888888888888889e14,5.00399958596721777777777777778e14,1.00079991719344355555555555556e15,
                    2.00159983438688711111111111111e15,4.00319966877377422222222222222e15,8.00639933754754844444444444444e15
                  },0)
                },0)
              ),
              // AgqQ_reg as^2 LM^2 NF^0 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^2 LM^2 NF^0 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^2 LM^2 NF^0 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    1.77777777777777777777777777778,1.77777777777777777777777777778,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556,
                    3.55555555555555555555555555556,3.55555555555555555555555555556,3.55555555555555555555555555556
                  },0)
                },0)
              )
            }
          )
        },0)
      },0),
      // AgqQ_reg as^3
      ome_logm<double>({
        // AgqQ_reg as^3 LM^0
        ome_nf<double>({
          // AgqQ_reg as^3 LM^0 NF^0
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^3 LM^0 NF^0 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^3 LM^0 NF^0 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^0 NF^0 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    -678.682381056504968950891436704,2103.42757982657281168223542486,-1916.26444303514054392991680763,
                    568.763937810286905662627670925,-18.9835029930838440938682125494,23.8344472143758771378872848339,
                    -2.53542718728455644986478212974,0.0706369575122060528331147570012,13.6819161911788480210892796215,
                    -7.55368155307619378105256064394,19.1773944705554683417285178851,-10.8506374612731415908055170946,
                    21.0393127341993378689416549984,-12.3932369653220750329764561944,21.4505908252045987438258402812,
                    -13.1026786528493470026165104157,21.2265875789630944267520776296,-13.3789661180132733138249570137,
                    20.7158473053028763153248502821,-13.4172244669149127635277328117,20.0810108477283626546283644989,
                    -13.3205169642629041923968506267,19.4022556909019218460824722587,-13.1464919683671184167294369635,
                    18.7202502899289058689871174992,-12.9287831637580940333880429646,18.0556520106354474902308198755,
                    -12.6876262433645845567020209287,17.4185772959177238588563500284,-12.4354588356833022011101400511,
                    16.813448037283408677569712055,-12.1800236938055301643038934901,16.2415740758111708223781134705,
                    -11.9261593633441890698683772061,15.7025724074973614960357272334,-11.6768688266317957275170147196,
                    15.1951647155272953172037466568,-11.4339754919910023238563861336,14.7176319327833927651416475025,
                    -11.1985355423989758747782774958,14.2680746984592346243590029099,-10.9711023570617892758352599701,
                    13.8445617525099761118873578075,-10.7518989227906214876890921081,13.4452126895143128071520939019,
                    -10.5409318029589153294368843493,13.0682419190915507261614991832,-10.3380673028244942031593464565,
                    12.7119796395380722908303486584,-10.143082791575164119284441266,12.3748792629237690159447972049,
                    -9.95570147350260822374343047556
                  },-1),
                  // AgqQ_reg as^3 LM^0 NF^0 [below 0.375] y=x H[0, y]^1
                  ome_x<double>({
                    -153.559798605251777318007964609,-38.2687033157460505611385150029,513.145919412345598438134279944,
                    -211.986934028416328155948585437,8.28640331276138859388030282212,15.3117402214865264872825582823,
                    11.46303724431937121242683753,10.0981944440093148992014540091,7.31000646077678105150620286396,
                    10.0469125086617212702404975559,4.58832384207009371293222631817,10.4799731208723161394537600442,
                    2.80189036918007619935905055372,10.8431577087071395517856499792,1.57570749581889056677194498804,
                    11.0796121597630652648608806641,0.700763101083061013481578361274,11.2095380304482070395954967699,
                    0.0569565289600401430084878301312,11.2602979459600806080450586385,-0.428390323224005703248552954499,
                    11.2543054161115801194635549756,-0.801384706854568464571322854599,11.2081599031490510818745708149,
                    -1.0924731752217223100975111447,11.1338188360581741762713666623,-1.32245026452381287363666077488,
                    11.0398505304862791779626364425,-1.50592933533619325985002325794,10.932411404684652028747208242,
                    -1.65343318813580318561851349893,10.81594853641545564806243081,-1.77269992358621028374499475581,
                    10.6936924833170301046047307643,-1.86952479602792261395526998694,10.5680015111012159069333578857,
                    -1.94831812141604323567282555487,10.4406030352681349617322968185,-2.01248409572877253153625024138,
                    10.3127642909544749953243540402,-2.06468360026501399731580157093,10.1854140924139962949644272736,
                    -2.10702004821495397225879603272,10.0592305353440748376761102027,-2.14117308416120282268127965773,
                    9.93470476571773761703771370844,-2.16849626794726637413797058879,9.81218776420459350086904984234,
                    -2.19008945151186657540296620451,9.69192495883355218134822078399,-2.20685309273976428889531642634,
                    9.57408203129528801328106620319
                  },-1),
                  // AgqQ_reg as^3 LM^0 NF^0 [below 0.375] y=x H[0, y]^2
                  ome_x<double>({
                    86.7301963893319292329527060593,-127.06974978959878935360255484,3.10562414266117969821673525377,
                    8.20027434842249657064471879287,-7.86853223593964334705075445816,1.98770919067215363511659807956,
                    -5.44901206573164245121917079589,0.448386103412558439013465468492,-4.26814416667444298249706832893,
                    -0.166743310035491105449953186579,-3.5425098579905912491022853911,-0.462561114323179418577847309609,
                    -3.04373782583227459339002117212,-0.617877623263595316234035539422,-2.67687922221635422015289061823,
                    -0.70236744182070460730072723018,-2.39435330375438880714629358313,-0.747925340494793838949220782419,
                    -2.16936099961237241620480036054,-0.770793677891943756888962083656,-1.98553618338242724572613047958,
                    -0.779895969878695558015664580353,-1.83226826314308517126039226981,-0.78042015073137040349651162255,
                    -1.70235264836929063968114700447,-0.775516209033306923725584099824,-1.59071335310709045483952508419,
                    -0.76716261643644400780094771521,-1.49366420519115364056225849244,-0.756635437944806846667499521095,
                    -1.40845973524605850632699597619,-0.744774829029594097279191539938,-1.33301070428440729705581269941,
                    -0.732142587351106614748601316414,-1.26569757730346096790504828128,-0.719118474686766424252984741637,
                    -1.20524455513154419400739451947,-0.705960811246152036300514076908,-1.15063230369619155927931325897,
                    -0.692845541009928055454856505435,-1.10103612837885947529448067915,-0.679891956944239914703402645735,
                    -1.05578130366575361878916122747,-0.667179956155565556229944498988,-1.0143102278059863142844177347,
                    -0.654761800520721428998008005464,-0.97615789057640065124494472333,-0.642670244742432089330918383642,
                    -0.940933289455357503847294160285,-0.630924222001821452352051847051,-0.908305170570037960677866628092
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^0 [below 0.375] y=x H[0, y]^3
                  ome_x<double>({
                    -14.2661179698216735253772290809,24.2853223593964334705075445816,2.28257887517146776406035665295,
                    0.79012345679012345679012345679,0.82962962962962962962962962963,0.434567901234567901234567901235,
                    0.541798941798941798941798941799,0.310405643738977072310405643739,0.409171075837742504409171075838,
                    0.243621399176954732510288065844,0.330415263748597081930415263749,0.201122334455667789001122334456,
                    0.277648277648277648277648277648,0.17148283814950481617148283815,0.239641839641839641839641839642,
                    0.149559082892416225749559082892,0.21089324618736383442265795207,0.132655531348341805858145727427,
                    0.188357604250277108894239957191,0.119211609270088802252544942603,0.170203285992759676970203285993,
                    0.108257174923841590508257174924,0.155258590041198736850910763954,0.099155809300736836968721026692,
                    0.142737520128824476650563607085,0.0914719848053181386514719848053,0.132092434314656536878759100981,
                    0.084897270082455267640452825638,0.122929716799448600214883739788,0.0792069573678769081067931642644,
                    0.114959007951221521855559675359,0.074233373158104340900039824771,0.107961333767785380688606495058,
                    0.069848814946854162540437050241,0.10176838412132529779588603118,0.0659542829477470000345817339281,
                    0.096248629581962915296248629582,0.0624718285536998987291384952204,0.0912977755083018240912977755083,
                    0.0593392312690558304593392312691,0.0868320478076575637551247307345,0.0565062158558093517443110939046,
                    0.0827833748990867538797204140369,0.0539317128464415286120712477302,0.0790958682431550648604912170804,
                    0.0515818409538216301501325655915,0.0757232101049554025193864848506,0.049428398487911284704377519672,
                    0.0726266866084495263826570817452,0.0474477198286722096245905769715,0.0697736872526788493175047796897
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^0 [below 0.375] y=x H[0, y]^4
                  ome_x<double>({
                    -1.11934156378600823045267489712,2.04115226337448559670781893004
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^0 NF^0 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^3 LM^0 NF^0 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^0 NF^0 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    37.9369032405002989208256115383,-299.032878449718564050021258822,-1154.2145399028593895671192994,
                    -2712.22488715861034218954101847,-4662.38336982631683340061703073,-8953.16875428124982061726120409,
                    -14759.8543826321147481646637954,-28128.892664908732019328925146,-45484.1098183307459332553068802,
                    -86196.8329225892122847827408374,-135634.109852270942090536926843,-254249.028736184961696619140064,
                    -381512.478431260564860253041808,-700572.709195139570822834556557,-956745.883181026744368464688011,
                    -1.68002200129686676544706916438e6,-1.78693246576806631348944943694e6,-2.690051614394892112476871593e6,
                    243058.328660290247693628447831,3.86144427964827367032446734661e6,2.79903850327419720307345934378e7,
                    6.89543368044102537409835034941e7,2.11477315325338142926252677272e8,4.73062025907242889473581203574e8,
                    1.21480069233242605910356537692e9,2.62385515576745282267731406791e9,6.234129391664084067626263691e9,
                    1.32235707632509625433006967445e10,3.00854486845043129961961974045e10,6.31148299603014934854136091825e10,
                    1.39703910916783951509453507412e11,2.90905001649678498449243018919e11,6.31891495044213849358794628797e11,
                    1.30875865263882040491856908497e12,2.80426202764649170353911858313e12,5.78471292324696199118540568681e12,
                    1.22677773604787077164865398408e13,2.52266032744815314006538474114e13,5.30715770561126405851885934973e13,
                    1.08856463404267173643948327374e14,2.27553327125748676876866612665e14,4.65770034948720419819145368265e14,
                    9.68601851601514321867139202849e14,1.97914573489783287459782967929e15,4.09815047919516055485586430077e15,
                    8.36139244868889559657616542046e15,1.72515145790078885227420614241e16,3.51532911716473655210903916925e16,
                    7.23087142458872542612441488651e16,1.4718021954573013919720268311e17,3.01953398355307705318629574201e17
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^0 NF^0 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^3 LM^0 NF^0 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    136.455201888561545206185019201,60.0884583652722593158130631992,120.096438521859572219773765984,
                    -85.1534138125251049423625027626,-136.715667758827868527556764903,-148.201280814762038939487546097,
                    -148.608959927986628614859154309,-144.818008382019880152328496396,-139.232648807990481623078556442,
                    -132.869341687913761181519892141,-126.207424100612383786813920241,-119.488395813935187417567655536,
                    -112.838389117511936461719360658,-106.32378275859956673730384452,-99.9784871259211353493738205555,
                    -93.8181608695656593011271717874,-87.8479812777582396279264066644,-82.0670564435308156272976182009,
                    -76.4710066681946775287723637108,-71.0535091401717385761423050324,-65.8072366137968079902494646206,
                    -60.7244324580800689699704144049,-55.7972628332601823384199065834,-51.0180299857460105028006108306,
                    -46.3792979482375223586982045295,-41.873962576238277979659619505,-37.4952861237158765905252590118,
                    -33.2369093068794679030634051129,-29.0928492387084893089105576151,-25.0574886984394852043099047375,
                    -21.1255603109425637309839710447,-17.2921279748115624557444772466,-13.5525670623234151703861304012,
                    -9.90254437273687728606984886601,-6.33799845916568358068379777478,-2.85512070795224240657790083974,
                    0.549662611539353341657242586959,3.87970721570911161814948363801,7.13816655340702859486269417443,
                    10.3280064353420380781876094228,13.4520185796912113616394527186,16.5128331246248433865983118811,
                    19.5129301680017936492634227728,22.4546503988787198831070567281,25.3402048864697867307058245687,
                    28.1716840909952969662211850875,30.951066158295502669833059245,33.6802245567246571012756628411,
                    36.3609351110633423004705960421,38.9948824842510291659402242693,41.5836661538146476200017609915
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    72.4459587795985427974020032719,101.196884321637650722037883041,161.519760841021005071499174588,
                    280.089622587672466050680803223,324.871944839022226687267375272,361.43137189249839694559185821,
                    390.902482643852763224747174887,415.221824326170847373992871999,435.809527711302765997713402666,
                    453.620199365310445659573813552,469.301782541676676076602005613,483.30648614291453731839857748,
                    495.95944482873232342075928922,507.500977900793050336676854262,518.113143334719362866271865056,
                    527.936881654235114074559301783,537.083393268111536651899914367,545.641888020386649449043171273,
                    553.684991721201492498750553179,561.272601241456269831051415775,568.454688264909322567271351544,
                    575.2733752563120970673956021,581.764497700736291404950655476,587.958797180260310013624721726,
                    593.882844804176783051140181692,599.559764713914018410368670337,605.009807308902846557335841747,
                    610.250808078421670840559396768,615.298558339024254840451260278,620.16710740036370968784102558,
                    624.869010824457769706796343966,629.415535916240301956865787143,633.816832991333844767583933808,
                    638.08207904098236398199379384,642.219598967925995108170734415,646.23696847050099873919594879,
                    650.141101813241533089757870842,653.938327074754780416728126949,657.634450959887825610012011147,
                    661.234814868307312911750180856,664.744343599828393935774829338,668.167587828998405445480774774,
                    671.508761283183639866381355867,674.771773398846676173485538292,677.960258101546115000962655251,
                    681.077599250075829799254455785,684.126953199168320421907986877,687.111268864490504451092790098,
                    690.033305615266060928129790538,692.895649271409171159312251727,695.700726441684307159062620032
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^2
                  ome_x<double>({
                    25.1431940783560578593195015748,48.5506014857634652667269089822,58.3686926834610869120546410096,
                    89.8529176491675340999833104198,102.675962916657246034139689021,113.013411476328027927143804247,
                    121.001340145738178818776177361,127.358743353935037809285961521,132.576405233501679280689337687,
                    136.970132546100244083839678756,140.747812342298558800672914108,144.051043416693119914276302309,
                    146.979544244542993779862850275,149.605676917399313359829153888,151.983350945487125194758069267,
                    154.153659225468448182457396676,156.148556486741370831041420921,157.993328791549986460281596276,
                    159.708294351668764002557223334,161.310002217115378981900399056,162.812094743221595239804020459,
                    164.225939795081395011695683585,165.561101993971898621216086426,166.825699339688834938216120007,
                    168.02667681393552191989916524,169.170018923298787439057214196,170.260916697289111262003557879,
                    171.303900275045668745619628676,172.302945184125863795053175728,173.261558286868277874974707824,
                    174.182847854103688644151536462,175.069581132065732327954154263,175.924231969164353624293362302,
                    176.749020478772967980150719013,177.54594627323246253696442412,178.316816471754700691952400957,
                    179.063269431835016258455178368,179.786794959506796491529029057,180.488751603418078287601226989,
                    181.170381520459045828490579392,181.832823308572838620369117901,182.477123129550740322182862871,
                    183.104244386645487432905259517,183.715076175417891049540981196,184.310440688846074607479325127,
                    184.891099727450051090810673099,185.457760440536865812567824367,186.011080404506536429471024696,
                    186.551672127582918778293360549,187.080107056645865817085971143,187.596919150489348439403618718
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^3
                  ome_x<double>({
                    -0.37311385459533607681755829904,8.51577503429355281207133058985,12.2578875171467764060356652949,
                    23.012345679012345679012345679,27.045267489711934156378600823,29.7887517146776406035665294925,
                    31.8244170096021947873799725652,33.4234763864393494023123652753,34.7305506564765824025083284343,
                    35.8305571885818799399046312627,36.7770592461950486641844666536,37.6057577539059020540502021984,
                    38.3415103415103415103415103415,39.0022487800265578043355821134,39.6012930827745642560457375272,
                    40.1487840006358524877043395562,40.6526052822349118645414941711,41.1189976004790819605634420449,
                    41.5529791492028238578710619379,41.9586410890924899257277042699,42.3393596894342221712936080554,
                    42.6979516425680856365410888975,43.0367897591839800302132603474,43.3578904989642689248338779954,
                    43.6629811261160877030616739366,43.9535518901168033945986215122,44.2308970415360288878981888858,
                    44.4961474137822859078177602457,44.7502965568203178347385760554,44.9942218865977832727687934244,
                    45.2287019422783498090185991505,45.4544305753731922722951249657,45.6720286991839803602636191468,
                    45.8820540824471175063159384369,46.085009563180376017352227251,46.2813499773857139699668683188,
                    46.4714880353523848632899039802,46.655799330774791396807548609,46.8346266311095641162456351686,
                    47.0082835688990555119748476619,47.1770578312405693154769552576,47.3412139267426449440205151471,
                    47.5009955951060657916328713413,47.6566279130879142455302317742,47.8083191414371959970514320465,
                    47.9562623499603643249857662972,48.1006368518227888927243436976,48.2416094732376163125083080619,
                    48.3793356806181218394488917058,48.5139605839022618679320475108,48.6456198319636263314129131081
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^0 [above 0.625] y=1-x H[0, y]^4
                  ome_x<double>({
                    -0.411522633744855967078189300412,-0.411522633744855967078189300412,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823,
                    -0.823045267489711934156378600823,-0.823045267489711934156378600823,-0.823045267489711934156378600823
                  },0)
                },0)
              )
            }
          ),
          // AgqQ_reg as^3 LM^0 NF^1
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^3 LM^0 NF^1 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^3 LM^0 NF^1 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^0 NF^1 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    10.2150416607980364864664192802,-11.7952885743782834000466661938,2.73715046002864787286283926973,
                    4.21399176954732510288065843621,-0.724279835390946502057613168724,-1.10617283950617283950617283951,
                    -1.07764060356652949245541838134,-0.985145992553399960807368214776,-0.888732118361747991377621007251,
                    -0.801149650532366581749297798681,-0.72438434907570710039845842315,-0.657688466083527811922873651269,
                    -0.599713299960213540460454040701,-0.549117951093259735235043877019,-0.504730658981099898207481981732,
                    -0.465571412844781451483391518665,-0.430833193718555270583489278375,-0.399853949587739399960584730789,
                    -0.372089686572170923936321453353,-0.347091549354796434004985521759,-0.324487186534468865204928311132,
                    -0.303965890284332411020643168603,-0.285266829451269108164413985957,-0.268169740630745414458076135043,
                    -0.252487546353959489827165710695,-0.238060477181627333408196665107,-0.224751367470872830310287682266,
                    -0.212441869511803825785076279724,-0.201029389223960866932803999835,-0.190424591581526066479982346155,
                    -0.180549358301696317567611765097,-0.171335106546543400442288412859,-0.162721397413563205429080673154,
                    -0.154654778332069155389301463975,-0.14708781528433463538084879544,-0.139978279891210382473022091762,
                    -0.133288463486188740562583074495,-0.126984595832840812789710699647,-0.121036350481432780366128711962,
                    -0.115416422185194104927720689121,-0.110100164512467292897121501247,-0.105065277955434467610424488277,
                    -0.100291540569694101800429903127,-0.0957605745740739031860002539488,-0.0914556434680040837302757740897,
                    -0.0873614751398000652474628137462,-0.0834641071863584558161075054329,-0.0797507512767421524662732891395,
                    -0.0762096738954081468205604553053,-0.0728300912163384310151223750232,-0.0696020762036716358620202858027,
                    -0.0665164763208299360026165941825
                  },-1)
                },0)
              ),
              // AgqQ_reg as^3 LM^0 NF^1 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^3 LM^0 NF^1 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^0 NF^1 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    10.835456667648905115819398851,35.7252248499033654507923032249,79.5238788001705598997621435777,
                    177.403785400899960289197268225,302.486519396403506077533193903,694.675025929552352516434276725,
                    1238.52799703539637245244558927,2732.71985778381520780601384367,5027.23211437931020722170710309,
                    10813.315563422989712891779657,20300.7442961629253661943500414,42935.7920734873823002663537683,
                    81735.2720488257463349026118035,170840.343743589181128312092553,328488.035690072319471377657219,
                    680686.766974184976933532086506,1.31861357300043511255447286848e6,2.71456472376267642888259134664e6,
                    5.28890916989390012530593405203e6,1.08325754723318746372795177851e7,2.12015093939558231655521186847e7,
                    4.32479806710354603725769764938e7,8.49544068876084196657388633149e7,1.72723454592316500518747827658e8,
                    3.40305175953752308600699439903e8,6.90005445187712972977794063553e8,1.36284444060486950486273700688e9,
                    2.75704545869932487212424793252e9,5.45684792300031289700136326501e9,1.10181033171136302033716596803e10,
                    2.18460039013284686704789447772e10,4.4037967265242592182338450772e10,8.74478949211835938445202595169e10,
                    1.76033180666395020899183235637e11,3.50012498911551188876522989247e11,7.03720823600554172748282496292e11,
                    1.40081919436081332693232337254e12,2.81344469980305935728676272026e12,5.60597113380376117706694978216e12,
                    1.12487257257692474796445086811e13,2.24333733014502822612974020004e13,4.49770631358047664252836056458e13,
                    8.9767072193310009704323993028e13,1.79844979252873631992735199009e14,3.59187639395115447717938093216e14,
                    7.19154566322124844677165896563e14,1.43717645289140886711744246321e15,2.87581263157336706852368348339e15,
                    5.75023107680966155965129914897e15,1.1500360450261110556019683112e16,2.3006404722203406165940007004e16
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^0 NF^1 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^3 LM^0 NF^1 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^0 NF^1 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    5.89764428718914170002333309689,4.31739737360889478644308618331,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802,
                    10.2150416607980364864664192802,10.2150416607980364864664192802,10.2150416607980364864664192802
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^1 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    -3.16049382716049382716049382716,4.74074074074074074074074074074,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358,
                    1.58024691358024691358024691358,1.58024691358024691358024691358,1.58024691358024691358024691358
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^1 [above 0.625] y=1-x H[0, y]^2
                  ome_x<double>({
                    -3.16049382716049382716049382716,-0.79012345679012345679012345679,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395
                  },0),
                  // AgqQ_reg as^3 LM^0 NF^1 [above 0.625] y=1-x H[0, y]^3
                  ome_x<double>({
                    -0.395061728395061728395061728395,-0.395061728395061728395061728395,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679,
                    -0.79012345679012345679012345679,-0.79012345679012345679012345679,-0.79012345679012345679012345679
                  },0)
                },0)
              )
            }
          )
        },0),
        // AgqQ_reg as^3 LM^1
        ome_nf<double>({
          // AgqQ_reg as^3 LM^1 NF^0
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^3 LM^1 NF^0 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^3 LM^1 NF^0 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^1 NF^0 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    -444.696026344289940347867233339,1493.07530203835698018429841674,-1334.69564815797393544277573073,
                    333.98865153319186883202850322,-89.370597711109463498989865459,41.9354374790831186801567529631,
                    -71.5592092623722076679286509682,32.8391298283297018080395970832,-51.1230449396324445596872227696,
                    27.1745178442789067559839764012,-39.4717252290444930075520949834,23.2078297801822189634013414746,
                    -32.0128277186106478585164511128,20.2884862648901313450389855172,-26.8377094650062288641429256092,
                    18.0518411472036825820199486617,-23.041831544120227199946975189,16.2827577759304817549570424477,
                    -20.1423681510275422278526219046,14.8474323336698847426552022287,-17.8580398142280494307704630269,
                    13.6586104606573998034360900903,-16.0138556528362536849284732516,12.6570525674197371806199544426,
                    -14.4952493548950942963068741342,11.8011476240898872106029735386,-13.2240889129498605894274892775,
                    11.0608123975189168016523287185,-12.1452775135061177588088339169,10.4137553228899613544618775418,
                    -11.2188615717286580728538503876,9.84310366679120881657584760616,-10.4151736800166219010089418711,
                    9.33584791353655823547196662525,-9.71173078054317281049504973316,8.88179316730451643093642601873,
                    -9.09118971108770120489335315458,8.47283480266480106517198166113,-8.53996237279598766240509792803,
                    8.10244712794758538482708606451,-8.04725504318268277802646042365,7.76531537408477554260093112952,
                    -7.60438775144968493979498406817,7.45706620630480190299919133496,-7.20430296459498335178972361773,
                    7.17406727733080980393761406141,-6.84120493940713317627512326873,6.91327601144277922131842448408,
                    -6.51029096644363816771728926086,6.67212405184304925611680603495,-6.2075483381348311624837453049,
                    6.44842791736585706665049926903
                  },-1),
                  // AgqQ_reg as^3 LM^1 NF^0 [below 0.375] y=x H[0, y]^1
                  ome_x<double>({
                    -133.633892818892358622155463901,89.8066970023528113773463008979,371.951790913138650562011352858,
                    -86.3868312757201646090534979424,-7.04526748971193415637860082305,8.83226337448559670781893004115,
                    3.99920987654320987654320987654,-7.25496733014193331653649113967,4.87245677812571236607517786656,
                    -11.5503834612426197208258180795,4.78103033074226489864349946654,-12.8847718155671765847598245653,
                    4.51687236202448935843487419857,-13.1790616077117454527922847482,4.2325754330826182051879431423,
                    -13.0507030198908528149494753099,3.96559547079274970431433016467,-12.7416350295439107434493824344,
                    3.72435714374143661019246829631,-12.3582586179417411539578817419,3.50885680860183418528397546754,
                    -11.9506932127716798099506268549,3.31667176716698662638122951754,-11.5433774275011702718944154627,
                    3.14489601019465283271238683899,-11.148239477305372888933953987,2.99077275032798965458340303743,
                    -10.7708341886812785138074436753,2.85187642081317626907116853081,-10.4133784363713340153023417296,
                    2.72613508020007549339516202991,-10.0763218102480253381606936102,2.61179907962263499347796732939,
                    -9.75918707380950293900013584726,2.50739561358590243462193491795,-9.46103055377419145999592102294,
                    2.41168372505076781167044948271,-9.18069777888646050218122944956,2.32361457735196751573201621898,
                    -8.91696578997180676907303887163,2.2422980085000483410280514607,-8.66862142741742543599905939308,
                    2.16697498010416920910352294737,-8.43450293964840144991062665014,2.09699509843255534470392084234,
                    -8.21352042617093707793132082132,2.03179833301432027344137911431,-8.00466407830257627272618796176,
                    1.9709001446478449930880110573,-7.80700546751521084874818571148,1.9138793590119980420192544266,
                    -7.61969498441460062522554595905
                  },-1),
                  // AgqQ_reg as^3 LM^1 NF^0 [below 0.375] y=x H[0, y]^2
                  ome_x<double>({
                    102.320987654320987654320987654,-85.2345679012345679012345679012,-3.16049382716049382716049382716,
                    9.08641975308641975308641975309,-3.5950617283950617283950617284,4.9975308641975308641975308642,
                    -2.34779541446208112874779541446,3.569664902998236331569664903,-1.77307466196355085243974132863,
                    2.8016460905349794238683127572,-1.43179947624392068836513280958,2.31290684624017957351290684624,
                    -1.20314253647586980920314253648,1.97205263871930538597205263872,-1.03844797178130511463844797178,
                    1.71992945326278659611992945326,-0.913870733478576615831517792302,1.52553861050593076736867586541,
                    -0.816216285084534138541706481163,1.37093350660602122590426683994,-0.737547572635291933537547572635,
                    1.24495751162417829084495751162,-0.672787223511861193020613310468,1.14029180695847362514029180696,
                    -0.61852925389157273215244229737,1.05192782526115859449192782526,-0.572400548696844993141289437586,
                    0.976318605948235577865207494837,-0.532695439464277267597829539081,0.910880009730584443228121389041,
                    -0.498155701121959928040758593224,0.853683791318199920350457984867,-0.467832446327069982983961478585,
                    0.803261371888822869215026077771,-0.440996331192409623782172801781,0.758474253899090500397689940174,
                    -0.417077394855172632950410728189,0.718426028367548835385092695034,-0.395623693869307904395623693869,
                    0.682401159594142050282401159594,-0.376272207166516109605540499849,0.649821482341807545059577579903,
                    -0.358727957896042600145455127493,0.620214697734077579038819348897,-0.342748762387005281062128607348,
                    0.593191170968948746726524504302,-0.328133910454806744250674767686,0.568426582610979774100341476228,
                    -0.314715641969947947658180687563,0.545648778029730410682791635173,-0.302352644761608347042520711988
                  },0),
                  // AgqQ_reg as^3 LM^1 NF^0 [below 0.375] y=x H[0, y]^3
                  ome_x<double>({
                    -5.92592592592592592592592592593,6.51851851851851851851851851852
                  },0),
                  // AgqQ_reg as^3 LM^1 NF^0 [below 0.375] y=x H[0, y]^4
                  ome_x<double>({
                    -2.37037037037037037037037037037,1.18518518518518518518518518519
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^1 NF^0 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^3 LM^1 NF^0 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^1 NF^0 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    47.7841282566074322120325584226,94.1499663820302961337763297519,-259.075647679071289451440418944,
                    -116.032227489075918345966074576,-418.243272984204433453959379096,986.057889248044181584500776125,
                    2426.24575141378897825656247154,10129.959858688133910208984441,23154.9437805977678130958290108,
                    62355.0928408830720617984894987,137496.671855117018607226994427,326054.981971757442774986639878,
                    703910.997334386215079055493118,1.57568788777995962385765744159e6,3.35475759955255246218093057132e6,
                    7.27474585241098219705270422143e6,1.5337584046925325370956900089e7,3.26135899275027251270145546204e7,
                    6.82637697149104843897200413099e7,1.43272218250874110439063705711e8,2.98219179530086707011399837689e8,
                    6.20183366393256025976966049225e8,1.28523196121934701216852666566e9,2.65487502676536601629564835942e9,
                    5.48223819061746754009112938893e9,1.12670089779153942841650133633e10,2.31977523762805272227286161808e10,
                    4.74873660175294089409186007233e10,9.75316335372936830830080649553e10,1.99027600403138770571618464466e11,
                    4.07917120716067229783581030413e11,8.3030344681064832067981301988e11,1.69869143333879092009316479227e12,
                    3.45043657956444159776944369041e12,7.04812813084642409936355301799e12,1.42915799193211718528612760659e13,
                    2.91532054969662438340343411504e13,5.90281327852396644509906678692e13,1.20265526061262209675683805667e14,
                    2.43206507558554360255640248287e14,4.949845489642723939198579338e14,9.99916256035159340770712394395e14,
                    2.03312301909599969092752058821e15,4.10332719494882706303200249284e15,8.33607653539907499418538110846e15,
                    1.68107187742592322003005897338e16,3.41250539919205527809639685846e16,6.87690451503021319831892318781e16,
                    1.39499517397283330420671983577e17,2.80946409851880251456570779886e17,5.69538294120585452381114315056e17
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^1 NF^0 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^3 LM^1 NF^0 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^1 NF^0 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    220.820485521892380902763117017,10.1869186856893023644638497898,161.946060748252515563375915317,
                    36.4126303125045800911072050841,-3.1562332613947228714873884136,-10.167857127019263558992708025,
                    -8.78418085073931348505184938126,-4.35115458277691895288914040689,1.33734520689208821270656406951,
                    7.53311136873048744745917200652,13.8849950783355525128261814425,20.2173346311565480383053594769,
                    26.4399133546955002452809416749,32.5069032156877252404549440733,38.39657206876763962062855967,
                    44.100642150472654279786104279,49.618453621223511047194471375,54.953651321368749941083980095,
                    60.1122542899218400693217471366,65.1015099816558258055654005143,69.929206790703194276844687773,
                    74.6032606320920523566788061347,79.131468515416786730576356883,83.5213653252077914731822092023,
                    87.7801449879135461065762173878,91.9146219702981883224493212108,95.9312179810593794520256549342,
                    99.8359642505575851814949519932,103.634513213742502637306646176,107.332155616054638344752712385,
                    110.933840475197255740095751931,114.444196250520920360057178905,117.867552173671805297042315115,
                    121.207959090377396386549870905,124.469209424605514050830277028,127.654856048520015526374066912,
                    130.768229954482423404606200811,133.812456698467966270669876419,136.790471630749688466874831441,
                    139.705033958315051627325041545,142.558739700121307287705932525,145.35403360496882951424943389,
                    148.093220105193457747102240304,150.77847337939926528480492307,153.411846595330159958073704651,
                    155.995280400565823648236000023,158.530610724600419357464601834,161.019575951405871291502008636,
                    163.463823517048450570481888135,165.864915982480509545391047861,168.224336627369101943996906792
                  },0),
                  // AgqQ_reg as^3 LM^1 NF^0 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    138.526275025578855121148833251,172.896645395949225491519203622,256.509340174614500365754456626,
                    332.62456651206306003653634963,363.241850462680343987153633581,389.624072684902566209375855803,
                    410.356912191075405715548695309,427.027872130606773818345369535,440.823723310584937923493601667,
                    452.525635351121317084211386723,462.652760571579870876098511944,471.560692372081817281524154182,
                    479.500904239027690961404568069,486.656556057927373108949963478,493.164581649062447353507317518,
                    499.129720990644605378481785309,504.633701783514287137052432769,509.741416007481234415546552221,
                    514.50518185496701867778431525,518.967752983837599815286017077,523.16448780431161002490114887,
                    527.124942722864766950045834955,530.874062152109643503032234751,534.433080936927318288804617947,
                    537.820218166661516200685532492,541.051217301160288039220656701,544.13977146813428214021647078,
                    547.097861850883711699818309811,549.936029514346488259565551669,552.663595695357205878140671538,
                    555.2888417880544788379238966,557.819157513340117616955874414,560.261163755627848658688427149,
                    562.62081506629201773679769356,564.903485723785542528202533488,567.11404240265122654467624729,
                    569.256905865200878509261935919,571.336103598847110672645197183,573.355314941718593711254116637,
                    575.317909942162614667031402603,577.226982964109310993179939878,579.085381865256160968927412798,
                    580.895733427569684128559673355,582.660465601354127879605100044,584.381827028767737557105876984,
                    586.061904235327052770628726532,587.702636814892554943903143231,589.305830881977782408639805586,
                    590.87317102270824726421495393,592.406230940602018649998485809,593.906482964150474730191533922
                  },0),
                  // AgqQ_reg as^3 LM^1 NF^0 [above 0.625] y=1-x H[0, y]^2
                  ome_x<double>({
                    10.5679012345679012345679012346,41.975308641975308641975308642,70.024691358024691358024691358,
                    115.851851851851851851851851852,134.37037037037037037037037037,147.091358024691358024691358025,
                    156.641975308641975308641975309,164.227160493827160493827160494,170.488183421516754850088183422,
                    175.802704291593180482069370958,180.410111699000587889476778366,184.470867404200737534070867404,
                    188.097354497354497354497354497,191.371117771117771117771117771,194.353074086407419740753074086,
                    197.089853356520023186689853357,199.617895684562351229017895685,201.966189075992997561625012605,
                    204.15815228102809802156207385,206.212962509212939209499237019,208.14651243340613825649945361,
                    209.97211347304226561192505465,211.701021601950394520053962778,213.342836577098703001695777754,
                    214.905809299491715249780489606,216.3970814411986395654004574,217.822874413658278691706250373,
                    219.188639951275668160947571466,220.499181261816978702258112776,221.758751359190396637617299732,
                    222.971133527894404422084854315,224.139707663765536956176231565,225.267505353928818517307255062,
                    226.35725590604496203130023142,227.411425048449398553383812327,228.432247653585729179910517285,
                    229.421755555511931759707868325,230.381801315557691805467914085,231.314078622104004199733524724,
                    232.220139879042454120639585981,233.101411433998219602720857536,233.959206815370836747695726088,
                    234.794738281796612116560525848,235.60912693509621266516844825,236.403411605349875166892965479,
                    237.178556682303737650471211331,237.935459039206094552828113688,238.674954172122956321575996116,
                    239.39782165878939457605001994,240.104790025332229203991030859,240.796541096130919050299924787
                  },0),
                  // AgqQ_reg as^3 LM^1 NF^0 [above 0.625] y=1-x H[0, y]^3
                  ome_x<double>({
                    -0.987654320987654320987654320988,-0.987654320987654320987654320988,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198,
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-1.97530864197530864197530864198
                  },0)
                },0)
              )
            }
          ),
          // AgqQ_reg as^3 LM^1 NF^1
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^3 LM^1 NF^1 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^3 LM^1 NF^1 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^1 NF^1 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    26.8641975308641975308641975309,-15.012345679012345679012345679,-1.97530864197530864197530864198,
                    7.50617283950617283950617283951,2.27160493827160493827160493827,1.08641975308641975308641975309,
                    0.602469135802469135802469135802,0.349911816578483245149911816578,0.199647266313932980599647266314,
                    0.102763080540858318636096413874,0.0368489124044679600235155790711,-0.00973117417561862006306450750895,
                    -0.0435765057987280209502431724654,-0.068683579794690905802016913128,-0.0875983510904145824780745415666,
                    -0.102009407406232803058199883597,-0.113074098153463232828312193392,-0.121607625006317816775333115202,
                    -0.128197587130049003687347913927,-0.133275967790358341900596229228,-0.137165813881786820669976271398,
                    -0.140112333814235627592549308604,-0.142304108506512913251685315535,-0.143887861666381437093965916604,
                    -0.144978935837446914022333039534,-0.145668845677819519072267311151,-0.146030804864032927305411092096,
                    -0.146123824547741905566922565449,-0.145995789590762990555343100421,-0.145685793193889359576234552428,
                    -0.145225926725857088947203430971,-0.14464266470832146126704757484,-0.143957945760789065436995586017,
                    -0.143190022968198333595318383347,-0.142354137793474443987579269298,-0.141463057812845792347563005406,
                    -0.140527508531165259220681120972,-0.139556522206686900857855156511,-0.138557721204373473901969249181,
                    -0.137537549366058216746921386872,-0.136501461857431051789130150596,-0.135454081658481212220061083972,
                    -0.134399329114304436866365922385,-0.133340529618788506480181870108,-0.132280503463840809261444216762,
                    -0.131221641077453065402514530507,-0.130165966240159592533985044519,-0.129115189370411455346007314592,
                    -0.128070752574310290494458154635,-0.127033867840764580620434593544,-0.126005549511744095941117237837,
                    -0.124986641955357871852078423135
                  },-1)
                },0)
              ),
              // AgqQ_reg as^3 LM^1 NF^1 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^3 LM^1 NF^1 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^1 NF^1 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    39.9830459958273331231679390341,99.3973264278632864200754566109,228.704248082027433776553137419,
                    422.480626863747398344230789323,864.293409743820432585925635742,1718.46441422035926008743363378,
                    3431.70952369634430313059614401,6900.69187734647292472490584126,13696.1821746331460970732440885,
                    27639.2162477623029818359182649,54740.8756793552104856265259753,110608.77281704667325391071493,
                    218905.53509712756844990067569,442490.655008915056950726028174,875592.084131815231734963999366,
                    1.76991246476667619043217628169e6,3.50261958009270016538504078672e6,7.07894087225476747508082417015e6,
                    1.40121763826183888670368909966e7,2.83119949778863880000719582506e7,5.60567191262668895555340000045e7,
                    1.13231385698686164891741396853e8,2.24260640081709219611210359759e8,4.5285767640290911891570619691e8,
                    8.97177795112767235387966724986e8,1.81116293159422191205566377425e9,3.58923890347433923759430619928e9,
                    7.24361536356437984762516633204e9,1.43589854489170889486229028338e10,2.89704938041437076392410664559e10,
                    5.74436853955799338397184091968e10,1.15866879508425362031180871142e11,2.29804144664367807141701671725e11,
                    4.63410283082059771562179605253e11,9.19327940978486268935627224345e11,1.85342451812015569202717900767e12,
                    3.67773303137890784307154557825e12,7.41287886883635804121213596617e12,1.47125251467852482366236498269e13,
                    2.96484175320433088838242559455e13,5.88561257516700528424561293436e13,1.18581950194481092589940864916e14,
                    2.35447304121782423585676294516e14,4.74283432471083638705741586289e14,9.41875570864566296865704937101e14,
                    1.89696561830264089600952552692e15,3.76782963786787884463460870661e15,7.58722487031437402449719703674e15,
                    1.50725607673541980546369131262e16,3.0346478670190453599987261867e16,6.02949619970606524449128130661e16
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^1 NF^1 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^3 LM^1 NF^1 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^1 NF^1 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    7.50617283950617283950617283951,19.3580246913580246913580246914,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309,
                    26.8641975308641975308641975309,26.8641975308641975308641975309,26.8641975308641975308641975309
                  },0),
                  // AgqQ_reg as^3 LM^1 NF^1 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    -9.48148148148148148148148148148,-2.37037037037037037037037037037,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519,
                    -11.8518518518518518518518518519,-11.8518518518518518518518518519,-11.8518518518518518518518518519
                  },0),
                  // AgqQ_reg as^3 LM^1 NF^1 [above 0.625] y=1-x H[0, y]^2
                  ome_x<double>({
                    -1.77777777777777777777777777778,-1.77777777777777777777777777778,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556,
                    -3.55555555555555555555555555556,-3.55555555555555555555555555556,-3.55555555555555555555555555556
                  },0)
                },0)
              )
            }
          )
        },0),
        // AgqQ_reg as^3 LM^2
        ome_nf<double>({
          // AgqQ_reg as^3 LM^2 NF^0
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^3 LM^2 NF^0 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^3 LM^2 NF^0 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^2 NF^0 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    -132.515734958508593854543208034,523.302959270470547090524467751,-414.651389464719115890823958124,
                    72.4938271604938271604938271605,10.2880658436213991769547325103,-0.441810699588477366255144032922,
                    1.94732510288065843621399176955,-0.48434030402284370538338792307,0.688821701520114218526916939615,
                    -0.468883849836230788611740992693,0.238794173362074596642497877066,-0.452893037045674531645350949504,
                    0.026317495947125576755206384836,-0.436509892402148294404186660079,-0.0882252404433926615448796970978,
                    -0.420255153279166303179327192351,-0.154677411555718434025312332191,-0.404560356887666896637805262693,
                    -0.194846639027296039556281343628,-0.389657800145917685765808673292,-0.21957596592808633919322577673,
                    -0.375639931654115721094313085126,-0.234755064070152146850130668858,-0.362518619632798655830763425718,
                    -0.243797949304530032774394136912,-0.350263331957599428790472116568,-0.248784753753548946331216782083,
                    -0.338823251634497246700068636517,-0.251033766847967151680524796589,-0.328139560081274669156535545217,
                    -0.251407380829107622846606284484,-0.318152060246558163785662747096,-0.250484556658896530193181908334,
                    -0.308802607976068668467925270222,-0.248662352192898828550224227244,-0.300036760425209815892687798772,
                    -0.246217904389905822294927816486,-0.291804438294152504269304795019,-0.243347459258596748027485475566,
                    -0.28406005239543118790567417702,-0.240191613410200228317919589918,-0.276762349562822013280778641877,
                    -0.23685202158094067963281571992,-0.269874122016206990643368916584,-0.233402682680464389781758812678,
                    -0.263361860998441698229686388981,-0.229897701629197660365744989729,-0.257195399179297340982584656924,
                    -0.226376713070193211412322132178,-0.251347565430777985555289255597,-0.222868725331168682336604817746,
                    -0.24579386358905801395942351096
                  },-1),
                  // AgqQ_reg as^3 LM^2 NF^0 [below 0.375] y=x H[0, y]^1
                  ome_x<double>({
                    -32.7901234567901234567901234568,52.9382716049382716049382716049,96.7901234567901234567901234568,
                    -3.16049382716049382716049382716,-2.37037037037037037037037037037,13.2740740740740740740740740741,
                    -1.3037037037037037037037037037,8.66878306878306878306878306878,-0.931216931216931216931216931217,
                    6.5467372134038800705467372134,-0.730864197530864197530864197531,5.28664421997755331088664421998,
                    -0.603367003367003367003367003367,4.44237244237244237244237244237,-0.514448514448514448514448514449,
                    3.83426943426943426943426943427,-0.448677248677248677248677248677,3.37429193899782135076252723312,
                    -0.39796659404502541757443718228,3.01372166800443374230783931506,-0.35763482781026640675763482781,
                    2.72325257588415483152325257588,-0.324771524771524771524771524772,2.48413744065917978961457222327,
                    -0.297467427902210510906163080076,2.28380032206119162640901771337,-0.274415954415954415954415954416,
                    2.1134789490345045900601456157,-0.254691810247365802921358476914,1.96687546879117760343813983661,
                    -0.237620872103630724320379492793,1.83934412721954434968895480575,-0.222700119474313022700119474313,
                    1.72738134028456609101770392093,-0.209546444840562487621311150723,1.62829414594120476473417649888,
                    -0.197862848843241000103745201784,1.53997807331140664473997807331,-0.187415485661099696187415485661,
                    1.46076440813282918546076440813,-0.178017693807167491378017693807,1.38931276492252102008199569175,
                    -0.169518647567428055232933281714,1.32453399838538806207552662459,-0.16179513853932458583621374319,
                    1.26553389189048103776785947329,-0.154745522861464890450397696775,1.21157136167928644031018375761,
                    -0.148285195463733854113132559016,1.16202698573519242212251330792,-0.142343159486016628873771730915,
                    1.11637899604286158908007647503
                  },-1),
                  // AgqQ_reg as^3 LM^2 NF^0 [below 0.375] y=x H[0, y]^2
                  ome_x<double>({
                    61.6296296296296296296296296296,-20.1481481481481481481481481481
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^2 NF^0 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^3 LM^2 NF^0 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^2 NF^0 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    70.4686061955346380342796214946,45.7443715690867519796886680881,-8.44124606244393266285056370632,
                    -143.708660795893502065330514109,-180.185113748821199274473357327,-263.975693338752718731174010651,
                    -247.77740763979561177853878812,78.8955075670443171929713072144,1008.58749795269912517487769124,
                    4365.93581054269595388182072862,11506.9690195435667024521441661,32025.7224027553784048744073524,
                    73427.7424615795843530409395126,180814.69701271124004430780473,394025.300559153751187972557151,
                    915420.990235979610185642188925,1.94459733656275822261799770594e6,4.3669793259885870029751956425e6,
                    9.1384335239703603441624677936e6,2.00726564686165974703483498936e7,4.15980690017553703586828989456e7,
                    8.99637490597104300919976181465e7,1.8518911455454536018356181797e8,3.95955329424587417682059204575e8,
                    8.11101963306685912137864148929e8,1.71914014401328829888263083436e9,3.50870259217564908360663696469e9,
                    7.38577241915912554203402181441e9,1.50313147001485135177682401947e10,3.1465964554958555232967883588e10,
                    6.38945583678084742320591716799e10,1.33147686173173681785534855204e11,2.69877081159645724110646842921e11,
                    5.60259501814416386090572480228e11,1.13389389487986665179020492278e12,2.34640976295229218831068156295e12,
                    4.74292638729487948626689250159e12,9.78785743388502335098655231783e12,1.97639815514502843461620539401e13,
                    4.0689988809243358774928024769e13,8.20892302949238428588899622356e13,1.68656577877141378119711516859e14,
                    3.39991569085848147457080957573e14,6.97265507006171101257156926089e14,1.40466812867210118584501542958e15,
                    2.876123624707773322525377427e15,5.79069166725560473785803751148e15,1.1839804491826474531827787476e16,
                    2.38256805341490260473641287476e16,4.86524386889902566109371111539e16,9.7860846247917968709591328654e16
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^2 NF^0 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^3 LM^2 NF^0 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^2 NF^0 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    39.0079934798676168609859392594,57.1808329860404563671587787656,60.8308017745500485367866933337,
                    15.2011721449204189071570637041,1.79376473751301149974965629666,-1.75388958347464282123799802432,
                    -2.96442456289851113399519967041,-3.22863217186802486541369299367,-3.01401052073843722788954753302,
                    -2.52868629502620446098305928614,-1.88138470772461715939575769884,-1.13485140397470285953757456344,
                    -0.327731949043800117187020765074,0.515106998534852201170037296724,1.37718094713708733169169610491,
                    2.24734078533925357618598292724,3.1179502493614160110468304865,3.98374644717350514698576778389,
                    4.84110253475194310348344490184,5.68753742284862545999128936925,6.52138245674580209792073811832,
                    7.34155062495128064150805553219,8.14737427863522676778996780313,8.93848951925053380860837083806,
                    9.71475290520025326611310855926,10.476180865153316369579742257,11.2229052684767378627287872846,
                    11.9551406228519894688913194369,12.6731597238156751001391387746,13.37727550354696913183714115,
                    14.0678274631852419012096586693,14.7451715182197012280280519649,15.4096724013001138836335976638,
                    16.0616979917313277867249868581,16.7016151031618668692067461148,17.3297863790688577782132400511,
                    17.946568032309718421111239014,18.5523082291076883892860273432,19.1473459655808367010761981303,
                    19.7320103207209816286057004168,20.3066199967335425178449327606,20.8714830781365657946969556784,
                    21.4268969566409953899280511919,21.973148380806197090857321929,22.5105135986784403520191671755,
                    23.0392585687390378577530515494,23.5596392200089213523864102972,24.0719017464502053430911073548,
                    24.5762829241543628833981599132,25.0730104424249512572139512564,25.5623032419137204996489998559
                  },0),
                  // AgqQ_reg as^3 LM^2 NF^0 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    -14.6172839506172839506172839506,17.3827160493827160493827160494,17.5802469135802469135802469136,
                    56.2962962962962962962962962963,70.8148148148148148148148148148,80.6913580246913580246913580247,
                    88.0197530864197530864197530864,93.776366843033509700176366843,98.4818342151675485008818342152,
                    102.441857730746619635508524397,105.849265138154027042915931805,108.832579765913099246432579766,
                    111.481289081289081289081289081,113.85994745994745994745994746,116.01650694984028317361650695,
                    117.987474254140920807587474254,119.801230867897534564201230868,121.480243213576546909880243214,
                    123.042576788982017740187674828,124.502959772584815584471587224,125.873546733815051668508840851,
                    127.164477765096960143399771883,128.384294984914179960619589103,129.540257648123219981253812635,
                    130.638583905869767582873878024,131.684638656272344072406889296,132.683081201381555848285331841,
                    133.637982541468081119995788736,134.552919456404996056910725651,135.43105064360387163381950818,
                    136.275178844053911164318808794,137.087801923195344032114301728,137.87115516891418114880088078,
                    138.627246548661474874589230225,139.357886279301205514319869955,140.064711770440422143732577799,
                    140.749208779120437359695506181,141.412729442641100880359026844,142.056507723846282670336138459,
                    142.681672699888451694961303435,143.289260044317901387568890779,143.880221988125373650325706382,
                    144.45543599423368870173018868,145.015712338968343135760686239,145.561800761025757441237007219,
                    146.094396311709163421800610522,146.614144518413891865659489163,147.121645955507270576881760875,
                    147.617460302077090473867861993,148.102109953899994576407222891,148.576083246920906644938339041
                  },0),
                  // AgqQ_reg as^3 LM^2 NF^0 [above 0.625] y=1-x H[0, y]^2
                  ome_x<double>({
                    -2.96296296296296296296296296296,-2.96296296296296296296296296296,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593,
                    -5.92592592592592592592592592593,-5.92592592592592592592592592593,-5.92592592592592592592592592593
                  },0)
                },0)
              )
            }
          )
        },0),
        // AgqQ_reg as^3 LM^3
        ome_nf<double>({
          // AgqQ_reg as^3 LM^3 NF^0
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^3 LM^3 NF^0 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^3 LM^3 NF^0 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^3 NF^0 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    -48.6584362139917695473251028807,96.7901234567901234567901234568,-64.,3.16049382716049382716049382716,
                    0.658436213991769547325102880658,0.460905349794238683127572016461,0.362139917695473251028806584362,
                    0.300999412110523221634332745444,0.258671369782480893592004703116,0.227317264354301391338428375465,
                    0.203017832647462277091906721536,0.183564035415887267739119590971,0.167601945379723157500935278713,
                    0.154249043137932026820915709805,0.142902365124587346809569031791,0.133134355356577578799801022023,
                    0.12463256907701352145796590241,0.117162914548535463568143306705,0.110546276123618171548454772856,
                    0.104643113472376171607911087329,0.0993430077250740018771207855029,0.0945573811070887094278907144404,
                    0.0902143124365346587568809791032,0.0862547722451104093616170910857,0.0826298410839473641406008555767,
                    0.0792986222937913759169797817141,0.0762266540044317822095599873378,0.0733846857303647427104217227674,
                    0.0707477250687127230337106880317,0.0682942871108047778971576332155,0.0660057978065640900889943035537,
                    0.0638661155284564010308664863108,0.0618611442984202840833665206425,0.0599785187598807670492258305878,
                    0.0582073457890451354503642085341,0.0565379911785140543310477951001,0.0549619024564558333621514449401,
                    0.0534714608788682862756936831011,0.052059857128083248940948746017,0.0507209863935010133840543197268,
                    0.0494493593908798587161160260575,0.0482400265598097576417359615192,0.0470885132131744597869259115872,
                    0.0459907638328259743776224522427,0.0449430940387012738433927064418,0.0439421490239750360336062317113,
                    0.0429848674615180251251104713263,0.0420684500583085569552147138059,0.0411903320732594039203145997267,
                    0.0403481592269164035459206009696,0.0395397665238935080204921474763,0.0387631595848215829541693220498
                  },-1),
                  // AgqQ_reg as^3 LM^3 NF^0 [below 0.375] y=x H[0, y]^1
                  ome_x<double>({
                    -5.53086419753086419753086419753,8.69135802469135802469135802469,20.5432098765432098765432098765
                  },-1),
                  // AgqQ_reg as^3 LM^3 NF^0 [below 0.375] y=x H[0, y]^2
                  ome_x<double>({
                    9.48148148148148148148148148148,-4.74074074074074074074074074074
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^3 NF^0 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^3 LM^3 NF^0 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^3 NF^0 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    -33.6660212426029158887063949052,-98.9792605811068344543696580361,-221.789529415096782272983800404,
                    -454.064838320071984182137194352,-880.258772717100264563139996493,-1762.92978736137923518999712803,
                    -3434.14621381431627764710147202,-6847.09986706516480107118515751,-13384.6705767002107917411629829,
                    -26649.9579578846928901111985698,-52234.4444854753238319554412076,-103957.626421739958105090069303,
                    -204181.375626627858488559781949,-406330.431253506559230408170846,-799380.034088410258143708077938,
                    -1.59090293164311019173931560042e6,-3.13399119476458947897088821933e6,-6.2379524505895558769993666351e6,
                    -1.23020308531677963487971479365e7,-2.4489813016661296237000018854e7,-4.83420674387314289177119876109e7,
                    -9.62500739936630812825714217697e7,-1.90145909903841445839180994772e8,-3.78642356781108101962948437555e8,
                    -7.4853934451716850064921996446e8,-1.49080425146849985486214752446e9,-2.94895381475321114342828341686e9,
                    -5.87400647099098322239893887196e9,-1.1625533646475710874399116003e10,
                    -2.31598762402224735896677254589e10,-4.58585498072237227922166477863e10,
                    -9.13686019167773576928125891021e10,-1.80994451073924822864089121118e11,
                    -3.60654893907667370054196675442e11,-7.1470392794627712186921539513e11,
                    -1.42429329011823573610244044538e12,-2.82347738384935997891123765183e12,
                    -5.6273155973455750992981236028e12,-1.11589499199476922620521525531e13,
                    -2.22423784182780787086763946089e13,-4.41193225381322767274409564062e13,
                    -8.79478299413157417327406252258e13,-1.74497067490430409178761458992e14,
                    -3.47872811288930017098670381424e14,-6.90382777944694475655939235088e14,
                    -1.37643662708473717986688869572e15,-2.73227389522597690243881507599e15,
                    -5.44781628185288969935487958234e15,-1.08163811702432412244584596123e16,
                    -2.15680232768961498087355718966e16,-4.28306937485246537305758431052e16
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^3 NF^0 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^3 LM^3 NF^0 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^3 NF^0 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    -16.,-14.024691358024691358024691358,-27.8518518518518518518518518519,-28.5102880658436213991769547325,
                    -28.5102880658436213991769547325,-28.4444444444444444444444444444,-28.3127572016460905349794238683,
                    -28.143445032333921222810111699,-27.9543797766019988242210464433,-27.7559082892416225749559082892,
                    -27.5540505584950029394473838918,-27.3523581493951864322234692605,-27.1529403381255233107084958937,
                    -26.9570407644481718555792629867,-26.7653713325670997628669586342,-26.5783114240257097399954542812,
                    -26.3960299782522004744226966449,-26.218562039912802440035555504,-26.0458576792636342382166637795,
                    -25.8778137393586692214525414082,-25.7142946808443884467276280142,-25.5551463782394049237628774707,
                    -25.4002052879460301759364925301,-25.2493045372764627500708189138,-25.1022779431567813434671543285,
                    -24.958962627050151915170386327,-24.819200672338007774441339829,-24.6828401263049903289246629052,
                    -24.5497355532960221483651874539,-24.4197482814157087140666771233,-24.2927464415329119601223888491,
                    -24.1686048674362111179513331735,-24.0472049052804901822405491461,-23.9284341660252583017183706629,
                    -23.8121862444299088138414707672,-23.6983604210188536075131800131,-23.5868613583612929539804478231,
                    -23.4775987994086023400402660918,-23.3704872730688040288430178733,-23.2654458103732045822196301448,
                    -23.1623976732944765304021256514,-23.0612700973590667305324091689,-22.9619940485534963368842812019,
                    -22.8645039945776834183952008691,-22.7687376901940131445108514342,-22.6746359762188496999649022365,
                    -22.5821425915747704609221993133,-22.4912039977460949450511225532,-22.4017692149420451312161737179,
                    -22.3137896692600021539782028986,-22.2272190501477179773942849669
                  },0),
                  // AgqQ_reg as^3 LM^3 NF^0 [above 0.625] y=1-x H[0, y]^1
                  ome_x<double>({
                    -1.97530864197530864197530864198,-1.97530864197530864197530864198,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395,
                    -3.95061728395061728395061728395,-3.95061728395061728395061728395,-3.95061728395061728395061728395
                  },0)
                },0)
              )
            }
          ),
          // AgqQ_reg as^3 LM^3 NF^1
          ome_piecewisex<double>(
            {0.375,0.625},
            {
              // AgqQ_reg as^3 LM^3 NF^1 [below 0.375]
              make_ome_genps(
                std::function<double(double)>(f_id<double>),
                // AgqQ_reg as^3 LM^3 NF^1 [below 0.375] y=x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^3 NF^1 [below 0.375] y=x H[0, y]^0
                  ome_x<double>({
                    2.37037037037037037037037037037,-2.37037037037037037037037037037,1.18518518518518518518518518519
                  },-1)
                },0)
              ),
              // AgqQ_reg as^3 LM^3 NF^1 [0.375,0.625]
              make_ome_genps(
                std::function<double(double)>(f_half<double>),
                // AgqQ_reg as^3 LM^3 NF^1 [0.375,0.625] y=1/2-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^3 NF^1 [0.375,0.625] y=1/2-x H[0, y]^0
                  ome_x<double>({
                    2.96296296296296296296296296296,8.2962962962962962962962962963,18.962962962962962962962962963,
                    37.9259259259259259259259259259,75.8518518518518518518518518519,151.703703703703703703703703704,
                    303.407407407407407407407407407,606.814814814814814814814814815,1213.62962962962962962962962963,
                    2427.25925925925925925925925926,4854.51851851851851851851851852,9709.03703703703703703703703704,
                    19418.0740740740740740740740741,38836.1481481481481481481481481,77672.2962962962962962962962963,
                    155344.592592592592592592592593,310689.185185185185185185185185,621378.37037037037037037037037,
                    1.24275674074074074074074074074e6,2.48551348148148148148148148148e6,4.97102696296296296296296296296e6,
                    9.94205392592592592592592592593e6,1.98841078518518518518518518519e7,3.97682157037037037037037037037e7,
                    7.95364314074074074074074074074e7,1.59072862814814814814814814815e8,3.1814572562962962962962962963e8,
                    6.36291451259259259259259259259e8,1.27258290251851851851851851852e9,2.54516580503703703703703703704e9,
                    5.09033161007407407407407407407e9,1.01806632201481481481481481481e10,2.03613264402962962962962962963e10,
                    4.07226528805925925925925925926e10,8.14453057611851851851851851852e10,1.6289061152237037037037037037e11,
                    3.25781223044740740740740740741e11,6.51562446089481481481481481481e11,1.30312489217896296296296296296e12,
                    2.60624978435792592592592592593e12,5.21249956871585185185185185185e12,1.04249991374317037037037037037e13,
                    2.08499982748634074074074074074e13,4.16999965497268148148148148148e13,8.33999930994536296296296296296e13,
                    1.66799986198907259259259259259e14,3.33599972397814518518518518519e14,6.67199944795629037037037037037e14,
                    1.33439988959125807407407407407e15,2.66879977918251614814814814815e15,5.3375995583650322962962962963e15
                  },0)
                },0)
              ),
              // AgqQ_reg as^3 LM^3 NF^1 [above 0.625]
              make_ome_genps(
                std::function<double(double)>(f_omx<double>),
                // AgqQ_reg as^3 LM^3 NF^1 [above 0.625] y=1-x
                ome_logx<double>({
                  // AgqQ_reg as^3 LM^3 NF^1 [above 0.625] y=1-x H[0, y]^0
                  ome_x<double>({
                    1.18518518518518518518518518519,1.18518518518518518518518518519,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037,
                    2.37037037037037037037037037037,2.37037037037037037037037037037,2.37037037037037037037037037037
                  },0)
                },0)
              )
            }
          )
        },0)
      },0)
    },2);


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


/* C interface */

double ome_AgqQ_reg(double as, double LM, double NF, double x)
{
  return(ome::AgqQ_reg(as, LM, NF, x));
}

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

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

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

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

int ome_AgqQ_reg_min_power()
{
  return(ome::AgqQ_reg.min_power());
}

int ome_AgqQ_reg_max_power()
{
  return(ome::AgqQ_reg.max_power());
}

int ome_AgqQ_reg_coeff_as_min_power(int order_as)
{
  return(ome::AgqQ_reg[order_as].min_power());
}

int ome_AgqQ_reg_coeff_as_max_power(int order_as)
{
  return(ome::AgqQ_reg[order_as].max_power());
}

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

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