Skip to content

Update Boost to 1.81 #3312

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

Merged
Stephan T. Lavavej (StephanTLavavej) merged 3 commits into from
Jan 12, 2023
Merged

Update Boost to 1.81 #3312

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
7b05d7a
Update Boost to 1.81
AZero13 Dec 29, 2022
981619b
switch to isclose
strega-nil Jan 5, 2023
fbc08d5
include utility for swap
strega-nil Jan 6, 2023
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
2 changes: 1 addition & 1 deletion boost-math
Open in desktop
Submodule boost-math updated 189 files
2 changes: 1 addition & 1 deletion docs/cgmanifest.json
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@
"type": "git",
"git": {
"repositoryUrl": "https://github.com/boostorg/math",
"commitHash": "db2a7cbb44a84c4c9cba62232c07f18b22bd965d"
"commitHash": "f395de082eca6ee993a15f2bdb90277eda518295"
}
}
},
Expand Down
31 changes: 24 additions & 7 deletions tests/std/tests/LWG3234_math_special_overloads/test.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,6 +5,7 @@
#include <cmath>
#include <limits>
#include <type_traits>
#include <utility>

using namespace std;

Expand Down Expand Up @@ -102,8 +103,24 @@ Ambiguous sph_neumann(unsigned int, Ambiguous) {
return Ambiguous{};
}

bool expect_epsilons(double expected, double calculated, unsigned int multiple) {
return abs((calculated - expected) / expected) <= multiple * numeric_limits<double>::epsilon();
template <class Float>
bool isclose(Float f, Float g, const int ulps = 1) {
if (f == g) {
return true;
}

if (f > g) {
swap(f, g);
}
// f < g
for (int i = 0; i < ulps; ++i) {
f = nextafter(f, numeric_limits<Float>::infinity());
if (f == g) {
return true;
}
}

return false;
}

void test_assoc_laguerre() {
Expand Down Expand Up @@ -197,7 +214,7 @@ void test_comp_ellint_2() {
static_assert(is_same_v<decltype(comp_ellint_2f(0)), float>);
static_assert(is_same_v<decltype(comp_ellint_2l(0)), long double>);
static_assert(is_same_v<decltype(comp_ellint_2(Ambiguous{})), Ambiguous>);
assert(comp_ellint_2(0) == acos(-1.0) / 2);
assert(isclose(comp_ellint_2(0), acos(-1.0) / 2));
}

void test_comp_ellint_3() {
Expand Down Expand Up @@ -285,7 +302,7 @@ void test_cyl_bessel_k() {
static_assert(is_same_v<decltype(cyl_bessel_kl(0, 0)), long double>);
static_assert(is_same_v<decltype(cyl_bessel_k(0, 0)), double>);
static_assert(is_same_v<decltype(cyl_bessel_k(Ambiguous{}, Ambiguous{})), Ambiguous>);
assert(expect_epsilons(cyl_bessel_k(0.5, 1), acos(-1.0) / 2 * (cyl_bessel_i(-0.5, 1) - cyl_bessel_i(0.5, 1)), 10));
assert(isclose(cyl_bessel_k(0.5, 1), acos(-1.0) / 2 * (cyl_bessel_i(-0.5, 1) - cyl_bessel_i(0.5, 1)), 15));
}

void test_cyl_neumann() {
Expand Down Expand Up @@ -507,7 +524,7 @@ void test_sph_bessel() {
static_assert(is_same_v<decltype(sph_besself(0u, 0)), float>);
static_assert(is_same_v<decltype(sph_bessell(0u, 0)), long double>);
static_assert(is_same_v<decltype(sph_bessel(0u, Ambiguous{})), Ambiguous>);
assert(expect_epsilons(sph_bessel(1, 2), sin(2) / (static_cast<double>(2) * 2) - cos(2) / 2, 2));
assert(isclose(sph_bessel(1, 2), sin(2) / (2.0 * 2) - cos(2) / 2, 2));
}

void test_sph_legendre() {
Expand All @@ -527,7 +544,7 @@ void test_sph_legendre() {
static_assert(is_same_v<decltype(sph_legendre(0u, 0u, Ambiguous{})), Ambiguous>);

const double pi = acos(-1.0);
assert(expect_epsilons(sph_legendre(3, 0, 1), 0.25 * sqrt(7 / pi) * (5 * pow(cos(1), 3) - 3 * cos(1)), 2));
assert(isclose(sph_legendre(3, 0, 1), 0.25 * sqrt(7 / pi) * (5 * pow(cos(1), 3) - 3 * cos(1))));
}

void test_sph_neumann() {
Expand All @@ -545,7 +562,7 @@ void test_sph_neumann() {
static_assert(is_same_v<decltype(sph_neumannf(0u, 0)), float>);
static_assert(is_same_v<decltype(sph_neumannl(0u, 0)), long double>);
static_assert(is_same_v<decltype(sph_neumann(0u, Ambiguous{})), Ambiguous>);
assert(expect_epsilons(sph_neumann(1, 1), -cos(1) - sin(1), 2));
assert(isclose(sph_neumann(1, 1), -cos(1) - sin(1), 2));
}

int main() {
Expand Down