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 bessel jy zero

Browse Source
This commit is contained in:
Matt Borland
2022-07-02 18:30:01 -07:00
parent 99adf50919
commit 06ce1cd511
1 changed files with 13 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
include/boost/math/special_functions/detail/bessel_jy_zero.hpp
View File

@@ -57,7 +57,9 @@
class equation_as_9_3_39_and_its_derivative
{
public:
equation_as_9_3_39_and_its_derivative(const T& zt) : zeta(zt) { }
explicit equation_as_9_3_39_and_its_derivative(const T& zt) : zeta(zt) { }
equation_as_9_3_39_and_its_derivative(const equation_as_9_3_39_and_its_derivative&) = default;
boost::math::tuple<T, T> operator()(const T& z) const
{
@@ -126,11 +128,11 @@
const T range_zmin = (std::max<T>)(z_estimate - T(1), T(1));
const T range_zmax = z_estimate + T(1);
const int my_digits10 = static_cast<int>(static_cast<float>(boost::math::tools::digits<T>() * 0.301F));
const auto my_digits10 = static_cast<int>(static_cast<float>(boost::math::tools::digits<T>() * 0.301F));
// Select the maximum allowed iterations based on the number
// of decimal digits in the numeric type T, being at least 12.
const std::uintmax_t iterations_allowed = static_cast<std::uintmax_t>((std::max)(12, my_digits10 * 2));
const auto iterations_allowed = static_cast<std::uintmax_t>((std::max)(12, my_digits10 * 2));
std::uintmax_t iterations_used = iterations_allowed;
@@ -187,6 +189,8 @@
const Policy& pol) : my_v(v),
my_pol(pol) { }
function_object_jv(const function_object_jv&) = default;
T operator()(const T& x) const
{
return boost::math::cyl_bessel_j(my_v, x, my_pol);
@@ -208,6 +212,8 @@
my_order_is_zero(order_is_zero),
my_pol(pol) { }
function_object_jv_and_jv_prime(const function_object_jv_and_jv_prime&) = default;
boost::math::tuple<T, T> operator()(const T& x) const
{
// Obtain Jv(x) and Jv'(x).
@@ -405,6 +411,8 @@
const Policy& pol) : my_v(v),
my_pol(pol) { }
function_object_yv(const function_object_yv&) = default;
T operator()(const T& x) const
{
return boost::math::cyl_neumann(my_v, x, my_pol);
@@ -424,6 +432,8 @@
const Policy& pol) : my_v(v),
my_pol(pol) { }
function_object_yv_and_yv_prime(const function_object_yv_and_yv_prime&) = default;
boost::math::tuple<T, T> operator()(const T& x) const
{
const T half_epsilon(boost::math::tools::epsilon<T>() / 2U);