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

No long double policy: Fix up lambert_w.

Browse Source
This commit is contained in:
jzmaddock
2020-07-20 09:29:18 +01:00
parent 0404e838e2
commit 43ee5254db
1 changed files with 6 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
include/boost/math/special_functions/lambert_w.hpp
View File

@@ -996,20 +996,18 @@ T lambert_w0_approx(T z)
//! and then refined using Halley iterations.
template <class T>
inline T get_near_singularity_param(T z)
inline T do_get_near_singularity_param(T z)
{
BOOST_MATH_STD_USING
const T p2 = 2 * (boost::math::constants::e<T>() * z + 1);
const T p = sqrt(p2);
return p;
}
inline float get_near_singularity_param(float z)
template <class T, class Policy>
inline T get_near_singularity_param(T z, const Policy)
{
return static_cast<float>(get_near_singularity_param((double)z));
}
inline double get_near_singularity_param(double z)
{
return static_cast<double>(get_near_singularity_param((long double)z));
typedef typename policies::precision<T, Policy>::type precision_type;
return static_cast<T>(get_near_singularity_param((precision_type)z));
}
// Forward declarations:
@@ -1220,7 +1218,7 @@ T lambert_w_negative_rational_float(T z, const Policy& pol)
else
{
// z is very close (within 0.01) of the singularity at e^-1.
return lambert_w_singularity_series(get_near_singularity_param(z));
return lambert_w_singularity_series(get_near_singularity_param(z, pol));
}
}