GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-04-23 10:09:36 +08:00
Code Issues Packages Projects Releases Wiki Activity

using std::frexp instead of frexp

Browse Source
This commit is contained in:
jdh8 2012-08-08 17:50:56 +08:00
parent 8cddcaf839
commit c5800a2452
1 changed files with 3 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h
View File

@ -283,9 +283,8 @@ EIGEN_STRONG_INLINE void MatrixExponential<MatrixType>::pade17(const MatrixType
template <typename MatrixType> template <typename MatrixType>
void MatrixExponential<MatrixType>::computeUV(float) void MatrixExponential<MatrixType>::computeUV(float)
{ {
using std::max; using std::frexp;
using std::pow; using std::pow;
using std::ceil;
if (m_l1norm < 4.258730016922831e-001) { if (m_l1norm < 4.258730016922831e-001) {
pade3(m_M); pade3(m_M);
} else if (m_l1norm < 1.880152677804762e+000) { } else if (m_l1norm < 1.880152677804762e+000) {
@ -302,9 +301,8 @@ void MatrixExponential<MatrixType>::computeUV(float)
template <typename MatrixType> template <typename MatrixType>
void MatrixExponential<MatrixType>::computeUV(double) void MatrixExponential<MatrixType>::computeUV(double)
{ {
using std::max; using std::frexp;
using std::pow; using std::pow;
using std::ceil;
if (m_l1norm < 1.495585217958292e-002) { if (m_l1norm < 1.495585217958292e-002) {
pade3(m_M); pade3(m_M);
} else if (m_l1norm < 2.539398330063230e-001) { } else if (m_l1norm < 2.539398330063230e-001) {
@ -325,9 +323,8 @@ void MatrixExponential<MatrixType>::computeUV(double)
template <typename MatrixType> template <typename MatrixType>
void MatrixExponential<MatrixType>::computeUV(long double) void MatrixExponential<MatrixType>::computeUV(long double)
{ {
using std::max; using std::frexp;
using std::pow; using std::pow;
using std::ceil;
#if LDBL_MANT_DIG == 53 // double precision #if LDBL_MANT_DIG == 53 // double precision
computeUV(double()); computeUV(double());
#elif LDBL_MANT_DIG <= 64 // extended precision #elif LDBL_MANT_DIG <= 64 // extended precision