Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

First batch changes double-fp-backend #1226

Merged
merged 7 commits into from
Dec 19, 2024
Merged

First batch changes double-fp-backend #1226

Show file tree
Hide file tree
Changes from 4 commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
a9bc5d6
First batch changes double-fp-backend
ckormanyos Dec 19, 2024
8a45d18
Try repairing table syntax errors
ckormanyos Dec 19, 2024
464922e
Attempt to restore 3 missing cover lines
ckormanyos Dec 19, 2024
dbc9e72
Repair typo on double min_exponent10
ckormanyos Dec 19, 2024
0308556
Increase tols for some updated MSVC runners
ckormanyos Dec 19, 2024
214e19c
Apply cyl_bessel_k0 logic change 113 bits
ckormanyos Dec 19, 2024
bf7fb37
Remove wrong constexpr and revert bessel_j change
ckormanyos Dec 19, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 3 additions & 1 deletion reporting/performance/test_kn.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -48,7 +48,9 @@ static const std::array<std::array<T, 3>, 9> kn_data = { {
{ { SC_(-10.0), SC_(1.0), SC_(1.80713289901029454691597861302340015908245782948536080022119e8) } },
{ { SC_(100.0), SC_(5.0), SC_(7.03986019306167654653386616796116726248616158936088056952477e115) } },
{ { SC_(100.0), SC_(80.0), SC_(8.39287107246490782848985384895907681748152272748337807033319e-12) } },
{ { SC_(-1000.0), SC_(700.0), SC_(6.51561979144735818903553852606383312984409361984128221539405e-31) } },
(std::numeric_limits<T>::is_specialized && (std::numeric_limits<T>::min_exponent10 > -300))
? std::array<T, 3> { { SC_(-1000.0), SC_(660.0), SC_(5.231804195291135637493087823300645230085362374177628598848978) } }
: std::array<T, 3> { { SC_(-1000.0), SC_(700.0), SC_(6.51561979144735818903553852606383312984409361984128221539405e-31) } },
} };

int main()
Expand Down
2 changes: 1 addition & 1 deletion test/test_bessel_i.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -214,7 +214,7 @@ void test_bessel(T, const char* name)
BOOST_CHECK_CLOSE_FRACTION(boost::math::cyl_bessel_i(T(0.5), T(11357)), SC_(7.173138695269929329584326974917488634629578339622112563648e4929), tolerance * mul);
}
#endif
BOOST_IF_CONSTEXPR (std::numeric_limits<T>::max_exponent > 1000)
BOOST_IF_CONSTEXPR (std::numeric_limits<T>::max_exponent10 > 304)
{
BOOST_IF_CONSTEXPR(std::is_floating_point<T>::value == false)
tolerance *= 4; // multiprecision type.
Expand Down
5 changes: 4 additions & 1 deletion test/test_bessel_j.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -281,7 +281,10 @@ void test_bessel(T, const char* name)
BOOST_CHECK_EQUAL(boost::math::cyl_bessel_j(T(0), T(2.5)), boost::math::cyl_bessel_j(T(0), T(-2.5)));
BOOST_CHECK_EQUAL(boost::math::cyl_bessel_j(T(1), T(2.5)), -boost::math::cyl_bessel_j(T(1), T(-2.5)));
#ifndef SYCL_LANGUAGE_VERSION
BOOST_CHECK_CLOSE_FRACTION(boost::math::cyl_bessel_j(364, T(38.5)), SC_(1.793940496519190500748409872348034004417458734118663909894e-309), tolerance);
BOOST_IF_CONSTEXPR (std::numeric_limits<T>::min_exponent10 <= -307)
{
BOOST_CHECK_CLOSE_FRACTION(boost::math::cyl_bessel_j(364, T(38.5)), SC_(1.793940496519190500748409872348034004417458734118663909894e-309), tolerance);
}
#endif
//
// Special cases at infinity:
Expand Down
4 changes: 3 additions & 1 deletion test/test_bessel_k.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -133,7 +133,9 @@ void test_bessel(T, const char* name)
{{ SC_(-10.0), SC_(1.0), SC_(1.80713289901029454691597861302340015908245782948536080022119e8) }},
{{ SC_(100.0), SC_(5.0), SC_(7.03986019306167654653386616796116726248616158936088056952477e115) }},
{{ SC_(100.0), SC_(80.0), SC_(8.39287107246490782848985384895907681748152272748337807033319e-12) }},
{{ SC_(-1000.0), SC_(700.0), SC_(6.51561979144735818903553852606383312984409361984128221539405e-31) }},
(std::numeric_limits<T>::is_specialized && (std::numeric_limits<T>::min_exponent10 > -300))
? std::array<typename table_type<T>::type, 3> {{ SC_(-1000.0), SC_(660.0), SC_(5.231804195291135637493087823300645230085362374177628598848978) }}
: std::array<typename table_type<T>::type, 3> {{ SC_(-1000.0), SC_(700.0), SC_(6.51561979144735818903553852606383312984409361984128221539405e-31) }},
}};
static const std::array<std::array<typename table_type<T>::type, 3>, 11> kv_data = {{
{{ SC_(0.5), SC_(0.875), SC_(0.558532231646608646115729767013630967055657943463362504577189) }},
Expand Down
Loading