comparison DEPENDENCIES/generic/include/boost/math/special_functions/expint.hpp @ 101:c530137014c0

Update Boost headers (1.58.0)

 #include <boost/math/tools/precision.hpp>
 #include <boost/math/tools/promotion.hpp>
 #include <boost/math/tools/fraction.hpp>
 #include <boost/math/tools/series.hpp>
 #include <boost/math/policies/error_handling.hpp>
+#include <boost/math/special_functions/math_fwd.hpp>
 #include <boost/math/special_functions/digamma.hpp>
 #include <boost/math/special_functions/log1p.hpp>
 #include <boost/math/special_functions/pow.hpp>
 
 namespace boost{ namespace math{

          BOOST_MATH_BIG_CONSTANT(T, 53, -0.0368031736257943745142),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.00399167106081113256961),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.000111507792921197858394)
       };
       static const T Q[6] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 53, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 53, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.37091387659397013215),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.056770677104207528384),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.00427347600017103698101),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.000131049900798434683324),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.528611029520217142048e-6)

          BOOST_MATH_BIG_CONSTANT(T, 53, -54844.4587226402067411),
          BOOST_MATH_BIG_CONSTANT(T, 53, -14751.4895786128450662),
          BOOST_MATH_BIG_CONSTANT(T, 53, -1185.45720315201027667)
       };
       static const T Q[12] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 53, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 53, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 53, 45.3058660811801465927),
          BOOST_MATH_BIG_CONSTANT(T, 53, 809.193214954550328455),
          BOOST_MATH_BIG_CONSTANT(T, 53, 7417.37624454689546708),
          BOOST_MATH_BIG_CONSTANT(T, 53, 38129.5594484818471461),
          BOOST_MATH_BIG_CONSTANT(T, 53, 113057.05869159631492),

          BOOST_MATH_BIG_CONSTANT(T, 64, -0.0237624819878732642231),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.00259113319641673986276),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.30853660894346057053e-4)
       };
       static const T Q[7] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 64, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 64, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.317978365797784100273),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.0393622602554758722511),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.00204062029115966323229),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.732512107100088047854e-5),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.202872781770207871975e-5),

          BOOST_MATH_BIG_CONSTANT(T, 64, -584705.537139793925189),
          BOOST_MATH_BIG_CONSTANT(T, 64, -65815.2605361889477244),
          BOOST_MATH_BIG_CONSTANT(T, 64, -2038.82870680427258038)
       };
       static const T Q[14] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 64, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 64, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 64, 64.1517806091379399478),
          BOOST_MATH_BIG_CONSTANT(T, 64, 1690.76044393722763785),
          BOOST_MATH_BIG_CONSTANT(T, 64, 24035.9534033068949426),
          BOOST_MATH_BIG_CONSTANT(T, 64, 203679.998633572361706),
          BOOST_MATH_BIG_CONSTANT(T, 64, 1074661.58459976978285),

          BOOST_MATH_BIG_CONSTANT(T, 113, -0.136528159460768830763009294683628406e-5),
          BOOST_MATH_BIG_CONSTANT(T, 113, -0.346839106212658259681029388908658618e-7),
          BOOST_MATH_BIG_CONSTANT(T, 113, -0.340500302777838063940402160594523429e-9)
       };
       static const T Q[10] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 113, 0.426568827778942588160423015589537302),
          BOOST_MATH_BIG_CONSTANT(T, 113, 0.0841384046470893490592450881447510148),
          BOOST_MATH_BIG_CONSTANT(T, 113, 0.0100557215850668029618957359471132995),
          BOOST_MATH_BIG_CONSTANT(T, 113, 0.000799334870474627021737357294799839363),
          BOOST_MATH_BIG_CONSTANT(T, 113, 0.434452090903862735242423068552687688e-4),

          BOOST_MATH_BIG_CONSTANT(T, 113, -2226.64251774578542836725386936102339),
          BOOST_MATH_BIG_CONSTANT(T, 113, -109.009437301400845902228611986479816),
          BOOST_MATH_BIG_CONSTANT(T, 113, -1.51492042209561411434644938098833499)
       };
       static const T Q[16] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 113, 46.734521442032505570517810766704587),
          BOOST_MATH_BIG_CONSTANT(T, 113, 908.694714348462269000247450058595655),
          BOOST_MATH_BIG_CONSTANT(T, 113, 9701.76053033673927362784882748513195),
          BOOST_MATH_BIG_CONSTANT(T, 113, 63254.2815292641314236625196594947774),
          BOOST_MATH_BIG_CONSTANT(T, 113, 265115.641285880437335106541757711092),

          BOOST_MATH_BIG_CONSTANT(T, 113, -5737221535080.30569711574295785864903),
          BOOST_MATH_BIG_CONSTANT(T, 113, -1077042281708.42654526404581272546244),
          BOOST_MATH_BIG_CONSTANT(T, 113, -68028222642.1941480871395695677675137)
       };
       static const T Q[20] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 113, 168.542326331163836642960118190147311),
          BOOST_MATH_BIG_CONSTANT(T, 113, 12535.7237814586576783518249115343619),
          BOOST_MATH_BIG_CONSTANT(T, 113, 544891.263372016404143120911148640627),
          BOOST_MATH_BIG_CONSTANT(T, 113, 15454474.7241010258634446523045237762),
          BOOST_MATH_BIG_CONSTANT(T, 113, 302495899.896629522673410325891717381),

          BOOST_MATH_BIG_CONSTANT(T, 53, 0.000116419523609765200999),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.798296365679269702435e-5),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.2777056254402008721e-6)
       };
       static const T Q[8] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 53, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 53, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 53, -1.17090412365413911947),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.62215109846016746276),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.195114782069495403315),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.0391523431392967238166),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.00504800158663705747345),

          / tools::evaluate_polynomial(Q, t);
       t = (z - r1) - r2;
       result *= t;
       if(fabs(t) < 0.1)
       {
-         result += boost::math::log1p(t / r);
+         result += boost::math::log1p(t / r, pol);
       }
       else
       {
          result += log(z / r);
       }

          BOOST_MATH_BIG_CONSTANT(T, 53, -0.000374885917942100256775),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.554086272024881826253e-4),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.396487648924804510056e-5)
       };
       static const T Q[8] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 53, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 53, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.744625566823272107711),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.329061095011767059236),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.100128624977313872323),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.0223851099128506347278),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.00365334190742316650106),

          BOOST_MATH_BIG_CONSTANT(T, 53, -0.00155941947035972031334),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.000209750022660200888349),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.138652200349182596186e-4)
       };
       static const T Q[9] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 53, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 53, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 53, 1.97017214039061194971),
          BOOST_MATH_BIG_CONSTANT(T, 53, 1.86232465043073157508),
          BOOST_MATH_BIG_CONSTANT(T, 53, 1.09601437090337519977),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.438873285773088870812),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.122537731979686102756),

          BOOST_MATH_BIG_CONSTANT(T, 53, -0.00207592267812291726961),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.000192178045857733706044),
          BOOST_MATH_BIG_CONSTANT(T, 53, -0.113161784705911400295e-9)
       };
       static const T Q[9] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 53, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 53, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 53, 2.84354408840148561131),
          BOOST_MATH_BIG_CONSTANT(T, 53, 3.6599610090072393012),
          BOOST_MATH_BIG_CONSTANT(T, 53, 2.75088464344293083595),
          BOOST_MATH_BIG_CONSTANT(T, 53, 1.2985244073998398643),
          BOOST_MATH_BIG_CONSTANT(T, 53, 0.383213198510794507409),

          BOOST_MATH_BIG_CONSTANT(T, 53, -2516.35323679844256203),
          BOOST_MATH_BIG_CONSTANT(T, 53, 18932.0850014925993025),
          BOOST_MATH_BIG_CONSTANT(T, 53, -38703.1431362056714134)
       };
       static const T Q[7] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 53, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 53, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 53, 61.9733592849439884145),
          BOOST_MATH_BIG_CONSTANT(T, 53, -2354.56211323420194283),
          BOOST_MATH_BIG_CONSTANT(T, 53, 22329.1459489893079041),
          BOOST_MATH_BIG_CONSTANT(T, 53, -70126.245140396567133),
          BOOST_MATH_BIG_CONSTANT(T, 53, 54738.2833147775537106),

          BOOST_MATH_BIG_CONSTANT(T, 64, 0.11293331317982763165e-4),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.629499283139417444244e-6),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.177833045143692498221e-7)
       };
       static const T Q[9] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 64, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 64, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 64, -1.20352377969742325748),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.66707904942606479811),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.223014531629140771914),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.0493340022262908008636),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.00741934273050807310677),

          / tools::evaluate_polynomial(Q, t);
       t = (z - r1) - r2;
       result *= t;
       if(fabs(t) < 0.1)
       {
-         result += boost::math::log1p(t / r);
+         result += boost::math::log1p(t / r, pol);
       }
       else
       {
          result += log(z / r);
       }

          BOOST_MATH_BIG_CONSTANT(T, 64, -0.877310384591205930343e-4),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.623067256376494930067e-5),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.377246883283337141444e-6)
       };
       static const T Q[10] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 64, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 64, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 64, 1.08073635708902053767),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.553681133533942532909),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.176763647137553797451),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.0387891748253869928121),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.0060603004848394727017),

          BOOST_MATH_BIG_CONSTANT(T, 64, -0.000392477245911296982776),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.44424044184395578775e-4),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.252788029251437017959e-5)
       };
       static const T Q[10] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 64, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 64, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 64, 2.00323265503572414261),
          BOOST_MATH_BIG_CONSTANT(T, 64, 1.94688958187256383178),
          BOOST_MATH_BIG_CONSTANT(T, 64, 1.19733638134417472296),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.513137726038353385661),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.159135395578007264547),

          BOOST_MATH_BIG_CONSTANT(T, 64, -0.504623302166487346677e-9),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.712662196671896837736e-10),
          BOOST_MATH_BIG_CONSTANT(T, 64, -0.533769629702262072175e-11)
       };
       static const T Q[9] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 64, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 64, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 64, 3.13286733695729715455),
          BOOST_MATH_BIG_CONSTANT(T, 64, 4.49281223045653491929),
          BOOST_MATH_BIG_CONSTANT(T, 64, 3.84900294427622911374),
          BOOST_MATH_BIG_CONSTANT(T, 64, 2.15205199043580378211),
          BOOST_MATH_BIG_CONSTANT(T, 64, 0.802912186540269232424),

          BOOST_MATH_BIG_CONSTANT(T, 64, 53732298.8764767916542),
          BOOST_MATH_BIG_CONSTANT(T, 64, -160695051.957997452509),
          BOOST_MATH_BIG_CONSTANT(T, 64, 137839271.592778020028)
       };
       static const T Q[9] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 64, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 64, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 64, 27.2103343964943718802),
          BOOST_MATH_BIG_CONSTANT(T, 64, -8785.48528692879413676),
          BOOST_MATH_BIG_CONSTANT(T, 64, 397530.290000322626766),
          BOOST_MATH_BIG_CONSTANT(T, 64, -7356441.34957799368252),
          BOOST_MATH_BIG_CONSTANT(T, 64, 63050914.5343400957524),

       result += z;
    }
    return result;
 }
 
-template <class T>
-void expint_i_imp_113a(T& result, const T& z)
+template <class T, class Policy>
+void expint_i_imp_113a(T& result, const T& z, const Policy& pol)
 {
    BOOST_MATH_STD_USING
    // Maximum Deviation Found:                     1.230e-36
    // Expected Error Term:                         -1.230e-36
    // Max Error found at long double precision =   Poly: 4.355299e-34   Cheb: 7.512581e-34

       BOOST_MATH_BIG_CONSTANT(T, 113, 0.154646053060262871360159325115980023e-9),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.143971277122049197323415503594302307e-11),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.306243138978114692252817805327426657e-13)
    };
    static const T Q[15] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, -1.40178870313943798705491944989231793),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.943810968269701047641218856758605284),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.405026631534345064600850391026113165),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.123924153524614086482627660399122762),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.0286364505373369439591132549624317707),

       / tools::evaluate_polynomial(Q, t);
    t = ((z - r1) - r2) - r3;
    result *= t;
    if(fabs(t) < 0.1)
    {
-      result += boost::math::log1p(t / r);
+      result += boost::math::log1p(t / r, pol);
    }
    else
    {
       result += log(z / r);
    }

       BOOST_MATH_BIG_CONSTANT(T, 113, -0.635187358949437991465353268374523944e-9),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.203015132965870311935118337194860863e-10),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.384276705503357655108096065452950822e-12)
    };
    static const T Q[15] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.58784732785354597996617046880946257),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.18550755302279446339364262338114098),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.55598993549661368604527040349702836),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.184290888380564236919107835030984453),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.0459658051803613282360464632326866113),

       BOOST_MATH_BIG_CONSTANT(T, 113, -0.502538501472133913417609379765434153e-10),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.326283053716799774936661568391296584e-13),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.869226483473172853557775877908693647e-15)
    };
    static const T Q[15] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, 2.23227220874479061894038229141871087),
       BOOST_MATH_BIG_CONSTANT(T, 113, 2.40221000361027971895657505660959863),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.65476320985936174728238416007084214),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.816828602963895720369875535001248227),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.306337922909446903672123418670921066),

       BOOST_MATH_BIG_CONSTANT(T, 113, -0.133291168587253145439184028259772437e-7),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.595500337089495614285777067722823397e-9),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.133141358866324100955927979606981328e-10)
    };
    static const T Q[14] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.72490783907582654629537013560044682),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.44524329516800613088375685659759765),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.778241785539308257585068744978050181),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.300520486589206605184097270225725584),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.0879346899691339661394537806057953957),

       BOOST_MATH_BIG_CONSTANT(T, 113, -0.565697065670893984610852937110819467e-6),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.350233957364028523971768887437839573e-7),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.105428907085424234504608142258423505e-8)
    };
    static const T Q[16] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, 3.17261315255467581204685605414005525),
       BOOST_MATH_BIG_CONSTANT(T, 113, 4.85267952971640525245338392887217426),
       BOOST_MATH_BIG_CONSTANT(T, 113, 4.74341914912439861451492872946725151),
       BOOST_MATH_BIG_CONSTANT(T, 113, 3.31108463283559911602405970817931801),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.74657006336994649386607925179848899),

       BOOST_MATH_BIG_CONSTANT(T, 113, -0.794712243338068631557849449519994144e-6),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.434084023639508143975983454830954835e-7),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.107839681938752337160494412638656696e-8)
    };
    static const T Q[12] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, 2.09913805456661084097134805151524958),
       BOOST_MATH_BIG_CONSTANT(T, 113, 2.07041755535439919593503171320431849),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.26406517226052371320416108604874734),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.529689923703770353961553223973435569),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.159578150879536711042269658656115746),

       BOOST_MATH_BIG_CONSTANT(T, 113, -0.29462181955852860250359064291292577e-4),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.21972450610957417963227028788460299e-5),
       BOOST_MATH_BIG_CONSTANT(T, 113, -0.720558173805289167524715527536874694e-7)
    };
    static const T Q[11] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, 2.95918362458402597039366979529287095),
       BOOST_MATH_BIG_CONSTANT(T, 113, 3.96472247520659077944638411856748924),
       BOOST_MATH_BIG_CONSTANT(T, 113, 3.15563251550528513747923714884142131),
       BOOST_MATH_BIG_CONSTANT(T, 113, 1.64674612007093983894215359287448334),
       BOOST_MATH_BIG_CONSTANT(T, 113, 0.58695020129846594405856226787156424),

       BOOST_MATH_BIG_CONSTANT(T, 113, -3661928.69330208734947103004900349266),
       BOOST_MATH_BIG_CONSTANT(T, 113, 8570600.041606912735872059184527855),
       BOOST_MATH_BIG_CONSTANT(T, 113, -6758379.93672362080947905580906028645)
    };
    static const T Q[10] = {    
-      BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+      BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
       BOOST_MATH_BIG_CONSTANT(T, 113, -99.4868026047611434569541483506091713),
       BOOST_MATH_BIG_CONSTANT(T, 113, 3879.67753690517114249705089803055473),
       BOOST_MATH_BIG_CONSTANT(T, 113, -76495.82413252517165830203774900806),
       BOOST_MATH_BIG_CONSTANT(T, 113, 820773.726408311894342553758526282667),
       BOOST_MATH_BIG_CONSTANT(T, 113, -4803087.64956923577571031564909646579),

 
    T result;
 
    if(z <= 6)
    {
-      expint_i_imp_113a(result, z);
+      expint_i_imp_113a(result, z, pol);
    }
    else if (z <= 10)
    {
       expint_i_113b(result, z);
    }

          BOOST_MATH_BIG_CONSTANT(T, 113, 93632.6567462673524739954389166550069),
          BOOST_MATH_BIG_CONSTANT(T, 113, -225025.189335919133214440347510936787),
          BOOST_MATH_BIG_CONSTANT(T, 113, 175864.614717440010942804684741336853)
       };
       static const T Q[9] = {    
-         BOOST_MATH_BIG_CONSTANT(T, 113, 1),
+         BOOST_MATH_BIG_CONSTANT(T, 113, 1.0),
          BOOST_MATH_BIG_CONSTANT(T, 113, -65.6998869881600212224652719706425129),
          BOOST_MATH_BIG_CONSTANT(T, 113, 1642.73850032324014781607859416890077),
          BOOST_MATH_BIG_CONSTANT(T, 113, -19937.2610222467322481947237312818575),
          BOOST_MATH_BIG_CONSTANT(T, 113, 124136.267326632742667972126625064538),
          BOOST_MATH_BIG_CONSTANT(T, 113, -384614.251466704550678760562965502293),