diff options
author | Chen-Pang He <jdh8@ms63.hinet.net> | 2013-07-05 23:28:57 +0800 |
---|---|---|
committer | Chen-Pang He <jdh8@ms63.hinet.net> | 2013-07-05 23:28:57 +0800 |
commit | 9e2b4eeac0dcf62f27e05bb0ae923d5c48ee98ed (patch) | |
tree | 56f6ca2b399c17f98b5a3fa39bd5501e6fc8bee7 | |
parent | c273a6c37c74949ee62136b73290f0f9568bb5e1 (diff) |
Const-correct the scaling functor.
-rw-r--r-- | unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h b/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h index f8bbd59c1..ca4532357 100644 --- a/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h +++ b/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h @@ -309,7 +309,7 @@ void MatrixExponential<MatrixType>::computeUV(float) const float maxnorm = 3.925724783138660f; std::frexp(m_l1norm / maxnorm, &m_squarings); 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); } } @@ -329,7 +329,7 @@ void MatrixExponential<MatrixType>::computeUV(double) const double maxnorm = 5.371920351148152; std::frexp(m_l1norm / maxnorm, &m_squarings); 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); } } @@ -352,7 +352,7 @@ void MatrixExponential<MatrixType>::computeUV(long double) const long double maxnorm = 4.0246098906697353063L; std::frexp(m_l1norm / maxnorm, &m_squarings); 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); } #elif LDBL_MANT_DIG <= 106 // double-double @@ -370,7 +370,7 @@ void MatrixExponential<MatrixType>::computeUV(long double) const long double maxnorm = 3.2579440895405400856599663723517L; std::frexp(m_l1norm / maxnorm, &m_squarings); 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); } #elif LDBL_MANT_DIG <= 112 // quadruple precison @@ -388,7 +388,7 @@ void MatrixExponential<MatrixType>::computeUV(long double) const long double maxnorm = 2.884233277829519311757165057717815L; std::frexp(m_l1norm / maxnorm, &m_squarings); 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); } #else |