cubicsuperancillary.hpp

Go to the documentation of this file.

#pragma once 
#include <vector>
 
namespace teqp {
 
namespace CubicSuperAncillary {
 
struct Chebyshev{
public:
    const std::vector<double> coeff;
    const double xmin, xmax;
    // Evaluate the expansion with Clenshaw's method
    double y(double x) const{
        // Scale to (-1, 1)
        double xscaled = (2*x - (xmax + xmin)) / (xmax - xmin);
        int Norder = static_cast<int>(coeff.size()) - 1;
        double u_k = 0, u_kp1 = coeff[Norder], u_kp2 = 0;
        for (int k = Norder-1; k > 0; k--){ // k must be signed!
            // Do the recurrent calculation
            u_k = 2.0*xscaled*u_kp1 - u_kp2 + coeff[k];
            // Update the values
            u_kp2 = u_kp1; u_kp1 = u_k;
        }
        return coeff[0] + xscaled*u_kp1 - u_kp2;
    };
};
 
// https://proquest.safaribooksonline.com/9780321637413
// https://web.stanford.edu/class/archive/cs/cs107/cs107.1202/lab1/
static int midpoint_Knuth(int x, int y) {
    return (x & y) + ((x ^ y) >> 1);
};
 
struct SuperAncillary{
public:
 
    const std::vector<Chebyshev> exps;
 
    int get_index(double x) const{
        int iL = 0, iR = static_cast<int>(exps.size()) - 1, iM;
        while (iR - iL > 1) {
            iM = midpoint_Knuth(iL, iR);
            if (x >= exps[iM].xmin) {
                iL = iM;
            }
            else {
                iR = iM;
            }
        }
        return (x < exps[iL].xmax) ? iL : iR;
    };
 
    double y(double x) const{
        // First check whether the input is possible
        if (x < exps[0].xmin) {
            throw std::invalid_argument("Ttilde (" + std::to_string(x) + ") is below the minimum of " + std::to_string(exps[0].xmin));
        }
        if (x > exps.back().xmax) {
            throw std::invalid_argument("Ttilde (" + std::to_string(x) + ") is above the maximum of " + std::to_string(exps.back().xmax));
        }
        // Bisection to find the expansion
        // we need
        auto i = get_index(x);
        // Evaluate the expansion
        return exps[i].y(x);
    }
};
 
const auto vdW_p = SuperAncillary{
{
  {
    {
      8.369016487003571e-13,
      1.3838182749488314e-12,
      8.006631877160537e-13,
      3.3527912014616007e-13,
      1.043509859614108e-13,
      2.4518354339464807e-14,
      4.362454910953413e-15,
      5.800290149785568e-16,
      5.5246311563342216e-17,
      3.341712352615597e-18,
      6.943507182585673e-20,
      -6.6172300100790814e-21,
      -4.799249815587726e-22,
      8.197853768324654e-24,
      1.6972506557505968e-24,
      -1.4268299756055458e-26,
      -5.664810160392086e-27,
      4.437342591868191e-31,
      -1.4685138788754898e-28
    },
    0.02962962962962963,
    0.03796296296296296
  },
  {
    {
      1.1901434626227093e-10,
      1.7361780173369773e-10,
      7.639416392753658e-11,
      2.219331080562664e-11,
      4.468985012784936e-12,
      6.345623163747038e-13,
      6.256832040916454e-14,
      3.9794258615342e-15,
      1.2126532982344795e-16,
      -2.510288446235168e-18,
      -3.1189177665897616e-19,
      -3.0550813286427195e-22,
      6.384476478410435e-22,
      2.6783405454063793e-25,
      -1.4724057241513977e-24,
      -5.563678192342752e-26,
      1.0097419586828951e-27,
      6.679443056687351e-26,
      5.4526065768876336e-27
    },
    0.03796296296296296,
    0.0462962962962963
  },
  {
    {
      3.341667610526653e-08,
      5.118389108542411e-08,
      2.4371256768154856e-08,
      7.57899054003236e-09,
      1.5635438033324503e-09,
      2.0715013642702e-10,
      1.5165995758982126e-11,
      1.6314454982752663e-13,
      -6.017223931292203e-14,
      -1.9704412673913643e-15,
      2.8051398344015177e-16,
      7.022929896522875e-18,
      -1.5199614910962848e-18,
      1.1958046052493347e-20,
      7.00911492046808e-21,
      -4.237329623030918e-22,
      -3.863191954564062e-23,
      -1.1166938269465874e-23,
      -1.2818068320304144e-23
    },
    0.0462962962962963,
    0.06296296296296297
  },
  {
    {
      8.410724206747044e-06,
      1.2481631400263709e-05,
      5.378159443205515e-06,
      1.3704527561658967e-06,
      1.9136168235190198e-07,
      8.957711663130537e-09,
      -9.462672897450059e-10,
      -6.968811539371602e-11,
      1.0243849900317855e-11,
      2.337729169280951e-13,
      -1.1494530386099005e-13,
      5.3160592525426036e-15,
      7.530730564512324e-16,
      -1.1431990843247566e-16,
      2.805445086019535e-18,
      7.988918627843267e-19,
      -1.0295338058018748e-19,
      7.210533399624742e-21,
      5.761312966431838e-22
    },
    0.06296296296296297,
    0.0962962962962963
  },
  {
    {
      0.0006802593682753608,
      0.0009003332567452364,
      0.0002880297751978097,
      4.118964226458961e-05,
      7.986981798077676e-07,
      -2.628431475282669e-07,
      1.4516582903911464e-08,
      1.5647152494426528e-09,
      -4.1929621787396065e-10,
      3.8557343395626107e-11,
      8.368753107866514e-13,
      -8.31125296022126e-13,
      1.3326079144488354e-13,
      -8.864161253017153e-15,
      -8.83270499045497e-16,
      3.3879136780332834e-16,
      -4.8188763160481214e-17,
      2.9076523496871995e-18,
      3.8719993601362204e-19
    },
    0.0962962962962963,
    0.16296296296296295
  },
  {
    {
      0.015871052250174933,
      0.017373046113870208,
      0.003617191914161977,
      0.00018958352208271672,
      -1.4572937633430827e-05,
      7.218425887058066e-07,
      2.6419947263015497e-08,
      -1.4409436469961143e-08,
      2.6380419456436925e-09,
      -3.4860393171849517e-10,
      3.335912880024321e-11,
      -1.2239621615490678e-12,
      -4.035014645542914e-13,
      1.3487363894996995e-13,
      -2.72259564068969e-14,
      4.280375966864147e-15,
      -5.359753439682091e-16,
      5.147249813874932e-17,
      -1.531435568635775e-18
    },
    0.16296296296296295,
    0.2962962962962963
  }
}
};
const auto vdW_rhoL = SuperAncillary{
{
  {
    {
      0.9316136201036852,
      -0.03865225088424542,
      -0.0008661148627486168,
      -3.8336595966463605e-05,
      -1.930441992932519e-06,
      -5.810841478132156e-08,
      6.379356029539984e-09,
      4.0565692471794534e-10,
      -3.1201431022198634e-10,
      -3.4029466744467385e-11,
      1.369238750159596e-11,
      1.18192261533423e-12,
      -7.751924102628038e-13,
      6.175615574477433e-15,
      4.091171845743702e-14,
      -6.980527267330672e-15,
      -1.2559397966072083e-15,
      1.4641066137244252e-15,
      -3.8163916471489756e-17
    },
    0.02962962962962963,
    0.0962962962962963
  },
  {
    {
      0.8455034879915954,
      -0.04802256898834967,
      -0.0015470962415359973,
      -7.972501988821878e-05,
      -4.017234579128948e-06,
      -2.921741757144769e-07,
      -2.6298690877779585e-08,
      -1.3988705679923719e-09,
      -1.5143918064008943e-10,
      -1.6808554548219945e-11,
      -1.6921186674068167e-13,
      -1.6256440638073855e-13,
      -9.825473767932635e-15,
      7.077671781985373e-16,
      -3.608224830031759e-16,
      -7.632783294297951e-16,
      -4.163336342344337e-17,
      8.257283745649602e-16,
      1.3877787807814457e-17
    },
    0.0962962962962963,
    0.16296296296296295
  },
  {
    {
      0.7327617276279949,
      -0.06612960981729517,
      -0.0032780016312522198,
      -0.00026019391083406673,
      -2.6595986463442967e-05,
      -3.184060496402852e-06,
      -4.0508839448810674e-07,
      -5.4371095223326726e-08,
      -7.554792333386295e-09,
      -1.076694275403689e-09,
      -1.5675695463901462e-10,
      -2.3183531483450537e-11,
      -3.4755739952707643e-12,
      -5.273004255457181e-13,
      -8.064382495120981e-14,
      -1.3183898417423734e-14,
      -1.9290125052862095e-15,
      4.0245584642661925e-16,
      -5.898059818321144e-17
    },
    0.16296296296296295,
    0.22962962962962963
  },
  {
    {
      0.6166811517612045,
      -0.04839066994648989,
      -0.0021672975061474767,
      -0.00017570976987404258,
      -1.8519040533800102e-05,
      -2.208561613849247e-06,
      -2.82949510460595e-07,
      -3.8035669329417043e-08,
      -5.291587154565569e-09,
      -7.554387740360546e-10,
      -1.1004187144836308e-10,
      -1.6289976312311438e-11,
      -2.4432955658681976e-12,
      -3.709324514211687e-13,
      -5.671157987663378e-14,
      -9.346690088563037e-15,
      -1.3600232051658168e-15,
      3.677613769070831e-16,
      -3.122502256758253e-17
    },
    0.22962962962962963,
    0.26296296296296295
  },
  {
    {
      0.5328306974405083,
      -0.034453741516764766,
      -0.0014682721939933072,
      -0.00012182552205808328,
      -1.2945521315733954e-05,
      -1.5482578362285837e-06,
      -1.9873352891802698e-07,
      -2.6746483347550587e-08,
      -3.724136381222376e-09,
      -5.319946913551199e-10,
      -7.753041264546567e-11,
      -1.1481385286948864e-11,
      -1.7225665338571616e-12,
      -2.6155466681387907e-13,
      -4.00304789316408e-14,
      -6.696032617270475e-15,
      -9.43689570931383e-16,
      3.885780586188048e-16,
      -1.734723475976807e-17
    },
    0.26296296296296295,
    0.2796296296296296
  },
  {
    {
      0.4734609830190833,
      -0.02424735785371843,
      -0.0010133644361935874,
      -8.529016931251368e-05,
      -9.098505499984672e-06,
      -1.089985510512731e-06,
      -1.4004471942871688e-07,
      -1.8859509249430362e-08,
      -2.627095759683007e-09,
      -3.7540060296437083e-10,
      -5.472236552783727e-11,
      -8.105231763533283e-12,
      -1.2161383011743965e-12,
      -1.8476886687324168e-13,
      -2.827946210537391e-14,
      -4.829470157119431e-15,
      -7.181755190543981e-16,
      4.2674197509029455e-16,
      6.938893903907228e-18
    },
    0.2796296296296296,
    0.287962962962963
  },
  {
    {
      0.4317586705206012,
      -0.016998875614382822,
      -0.0007071969216258488,
      -6.0002743216675675e-05,
      -6.413721659864063e-06,
      -7.690254281537967e-07,
      -9.885537948645107e-08,
      -1.3316815754282896e-08,
      -1.855412701950998e-09,
      -2.6517297732620015e-10,
      -3.8659075940472576e-11,
      -5.726488727653134e-12,
      -8.592571099086399e-13,
      -1.30593452718486e-13,
      -1.9970136655445003e-14,
      -3.4555691641457997e-15,
      -4.891920202254596e-16,
      3.5041414214731503e-16,
      -6.938893903907228e-18
    },
    0.287962962962963,
    0.29212962962962963
  },
  {
    {
      0.4025301670915846,
      -0.011913454752168377,
      -0.0004965958130950694,
      -4.231862747736473e-05,
      -4.528099139627234e-06,
      -5.431755606445654e-07,
      -6.984058851383645e-08,
      -9.409722008524302e-09,
      -1.311186739044734e-09,
      -1.8740799190286594e-10,
      -2.732347334499785e-11,
      -4.0475123253003176e-12,
      -6.073752611968075e-13,
      -9.242259735309233e-14,
      -1.4089424071883627e-14,
      -2.5500435096859064e-15,
      -3.0878077872387166e-16,
      3.9898639947466563e-16,
      -3.122502256758253e-17
    },
    0.29212962962962963,
    0.29421296296296295
  },
  {
    {
      0.3820369525637211,
      -0.008358544107590068,
      -0.00034987653550601003,
      -2.988469717726469e-05,
      -3.1993345159762876e-06,
      -3.8386783924454493e-07,
      -4.936325082721682e-08,
      -6.651310394817367e-09,
      -9.268700876252645e-10,
      -1.3248298638690592e-10,
      -1.93161493744487e-11,
      -2.861148817867587e-12,
      -4.2926079357741287e-13,
      -6.528458329491116e-14,
      -9.988537774674455e-15,
      -2.067790383364354e-15,
      -2.740863092043355e-16,
      2.3592239273284576e-16,
      -1.734723475976807e-17
    },
    0.29421296296296295,
    0.2952546296296296
  },
  {
    {
      0.3676497053323325,
      -0.005873143997041752,
      -0.0002469394057881408,
      -2.11177904382992e-05,
      -2.261380224737042e-06,
      -2.7135935164526725e-07,
      -3.489748082521893e-08,
      -4.70234896160493e-09,
      -6.552974517182175e-10,
      -9.366743838890024e-11,
      -1.3656992203792129e-11,
      -2.0231802344561345e-12,
      -3.0342742207700724e-13,
      -4.5678738569421284e-14,
      -7.029099524658022e-15,
      -1.4190038033490282e-15,
      -1.9081958235744878e-16,
      2.706168622523819e-16,
      0.0
    },
    0.2952546296296296,
    0.295775462962963
  },
  {
    {
      0.3575341743707003,
      -0.004132704229989258,
      -0.0001744465560134939,
      -1.4927613207975365e-05,
      -1.5987220087325393e-06,
      -1.918530817199282e-07,
      -2.4673554255283392e-08,
      -3.324766725820716e-09,
      -4.633303775369857e-10,
      -6.62290049191494e-11,
      -9.656359045706608e-12,
      -1.4302690976020216e-12,
      -2.1446733278196461e-13,
      -3.235953172087136e-14,
      -4.947431353485854e-15,
      -8.743006318923108e-16,
      -1.249000902703301e-16,
      2.706168622523819e-16,
      -1.3877787807814457e-17
    },
    0.295775462962963,
    0.29603587962962963
  },
  {
    {
      0.35041246164734846,
      -0.0029115644765456902,
      -0.0001232928553679749,
      -1.0553674979200955e-05,
      -1.130355650916931e-06,
      -1.3565108163485218e-07,
      -1.744588612809239e-08,
      -2.3508600899280196e-09,
      -3.2761145879467435e-10,
      -4.682861390326032e-11,
      -6.827677312415403e-12,
      -1.011458278243893e-12,
      -1.5187850976872141e-13,
      -2.3189783426857957e-14,
      -3.4139358007223564e-15,
      -1.1136924715771102e-15,
      2.7755575615628914e-17,
      5.134781488891349e-16,
      -1.231653667943533e-16
    },
    0.29603587962962963,
    0.29616608796296295
  },
  {
    {
      0.3453929468556045,
      -0.002053229816890382,
      -8.715997446502935e-05,
      -7.461959008286251e-06,
      -7.992427072806008e-07,
      -9.591642778292164e-08,
      -1.2335767615723192e-08,
      -1.662271872232557e-09,
      -2.3165206211595013e-10,
      -3.311351193246992e-11,
      -4.8279401310136194e-12,
      -7.147511749128199e-13,
      -1.0711917464156784e-13,
      -1.676089822488791e-14,
      -2.4980018054066022e-15,
      -2.983724378680108e-16,
      -5.551115123125783e-17,
      2.393918396847994e-16,
      -3.469446951953614e-18
    },
    0.29616608796296295,
    0.2962311921296296
  },
  {
    {
      0.34185203830414934,
      -0.0014490000408041367,
      -6.162384665049647e-05,
      -5.276183904001641e-06,
      -5.651359926590882e-07,
      -6.782196396853957e-08,
      -8.722585896625557e-09,
      -1.1753904204103716e-09,
      -1.6380120870795878e-10,
      -2.3415179517538576e-11,
      -3.413873350677221e-12,
      -5.057065877167588e-13,
      -7.568598525686809e-14,
      -1.1497747198774277e-14,
      -1.7867651802561113e-15,
      -9.783840404509192e-16,
      -3.8163916471489756e-17,
      8.673617379884035e-16,
      4.5102810375396984e-17
    },
    0.2962311921296296,
    0.296263744212963
  },
  {
    {
      0.33935254842481155,
      -0.001023145052791774,
      -4.35719527140041e-05,
      -3.7307483613316372e-06,
      -3.99606561996696e-07,
      -4.7956948610178296e-08,
      -6.167757610198166e-09,
      -8.311229968138711e-10,
      -1.1582452358327444e-10,
      -1.6557800269767498e-11,
      -2.4139960863589494e-12,
      -3.5761671401957074e-13,
      -5.354744425645208e-14,
      -9.089951014118469e-15,
      -1.2524703496552547e-15,
      3.0531133177191805e-16,
      1.0408340855860843e-16,
      1.1796119636642288e-15,
      7.45931094670027e-17
    },
    0.296263744212963,
    0.29628002025462963
  },
  {
    {
      0.3375873324362272,
      -0.0007227370941488807,
      -3.080906979128009e-05,
      -2.638010219744441e-06,
      -2.8256276612612097e-07,
      -3.391053622625595e-08,
      -4.361249632267583e-09,
      -5.876913226898761e-10,
      -8.18994629059322e-11,
      -1.1706968727764888e-11,
      -1.7064231971897215e-12,
      -2.5467475350815505e-13,
      -3.8767600241129685e-14,
      -3.230055112268815e-15,
      -1.0408340855860843e-16,
      -3.2959746043559335e-16,
      4.718447854656915e-16,
      -2.255140518769849e-16,
      -5.828670879282072e-16
    },
    0.29628002025462963,
    0.29628815827546295
  },
  {
    {
      0.33634024377980126,
      -0.0005106811254640582,
      -2.1784964224096787e-05,
      -1.8653452834756223e-06,
      -1.9980143188894073e-07,
      -2.3978314021844138e-08,
      -3.083862751968036e-09,
      -4.15558362026891e-10,
      -5.791153667722426e-11,
      -8.274964047316757e-12,
      -1.206382216345503e-12,
      -1.7613341340982913e-13,
      -2.653433028854124e-14,
      -5.748873599387139e-15,
      -8.881784197001252e-16,
      -1.9359513991901167e-15,
      3.8163916471489756e-17,
      2.723515857283587e-15,
      2.0122792321330962e-16
    },
    0.29628815827546295,
    0.2962922272858796
  },
  {
    {
      0.3354589771042352,
      -0.0003609193232701237,
      -1.540417623366666e-05,
      -1.318994892206965e-06,
      -1.412807294452756e-07,
      -1.695521064867811e-08,
      -2.180618217245689e-09,
      -2.9384074340388366e-10,
      -4.09495146125316e-11,
      -5.853536405586723e-12,
      -8.535151752031567e-13,
      -1.2840770113875521e-13,
      -1.8620521791135047e-14,
      8.777700788442644e-16,
      -1.1449174941446927e-16,
      -3.2612801348363973e-16,
      -7.28583859910259e-17,
      -3.642919299551295e-16,
      -6.418476861114186e-17
    },
    0.2962922272858796,
    0.296294261791088
  },
  {
    {
      0.33483610810587444,
      -0.00025511470266906255,
      -1.0892355099323403e-05,
      -9.326690283625871e-07,
      -9.990048497124371e-08,
      -1.198913999328477e-08,
      -1.5419303056596334e-09,
      -2.0778637838114733e-10,
      -2.8956462228002522e-11,
      -4.139157766536172e-12,
      -6.029898802495381e-13,
      -8.481410018745805e-14,
      -1.3000017728970192e-14,
      -4.253541963095131e-15,
      -4.961309141293668e-16,
      -3.400058012914542e-16,
      2.42861286636753e-16,
      3.427813588530171e-15,
      2.5847379792054426e-16
    },
    0.296294261791088,
    0.29629527904369213
  },
  {
    {
      0.3341194565621786,
      -0.0005254907183565548,
      -0.00010585696826935151,
      -4.5948894997340645e-05,
      -2.5987955723658085e-05,
      -1.6924206122554547e-05,
      -1.2054006867302791e-05,
      -9.1430248983157e-06,
      -7.271484405917822e-06,
      -6.0039811285959865e-06,
      -5.112767081889152e-06,
      -4.469487214352197e-06,
      -3.9975640410763735e-06,
      -3.6492975527165072e-06,
      -3.3940536075909022e-06,
      -3.2118298490882435e-06,
      -3.0896085961677755e-06,
      -3.0192461827061035e-06,
      -1.498133949695904e-06
    },
    0.29629527904369213,
    0.2962962962962963
  }
}
};
const auto vdW_rhoV = SuperAncillary{
{
  {
    {
      2.2489373895047762e-11,
      3.688517107894903e-11,
      2.0894583452709746e-11,
      8.476583933451082e-12,
      2.530366544413405e-12,
      5.637212489418141e-13,
      9.363230492731065e-14,
      1.1333103251162971e-14,
      9.33519994749376e-16,
      4.1362278059373704e-17,
      -4.84986420418459e-19,
      -1.65650780002939e-19,
      -4.018942009006866e-21,
      4.891126970541962e-22,
      2.0177374582077214e-23,
      -1.6609931787362484e-24,
      -6.522933053091502e-26,
      3.332148463653554e-27,
      -4.030290364774149e-27
    },
    0.02962962962962963,
    0.03796296296296296
  },
  {
    {
      2.6379937128657115e-09,
      3.780798221053315e-09,
      1.6044278594036214e-09,
      4.4331172895958845e-10,
      8.361702765680484e-11,
      1.0880537367567157e-11,
      9.447363806147374e-13,
      4.7714402398820264e-14,
      5.233241126040141e-16,
      -8.52742988949115e-17,
      -3.190159512313497e-18,
      1.504383250318862e-19,
      7.755246373275116e-21,
      -3.752718106348484e-22,
      -1.6549266806029178e-23,
      -6.308867757850729e-25,
      -7.270142102516845e-26,
      1.3231658626580658e-24,
      1.0339757656912846e-25
    },
    0.03796296296296296,
    0.0462962962962963
  },
  {
    {
      5.488117599580858e-07,
      8.246318599119535e-07,
      3.7495338087883065e-07,
      1.0860523684389838e-07,
      2.0169156352224874e-08,
      2.2464415506045487e-09,
      1.0692430180323234e-10,
      -5.062128204483255e-12,
      -7.086734785957356e-13,
      1.8598468871247458e-14,
      3.791438263524812e-15,
      -1.5845682631384665e-16,
      -1.6139528628166535e-17,
      1.4249482656836079e-18,
      2.1099309474696353e-20,
      -9.177568896275842e-21,
      2.61595868719895e-23,
      -1.4175807747627512e-22,
      -2.0384832220603676e-22
    },
    0.0462962962962963,
    0.06296296296296297
  },
  {
    {
      9.198583413591418e-05,
      0.0001323418019880509,
      5.269589299645468e-05,
      1.1673085614775962e-05,
      1.217520790883859e-06,
      -5.281166426031512e-09,
      -9.943118098154909e-09,
      2.878024017224007e-10,
      1.005351237026266e-10,
      -8.653088852547337e-12,
      -5.972466878192423e-13,
      1.4256395010153724e-13,
      -4.811359136869256e-15,
      -1.1177042446084331e-15,
      1.499644436891851e-16,
      -3.257172554437806e-18,
      -1.0356698315857932e-18,
      1.720668023008306e-19,
      -4.5858893159939854e-21
    },
    0.06296296296296297,
    0.0962962962962963
  },
  {
    {
      0.004754326043670642,
      0.005941275521721304,
      0.001655883437830305,
      0.00017737301927121552,
      -2.1044306045037264e-06,
      -2.446112805967309e-07,
      1.8602953134177618e-07,
      -1.3107585751750864e-08,
      -4.852484754108938e-10,
      4.195957567265955e-10,
      -6.078596600608724e-11,
      3.1842478720250315e-12,
      7.401045989713535e-13,
      -2.1153502669131463e-13,
      2.655919239960905e-14,
      -7.025291264527167e-16,
      -4.339586958009012e-16,
      1.0583507269353032e-16,
      -1.1385816876992305e-17
    },
    0.0962962962962963,
    0.16296296296296295
  },
  {
    {
      0.036691881497644734,
      0.027888941391814935,
      0.003974549212687695,
      0.0002715266512105239,
      2.3797626951819117e-05,
      3.422813454436401e-06,
      3.9210812700050945e-07,
      5.451342389246796e-08,
      7.630582787377174e-09,
      1.0642235220240798e-09,
      1.580515761670509e-10,
      2.3087444716435312e-11,
      3.480158435736902e-12,
      5.272167251380022e-13,
      8.05698823630463e-14,
      1.2436666280146724e-14,
      1.932265111803666e-15,
      3.380542373809803e-16,
      4.8030156241107846e-17
    },
    0.16296296296296295,
    0.22962962962962963
  },
  {
    {
      0.09807361548934067,
      0.03140331546583538,
      0.0023425610291630955,
      0.0001748056758866679,
      1.8480646404825905e-05,
      2.210790222235623e-06,
      2.8286844096486047e-07,
      3.8038013492475153e-08,
      5.2915361961961005e-09,
      7.554389618198709e-10,
      1.1004225222016606e-10,
      1.6290535326951572e-11,
      2.443793865186672e-12,
      3.706427526006806e-13,
      5.67323965583455e-14,
      8.677086826835989e-15,
      1.353951672999898e-15,
      3.0357660829594124e-16,
      3.209238430557093e-17
    },
    0.22962962962962963,
    0.26296296296296295
  },
  {
    {
      0.15708895435350148,
      0.026473214602337546,
      0.0015096794373112487,
      0.00012167462420701724,
      1.2944662426067101e-05,
      1.5482932451916717e-06,
      1.987328542198824e-07,
      2.6746493632726076e-08,
      3.7241364748974437e-09,
      5.31994449361195e-10,
      7.753074050820263e-11,
      1.1481976827654172e-11,
      1.7230132251522257e-12,
      2.6130746871855237e-13,
      4.001920322904695e-14,
      6.041174505089231e-15,
      9.497611030973019e-16,
      2.697495005143935e-16,
      2.2985086056692694e-17
    },
    0.26296296296296295,
    0.2796296296296296
  },
  {
    {
      0.20464143088991663,
      0.02037941444367245,
      0.0010233677370924937,
      8.527036209769225e-05,
      9.098480032524933e-06,
      1.0899862908780866e-06,
      1.4004471174042243e-07,
      1.8859509351779047e-08,
      2.627095993870676e-09,
      3.754003670419781e-10,
      5.472267257389252e-11,
      8.105812895897735e-12,
      1.2166136154068141e-12,
      1.8451386252227309e-13,
      2.827425793494598e-14,
      4.210173876195711e-15,
      6.886852199627924e-16,
      2.3418766925686896e-16,
      -3.469446951953614e-18
    },
    0.2796296296296296,
    0.287962962962963
  },
  {
    {
      0.2405791169899757,
      0.015094667641419904,
      0.0007096529008297752,
      6.0000237077412547e-05,
      6.413720740528275e-06,
      7.690254484101627e-07,
      9.885537929910093e-08,
      1.331681581152877e-08,
      1.8554129257303265e-09,
      2.651727292607431e-10,
      3.8659369108740016e-11,
      5.727075064188014e-12,
      8.597150769062978e-13,
      1.303020191745219e-13,
      1.9978810272824887e-14,
      2.831068712794149e-15,
      4.666406150377611e-16,
      2.8622937353617317e-16,
      1.1275702593849246e-17
    },
    0.287962962962963,
    0.29212962962962963
  },
  {
    {
      0.26696049015665074,
      0.010968687650483065,
      0.0004972041585612832,
      4.231831304510672e-05,
      4.5280991004606475e-06,
      5.431755610851852e-07,
      6.98405883993447e-08,
      9.409722064035453e-09,
      1.3111869472115512e-09,
      1.874077785318784e-10,
      2.732374916603053e-11,
      4.048053559024822e-12,
      6.078245545770855e-13,
      9.212769436217627e-14,
      1.4092893518835581e-14,
      1.93421667571414e-15,
      3.0184188481996443e-16,
      2.393918396847994e-16,
      3.642919299551295e-17
    },
    0.29212962962962963,
    0.29421296296296295
  },
  {
    {
      0.2860388317887753,
      0.007887981700325456,
      0.00035002791393559315,
      2.9884657820819505e-05,
      3.199334514061153e-06,
      3.8386783914393097e-07,
      4.936325074221537e-08,
      6.651310417368772e-09,
      9.268703131393163e-10,
      1.3248275219923666e-10,
      1.9316418256587475e-11,
      2.8617264807850873e-12,
      4.2971702585159477e-13,
      6.502437477351464e-14,
      9.985068327722502e-15,
      1.4536982728685643e-15,
      2.463307335887066e-16,
      3.937822290467352e-16,
      1.3877787807814457e-17
    },
    0.29421296296296295,
    0.2952546296296296
  },
  {
    {
      0.2997208026355501,
      0.005638316441575676,
      0.00024697716181993484,
      2.1117785516101134e-05,
      2.2613802245774473e-06,
      2.713593515168977e-07,
      3.4897480762768884e-08,
      4.702349017116081e-09,
      6.552976390683529e-10,
      9.366720593595446e-11,
      1.3657262820654381e-11,
      2.023738815415399e-12,
      3.0388191962771316e-13,
      4.5428938388880624e-14,
      7.004813395994347e-15,
      8.222589276130066e-16,
      1.6653345369377348e-16,
      3.5041414214731503e-16,
      1.214306433183765e-17
    },
    0.2952546296296296,
    0.295775462962963
  },
  {
    {
      0.3094842337937238,
      0.004015403659434553,
      0.00017445598394503006,
      1.4927612592627315e-05,
      1.598722008659681e-06,
      1.918530815638031e-07,
      2.4673554175486112e-08,
      3.3247667709235262e-09,
      4.633306065204845e-10,
      6.622874818007496e-11,
      9.656650479250573e-12,
      1.4308033924326224e-12,
      2.1489754420400686e-13,
      3.2085445411667024e-14,
      4.933553565678039e-15,
      3.157196726277789e-16,
      1.0408340855860843e-16,
      3.5735303605122226e-16,
      -3.469446951953614e-18
    },
    0.295775462962963,
    0.29603587962962963
  },
  {
    {
      0.3164300310090633,
      0.002852942467447068,
      0.00012329521096597976,
      1.0553674902311072e-05,
      1.1303556508510115e-06,
      1.356510815030132e-07,
      1.744588605176456e-08,
      2.35086013503083e-09,
      3.276116843087262e-10,
      4.6828395328102346e-11,
      6.8279583376185116e-12,
      1.0120272675440134e-12,
      1.523364767663793e-13,
      2.2915697117653622e-14,
      3.4208746946262636e-15,
      5.169475958410885e-16,
      -5.898059818321144e-17,
      1.0755285551056204e-16,
      1.1275702593849246e-16
    },
    0.29603587962962963,
    0.29616608796296295
  },
  {
    {
      0.3213616216197866,
      0.002023925878183369,
      8.716056319137827e-05,
      7.461958998724455e-06,
      7.992427072528452e-07,
      9.5916427661491e-08,
      1.2335767528987018e-08,
      1.662271934682602e-09,
      2.316522806911081e-10,
      3.3113272540630234e-11,
      4.828210747875872e-12,
      7.152750614025649e-13,
      1.0757367219227376e-13,
      1.6459056340067946e-14,
      2.4875934645507414e-15,
      -2.8796409701215e-16,
      1.734723475976807e-17,
      4.0245584642661925e-16,
      8.673617379884035e-18
    },
    0.29616608796296295,
    0.2962311921296296
  },
  {
    {
      0.32485857647079686,
      0.0014343498375109544,
      6.162399381036757e-05,
      5.27618390287754e-06,
      5.651359926105159e-07,
      6.782196385404782e-08,
      8.722585830706064e-09,
      1.1753904828604167e-09,
      1.6380142034422285e-10,
      2.3414967881274507e-11,
      3.4141717231150892e-12,
      5.062547603351675e-13,
      7.616823838318965e-14,
      1.1230599783473849e-14,
      1.7832957333041577e-15,
      3.7470027081099033e-16,
      6.938893903907228e-18,
      -2.185751579730777e-16,
      -4.336808689942018e-17
    },
    0.2962311921296296,
    0.296263744212963
  },
  {
    {
      0.3273360915970311,
      0.0010158203926102366,
      4.357198950115382e-05,
      3.7307483612414316e-06,
      3.9960656192036814e-07,
      4.795694848874765e-08,
      6.167757533870333e-09,
      8.311230384472346e-10,
      1.1582474909732632e-10,
      1.6557578225162572e-11,
      2.414311806031577e-12,
      3.581648866379794e-13,
      5.39950029132541e-14,
      8.791578576250458e-15,
      1.2351231148954867e-15,
      -8.847089727481716e-16,
      -1.1102230246251565e-16,
      -5.551115123125783e-16,
      -7.112366251504909e-17
    },
    0.296263744212963,
    0.29628002025462963
  },
  {
    {
      0.3290903207332936,
      0.0007190748744181456,
      3.080907898764598e-05,
      2.6380102197860744e-06,
      2.8256276604285424e-07,
      3.391053608747807e-08,
      4.3612495316536215e-09,
      5.876913573843456e-10,
      8.189969188943103e-11,
      1.1706732805372155e-11,
      1.7067146307336856e-12,
      2.552090483387559e-13,
      3.920475055707584e-14,
      2.9976021664879227e-15,
      1.0408340855860843e-16,
      -2.7755575615628914e-16,
      -4.891920202254596e-16,
      8.708311849403572e-16,
      5.863365348801608e-16
    },
    0.29628002025462963,
    0.29628815827546295
  },
  {
    {
      0.3303319160946099,
      0.0005088500431879688,
      2.178496652304948e-05,
      1.8653452835380724e-06,
      1.998014318195518e-07,
      2.39783138934746e-08,
      3.0838626687013093e-09,
      4.1555842447693614e-10,
      5.791174831348833e-11,
      8.2747732277344e-12,
      1.2066458943138514e-12,
      1.7669546381604562e-13,
      2.7006175074006933e-14,
      5.43315392675936e-15,
      8.916478666520788e-16,
      1.3530843112619095e-15,
      -5.898059818321144e-17,
      -2.0643209364124004e-15,
      -2.0990154059319366e-16
    },
    0.29628815827546295,
    0.2962922272858796
  },
  {
    {
      0.3312104361553831,
      0.0003600037890292736,
      1.5404176808297282e-05,
      1.3189948922485983e-06,
      1.412807293793561e-07,
      1.695521050990023e-08,
      2.1806181235706212e-09,
      2.9384080932337575e-10,
      4.094972624879567e-11,
      5.85331783042875e-12,
      8.538031393001688e-13,
      1.289766904388756e-13,
      1.903685542536948e-14,
      -1.1622647289044608e-15,
      9.367506770274758e-17,
      -2.6020852139652106e-16,
      5.551115123125783e-17,
      9.957312752106873e-16,
      5.551115123125783e-17
    },
    0.2962922272858796,
    0.296294261791088
  },
  {
    {
      0.331831931854538,
      0.00025465693727301855,
      1.0892355242906465e-05,
      9.326690283938122e-07,
      9.990048491226311e-08,
      1.1989139896140255e-08,
      1.541930194637331e-09,
      2.0778642001451075e-10,
      2.8956680803160495e-11,
      4.138932252484295e-12,
      6.032882526874062e-13,
      8.537268114672258e-14,
      1.3461454173580023e-14,
      3.986394547794703e-15,
      4.787836793695988e-16,
      -2.6020852139652106e-16,
      -2.7755575615628914e-16,
      -2.789435349370706e-15,
      -2.7755575615628914e-16
    },
    0.296294261791088,
    0.29629527904369213
  },
  {
    {
      0.33254766786859097,
      0.0005250329541100184,
      0.00010585696841291029,
      4.5948894997410034e-05,
      2.5987955723599104e-05,
      1.6924206122436586e-05,
      1.2054006867219524e-05,
      9.143024898339985e-06,
      7.271484406146805e-06,
      6.0039811284016975e-06,
      5.112767082184055e-06,
      4.469487214910778e-06,
      3.997564041510054e-06,
      3.6492975524389515e-06,
      3.394053607587433e-06,
      3.2118298484914987e-06,
      3.0896085961538977e-06,
      3.0192461833444817e-06,
      1.4981339496820262e-06
    },
    0.29629527904369213,
    0.2962962962962963
  }
}
};
 
const auto SRK_p = SuperAncillary{
{
  {
    {
      2.9225473692726744e-13,
      4.846211782171669e-13,
      2.824570002653489e-13,
      1.1955110559200177e-13,
      3.771710795961595e-14,
      9.009392578047045e-15,
      1.6354246747926244e-15,
      2.2297564454973646e-16,
      2.1977303354380907e-17,
      1.4079420760456038e-18,
      3.645901247485968e-20,
      -2.3224668904548756e-21,
      -2.1119277731077838e-22,
      1.3269831236005753e-24,
      7.111685153229056e-25,
      2.0749322108366765e-27,
      -2.355803670950289e-27,
      1.89511506527704e-29,
      -1.203321029253145e-29
    },
    0.020267685653535945,
    0.02596797224359293
  },
  {
    {
      4.42411833611544e-11,
      6.486786260066829e-11,
      2.8834316573742172e-11,
      8.49036232570771e-12,
      1.7382606171711586e-12,
      2.519225194793423e-13,
      2.5510836572731127e-14,
      1.688802065855514e-15,
      5.658115744876586e-17,
      -7.526264008617426e-19,
      -1.3472892665001768e-19,
      -9.922698129701787e-22,
      2.6795301476514526e-22,
      2.2940832430296035e-24,
      -5.789860391087721e-25,
      1.1649897848303902e-26,
      2.589988124021626e-26,
      5.472801416061292e-26,
      2.268132874691453e-26
    },
    0.02596797224359293,
    0.03166825883364991
  },
  {
    {
      1.3403609277725988e-08,
      2.0613817317201788e-08,
      9.908855672700315e-09,
      3.1242347660239074e-09,
      6.566525761190694e-10,
      8.934527720041564e-11,
      6.864947922093292e-12,
      1.0961157691701013e-13,
      -2.552558822185046e-14,
      -1.0755067059569204e-15,
      1.1595758568473605e-16,
      4.364935105515505e-18,
      -6.518267188981206e-19,
      -2.420947626615306e-21,
      3.3222352209999887e-21,
      -1.4160136552428753e-22,
      -1.0359144702519558e-23,
      3.0017608947725106e-24,
      -1.1502980393315541e-24
    },
    0.03166825883364991,
    0.04306883201376388
  },
  {
    {
      3.6266120105336493e-06,
      5.40530078063072e-06,
      2.3535707621866956e-06,
      6.09797026711632e-07,
      8.755959081374292e-08,
      4.4744907993962915e-09,
      -4.067655714209484e-10,
      -3.591728379118743e-11,
      4.465366516258063e-12,
      1.6493814952480902e-13,
      -5.3226243319885286e-14,
      1.7922687946152793e-15,
      4.005718486965302e-16,
      -4.988684915276682e-17,
      4.789814152406683e-19,
      4.147268524381617e-19,
      -4.591183271914325e-20,
      4.350970022028926e-22,
      -1.0629270871306406e-21
    },
    0.04306883201376388,
    0.06586997837399182
  },
  {
    {
      0.0003094297409003984,
      0.0004117661948641284,
      0.0001333784156235523,
      1.953127675250332e-05,
      4.4650824621707415e-07,
      -1.245752753612688e-07,
      6.1713045177982556e-09,
      8.457787641085228e-10,
      -2.019178315130611e-10,
      1.671511339141844e-11,
      7.461694230598466e-13,
      -4.260132318823111e-13,
      6.171897036311516e-14,
      -3.212145765269425e-15,
      -5.836180885150665e-16,
      1.719715110942645e-16,
      -2.147440279526465e-17,
      8.273394312306379e-19,
      2.34734005507848e-19
    },
    0.06586997837399182,
    0.1114722710944477
  },
  {
    {
      0.0074771598568592074,
      0.008244590094104293,
      0.0017465545422224592,
      9.786738241774412e-05,
      -6.516977077872878e-06,
      3.2958650511163794e-07,
      1.5363696321735632e-08,
      -7.20353161771562e-09,
      1.2861343701874766e-09,
      -1.6587511786082245e-10,
      1.5108394000934097e-11,
      -3.6747905743763146e-13,
      -2.33340277167645e-13,
      7.039936125372015e-14,
      -1.3651690997269761e-14,
      2.071300487897776e-15,
      -2.45029690981724e-16,
      2.15485181781494e-17,
      1.4026865606531214e-18
    },
    0.1114722710944477,
    0.20267685653535944
  }
}
};
const auto SRK_rhoL = SuperAncillary{
{
  {
    {
      0.9087646203247287,
      -0.050607640245558645,
      -0.0008181930140446145,
      -3.676466126572436e-05,
      -1.8153449552124612e-06,
      -4.2148973304623194e-08,
      8.967883102783958e-09,
      5.813281384048352e-10,
      -3.659740635542086e-10,
      -4.441783746367278e-11,
      1.5951594212193498e-11,
      1.7195689316906737e-12,
      -9.141229440068344e-13,
      -1.4432899320127035e-14,
      5.0730253331465747e-14,
      -6.987466161234579e-15,
      -1.8318679906315083e-15,
      1.5126788710517758e-15,
      -1.734723475976807e-17
    },
    0.020267685653535945,
    0.06586997837399182
  },
  {
    {
      0.7992363370107459,
      -0.059445681918569984,
      -0.0014502328603496623,
      -7.138164945971748e-05,
      -3.2189087307149533e-06,
      -2.486464586864967e-07,
      -2.4319133658246006e-08,
      -1.044768556390796e-09,
      -1.3343107174712543e-10,
      -1.679712618996021e-11,
      2.2661039711380226e-13,
      -1.6726897644758765e-13,
      -1.0783041126671833e-14,
      1.4363510381087963e-15,
      -4.0245584642661925e-16,
      -7.28583859910259e-16,
      -4.85722573273506e-17,
      7.632783294297951e-16,
      1.734723475976807e-17
    },
    0.06586997837399182,
    0.1114722710944477
  },
  {
    {
      0.6649628598972782,
      -0.0760521559973085,
      -0.002952854145734922,
      -0.0002237665443423284,
      -2.2837440018151633e-05,
      -2.7674939431024392e-06,
      -3.5249351108129767e-07,
      -4.743986068561634e-08,
      -6.602545343792343e-09,
      -9.417646246179245e-10,
      -1.37257025190074e-10,
      -2.031076695718781e-11,
      -3.0466879019641624e-12,
      -4.624911564832246e-13,
      -7.074202335033419e-14,
      -1.1532441668293814e-14,
      -1.7277845820728999e-15,
      3.885780586188048e-16,
      -2.42861286636753e-17
    },
    0.1114722710944477,
    0.15707456381490356
  },
  {
    {
      0.5358548459176075,
      -0.051616853474821534,
      -0.0019014483430589563,
      -0.0001522452155350712,
      -1.6138121186740073e-05,
      -1.9310179919518333e-06,
      -2.4780931559553743e-07,
      -3.3348723178783235e-08,
      -4.64304821673478e-09,
      -6.632207977252946e-10,
      -9.664995193059411e-11,
      -1.4312308282971031e-11,
      -2.1472545963518996e-12,
      -3.260794412263124e-13,
      -4.987676938128516e-14,
      -8.229528170033973e-15,
      -1.2004286453759505e-15,
      3.469446951953614e-16,
      -2.0816681711721685e-17
    },
    0.15707456381490356,
    0.1798757101751315
  },
  {
    {
      0.4485732937583308,
      -0.03477837295500611,
      -0.00128197185101675,
      -0.00010638749439246392,
      -1.1344133512683485e-05,
      -1.3588779068512968e-06,
      -1.7457610832632975e-07,
      -2.350823218033593e-08,
      -3.2744967258524493e-09,
      -4.678942372438044e-10,
      -6.820332146273422e-11,
      -1.0101888076041732e-11,
      -1.5157701482859665e-12,
      -2.302880108828731e-13,
      -3.522182545623309e-14,
      -5.863365348801608e-15,
      -8.500145032286355e-16,
      3.0878077872387166e-16,
      -1.9081958235744878e-17
    },
    0.1798757101751315,
    0.19127628335524546
  },
  {
    {
      0.3896887959074917,
      -0.02352363480040463,
      -0.0008864377996999455,
      -7.48085471366268e-05,
      -7.994798847377887e-06,
      -9.585173595394458e-07,
      -1.2320694868853121e-07,
      -1.6596597669088498e-08,
      -2.312318735792074e-09,
      -3.3046667832770105e-10,
      -4.817745855034694e-11,
      -7.136409518881948e-12,
      -1.0708690878491467e-12,
      -1.627031842588167e-13,
      -2.490022077417109e-14,
      -4.173744683200198e-15,
      -5.93275428784068e-16,
      3.0184188481996443e-16,
      -1.214306433183765e-17
    },
    0.19127628335524546,
    0.19697656994530244
  },
  {
    {
      0.3497350745917645,
      -0.016032294782450524,
      -0.0006202549276113981,
      -5.274770800074019e-05,
      -5.643402364315098e-06,
      -6.769284142359167e-07,
      -8.70356488269497e-08,
      -1.17262048994482e-08,
      -1.6339517863728226e-09,
      -2.3353845593487854e-10,
      -3.404898932246603e-11,
      -5.04390973432578e-12,
      -7.568771998034407e-13,
      -1.1509196373715724e-13,
      -1.7600504387260685e-14,
      -3.084338340286763e-15,
      -4.440892098500626e-16,
      2.8449465006019636e-16,
      -3.469446951953614e-18
    },
    0.19697656994530244,
    0.19982671324033094
  },
  {
    {
      0.3224124351709759,
      -0.011013634884520482,
      -0.00043639529901024973,
      -3.724427437728703e-05,
      -3.986980578418037e-06,
      -4.783591501592022e-07,
      -6.151332742135307e-08,
      -8.288353604141152e-09,
      -1.1549863634940927e-09,
      -1.6508804739912009e-10,
      -2.4069971710227733e-11,
      -3.5656928798477594e-12,
      -5.350303533546708e-13,
      -8.137934770502397e-14,
      -1.2437967322753707e-14,
      -2.203098814490545e-15,
      -3.157196726277789e-16,
      2.5673907444456745e-16,
      0.0
    },
    0.19982671324033094,
    0.20125178488784518
  },
  {
    {
      0.30358578390142943,
      -0.0076188736955520135,
      -0.0003078315487315582,
      -2.6316343990171603e-05,
      -2.8179719795684566e-06,
      -3.381439433146338e-07,
      -4.348579207141823e-08,
      -5.8595724626575585e-09,
      -8.165597746578879e-10,
      -1.1671785066225127e-10,
      -1.701778301610446e-11,
      -2.5209348497590156e-12,
      -3.7829808730016623e-13,
      -5.745404152435185e-14,
      -8.739536871971154e-15,
      -1.5508427875232655e-15,
      -1.9081958235744878e-16,
      2.671474153004283e-16,
      -3.642919299551295e-17
    },
    0.20125178488784518,
    0.2019643207116023
  },
  {
    {
      0.29052982685479933,
      -0.005300319842511697,
      -0.00021741216965811758,
      -1.8601576619783206e-05,
      -1.99216392636134e-06,
      -2.390659452899957e-07,
      -3.0745310054813846e-08,
      -4.142926245120915e-09,
      -5.773458308655499e-10,
      -8.252584726697876e-11,
      -1.2032597140887447e-11,
      -1.7825428633155838e-12,
      -2.674006849279209e-13,
      -4.087702398791748e-14,
      -6.238065619612598e-15,
      -1.2177758801357186e-15,
      -1.457167719820518e-16,
      3.5735303605122226e-16,
      1.214306433183765e-17
    },
    0.2019643207116023,
    0.20232058862348087
  },
  {
    {
      0.281429493928621,
      -0.00370348431999927,
      -0.00015364397898540855,
      -1.3150855018025098e-05,
      -1.408515659318682e-06,
      -1.6903172642135367e-07,
      -2.1738876799376472e-08,
      -2.929343496937964e-09,
      -4.0822773633708564e-10,
      -5.835245481256024e-11,
      -8.508027615761193e-12,
      -1.2602800747441023e-12,
      -1.8902587828328876e-13,
      -2.898029038966854e-14,
      -4.3923198411732756e-15,
      -6.83481049534862e-16,
      -7.28583859910259e-17,
      6.210310043996969e-16,
      3.642919299551295e-17
    },
    0.20232058862348087,
    0.20249872257942014
  },
  {
    {
      0.27506157036291856,
      -0.0025962448144549836,
      -0.00010861133767400138,
      -9.298192002037131e-06,
      -9.959154303774975e-07,
      -1.1951873122098555e-07,
      -1.5371233163044562e-08,
      -2.0713067025446286e-09,
      -2.886544493013732e-10,
      -4.1260670227694085e-11,
      -6.015958564642432e-12,
      -8.912228594004645e-13,
      -1.3355809513893036e-13,
      -1.9872992140790302e-14,
      -3.0531133177191805e-15,
      -7.181755190543981e-16,
      -8.500145032286355e-17,
      2.96637714392034e-16,
      2.6020852139652106e-18
    },
    0.20249872257942014,
    0.2025877895573898
  },
  {
    {
      0.2705926697847878,
      -0.0018244565862979688,
      -7.678880274703502e-05,
      -6.574507803324206e-06,
      -7.04198907909917e-07,
      -8.451082581135971e-08,
      -1.0868935571006766e-08,
      -1.4646166258958093e-09,
      -2.0410736989440181e-10,
      -2.9175961993588295e-11,
      -4.253862886938187e-12,
      -6.296716620335374e-13,
      -9.435334458185451e-14,
      -1.411544492402328e-14,
      -2.1076890233118206e-15,
      1.1622647289044608e-16,
      -1.3877787807814457e-17,
      2.0816681711721685e-16,
      -6.071532165918825e-18
    },
    0.2025877895573898,
    0.20263232304637463
  },
  {
    {
      0.26744977424244143,
      -0.0012843645575066134,
      -5.4294008668638594e-05,
      -4.6487705078930575e-06,
      -4.979368752669466e-07,
      -5.97575837875397e-08,
      -7.685439108034986e-09,
      -1.035633939716618e-09,
      -1.4432471884262554e-10,
      -2.0630074037963908e-11,
      -3.0077728502275747e-12,
      -4.461934294264225e-13,
      -6.715808464896611e-14,
      -1.0720591081536668e-14,
      -1.4190038033490282e-15,
      -4.440892098500626e-16,
      1.8735013540549517e-16,
      8.378714388967978e-16,
      -1.214306433183765e-16
    },
    0.20263232304637463,
    0.20265458979086703
  },
  {
    {
      0.2652360092057196,
      -0.0009053089022688643,
      -3.839029629034568e-05,
      -3.2871387625962867e-06,
      -3.5209208477243736e-07,
      -4.2254784272249823e-08,
      -5.434404972901041e-09,
      -7.32300850167511e-10,
      -1.0205274299235345e-10,
      -1.458767828754759e-11,
      -2.126882003850028e-12,
      -3.151055805172831e-13,
      -4.713764101271778e-14,
      -7.455841499748317e-15,
      -1.1188966420050406e-15,
      -4.2500725161431774e-16,
      -1.9081958235744878e-17,
      4.0072112295064244e-16,
      1.734723475976807e-17
    },
    0.20265458979086703,
    0.20266572316311324
  },
  {
    {
      0.2636749575957985,
      -0.0006387083278175353,
      -2.7145557162256556e-05,
      -2.324344536088624e-06,
      -2.4896583217942636e-07,
      -2.9878569729469007e-08,
      -3.8426970945526495e-09,
      -5.178131846916845e-10,
      -7.21618736987395e-11,
      -1.0313611076284701e-11,
      -1.5037519840443991e-12,
      -2.2252685805135286e-13,
      -3.314015728506092e-14,
      -3.944761184371259e-15,
      -7.4593109467002705e-16,
      -1.061650767297806e-15,
      -7.28583859910259e-17,
      3.1051550219984847e-16,
      -8.673617379884035e-19
    },
    0.20266572316311324,
    0.20267128984923632
  },
  {
    {
      0.26257329202514273,
      -0.0004509124480249222,
      -1.9194637512940163e-05,
      -1.6435549859191573e-06,
      -1.7604512139550443e-07,
      -2.11273140666296e-08,
      -2.7171948059623308e-09,
      -3.6615067228185083e-10,
      -5.1026160727274394e-11,
      -7.291847681223373e-12,
      -1.0631721197862376e-12,
      -1.570323732158485e-13,
      -2.3507237822961713e-14,
      -4.498137973207861e-15,
      -5.672545766444159e-16,
      3.8163916471489756e-16,
      -1.1622647289044608e-16,
      -1.5334955527634975e-15,
      -1.6046192152785466e-16
    },
    0.20267128984923632,
    0.20267407319229788
  },
  {
    {
      0.261795381349911,
      -0.00031848138784485065,
      -1.3572598294689722e-05,
      -1.1621671743376055e-06,
      -1.2448259038670695e-07,
      -1.493926039992932e-08,
      -1.921346123867629e-09,
      -2.5890705378228684e-10,
      -3.6080990625797504e-11,
      -5.1575497345135446e-12,
      -7.518204808709683e-13,
      -1.1052617154838629e-13,
      -1.6753959330984003e-14,
      -5.920611223508843e-15,
      -8.795048023202412e-16,
      -2.8033131371785203e-15,
      -1.5265566588595902e-16,
      1.177877240188252e-15,
      8.153200337090993e-17
    },
    0.20267407319229788,
    0.20267546486382865
  },
  {
    {
      0.2612458584790565,
      -0.00022501920934304802,
      -9.597255076118291e-06,
      -8.217756914414426e-07,
      -8.802244447041196e-08,
      -1.0563647579939217e-08,
      -1.358598145442813e-09,
      -1.83071792100864e-10,
      -2.551181141341452e-11,
      -3.646784263455771e-12,
      -5.306519113013053e-13,
      -7.894206122127656e-14,
      -1.3612375115990005e-14,
      -1.3461454173580023e-15,
      7.650130529057719e-16,
      -1.474514954580286e-16,
      6.262351748276274e-16,
      -3.9517000782751666e-15,
      -1.343543332144037e-15
    },
    0.20267546486382865,
    0.20267616069959404
  },
  {
    {
      0.26061394134803445,
      -0.0004632477727018712,
      -9.327071115945024e-05,
      -4.048564047778869e-05,
      -2.2898027861414927e-05,
      -1.4911944355367607e-05,
      -1.0620804294449865e-05,
      -8.055933548379093e-06,
      -6.4069163156280085e-06,
      -5.290117203090791e-06,
      -4.504867114365346e-06,
      -3.938072224525177e-06,
      -3.5222599736994126e-06,
      -3.2154018167716847e-06,
      -2.990505978434088e-06,
      -2.8299483429503752e-06,
      -2.7222590090695414e-06,
      -2.660262590837778e-06,
      -1.3200081928170515e-06
    },
    0.20267616069959404,
    0.20267685653535944
  }
}
};
const auto SRK_rhoV = SuperAncillary{
{
  {
    {
      1.147791769634607e-11,
      1.8883400952870644e-11,
      1.0781377594967363e-11,
      4.424169417038083e-12,
      1.340057930622405e-12,
      3.0392837657011474e-13,
      5.1611350145284516e-14,
      6.429591419258869e-15,
      5.526656713208219e-16,
      2.6838869757647005e-17,
      -7.144055741834159e-20,
      -9.585729650435423e-20,
      -3.2421103573493004e-21,
      2.5635377164771786e-22,
      1.5013435876237105e-23,
      -8.328903877850174e-25,
      -5.167748904012326e-26,
      3.2366962941218114e-27,
      -1.163964265653603e-27
    },
    0.020267685653535945,
    0.02596797224359293
  },
  {
    {
      1.4330520293934438e-09,
      2.0653900997266785e-09,
      8.861422337329386e-10,
      2.484390807063457e-10,
      4.771871094962088e-11,
      6.3539423424235915e-12,
      5.695290133656733e-13,
      3.041636735038608e-14,
      4.618066754497384e-16,
      -4.882537910833838e-17,
      -2.2626929302787274e-18,
      7.735800358738017e-20,
      5.494308515884454e-21,
      -1.9098299796206625e-22,
      -1.0551803468236254e-23,
      9.358288473073072e-25,
      8.243533350687156e-25,
      1.8130926610110065e-24,
      6.625926732877158e-25
    },
    0.02596797224359293,
    0.03166825883364991
  },
  {
    {
      3.216589566610815e-07,
      4.855615137205529e-07,
      2.2313237230527845e-07,
      6.564126429055284e-08,
      1.2457975717948848e-08,
      1.4356027589196413e-09,
      7.456366737354983e-11,
      -2.6879017209951425e-12,
      -4.739337426481451e-13,
      7.722229482544613e-15,
      2.5616030504147617e-15,
      -7.662739526444225e-17,
      -1.1817890172164798e-17,
      8.184648689324978e-19,
      2.707129500338779e-20,
      -5.7904193842360474e-21,
      1.4315394475995835e-22,
      8.504450672810816e-23,
      -2.4996364135586805e-23
    },
    0.03166825883364991,
    0.04306883201376388
  },
  {
    {
      5.7939913777071644e-05,
      8.379335832127809e-05,
      3.377936608738917e-05,
      7.637987852998414e-06,
      8.302628510075688e-07,
      1.4927574593586637e-09,
      -6.691320992196064e-09,
      1.2878253547929106e-10,
      7.094465860674562e-11,
      -5.0660797908678295e-12,
      -4.994533593074178e-13,
      9.317691594005963e-14,
      -1.768263545468354e-15,
      -8.453117101310306e-16,
      9.273541699666695e-17,
      -5.0145673966434665e-19,
      -7.953374676921863e-19,
      8.536503921547246e-20,
      -1.495542547495874e-20
    },
    0.04306883201376388,
    0.06586997837399182
  },
  {
    {
      0.0031631846269420103,
      0.003982211787182989,
      0.0011297528144810199,
      0.00012650572507611033,
      -5.2717473574936e-07,
      -1.4260061216196563e-07,
      1.3326423389224526e-07,
      -8.09481292393274e-09,
      -4.78981067790718e-10,
      3.091240719345703e-10,
      -4.0682128710022676e-11,
      1.5201595949913826e-12,
      6.298290797186528e-13,
      -1.5266145620318645e-13,
      1.6769319426449512e-14,
      2.1858532236844475e-17,
      -3.670075175840268e-16,
      7.421194464073827e-17,
      -6.8516495103400354e-18
    },
    0.06586997837399182,
    0.1114722710944477
  },
  {
    {
      0.02527054526018694,
      0.01958885465789757,
      0.0029253705571776014,
      0.00022445447854577003,
      2.1127284490167714e-05,
      2.925789705565077e-06,
      3.436409961475248e-07,
      4.7523257336800576e-08,
      6.65800338960857e-09,
      9.32793360239846e-10,
      1.3817232769183438e-10,
      2.024545819618445e-11,
      3.049514308607615e-12,
      4.623933614472664e-13,
      7.069220426050848e-14,
      1.0914066959216706e-14,
      1.6929274822274909e-15,
      2.8888566885876266e-16,
      4.003416521902725e-17
    },
    0.1114722710944477,
    0.15707456381490356
  },
  {
    {
      0.06954212908713689,
      0.023253617876107757,
      0.0018817558019952863,
      0.00015099612130791773,
      1.6119981407421245e-05,
      1.9323965214038813e-06,
      2.477548253195308e-07,
      3.3350354563480705e-08,
      4.6430122303299515e-09,
      6.632208966045328e-10,
      9.665026374713892e-11,
      1.4312801378119078e-11,
      2.1476913129869768e-12,
      3.25807089640584e-13,
      4.9888478764748e-14,
      7.651431571664702e-15,
      1.1947907940790259e-15,
      2.40692882291782e-16,
      2.3635607360183997e-17
    },
    0.15707456381490356,
    0.1798757101751315
  },
  {
    {
      0.11419041805449423,
      0.020521580327712844,
      0.0012740772212567466,
      0.0001062161198180098,
      1.1343925939054193e-05,
      1.3588982222490983e-06,
      1.7457566693114546e-07,
      2.3508239255405627e-08,
      3.274496806517091e-09,
      4.678940377506047e-10,
      6.82035469767861e-11,
      1.0102339971507224e-11,
      1.5161908187288908e-12,
      2.3003040444669054e-13,
      3.523483588230292e-14,
      5.3507545616504615e-15,
      8.448103328007051e-16,
      2.393918396847994e-16,
      2.0816681711721685e-17
    },
    0.1798757101751315,
    0.19127628335524546
  },
  {
    {
      0.15165818890499952,
      0.016369427376503674,
      0.0008840770152282037,
      7.478700795372859e-05,
      7.99480167670754e-06,
      9.58517782677533e-07,
      1.2320694370814012e-07,
      1.659659774975314e-08,
      2.3123189231422092e-09,
      3.3046648403867174e-10,
      4.817769534010141e-11,
      7.1368752921352474e-12,
      1.0712420533964817e-12,
      1.6248200701562965e-13,
      2.4886342986363275e-14,
      3.694961003830599e-15,
      5.924080670460796e-16,
      2.1337098754514727e-16,
      1.7780915628762273e-17
    },
    0.19127628335524546,
    0.19697656994530244
  },
  {
    {
      0.1808715039787546,
      0.0124478423701576,
      0.0006196163285778558,
      5.2745024120972456e-05,
      5.643402899935257e-06,
      6.769284248055868e-07,
      8.703564869511071e-08,
      1.1726204920264882e-08,
      1.6339519685187875e-09,
      2.335382477680614e-10,
      3.40492009587301e-11,
      5.044305251278303e-12,
      7.572241444986361e-13,
      1.1483695938618865e-13,
      1.7590096046404824e-14,
      2.5847379792054426e-15,
      4.198030811863873e-16,
      2.3765711620882257e-16,
      3.469446951953614e-18
    },
    0.19697656994530244,
    0.19982671324033094
  },
  {
    {
      0.202815039861727,
      0.00921945969494375,
      0.00043622961797678585,
      3.7243939807177565e-05,
      3.986980621388872e-06,
      4.783591503864509e-07,
      6.151332734329051e-08,
      8.28835361281477e-09,
      1.154986554313675e-09,
      1.6508787566149596e-10,
      2.4070209367343942e-11,
      3.5661248259932776e-12,
      5.35387706390722e-13,
      8.11677114409548e-14,
      1.2444906216657614e-14,
      1.7468665403086447e-15,
      2.8102520310824275e-16,
      2.3071822230491534e-16,
      6.938893903907228e-18
    },
    0.19982671324033094,
    0.20125178488784518
  },
  {
    {
      0.2189498032807097,
      0.006721284919406446,
      0.0003077893762559753,
      2.6316302236593578e-05,
      2.817971982508813e-06,
      3.381439431966726e-07,
      4.3485791989886224e-08,
      5.859572497352028e-09,
      8.165599741510876e-10,
      1.16717652903775e-10,
      1.701801546905024e-11,
      2.5213789389688657e-12,
      3.7863462365450573e-13,
      5.725107887766256e-14,
      8.748210489351038e-15,
      1.0911410663894117e-15,
      1.682681771697503e-16,
      2.237793284010081e-16,
      3.209238430557093e-17
    },
    0.20125178488784518,
    0.2019643207116023
  },
  {
    {
      0.23065921854671798,
      0.004851398420450417,
      0.00021740153263174328,
      1.8601571405155864e-05,
      1.9921639265139957e-06,
      2.3906594516856505e-07,
      3.074531000103742e-08,
      4.142926290223725e-09,
      5.773459800517688e-10,
      8.252565471267292e-11,
      1.2032836532727131e-11,
      1.7829522580559143e-12,
      2.677597726874481e-13,
      4.066538772384831e-14,
      6.241535066564552e-15,
      7.580741590018647e-16,
      1.3530843112619095e-16,
      1.5612511283791264e-16,
      -9.540979117872439e-18
    },
    0.2019643207116023,
    0.20232058862348087
  },
  {
    {
      0.23908612939229767,
      0.0034789916333673055,
      0.00015364130799821497,
      1.3150854366506329e-05,
      1.4085156592926612e-06,
      1.6903172630859664e-07,
      2.173887675774311e-08,
      2.9293435507143917e-09,
      4.0822789593164543e-10,
      5.835226572770136e-11,
      8.508240986748739e-12,
      1.2606964083783367e-12,
      1.8938843548976791e-13,
      2.8749572167363624e-14,
      4.397524011601206e-15,
      2.5847379792054426e-16,
      4.5102810375396984e-17,
      -9.540979117872439e-17,
      -2.862293735361732e-17
    },
    0.20232058862348087,
    0.20249872257942014
  },
  {
    {
      0.24511730390452105,
      0.002483990450115506,
      0.0001086106684613404,
      9.298191920630028e-06,
      9.959154303011697e-07,
      1.1951873112384104e-07,
      1.537123309712507e-08,
      2.0713067459127155e-09,
      2.8865461930427383e-10,
      4.126050022479344e-11,
      6.016184078694309e-12,
      8.916461319286029e-13,
      1.339206523454095e-13,
      1.9670029494101016e-14,
      3.0687258290029718e-15,
      2.6020852139652106e-16,
      7.112366251504909e-17,
      2.0296264668928643e-16,
      -8.673617379884035e-19
    },
    0.20249872257942014,
    0.2025877895573898
  },
  {
    {
      0.24941782042567914,
      0.0017683273954829087,
      7.678863526062486e-05,
      6.57450779313444e-06,
      7.041989078648142e-07,
      8.45108257072763e-08,
      1.086893550855672e-08,
      1.4646166571208319e-09,
      2.0410752775423813e-10,
      2.9175781582346794e-11,
      4.254084931543112e-12,
      6.301330984781472e-13,
      9.470722817095378e-14,
      1.391421700080997e-14,
      2.114627917215728e-15,
      -5.585809592645319e-16,
      -1.214306433183765e-17,
      2.8622937353617317e-16,
      3.469446951953614e-18
    },
    0.2025877895573898,
    0.20263232304637463
  },
  {
    {
      0.2524765215536145,
      0.0012562994595141843,
      5.429396677408306e-05,
      4.648770506642322e-06,
      4.979368752443952e-07,
      5.975758369559936e-08,
      7.685439031707153e-09,
      1.0356339657374702e-09,
      1.4432488711080271e-10,
      2.0629898830892834e-11,
      3.0079948948324997e-12,
      4.4662190612498875e-13,
      6.752064185544526e-14,
      1.0500281200087613e-14,
      1.4103301859691442e-15,
      -1.734723475976807e-17,
      -1.8908485888147197e-16,
      -3.4867941867133823e-16,
      1.2663481374630692e-16
    },
    0.20263232304637463,
    0.20265458979086703
  },
  {
    {
      0.2546481887862289,
      0.0008912762275735033,
      3.839028581376948e-05,
      3.287138762469652e-06,
      3.5209208471866094e-07,
      4.225478415602335e-08,
      5.434404913920443e-09,
      7.32300886596704e-10,
      1.02052919934148e-10,
      1.4587508284646944e-11,
      2.1270849664967173e-12,
      3.1551150581066167e-13,
      4.750713711310084e-14,
      7.212980213111564e-15,
      1.1084883011491797e-15,
      -6.591949208711867e-17,
      1.734723475976807e-18,
      9.367506770274758e-17,
      -9.540979117872439e-18
    },
    0.20265458979086703,
    0.20266572316311324
  },
  {
    {
      0.2561881913448183,
      0.0006316919590384738,
      2.7145554542667982e-05,
      2.3243445360920933e-06,
      2.4896583211350687e-07,
      2.9878569618446704e-08,
      3.842697030367881e-09,
      5.178132280597714e-10,
      7.216204023219319e-11,
      1.0313451481724911e-11,
      1.5039688244788962e-12,
      2.2294839585601522e-13,
      3.350444921501605e-14,
      3.740063814205996e-15,
      7.407269242420966e-16,
      5.967448757360216e-16,
      6.591949208711867e-17,
      1.8908485888147197e-16,
      -7.806255641895632e-18
    },
    0.20266572316311324,
    0.20267128984923632
  },
  {
    {
      0.2572793323524821,
      0.00044740425577668823,
      1.919463685796409e-05,
      1.6435549859573212e-06,
      1.7604512135387107e-07,
      2.1127313966015637e-08,
      2.717194729634498e-09,
      3.661506878943621e-10,
      5.1026320321834184e-11,
      7.29169849500444e-12,
      1.0634045727320185e-12,
      1.5748687076655443e-13,
      2.3859386688585005e-14,
      4.270889197854899e-15,
      5.672545766444159e-16,
      -8.448103328007051e-16,
      9.71445146547012e-17,
      2.0105445086571194e-15,
      1.448494102440634e-16
    },
    0.20267128984923632,
    0.20267407319229788
  },
  {
    {
      0.2580519807368856,
      0.0003167272897560327,
      1.3572598130872846e-05,
      1.1621671743965861e-06,
      1.2448259037282916e-07,
      1.4939260268090337e-08,
      1.921346066621754e-09,
      2.589070850073094e-10,
      3.6081145016186866e-11,
      5.157388405230279e-12,
      7.520546685402252e-13,
      1.1097546492866428e-13,
      1.7126924878319016e-14,
      5.67774993687209e-15,
      8.673617379884035e-16,
      2.3071822230491534e-15,
      1.43982048506075e-16,
      -6.661338147750939e-16,
      -7.892991815694472e-17
    },
    0.20267407319229788,
    0.20267546486382865
  },
  {
    {
      0.2585988724599927,
      0.00022414215980744862,
      9.597255035097285e-06,
      8.217756915021579e-07,
      8.802244440969664e-08,
      1.0563647474121085e-08,
      1.358598076053874e-09,
      1.8307183026478047e-10,
      2.5511984885762118e-11,
      3.646605586937746e-12,
      5.308739559062303e-13,
      7.936012957898697e-14,
      1.395758508770939e-14,
      1.1587952819525071e-15,
      -7.546047120499111e-16,
      -3.0704605524789486e-16,
      -6.401129626354418e-16,
      4.4374226515486725e-15,
      1.3426759704060487e-15
    },
    0.20267546486382865,
    0.20267616069959404
  },
  {
    {
      0.2592290354916159,
      0.00046237072283881325,
      9.327071111846566e-05,
      4.04856404778442e-05,
      2.2898027861355946e-05,
      1.4911944355256584e-05,
      1.0620804294380476e-05,
      8.055933548441543e-06,
      6.4069163158222975e-06,
      5.29011720291038e-06,
      4.504867114601269e-06,
      3.9380722249553884e-06,
      3.522259974053296e-06,
      3.215401816561783e-06,
      2.9905059784306187e-06,
      2.829948342482e-06,
      2.7222590090452553e-06,
      2.66026259132697e-06,
      1.3200081928118473e-06
    },
    0.20267616069959404,
    0.20267685653535944
  }
}
};
 
const auto PR_p = SuperAncillary{
{
  {
    {
      4.2115936784941e-14,
      7.082139279778615e-14,
      4.281441668412217e-14,
      1.9123012372347713e-14,
      6.460328562307523e-15,
      1.6765603431288853e-15,
      3.36292637163962e-16,
      5.1860921648633566e-17,
      6.0107675456840275e-18,
      4.93726813515645e-19,
      2.4083621285564296e-20,
      9.47616424537102e-23,
      -7.1051225501406e-23,
      -3.0677436333050444e-24,
      1.4628333292452202e-25,
      1.2639227645878215e-26,
      -3.8533659659909825e-28,
      -3.5261080979808146e-29,
      3.32550323497344e-30
    },
    0.01701444200703503,
    0.021799753821513633
  },
  {
    {
      8.938956773181632e-12,
      1.3439327305351806e-11,
      6.288329039800055e-12,
      1.9826748553822934e-12,
      4.414511850216345e-13,
      7.088859014953733e-14,
      8.181698505611952e-15,
      6.527653757029151e-16,
      3.162152136162365e-17,
      4.279990806191533e-19,
      -4.6206272928403434e-20,
      -2.0797341030133568e-21,
      6.177896652744867e-23,
      4.7788231227984526e-24,
      -1.3861706054362575e-25,
      -1.4215273512082633e-26,
      -1.669229675447661e-27,
      4.487040828897115e-27,
      7.7150596530614955e-28
    },
    0.021799753821513633,
    0.026585065635992236
  },
  {
    {
      4.085507279084059e-09,
      6.410401559266648e-09,
      3.231824978266208e-09,
      1.0923069013593605e-09,
      2.52024445187971e-10,
      3.905509723141962e-11,
      3.73471874979177e-12,
      1.5022634999916855e-13,
      -8.12601770161866e-15,
      -9.644616840775853e-16,
      2.3559108130823366e-17,
      4.723565049544365e-18,
      -1.6713754821230258e-19,
      -1.9677138009092646e-20,
      1.4426017766021298e-21,
      3.847561149043651e-23,
      -1.0404785039052024e-23,
      1.3764802380765226e-24,
      4.620579202932928e-25
    },
    0.026585065635992236,
    0.03615568926494944
  },
  {
    {
      1.6349625391332918e-06,
      2.489669257624972e-06,
      1.1427314989998866e-06,
      3.2189532774728283e-07,
      5.2973211090187104e-08,
      3.881322341579814e-09,
      -1.402363237723129e-10,
      -3.26068965484856e-11,
      1.5764214067874998e-12,
      2.669115963311214e-13,
      -2.6973487736840458e-14,
      -1.1725895938982276e-15,
      3.409623244140794e-16,
      -1.3374770488526952e-17,
      -2.1340688015269687e-18,
      3.0499875924358835e-19,
      -9.006032316747658e-21,
      5.945360652724886e-23,
      7.708806321111372e-22
    },
    0.03615568926494944,
    0.05529693652286385
  },
  {
    {
      0.00018641763706717313,
      0.0002549265195801485,
      8.796517801636943e-05,
      1.4493063101704085e-05,
      5.93659937733161e-07,
      -8.865727397666071e-08,
      1.5453090824947498e-09,
      1.0034143700587606e-09,
      -1.473603817532662e-10,
      4.333381441536222e-12,
      1.903659250210749e-12,
      -3.854457196096114e-13,
      2.9156300395036635e-14,
      2.2880752244019054e-15,
      -9.628463327642823e-16,
      1.2770440460045229e-16,
      -4.453540417985501e-18,
      -1.6186270226437646e-18,
      3.6149778002037437e-19
    },
    0.05529693652286385,
    0.09357943103869266
  },
  {
    {
      0.005448712641404622,
      0.0062362159058999755,
      0.001443947803998351,
      0.00010960884095892343,
      -2.254571049831034e-06,
      3.108042103893788e-07,
      3.183233642737454e-08,
      -6.3983755658182065e-09,
      1.1291747205506315e-09,
      -1.183401464725808e-10,
      7.412891426996499e-12,
      7.92415429334923e-13,
      -3.6369037393225306e-13,
      7.890138331101762e-14,
      -1.258085161656014e-14,
      1.5166023276634577e-15,
      -1.0518116325825e-16,
      -6.600080725005508e-18,
      5.5294310796760726e-18
    },
    0.09357943103869266,
    0.1701444200703503
  }
}
};
const auto PR_rhoL = SuperAncillary{
{
  {
    {
      0.9207205305399191,
      -0.045064618991104796,
      -0.0010750229061558397,
      -4.150471723450566e-05,
      -2.0449656082804912e-06,
      -8.411570305566496e-08,
      2.630701059769258e-09,
      7.618859571012493e-10,
      -1.2385135972348138e-10,
      -4.380963647410141e-11,
      3.6138314563061158e-12,
      2.0083795737591004e-12,
      -2.442004931602071e-13,
      -8.769374115757955e-14,
      2.0539125955565396e-14,
      1.4155343563970746e-15,
      -1.4502288259166107e-15,
      9.645062526431047e-16,
      4.85722573273506e-17
    },
    0.01701444200703503,
    0.05529693652286385
  },
  {
    {
      0.8199443180127295,
      -0.056330799734005496,
      -0.0018219491199177318,
      -8.821642822268161e-05,
      -4.3808005871781575e-06,
      -2.7804292488525784e-07,
      -2.5617864357618814e-08,
      -1.4936050521385802e-09,
      -1.1656688808647786e-10,
      -1.80864212495635e-11,
      -3.4572344986827375e-13,
      -7.897849041427207e-14,
      -1.9761969838327786e-14,
      1.0269562977782698e-15,
      -1.5265566588595902e-16,
      -7.91033905045424e-16,
      -2.7755575615628914e-17,
      8.049116928532385e-16,
      2.0816681711721685e-17
    },
    0.05529693652286385,
    0.09357943103869266
  },
  {
    {
      0.6880146598953387,
      -0.07711229066058617,
      -0.0036751053290779087,
      -0.00027023163542316125,
      -2.612667202559621e-05,
      -3.0579327836427472e-06,
      -3.869841956188891e-07,
      -5.180305783641925e-08,
      -7.196452149471622e-09,
      -1.0250969587066727e-09,
      -1.492256387902735e-10,
      -2.2068007143882795e-11,
      -3.3078817462950383e-12,
      -5.018416238122825e-13,
      -7.670947210769441e-14,
      -1.2490009027033011e-14,
      -1.887379141862766e-15,
      3.5388358909926865e-16,
      -2.7755575615628914e-17
    },
    0.09357943103869266,
    0.13186192555452148
  },
  {
    {
      0.5539049274549523,
      -0.05513150079759916,
      -0.00226486529446273,
      -0.00017169038790523783,
      -1.771968913242411e-05,
      -2.1044071880221837e-06,
      -2.6937613983174513e-07,
      -3.620174487267702e-08,
      -5.035860249635871e-09,
      -7.188804967972473e-10,
      -1.0471176703497065e-10,
      -1.5500455086137066e-11,
      -2.3248417080345973e-12,
      -3.5293989952833726e-13,
      -5.399847236020605e-14,
      -8.854028621385623e-15,
      -1.2975731600306517e-15,
      2.8449465006019636e-16,
      -4.5102810375396984e-17
    },
    0.13186192555452148,
    0.1510031728124359
  },
  {
    {
      0.459674826377972,
      -0.03799561202969537,
      -0.0014593552334729865,
      -0.00011667442436243125,
      -1.2328894690105674e-05,
      -1.4734787574105512e-06,
      -1.89107481913392e-07,
      -2.5449463390836424e-08,
      -3.5434242970366903e-09,
      -5.061691760177567e-10,
      -7.376567759398434e-11,
      -1.0923841692322966e-11,
      -1.6389077595047041e-12,
      -2.489328188026718e-13,
      -3.807371085073896e-14,
      -6.283168429987995e-15,
      -9.055256544598933e-16,
      3.0531133177191805e-16,
      -2.2551405187698492e-17
    },
    0.1510031728124359,
    0.16057379644139308
  },
  {
    {
      0.39512265522451157,
      -0.025874636107352526,
      -0.000981657847521144,
      -8.126867623134762e-05,
      -8.658142312312905e-06,
      -1.0370940123176353e-06,
      -1.3324415930673905e-07,
      -1.794336704219468e-08,
      -2.4994516144294376e-09,
      -3.5715819884929445e-10,
      -5.20629060696276e-11,
      -7.711369737206653e-12,
      -1.157088314052146e-12,
      -1.7580381594939354e-13,
      -2.6891683324592464e-14,
      -4.472117121068209e-15,
      -6.314393452555578e-16,
      2.636779683484747e-16,
      -1.734723475976807e-17
    },
    0.16057379644139308,
    0.1653591082558717
  },
  {
    {
      0.3511712698170222,
      -0.01762984199424723,
      -0.0006775212105428925,
      -5.710620124233304e-05,
      -6.102283057512342e-06,
      -7.316593272181648e-07,
      -9.405098579107207e-08,
      -1.2669533731163307e-08,
      -1.7652206779628088e-09,
      -2.5228212227612623e-10,
      -3.67796765265993e-11,
      -5.448200918189983e-12,
      -8.175543575461575e-13,
      -1.2424436479641088e-13,
      -1.8981344274138223e-14,
      -3.202299536653186e-15,
      -4.440892098500626e-16,
      3.191891195797325e-16,
      -3.122502256758253e-17
    },
    0.1653591082558717,
    0.167751764163111
  },
  {
    {
      0.3211597705559426,
      -0.012077354670507117,
      -0.00047368246685701096,
      -4.0264768868699535e-05,
      -4.308027496135319e-06,
      -5.167712383395695e-07,
      -6.644526837421005e-08,
      -8.952242693677226e-09,
      -1.247437774604121e-09,
      -1.7829611051456418e-10,
      -2.5995025576541764e-11,
      -3.8507357025263644e-12,
      -5.778121037192108e-13,
      -8.784639682346551e-14,
      -1.3433698597964394e-14,
      -2.373101715136272e-15,
      -3.365363543395006e-16,
      1.2836953722228372e-16,
      -2.2551405187698492e-17
    },
    0.167751764163111,
    0.16894809211673065
  },
  {
    {
      0.30054357607263954,
      -0.008327461617588218,
      -0.00033318306039584883,
      -2.8432365464115678e-05,
      -3.043795371863306e-06,
      -3.652042762300467e-07,
      -4.696313324337176e-08,
      -6.3279062684218346e-09,
      -8.818025996892853e-10,
      -1.260411150449947e-10,
      -1.8376970084554856e-11,
      -2.7223640008955385e-12,
      -4.0846839799435486e-13,
      -6.197126145579546e-14,
      -9.492406860545088e-15,
      -1.7208456881689926e-15,
      -2.393918396847994e-16,
      2.5673907444456745e-16,
      -8.673617379884035e-18
    },
    0.16894809211673065,
    0.16954625609354046
  },
  {
    {
      0.2862916733850893,
      -0.005776207912322414,
      -0.00023501417597186075,
      -2.0091122712718318e-05,
      -2.151427414113366e-06,
      -2.581648955793381e-07,
      -3.320061314066036e-08,
      -4.473697898244033e-09,
      -6.234334853916224e-10,
      -8.9112786594292e-11,
      -1.2992894954377832e-11,
      -1.9243946713931592e-12,
      -2.8871002810682e-13,
      -4.3874626154405405e-14,
      -6.7133798520302435e-15,
      -1.2906342661267445e-15,
      -1.5612511283791264e-16,
      3.0878077872387166e-16,
      1.734723475976807e-18
    },
    0.16954625609354046,
    0.16984533808194538
  },
  {
    {
      0.2763847311775531,
      -0.004026323535076086,
      -0.00016598517183828362,
      -1.4201795987777571e-05,
      -1.520984803850961e-06,
      -1.8252416650321734e-07,
      -2.3473785675659498e-08,
      -3.1630968764378986e-09,
      -4.408004100620033e-10,
      -6.300811088155722e-11,
      -9.186814503570062e-12,
      -1.3607726057074387e-12,
      -2.0409195167214733e-13,
      -3.113134749987978e-14,
      -4.758346494604382e-15,
      -9.957312752106873e-16,
      -1.682681771697503e-16,
      -2.445960101127298e-16,
      -4.2500725161431774e-17
    },
    0.16984533808194538,
    0.16999487907614785
  },
  {
    {
      0.2694673741877153,
      -0.0028173393805495644,
      -0.00011730308752334363,
      -1.004050356021198e-05,
      -1.0753911961810386e-06,
      -1.290548921391993e-07,
      -1.6597556396977242e-08,
      -2.236546208650436e-09,
      -3.1168101606582344e-10,
      -4.4551562092265407e-11,
      -6.495835119801896e-12,
      -9.626032609899582e-13,
      -1.4432725847779437e-13,
      -2.216109240560371e-14,
      -3.3948538424866115e-15,
      -1.0026701691145945e-15,
      -1.5092094240998222e-16,
      -2.7582103268031233e-16,
      -4.7704895589362195e-17
    },
    0.16999487907614785,
    0.17006964957324908
  },
  {
    {
      0.264620888720384,
      -0.001977082335539313,
      -8.2923124407136e-05,
      -7.099113954002478e-06,
      -7.603784178723816e-07,
      -9.125234246010194e-08,
      -1.1735920570424274e-08,
      -1.581440811637691e-09,
      -2.2038751899822184e-10,
      -3.1502528016758014e-11,
      -4.593121022411495e-12,
      -6.799057844508738e-13,
      -1.0198439315267649e-13,
      -1.6341095143701523e-14,
      -2.426878142891553e-15,
      -6.210310043996969e-16,
      -5.898059818321144e-17,
      4.527628272299467e-16,
      1.3010426069826053e-17
    },
    0.17006964957324908,
    0.17010703482179967
  },
  {
    {
      0.2612166158875983,
      -0.001390391541082515,
      -5.862766135761137e-05,
      -5.019621645544273e-06,
      -5.376553057923755e-07,
      -6.452400296876049e-08,
      -8.298434496858964e-09,
      -1.1182347825489103e-09,
      -1.558359847414481e-10,
      -2.227549741884438e-11,
      -3.247855109855813e-12,
      -4.816719939571001e-13,
      -7.213847574849552e-14,
      -1.1423154089307275e-14,
      -1.7156415177410622e-15,
      -5.048045315092509e-16,
      -7.632783294297951e-17,
      -9.367506770274758e-17,
      -3.209238430557093e-17
    },
    0.17010703482179967,
    0.170125727446075
  },
  {
    {
      0.2588208891336881,
      -0.0009793198302267478,
      -4.145328900037257e-05,
      -3.5493342854138465e-06,
      -3.8017496464121114e-07,
      -4.562495321824844e-08,
      -5.867838735740261e-09,
      -7.907071299800839e-10,
      -1.1019215748198086e-10,
      -1.5751122628415715e-11,
      -2.296527551459704e-12,
      -3.4033192930493783e-13,
      -5.1030360492809734e-14,
      -9.162809400109495e-15,
      -1.2975731600306517e-15,
      -4.354155924701786e-16,
      -1.3010426069826053e-16,
      -8.690964614643804e-16,
      -9.107298248878237e-17
    },
    0.170125727446075,
    0.17013507375821263
  },
  {
    {
      0.25713261899860995,
      -0.0006905565301747237,
      -2.931094851913764e-05,
      -2.509732105073978e-06,
      -2.688226164813262e-07,
      -3.226156955644932e-08,
      -4.149175109394165e-09,
      -5.591129538518036e-10,
      -7.791681362612213e-11,
      -1.1135805819129097e-11,
      -1.6232518801340134e-12,
      -2.40765740677773e-13,
      -3.687848637579094e-14,
      -5.707240235963695e-15,
      -4.597017211338539e-16,
      -1.2160411566597418e-15,
      2.203098814490545e-16,
      -1.2923689896027213e-15,
      -5.828670879282072e-16
    },
    0.17013507375821263,
    0.17013974691428146
  },
  {
    {
      0.2559417264393482,
      -0.0004873303221206906,
      -2.0725636255014904e-05,
      -1.7746393138461797e-06,
      -1.900857017015506e-07,
      -2.2812321970669402e-08,
      -2.933903902810342e-09,
      -3.9535087530129243e-10,
      -5.5095764756041277e-11,
      -7.875578661442617e-12,
      -1.1483054090932754e-12,
      -1.7132302521094545e-13,
      -2.5521251778570786e-14,
      -4.4114017994090204e-15,
      -6.29704621779581e-16,
      -3.0184188481996443e-16,
      1.8561541192951836e-16,
      2.534430998402115e-15,
      1.951563910473908e-16
    },
    0.17013974691428146,
    0.17014208349231588
  },
  {
    {
      0.2551010899484165,
      -0.0003441099749644915,
      -1.465512048574695e-05,
      -1.254856211765562e-06,
      -1.3441067886650615e-07,
      -1.613073148838684e-08,
      -2.0745819472034466e-09,
      -2.7955775715182707e-10,
      -3.895907554896105e-11,
      -5.572540692777572e-12,
      -8.118037492232943e-13,
      -1.156644224842296e-13,
      -1.7319479184152442e-14,
      -6.123573870198129e-16,
      -6.765421556309548e-17,
      1.5838025335668249e-15,
      2.5326962749261384e-16,
      1.56472057533108e-15,
      1.1622647289044608e-16
    },
    0.17014208349231588,
    0.1701432517813331
  },
  {
    {
      0.25450739757424595,
      -0.0002430797888991077,
      -1.0362693722733643e-05,
      -8.873161782604194e-07,
      -9.50426266178167e-08,
      -1.1406137241307124e-08,
      -1.4669526513660935e-09,
      -1.976803373771041e-10,
      -2.754695430096099e-11,
      -3.937652287566706e-12,
      -5.726478319312278e-13,
      -8.348009783443189e-14,
      -1.5083420623618338e-14,
      -2.8484159475539172e-15,
      9.280770596475918e-16,
      -8.326672684688674e-17,
      1.0928757898653885e-15,
      -1.9723805921856297e-15,
      -1.4632392519864368e-15
    },
    0.1701432517813331,
    0.1701438359258417
  },
  {
    {
      0.25382485353286033,
      -0.0005003082472470298,
      -0.00010070946646463866,
      -4.371453981542178e-05,
      -2.472424034292986e-05,
      -1.6101233343186055e-05,
      -1.1467857101349344e-05,
      -8.698427342242881e-06,
      -6.917894183581727e-06,
      -5.7120257463023905e-06,
      -4.864148734601281e-06,
      -4.252149622345491e-06,
      -3.803174641357071e-06,
      -3.471843289028567e-06,
      -3.2290110822119633e-06,
      -3.0556483156078673e-06,
      -2.9393703138338306e-06,
      -2.8724294233051922e-06,
      -1.425284251615512e-06
    },
    0.1701438359258417,
    0.1701444200703503
  }
}
};
const auto PR_rhoV = SuperAncillary{
{
  {
    {
      1.967544593078315e-12,
      3.286395018856467e-12,
      1.951259019068106e-12,
      8.480414006053818e-13,
      2.7638255525455143e-13,
      6.855112175816157e-14,
      1.2987941190221208e-14,
      1.8596645822713178e-15,
      1.9425469025161906e-16,
      1.3440422453482796e-17,
      4.1371292677837237e-19,
      -2.0472087922962698e-20,
      -2.3984414514408427e-21,
      -1.0726839771583605e-23,
      8.190099610752978e-24,
      1.4834694490433158e-25,
      -2.890640271333655e-26,
      -1.2318556073091862e-28,
      2.6437933681383566e-28
    },
    0.01701444200703503,
    0.021799753821513633
  },
  {
    {
      3.4427283919277603e-10,
      5.10033605480091e-10,
      2.3125713682061785e-10,
      6.97923116109859e-11,
      1.468781950176213e-11,
      2.192555878654968e-12,
      2.289969778559488e-13,
      1.563184655125417e-14,
      5.345224322369889e-16,
      -8.665670712890633e-18,
      -1.4805111905363385e-18,
      -1.2881526260516941e-20,
      3.1913825586258806e-21,
      4.0303951255023627e-23,
      -7.672322324660242e-24,
      -2.068961273341252e-25,
      -4.90734591919887e-26,
      1.348005514841665e-25,
      2.4183319910455338e-26
    },
    0.021799753821513633,
    0.026585065635992236
  },
  {
    {
      1.1651253109754791e-07,
      1.799247246805057e-07,
      8.718524363492525e-08,
      2.7723805947266338e-08,
      5.859084478366607e-09,
      7.937371730144063e-10,
      5.851825200412916e-11,
      3.5459917618408365e-13,
      -2.840460099622861e-13,
      -9.001110583611441e-15,
      1.453688138264624e-15,
      4.153752300772462e-17,
      -8.493292717205814e-18,
      1.0460073661948098e-20,
      4.4932489647971677e-20,
      -2.1592128174018962e-21,
      -1.6723265540349414e-22,
      6.45330124762073e-23,
      1.344814730252227e-23
    },
    0.026585065635992236,
    0.03615568926494944
  },
  {
    {
      3.098979961669226e-05,
      4.5977523083262384e-05,
      1.9714735728750365e-05,
      4.9340743435473964e-06,
      6.495334627828851e-07,
      2.0013065084957464e-08,
      -4.62369606718499e-09,
      -1.6742776075828521e-10,
      5.6102600385966035e-11,
      -2.1810195956609656e-13,
      -6.246733592709029e-13,
      4.045717710917202e-14,
      3.97102187587516e-15,
      -7.21043417103846e-16,
      1.9805582453051167e-17,
      4.946315069940491e-18,
      -6.192935810062008e-19,
      5.64732747802203e-20,
      1.2579762755706445e-20
    },
    0.03615568926494944,
    0.05529693652286385
  },
  {
    {
      0.0005902854802776641,
      0.0006119857879876187,
      0.000137226493734894,
      1.3534511317666194e-05,
      2.5757985750152563e-07,
      -3.0381140235968494e-08,
      1.4660600077037917e-09,
      1.4064389556322658e-10,
      -1.7247649444289924e-11,
      4.829482286631236e-13,
      7.904912957394129e-14,
      -9.965950794102972e-15,
      3.496086138144762e-16,
      3.581424707580633e-17,
      -5.7462715141731735e-18,
      5.929230630780102e-20,
      6.776263578034403e-21,
      5.166900978251232e-19,
      1.1858461261560205e-20
    },
    0.05529693652286385,
    0.07443818378077825
  },
  {
    {
      0.0034551514866857118,
      0.002398865145025046,
      0.0003121082643740536,
      1.5218984619397858e-05,
      9.690936261096291e-08,
      1.2913919961681083e-08,
      1.4096536882115922e-09,
      -3.3114702454217944e-11,
      2.594155481248396e-12,
      1.6231949866250123e-13,
      -1.3332379904945624e-14,
      1.0818033751788803e-15,
      -1.3362791775883842e-17,
      -2.2497195079074217e-18,
      -9.486769009248164e-20,
      -1.2197274440461925e-18,
      -2.303929616531697e-19,
      3.6591823321385775e-18,
      3.9979955110402976e-19
    },
    0.07443818378077825,
    0.09357943103869266
  },
  {
    {
      0.020288628115803765,
      0.016610099415744563,
      0.00272192651625002,
      0.0002394799179501481,
      2.3315294231581688e-05,
      3.169727025417485e-06,
      3.793923939678003e-07,
      5.156601162595916e-08,
      7.275292876224738e-09,
      1.0148922060279446e-09,
      1.5007760583319457e-10,
      2.202715483970441e-11,
      3.3070396464676688e-12,
      5.021763712330374e-13,
      7.666279720416891e-14,
      1.1847890290378471e-14,
      1.843794214528849e-15,
      3.1604493327952454e-16,
      4.651227319962814e-17
    },
    0.09357943103869266,
    0.13186192555452148
  },
  {
    {
      0.059633242372117805,
      0.021523786193835045,
      0.001926244019335927,
      0.00016436263238290855,
      1.7587241022322224e-05,
      2.103887459802955e-06,
      2.692872191701126e-07,
      3.6202589339100893e-08,
      5.035822049290206e-09,
      7.188798449749012e-10,
      1.0471208969353718e-10,
      1.5500948181285112e-11,
      2.3252866646061854e-12,
      3.5267275211303684e-13,
      5.3994569232385103e-14,
      8.295447662121092e-15,
      1.2962721174236691e-15,
      2.658463726934457e-16,
      3.230922474006803e-17
    },
    0.13186192555452148,
    0.1510031728124359
  },
  {
    {
      0.10209461339960167,
      0.020078909930632285,
      0.0013564703941337276,
      0.00011554892995170075,
      1.2319437775788929e-05,
      1.4734155808103153e-06,
      1.8910621422518065e-07,
      2.5449462346532892e-08,
      3.5434243221901807e-09,
      5.06168969585663e-10,
      7.376593953722921e-11,
      1.0924287516256292e-11,
      1.6393041438189648e-12,
      2.4869169223951104e-13,
      3.807978238290488e-14,
      5.776629175002768e-15,
      9.063930161978817e-16,
      2.411265631607762e-16,
      2.0383000842727483e-17
    },
    0.1510031728124359,
    0.16057379644139308
  },
  {
    {
      0.1393729010835231,
      0.016593251298797226,
      0.0009532795808090732,
      8.111304465426971e-05,
      8.657481578122454e-06,
      1.037091311451195e-06,
      1.332441384432198e-07,
      1.794336703699051e-08,
      2.4994517905038705e-09,
      3.571580193054147e-10,
      5.2063129848956e-11,
      7.711812959054765e-12,
      1.157456075429053e-12,
      1.7556355674797075e-13,
      2.688300970721258e-14,
      3.9968028886505635e-15,
      6.323067069935462e-16,
      2.2985086056692694e-16,
      1.5178830414797062e-17
    },
    0.16057379644139308,
    0.1653591082558717
  },
  {
    {
      0.1692940156782449,
      0.012902057653797903,
      0.0006700678810102369,
      5.7085725788112024e-05,
      6.102239057614867e-06,
      7.316592307484576e-07,
      9.405098538861623e-08,
      1.266953380749114e-08,
      1.7652208635782207e-09,
      2.5228192451764997e-10,
      3.677991071426856e-11,
      5.448613782377265e-12,
      8.179099758587327e-13,
      1.240327285323417e-13,
      1.8984813721090177e-14,
      2.7564756033271465e-15,
      4.215378046623641e-16,
      2.0122792321330962e-16,
      2.688821387764051e-17
    },
    0.1653591082558717,
    0.167751764163111
  },
  {
    {
      0.1921856564900402,
      0.00969084706538421,
      0.0004717724978338892,
      4.026214237729843e-05,
      4.308024652975584e-06,
      5.167712350106352e-07,
      6.644526831522946e-08,
      8.952242742249483e-09,
      1.2474379341986808e-09,
      1.7829594398111048e-10,
      2.599524762114669e-11,
      3.851164179224931e-12,
      5.781538442439782e-13,
      8.762435221854048e-14,
      1.3431963874488417e-14,
      1.915134717478395e-15,
      3.2612801348363973e-16,
      3.2439329000766293e-16,
      2.3418766925686896e-17
    },
    0.167751764163111,
    0.16894809211673065
  },
  {
    {
      0.2092148975305195,
      0.007128445272373621,
      0.0003326996220906972,
      2.8432032853142192e-05,
      3.043795191084303e-06,
      3.652042760392271e-07,
      4.6963133177452265e-08,
      6.327906323932986e-09,
      8.818027766310799e-10,
      1.2604095198098797e-10,
      1.8377181720818925e-11,
      2.7227942123175808e-12,
      4.087806482200307e-13,
      6.178738076734192e-14,
      9.485467966641181e-15,
      1.2732870313669764e-15,
      2.237793284010081e-16,
      2.3418766925686896e-16,
      7.806255641895632e-18
    },
    0.16894809211673065,
    0.16954625609354046
  },
  {
    {
      0.22166645905817658,
      0.005175245300316345,
      0.0002348925657099949,
      2.0091080864132957e-05,
      2.1514274026867425e-06,
      2.581648954613769e-07,
      3.3200613073006147e-08,
      4.473697952020461e-09,
      6.234336501903526e-10,
      8.911260965249745e-11,
      1.299311526425928e-11,
      1.9248335564325814e-12,
      2.8905870752549134e-13,
      4.368554129552393e-14,
      6.711645128554267e-15,
      8.378714388967978e-16,
      1.5265566588595902e-16,
      1.8214596497756474e-16,
      -2.6020852139652106e-18
    },
    0.16954625609354046,
    0.16984533808194538
  },
  {
    {
      0.23067150083801227,
      0.003725476879880729,
      0.00016595467496789695,
      1.4201790739628434e-05,
      1.5209848030911521e-06,
      1.8252416639913394e-07,
      2.3473785595862218e-08,
      3.1630968937851334e-09,
      4.4080058700379787e-10,
      6.3007946082827e-11,
      9.187020935663703e-12,
      1.3611941435121011e-12,
      2.0445450887862648e-13,
      3.09127723419067e-14,
      4.751407600700475e-15,
      5.238864897449957e-16,
      1.491862189340054e-16,
      7.042977312465837e-16,
      4.163336342344337e-17
    },
    0.16984533808194538,
    0.16999487907614785
  },
  {
    {
      0.237137473422725,
      0.0026668244973612407,
      0.00011729545148811069,
      1.004050290311434e-05,
      1.0753911961168539e-06,
      1.290548920437895e-07,
      1.659755634146609e-08,
      2.2365462572226935e-09,
      3.116811773951067e-10,
      4.455141290604647e-11,
      6.496051960236393e-12,
      9.630022473894329e-13,
      1.446689990025618e-13,
      2.1937313077202703e-14,
      3.3931191190106347e-15,
      5.412337245047638e-16,
      1.3357370765021415e-16,
      7.268491364342822e-16,
      4.597017211338539e-17
    },
    0.16999487907614785,
    0.17006964957324908
  },
  {
    {
      0.24175815792274802,
      0.0019018019777060826,
      8.292121391923792e-05,
      7.0991138717956676e-06,
      7.603784178359524e-07,
      9.125234237510049e-08,
      1.17359205131784e-08,
      1.5814408688835657e-09,
      2.2038768553167554e-10,
      3.150234934023999e-11,
      4.593327454505136e-12,
      6.803307917024881e-13,
      1.023348072948238e-13,
      1.6110376921396607e-14,
      2.4338170367954604e-15,
      1.5265566588595902e-16,
      3.122502256758253e-17,
      3.642919299551295e-17,
      -1.3877787807814457e-17
    },
    0.17006964957324908,
    0.17010703482179967
  },
  {
    {
      0.24504950305352938,
      0.001352745629684961,
      5.862718355055703e-05,
      5.019621635323282e-06,
      5.376553057160477e-07,
      6.452400284732984e-08,
      8.298434437878366e-09,
      1.118234846733679e-09,
      1.5583615474434875e-10,
      2.227533782428459e-11,
      3.248068480843358e-12,
      4.820623067391949e-13,
      7.248195099673893e-14,
      1.1211517825238104e-14,
      1.7034984534092246e-15,
      6.418476861114186e-17,
      6.591949208711867e-17,
      5.533767888366015e-16,
      2.949029909160572e-17
    },
    0.17010703482179967,
    0.170125727446075
  },
  {
    {
      0.2473887591484516,
      0.0009604954409799197,
      4.1453169525423536e-05,
      3.5493342841735193e-06,
      3.8017496460478195e-07,
      4.5624953119369205e-08,
      5.867838671555492e-09,
      7.907071629398299e-10,
      1.1019232575015803e-10,
      1.5750949156068117e-11,
      2.296747861341153e-12,
      3.407586712800281e-13,
      5.137903991148107e-14,
      8.93209117780458e-15,
      1.2923689896027213e-15,
      -6.938893903907228e-18,
      1.1796119636642288e-16,
      1.3270634591222574e-15,
      9.627715291671279e-17
    },
    0.170125727446075,
    0.17013507375821263
  },
  {
    {
      0.2490487922516141,
      0.0006811439771107868,
      2.9310918647461326e-05,
      2.509732104961221e-06,
      2.6882261645530536e-07,
      3.226156947144787e-08,
      4.149175027862162e-09,
      5.591129972198905e-10,
      7.791697148595844e-11,
      1.113566183708059e-11,
      1.6234808636328424e-12,
      2.411647270772477e-13,
      3.721675745360642e-14,
      5.479991460610734e-15,
      4.753142324176451e-16,
      7.563394355258879e-16,
      -2.2724877535296173e-16,
      1.7572748811645056e-15,
      5.733261088103347e-16
    },
    0.17013507375821263,
    0.17013974691428146
  },
  {
    {
      0.2502255658692596,
      0.0004826239559720131,
      2.0725628786679925e-05,
      1.774639313863527e-06,
      1.9008570164603944e-07,
      2.2812321890872123e-08,
      2.933903833421403e-09,
      3.953508978526976e-10,
      5.5095907003366307e-11,
      7.875412127988923e-12,
      1.1485309231451524e-12,
      1.7173241995127597e-13,
      2.5854318685958333e-14,
      4.182418300580082e-15,
      6.245004513516506e-16,
      -1.4224732503009818e-16,
      -1.9255430583342559e-16,
      -2.048708425128609e-15,
      -2.0036056147532122e-16
    },
    0.17013974691428146,
    0.17014208349231588
  },
  {
    {
      0.25105914278296176,
      0.00034175676948506976,
      1.4655118618528765e-05,
      1.2548562118210732e-06,
      1.3441067884568947e-07,
      1.6130731377364538e-08,
      2.0745818743450606e-09,
      2.7955778143795573e-10,
      3.895922473517999e-11,
      5.572348138471739e-12,
      8.120119160404116e-13,
      1.1612238948188747e-13,
      1.7666423879347803e-14,
      4.0072112295064244e-16,
      5.724587470723463e-17,
      -2.0209528495129803e-15,
      -2.723515857283587e-16,
      -1.0651202142497596e-15,
      -1.1188966420050406e-16
    },
    0.17014208349231588,
    0.1701432517813331
  },
  {
    {
      0.25164930534191154,
      0.0002419031805581083,
      1.0362693255834554e-05,
      8.87316178307257e-07,
      9.50426266022042e-08,
      1.140613715977512e-08,
      1.4669525750382606e-09,
      1.976803529896154e-10,
      2.7547124303861636e-11,
      3.937477080495633e-12,
      5.728750807065808e-13,
      8.391204397995011e-14,
      1.5444243106621514e-14,
      2.6229018956769323e-15,
      -9.194034422677078e-16,
      -3.5735303605122226e-16,
      -1.0928757898653885e-15,
      2.4494295480792516e-15,
      1.4641066137244252e-15
    },
    0.1701432517813331,
    0.1701438359258417
  },
  {
    {
      0.2523294961628807,
      0.0004991316351717207,
      0.000100709465997776,
      4.37145398154426e-05,
      2.472424034286394e-05,
      1.6101233343113197e-05,
      1.1467857101271281e-05,
      8.698427342246351e-06,
      6.917894183770812e-06,
      5.712025746135857e-06,
      4.864148734845877e-06,
      4.252149622756621e-06,
      3.80317464170922e-06,
      3.4718432888186657e-06,
      3.229011082201555e-06,
      3.055648315165513e-06,
      2.9393703138234223e-06,
      2.872429423783976e-06,
      1.4252842516137773e-06
    },
    0.1701438359258417,
    0.1701444200703503
  }
}
};
 
static inline double supercubic(int EOS, int prop, double Ttilde){
    switch(EOS){
        case 0:{
            switch(prop){
                case 100:
                    return vdW_p.y(Ttilde);
                case 101:
                    return vdW_rhoL.y(Ttilde);
                case 102:
                    return vdW_rhoV.y(Ttilde);
                default:
                    return -2;
            }
        }
        case 1:{
            switch(prop){
                case 100:
                    return SRK_p.y(Ttilde);
                case 101:
                    return SRK_rhoL.y(Ttilde);
                case 102:
                    return SRK_rhoV.y(Ttilde);
                default:
                    return -2;
            }
        }
        case 2:{
            switch(prop){
                case 100:
                    return PR_p.y(Ttilde);
                case 101:
                    return PR_rhoL.y(Ttilde);
                case 102:
                    return PR_rhoV.y(Ttilde);
                default:
                    return -2;
            }
        }
        default:
            return -1;
    }
}
 
const int VDW_CODE = 0, SRK_CODE = 1, PR_CODE = 2, UNKNOWN_CODE = -1;
const int P_CODE = 100, RHOL_CODE = 101, RHOV_CODE = 102;
 
}; // namespace CubicSuperAncillary
 
}; // namespace teqp