boost/math
2
0
Fork 0
mirror of https://github.com/boostorg/math.git synced 2026-02-27 17:12:22 +00:00
Code Issues Projects Releases Wiki Activity

Try correction mentioned in the issue thread

Browse Source
This commit is contained in:
Christopher Kormanyos
2021-01-31 14:59:56 +01:00
parent 88c7cc8191
commit 2eac693e12
2 changed files with 38 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
include/boost/math/cstdfloat/cstdfloat_complex_std.hpp
View File

@@ -509,7 +509,37 @@
inline complex<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE> pow(const BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE& x,
const complex<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE>& a)
{
return std::exp(a * std::log(x));
complex<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE> result;
if(x > BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE(0))
{
result = std::exp(a * std::log(x));
}
else if(x < BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE(0))
{
using std::atan2;
using std::log;
const complex<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE>
cpx_lg_x
(
log(-x),
atan2(BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE(0), BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE(-1))
);
result = std::exp(a * cpx_lg_x);
}
else
{
result =
complex<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE>
(
-std::numeric_limits<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE>::infinity(),
BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE(0)
);
}
return result;
}
inline complex<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE> sinh(const complex<BOOST_CSTDFLOAT_EXTENDED_COMPLEX_FLOAT_TYPE>& x)