@@ -562,6 +562,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
562562 // Maximum Deviation Found: 2.852e-18
563563 // Expected Error Term: 2.852e-18
564564 // Max Error found at double precision = Poly: 2.636335e-16 Cheb: 4.187027e-16
565+ // LCOV_EXCL_START
565566 BOOST_MATH_STATIC const T P[10 ] = {
566567 BOOST_MATH_BIG_CONSTANT (T, 53 , 2.98677224343598593013 ),
567568 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.356343618769377415068 ),
@@ -590,6 +591,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
590591 BOOST_MATH_STATIC_LOCAL_VARIABLE const T r1 = static_cast <T>(c1 / c2);
591592 BOOST_MATH_STATIC_LOCAL_VARIABLE const T r2 = BOOST_MATH_BIG_CONSTANT (T, 53 , 0.131401834143860282009280387409357165515556574352422001206362e-16 );
592593 BOOST_MATH_STATIC_LOCAL_VARIABLE const T r = static_cast <T>(BOOST_MATH_BIG_CONSTANT (T, 53 , 0.372507410781366634461991866580119133535689497771654051555657435242200120636201854384926049951548942392 ));
594+ // LCOV_EXCL_STOP
593595 T t = (z / 3 ) - 1 ;
594596 result = tools::evaluate_polynomial (P, t)
595597 / tools::evaluate_polynomial (Q, t);
@@ -609,6 +611,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
609611 // Maximum Deviation Found: 6.546e-17
610612 // Expected Error Term: 6.546e-17
611613 // Max Error found at double precision = Poly: 6.890169e-17 Cheb: 6.772128e-17
614+ // LCOV_EXCL_START
612615 BOOST_MATH_STATIC_LOCAL_VARIABLE const T Y = 1 .158985137939453125F ;
613616 BOOST_MATH_STATIC const T P[8 ] = {
614617 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.00139324086199402804173 ),
@@ -630,6 +633,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
630633 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.000402453408512476836472 ),
631634 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.263649630720255691787e-4 )
632635 };
636+ // LCOV_EXCL_STOP
633637 T t = z / 2 - 4 ;
634638 result = Y + tools::evaluate_polynomial (P, t)
635639 / tools::evaluate_polynomial (Q, t);
@@ -641,7 +645,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
641645 // Maximum Deviation Found: 1.843e-17
642646 // Expected Error Term: -1.842e-17
643647 // Max Error found at double precision = Poly: 4.375868e-17 Cheb: 5.860967e-17
644-
648+ // LCOV_EXCL_START
645649 BOOST_MATH_STATIC_LOCAL_VARIABLE const T Y = 1 .0869731903076171875F ;
646650 BOOST_MATH_STATIC const T P[9 ] = {
647651 BOOST_MATH_BIG_CONSTANT (T, 53 , -0.00893891094356945667451 ),
@@ -665,6 +669,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
665669 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.00278170769163303669021 ),
666670 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.000159150281166108755531 )
667671 };
672+ // LCOV_EXCL_STOP
668673 T t = z / 5 - 3 ;
669674 result = Y + tools::evaluate_polynomial (P, t)
670675 / tools::evaluate_polynomial (Q, t);
@@ -676,8 +681,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
676681 // Maximum Deviation Found: 5.102e-18
677682 // Expected Error Term: 5.101e-18
678683 // Max Error found at double precision = Poly: 1.441088e-16 Cheb: 1.864792e-16
679-
680-
684+ // LCOV_EXCL_START
681685 BOOST_MATH_STATIC_LOCAL_VARIABLE const T Y = 1 .03937530517578125F ;
682686 BOOST_MATH_STATIC const T P[9 ] = {
683687 BOOST_MATH_BIG_CONSTANT (T, 53 , -0.00356165148914447597995 ),
@@ -700,6 +704,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
700704 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.0651165455496281337831 ),
701705 BOOST_MATH_BIG_CONSTANT (T, 53 , 0.00488071077519227853585 )
702706 };
707+ // LCOV_EXCL_STOP
703708 T t = z / 10 - 3 ;
704709 result = Y + tools::evaluate_polynomial (P, t)
705710 / tools::evaluate_polynomial (Q, t);
@@ -709,6 +714,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
709714 else
710715 {
711716 // Max Error found at double precision = 3.381886e-17
717+ // LCOV_EXCL_START
712718 BOOST_MATH_STATIC_LOCAL_VARIABLE const T exp40 = static_cast <T>(BOOST_MATH_BIG_CONSTANT (T, 53 , 2.35385266837019985407899910749034804508871617254555467236651e17 ));
713719 BOOST_MATH_STATIC_LOCAL_VARIABLE const T Y= 1 .013065338134765625F ;
714720 BOOST_MATH_STATIC const T P[6 ] = {
@@ -728,6 +734,7 @@ BOOST_MATH_GPU_ENABLED T expint_i_imp(T z, const Policy& pol, const boost::math:
728734 BOOST_MATH_BIG_CONSTANT (T, 53 , 54738.2833147775537106 ),
729735 BOOST_MATH_BIG_CONSTANT (T, 53 , 8297.16296356518409347 )
730736 };
737+ // LCOV_EXCL_STOP
731738 T t = 1 / z;
732739 result = Y + tools::evaluate_polynomial (P, t)
733740 / tools::evaluate_polynomial (Q, t);
0 commit comments