boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-01-23 17:52:09 +00:00
Code Issues Projects Releases Wiki Activity

Fix deprecated implicit copy constructor in airy ai bi zero

Browse Source
This commit is contained in:
Matt Borland
2022-07-02 18:48:15 -07:00
parent 06ce1cd511
commit 0ed3134719
1 changed files with 73 additions and 29 deletions
Download Patch File Download Diff File Expand all files Collapse all files

102
include/boost/math/special_functions/detail/airy_ai_bi_zero.hpp
View File

@@ -46,7 +46,7 @@
* one_over_z_squared + ( T(77125UL) / 82944UL))
* one_over_z_squared - ( T(5U) / 36U))
* one_over_z_squared + ( T(5U) / 48U))
* one_over_z_squared + (1)));
* one_over_z_squared + 1));
return fz;
}
@@ -60,23 +60,43 @@
switch(m)
{
case 0: { guess = T(0); break; }
case 1: { guess = T(-2.33810741045976703849); break; }
case 2: { guess = T(-4.08794944413097061664); break; }
case 3: { guess = T(-5.52055982809555105913); break; }
case 4: { guess = T(-6.78670809007175899878); break; }
case 5: { guess = T(-7.94413358712085312314); break; }
case 6: { guess = T(-9.02265085334098038016); break; }
case 7: { guess = T(-10.0401743415580859306); break; }
case 8: { guess = T(-11.0085243037332628932); break; }
case 9: { guess = T(-11.9360155632362625170); break; }
case 10:{ guess = T(-12.8287767528657572004); break; }
case 0:
guess = T(0);
break;
case 1:
guess = T(-2.33810741045976703849);
break;
case 2:
guess = T(-4.08794944413097061664);
break;
case 3:
guess = T(-5.52055982809555105913);
break;
case 4:
guess = T(-6.78670809007175899878);
break;
case 5:
guess = T(-7.94413358712085312314);
break;
case 6:
guess = T(-9.02265085334098038016);
break;
case 7:
guess = T(-10.0401743415580859306);
break;
case 8:
guess = T(-11.0085243037332628932);
break;
case 9:
guess = T(-11.9360155632362625170);
break;
case 10:
guess = T(-12.8287767528657572004);
break;
default:
{
const T t(((boost::math::constants::pi<T>() * 3) * ((T(m) * 4) - 1)) / 8);
guess = -boost::math::detail::airy_zero::equation_as_10_4_105(t, pol);
break;
}
}
return guess;
@@ -86,7 +106,9 @@
class function_object_ai_and_ai_prime
{
public:
function_object_ai_and_ai_prime(const Policy& pol) : my_pol(pol) { }
explicit function_object_ai_and_ai_prime(const Policy& pol) : my_pol(pol) { }
function_object_ai_and_ai_prime(const function_object_ai_and_ai_prime&) = default;
boost::math::tuple<T, T> operator()(const T& x) const
{
@@ -111,23 +133,43 @@
switch(m)
{
case 0: { guess = T(0); break; }
case 1: { guess = T(-1.17371322270912792492); break; }
case 2: { guess = T(-3.27109330283635271568); break; }
case 3: { guess = T(-4.83073784166201593267); break; }
case 4: { guess = T(-6.16985212831025125983); break; }
case 5: { guess = T(-7.37676207936776371360); break; }
case 6: { guess = T(-8.49194884650938801345); break; }
case 7: { guess = T(-9.53819437934623888663); break; }
case 8: { guess = T(-10.5299135067053579244); break; }
case 9: { guess = T(-11.4769535512787794379); break; }
case 10: { guess = T(-12.3864171385827387456); break; }
case 0:
guess = T(0);
break;
case 1:
guess = T(-1.17371322270912792492);
break;
case 2:
guess = T(-3.27109330283635271568);
break;
case 3:
guess = T(-4.83073784166201593267);
break;
case 4:
guess = T(-6.16985212831025125983);
break;
case 5:
guess = T(-7.37676207936776371360);
break;
case 6:
guess = T(-8.49194884650938801345);
break;
case 7:
guess = T(-9.53819437934623888663);
break;
case 8:
guess = T(-10.5299135067053579244);
break;
case 9:
guess = T(-11.4769535512787794379);
break;
case 10:
guess = T(-12.3864171385827387456);
break;
default:
{
const T t(((boost::math::constants::pi<T>() * 3) * ((T(m) * 4) - 3)) / 8);
guess = -boost::math::detail::airy_zero::equation_as_10_4_105(t, pol);
break;
}
}
return guess;
@@ -137,7 +179,9 @@
class function_object_bi_and_bi_prime
{
public:
function_object_bi_and_bi_prime(const Policy& pol) : my_pol(pol) { }
explicit function_object_bi_and_bi_prime(const Policy& pol) : my_pol(pol) { }
function_object_bi_and_bi_prime(const function_object_bi_and_bi_prime&) = default;
boost::math::tuple<T, T> operator()(const T& x) const
{