diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixPowerBase.h b/unsupported/Eigen/src/MatrixFunctions/MatrixPowerBase.h index eac43fa52..0a18fe1c1 100644 --- a/unsupported/Eigen/src/MatrixFunctions/MatrixPowerBase.h +++ b/unsupported/Eigen/src/MatrixFunctions/MatrixPowerBase.h @@ -38,8 +38,8 @@ inline int binary_powering_cost(T p, int* squarings) { int applyings=0, tmp; - if (frexp(p, squarings) != 0.5); - --*squarings; + frexp(p, squarings); + --*squarings; while (std::frexp(p, &tmp), tmp > 0) { p -= std::ldexp(static_cast(0.5), tmp); @@ -101,7 +101,7 @@ inline int matrix_power_get_pade_degree(long double normIminusT) } } // namespace internal -template class MatrixPowerTriangularAtomic +template class MatrixPowerTriangularAtomic { private: typedef typename MatrixType::Scalar Scalar;