boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-02-26 16:52:27 +00:00
Code Issues Projects Releases Wiki Activity

Mark up zeta constants for coverage.

Browse Source
This commit is contained in:
jzmaddock
2024-07-01 12:25:15 +01:00
parent d441001a5e
commit d8a4900e59
1 changed files with 44 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

52
include/boost/math/special_functions/zeta.hpp
View File

@@ -209,6 +209,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
// Maximum Deviation Found: 2.020e-18
// Expected Error Term: -2.020e-18
// Max error found at double precision: 3.994987e-17
// LCOV_EXCL_START
static const T P[6] = {
static_cast<T>(0.24339294433593750202L),
static_cast<T>(-0.49092470516353571651L),
@@ -225,6 +226,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
static_cast<T>(0.00024978985622317935355L),
static_cast<T>(-0.101855788418564031874e-4L),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, sc) / tools::evaluate_polynomial(Q, sc);
result -= 1.2433929443359375F;
result += (sc);
@@ -234,6 +236,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
{
// Maximum Deviation Found: 9.007e-20
// Expected Error Term: 9.007e-20
// LCOV_EXCL_START
static const T P[6] = {
static_cast<T>(0.577215664901532860516L),
static_cast<T>(0.243210646940107164097L),
@@ -250,6 +253,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
static_cast<T>(0.000255784226140488490982L),
static_cast<T>(0.10991819782396112081e-4L),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(-sc)) / tools::evaluate_polynomial(Q, T(-sc));
result += 1 / (-sc);
}
@@ -257,6 +261,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
{
// Maximum Deviation Found: 5.946e-22
// Expected Error Term: -5.946e-22
// LCOV_EXCL_START
static const float Y = 0.6986598968505859375;
static const T P[6] = {
static_cast<T>(-0.0537258300023595030676L),
@@ -275,6 +280,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
static_cast<T>(0.106951867532057341359e-4L),
static_cast<T>(0.236276623974978646399e-7L),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 2)) / tools::evaluate_polynomial(Q, T(s - 2));
result += Y + 1 / (-sc);
}
@@ -283,7 +289,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
// Maximum Deviation Found: 2.955e-17
// Expected Error Term: 2.955e-17
// Max error found at double precision: 2.009135e-16
// LCOV_EXCL_START
static const T P[6] = {
static_cast<T>(-2.49710190602259410021L),
static_cast<T>(-2.60013301809475665334L),
@@ -303,6 +309,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
static_cast<T>(0.718833729365459760664e-8L),
static_cast<T>(-0.1129200113474947419e-9L),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 4)) / tools::evaluate_polynomial(Q, T(s - 4));
result = 1 + exp(result);
}
@@ -311,6 +318,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
// Maximum Deviation Found: 7.117e-16
// Expected Error Term: 7.117e-16
// Max error found at double precision: 9.387771e-16
// LCOV_EXCL_START
static const T P[7] = {
static_cast<T>(-4.78558028495135619286L),
static_cast<T>(-1.89197364881972536382L),
@@ -331,6 +339,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
static_cast<T>(-0.833378440625385520576e-10L),
static_cast<T>(0.699841545204845636531e-12L),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 7)) / tools::evaluate_polynomial(Q, T(s - 7));
result = 1 + exp(result);
}
@@ -338,6 +347,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
{
// Max error in interpolated form: 1.668e-17
// Max error found at long double precision: 1.669714e-17
// LCOV_EXCL_START
static const T P[8] = {
static_cast<T>(-10.3948950573308896825L),
static_cast<T>(-2.85827219671106697179L),
@@ -358,6 +368,7 @@ inline T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<in
static_cast<T>(0.118507153474022900583e-7L),
static_cast<T>(0.222609483627352615142e-14L),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 15)) / tools::evaluate_polynomial(Q, T(s - 15));
result = 1 + exp(result);
}
@@ -383,6 +394,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
// Maximum Deviation Found: 3.099e-20
// Expected Error Term: 3.099e-20
// Max error found at long double precision: 5.890498e-20
// LCOV_EXCL_START
static const T P[6] = {
BOOST_MATH_BIG_CONSTANT(T, 64, 0.243392944335937499969),
BOOST_MATH_BIG_CONSTANT(T, 64, -0.496837806864865688082),
@@ -400,6 +412,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
BOOST_MATH_BIG_CONSTANT(T, 64, -0.159600883054550987633e-4),
BOOST_MATH_BIG_CONSTANT(T, 64, 0.339770279812410586032e-6),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, sc) / tools::evaluate_polynomial(Q, sc);
result -= 1.2433929443359375F;
result += (sc);
@@ -410,7 +423,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
// Maximum Deviation Found: 1.059e-21
// Expected Error Term: 1.059e-21
// Max error found at long double precision: 1.626303e-19
// LCOV_EXCL_START
static const T P[6] = {
BOOST_MATH_BIG_CONSTANT(T, 64, 0.577215664901532860605),
BOOST_MATH_BIG_CONSTANT(T, 64, 0.222537368917162139445),
@@ -428,6 +441,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
BOOST_MATH_BIG_CONSTANT(T, 64, 0.635994377921861930071e-5),
BOOST_MATH_BIG_CONSTANT(T, 64, 0.226583954978371199405e-7),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(-sc)) / tools::evaluate_polynomial(Q, T(-sc));
result += 1 / (-sc);
}
@@ -435,6 +449,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
{
// Maximum Deviation Found: 5.946e-22
// Expected Error Term: -5.946e-22
// LCOV_EXCL_START
static const float Y = 0.6986598968505859375;
static const T P[7] = {
BOOST_MATH_BIG_CONSTANT(T, 64, -0.053725830002359501027),
@@ -455,12 +470,14 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
BOOST_MATH_BIG_CONSTANT(T, 64, 0.278090318191657278204e-6),
BOOST_MATH_BIG_CONSTANT(T, 64, -0.19683620233222028478e-8),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 2)) / tools::evaluate_polynomial(Q, T(s - 2));
result += Y + 1 / (-sc);
}
else if(s <= 7)
{
// Max error found at long double precision: 8.132216e-19
// LCOV_EXCL_START
static const T P[8] = {
BOOST_MATH_BIG_CONSTANT(T, 64, -2.49710190602259407065),
BOOST_MATH_BIG_CONSTANT(T, 64, -3.36664913245960625334),
@@ -482,6 +499,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
BOOST_MATH_BIG_CONSTANT(T, 64, -0.927884739284359700764e-8),
BOOST_MATH_BIG_CONSTANT(T, 64, 0.119810501805618894381e-9),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 4)) / tools::evaluate_polynomial(Q, T(s - 4));
result = 1 + exp(result);
}
@@ -489,6 +507,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
{
// Max error in interpolated form: 1.133e-18
// Max error found at long double precision: 2.183198e-18
// LCOV_EXCL_START
static const T P[9] = {
BOOST_MATH_BIG_CONSTANT(T, 64, -4.78558028495135548083),
BOOST_MATH_BIG_CONSTANT(T, 64, -3.23873322238609358947),
@@ -511,6 +530,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
BOOST_MATH_BIG_CONSTANT(T, 64, 0.117957556472335968146e-7),
BOOST_MATH_BIG_CONSTANT(T, 64, -0.193432300973017671137e-12),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 7)) / tools::evaluate_polynomial(Q, T(s - 7));
result = 1 + exp(result);
}
@@ -518,6 +538,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
{
// Max error in interpolated form: 1.668e-17
// Max error found at long double precision: 1.669714e-17
// LCOV_EXCL_START
static const T P[9] = {
BOOST_MATH_BIG_CONSTANT(T, 64, -10.3948950573308861781),
BOOST_MATH_BIG_CONSTANT(T, 64, -2.82646012777913950108),
@@ -541,6 +562,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 64>&
BOOST_MATH_BIG_CONSTANT(T, 64, -0.939798249922234703384e-16),
BOOST_MATH_BIG_CONSTANT(T, 64, 0.264584017421245080294e-18),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 15)) / tools::evaluate_polynomial(Q, T(s - 15));
result = 1 + exp(result);
}
@@ -566,7 +588,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
// Maximum Deviation Found: 9.493e-37
// Expected Error Term: 9.492e-37
// Max error found at long double precision: 7.281332e-31
// LCOV_EXCL_START
static const T P[10] = {
BOOST_MATH_BIG_CONSTANT(T, 113, -1.0),
BOOST_MATH_BIG_CONSTANT(T, 113, -0.0353008629988648122808504280990313668),
@@ -592,6 +614,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, -0.217062446168217797598596496310953025e-9),
BOOST_MATH_BIG_CONSTANT(T, 113, 0.315823200002384492377987848307151168e-11),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, sc) / tools::evaluate_polynomial(Q, sc);
result += (sc);
result /= (sc);
@@ -600,7 +623,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
{
// Maximum Deviation Found: 1.616e-37
// Expected Error Term: -1.615e-37
// LCOV_EXCL_START
static const T P[10] = {
BOOST_MATH_BIG_CONSTANT(T, 113, 0.577215664901532860606512090082402431),
BOOST_MATH_BIG_CONSTANT(T, 113, 0.255597968739771510415479842335906308),
@@ -626,6 +649,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, -0.377105263588822468076813329270698909e-11),
BOOST_MATH_BIG_CONSTANT(T, 113, -0.581926559304525152432462127383600681e-13),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(-sc)) / tools::evaluate_polynomial(Q, T(-sc));
result += 1 / (-sc);
}
@@ -634,7 +658,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
// Maximum Deviation Found: 1.891e-36
// Expected Error Term: -1.891e-36
// Max error found: 2.171527e-35
// LCOV_EXCL_START
static const float Y = 0.6986598968505859375;
static const T P[11] = {
BOOST_MATH_BIG_CONSTANT(T, 113, -0.0537258300023595010275848333539748089),
@@ -663,6 +687,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, 0.105677416606909614301995218444080615e-11),
BOOST_MATH_BIG_CONSTANT(T, 113, 0.547223964564003701979951154093005354e-15),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 2)) / tools::evaluate_polynomial(Q, T(s - 2));
result += Y + 1 / (-sc);
}
@@ -670,7 +695,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
{
// Max error in interpolated form: 1.510e-37
// Max error found at long double precision: 2.769266e-34
// LCOV_EXCL_START
static const T Y = 3.28348541259765625F;
static const T P[13] = {
@@ -704,6 +729,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, -0.275363878344548055574209713637734269e-13),
BOOST_MATH_BIG_CONSTANT(T, 113, 0.221564186807357535475441900517843892e-15),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 4)) / tools::evaluate_polynomial(Q, T(s - 4));
result -= Y;
result = 1 + exp(result);
@@ -712,7 +738,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
{
// Max error in interpolated form: 1.999e-34
// Max error found at long double precision: 2.156186e-33
// LCOV_EXCL_START
static const T P[13] = {
BOOST_MATH_BIG_CONSTANT(T, 113, -4.0545627381873738086704293881227365),
BOOST_MATH_BIG_CONSTANT(T, 113, -4.70088348734699134347906176097717782),
@@ -744,6 +770,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, 0.294670713571839023181857795866134957e-16),
BOOST_MATH_BIG_CONSTANT(T, 113, -0.147003914536437243143096875069813451e-18),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 6)) / tools::evaluate_polynomial(Q, T(s - 6));
result = 1 + exp(result);
}
@@ -751,6 +778,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
{
// Max error in interpolated form: 1.641e-32
// Max error found at long double precision: 1.696121e-32
// LCOV_EXCL_START
static const T P[13] = {
BOOST_MATH_BIG_CONSTANT(T, 113, -6.91319491921722925920883787894829678),
BOOST_MATH_BIG_CONSTANT(T, 113, -3.65491257639481960248690596951049048),
@@ -782,6 +810,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, 0.887948682401000153828241615760146728e-19),
BOOST_MATH_BIG_CONSTANT(T, 113, -0.34980761098820347103967203948619072e-21),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 10)) / tools::evaluate_polynomial(Q, T(s - 10));
result = 1 + exp(result);
}
@@ -789,7 +818,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
{
// Max error in interpolated form: 1.563e-31
// Max error found at long double precision: 1.562725e-31
// LCOV_EXCL_START
static const T P[13] = {
BOOST_MATH_BIG_CONSTANT(T, 113, -11.7824798233959252791987402769438322),
BOOST_MATH_BIG_CONSTANT(T, 113, -4.36131215284987731928174218354118102),
@@ -821,6 +850,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, -0.291354445847552426900293580511392459e-22),
BOOST_MATH_BIG_CONSTANT(T, 113, 0.73614324724785855925025452085443636e-25),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 17)) / tools::evaluate_polynomial(Q, T(s - 17));
result = 1 + exp(result);
}
@@ -828,6 +858,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
{
// Max error in interpolated form: 2.311e-27
// Max error found at long double precision: 2.297544e-27
// LCOV_EXCL_START
static const T P[14] = {
BOOST_MATH_BIG_CONSTANT(T, 113, -20.7944102007844314586649688802236072),
BOOST_MATH_BIG_CONSTANT(T, 113, -4.95759941987499442499908748130192187),
@@ -862,6 +893,7 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
BOOST_MATH_BIG_CONSTANT(T, 113, -0.557103423021951053707162364713587374e-31),
BOOST_MATH_BIG_CONSTANT(T, 113, 0.618708773442584843384712258199645166e-34),
};
// LCOV_EXCL_END
result = tools::evaluate_polynomial(P, T(s - 30)) / tools::evaluate_polynomial(Q, T(s - 30));
result = 1 + exp(result);
}
@@ -879,9 +911,11 @@ T zeta_imp_prec(T s, T sc, const Policy&, const std::integral_constant<int, 113>
template <class T, class Policy>
T zeta_imp_odd_integer(int s, const T&, const Policy&, const std::true_type&)
{
// LCOV_EXCL_START
static const T results[] = {
BOOST_MATH_BIG_CONSTANT(T, 113, 1.2020569031595942853997381615114500), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0369277551433699263313654864570342), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0083492773819228268397975498497968), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0020083928260822144178527692324121), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0004941886041194645587022825264699), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0001227133475784891467518365263574), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000305882363070204935517285106451), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000076371976378997622736002935630), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000019082127165539389256569577951), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000004769329867878064631167196044), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000001192199259653110730677887189), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000298035035146522801860637051), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000074507117898354294919810042), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000018626597235130490064039099), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000004656629065033784072989233), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000001164155017270051977592974), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000291038504449709968692943), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000072759598350574810145209), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000018189896503070659475848), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000004547473783042154026799), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000001136868407680227849349), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000284217097688930185546), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000071054273952108527129), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000017763568435791203275), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000004440892103143813364), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000001110223025141066134), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000277555756213612417), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000069388939045441537), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000017347234760475766), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000004336808690020650), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000001084202172494241), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000271050543122347), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000067762635780452), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000016940658945098), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000004235164736273), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000001058791184068), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000264697796017), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000066174449004), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000016543612251), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000004135903063), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000001033975766), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000258493941), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000064623485), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000016155871), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000004038968), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000001009742), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000252435), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000063109), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000015777), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000003944), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000000986), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000000247), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000000062), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000000015), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000000004), BOOST_MATH_BIG_CONSTANT(T, 113, 1.0000000000000000000000000000000001),
};
// LCOV_EXCL_END
return s > 113 ? 1 : results[(s - 3) / 2];
}
@@ -891,9 +925,11 @@ T zeta_imp_odd_integer(int s, const T& sc, const Policy& pol, const std::false_t
#ifdef BOOST_MATH_NO_THREAD_LOCAL_WITH_NON_TRIVIAL_TYPES
static_assert(std::is_trivially_destructible<T>::value, "Your platform does not support thread_local with non-trivial types, last checked with Mingw-x64-8.1, Jan 2021. Please try a Mingw build with the POSIX threading model, see https://sourceforge.net/p/mingw-w64/bugs/527/");
#endif
// LCOV_EXCL_START
static BOOST_MATH_THREAD_LOCAL bool is_init = false;
static BOOST_MATH_THREAD_LOCAL T results[50] = {};
static BOOST_MATH_THREAD_LOCAL int digits = tools::digits<T>();
// LCOV_EXCL_END
int current_digits = tools::digits<T>();
if(digits != current_digits)
{