diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2013-08-11 17:52:43 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2013-08-11 17:52:43 +0200 |
| commit | 6719e56b5bfe9ae4badc9a6e894c5824f663d62e (patch) | |
| tree | 975a7bb2e9f1bf605814309fc7120cdb097b29ca /unsupported/Eigen/src/MatrixFunctions | |
| parent | c13e9bbabf885abad49aafb77535192262355448 (diff) | |
Ref<> objects must be nested by reference because they potentially store a temporary object
Diffstat (limited to 'unsupported/Eigen/src/MatrixFunctions')
| -rw-r--r-- | unsupported/Eigen/src/MatrixFunctions/MatrixPower.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h b/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h index ee665c18e..b91d8632a 100644 --- a/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h +++ b/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h @@ -299,7 +299,7 @@ MatrixPowerAtomic<MatrixType>::computeSuperDiag(const ComplexScalar& curr, const ComplexScalar logCurr = log(curr); ComplexScalar logPrev = log(prev); int unwindingNumber = ceil((numext::imag(logCurr - logPrev) - M_PI) / (2*M_PI)); - ComplexScalar w = numext::atanh2(curr - prev, curr + prev) + ComplexScalar(0, M_PI*unwindingNumber); + ComplexScalar w = numext::atanh2(curr - prev, curr + prev) + ComplexScalar(RealScalar(0), RealScalar(M_PI*unwindingNumber)); return RealScalar(2) * exp(RealScalar(0.5) * p * (logCurr + logPrev)) * sinh(p * w) / (curr - prev); } |