GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-15 13:15:57 +08:00
Code Issues Packages Projects Releases Wiki Activity

Const-correct the scaling functor.

Browse Source
This commit is contained in:
Chen-Pang He 2013-07-05 23:28:57 +08:00
parent c273a6c37c
commit 9e2b4eeac0
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -309,7 +309,7 @@ void MatrixExponential<MatrixType>::computeUV(float)
const float maxnorm = 3.925724783138660f; const float maxnorm = 3.925724783138660f;
std::frexp(m_l1norm / maxnorm, &m_squarings); std::frexp(m_l1norm / maxnorm, &m_squarings);
if (m_squarings < 0) m_squarings = 0; if (m_squarings < 0) m_squarings = 0;
MatrixType A = CwiseUnaryOp<ScalingOp, MatrixType>(m_M, m_squarings); MatrixType A = CwiseUnaryOp<ScalingOp, const MatrixType>(m_M, m_squarings);
pade7(A); pade7(A);
} }
} }
@ -329,7 +329,7 @@ void MatrixExponential<MatrixType>::computeUV(double)
const double maxnorm = 5.371920351148152; const double maxnorm = 5.371920351148152;
std::frexp(m_l1norm / maxnorm, &m_squarings); std::frexp(m_l1norm / maxnorm, &m_squarings);
if (m_squarings < 0) m_squarings = 0; if (m_squarings < 0) m_squarings = 0;
MatrixType A = CwiseUnaryOp<ScalingOp, MatrixType>(m_M, m_squarings); MatrixType A = CwiseUnaryOp<ScalingOp, const MatrixType>(m_M, m_squarings);
pade13(A); pade13(A);
} }
} }
@ -352,7 +352,7 @@ void MatrixExponential<MatrixType>::computeUV(long double)
const long double maxnorm = 4.0246098906697353063L; const long double maxnorm = 4.0246098906697353063L;
std::frexp(m_l1norm / maxnorm, &m_squarings); std::frexp(m_l1norm / maxnorm, &m_squarings);
if (m_squarings < 0) m_squarings = 0; if (m_squarings < 0) m_squarings = 0;
MatrixType A = CwiseUnaryOp<ScalingOp, MatrixType>(m_M, m_squarings); MatrixType A = CwiseUnaryOp<ScalingOp, const MatrixType>(m_M, m_squarings);
pade13(A); pade13(A);
} }
#elif LDBL_MANT_DIG <= 106 // double-double #elif LDBL_MANT_DIG <= 106 // double-double
@ -370,7 +370,7 @@ void MatrixExponential<MatrixType>::computeUV(long double)
const long double maxnorm = 3.2579440895405400856599663723517L; const long double maxnorm = 3.2579440895405400856599663723517L;
std::frexp(m_l1norm / maxnorm, &m_squarings); std::frexp(m_l1norm / maxnorm, &m_squarings);
if (m_squarings < 0) m_squarings = 0; if (m_squarings < 0) m_squarings = 0;
MatrixType A = CwiseUnaryOp<ScalingOp, MatrixType>(m_M, m_squarings); MatrixType A = CwiseUnaryOp<ScalingOp, const MatrixType>(m_M, m_squarings);
pade17(A); pade17(A);
} }
#elif LDBL_MANT_DIG <= 112 // quadruple precison #elif LDBL_MANT_DIG <= 112 // quadruple precison
@ -388,7 +388,7 @@ void MatrixExponential<MatrixType>::computeUV(long double)
const long double maxnorm = 2.884233277829519311757165057717815L; const long double maxnorm = 2.884233277829519311757165057717815L;
std::frexp(m_l1norm / maxnorm, &m_squarings); std::frexp(m_l1norm / maxnorm, &m_squarings);
if (m_squarings < 0) m_squarings = 0; if (m_squarings < 0) m_squarings = 0;
MatrixType A = CwiseUnaryOp<ScalingOp, MatrixType>(m_M, m_squarings); MatrixType A = CwiseUnaryOp<ScalingOp, const MatrixType>(m_M, m_squarings);
pade17(A); pade17(A);
} }
#else #else