diff --git a/src/math/fwd.hpp b/src/math/fwd.hpp
index 0aa74e6345093b74d3f6bb818d9cd13c86cb4a41..97fb2a78d62786d0bdd368aeeaebd9ea6ec93317 100644
--- a/src/math/fwd.hpp
+++ b/src/math/fwd.hpp
@@ -57,10 +57,14 @@ namespace se3
   {
     using std::fabs;
     using std::sqrt;
+    using std::acos;
+    using std::pow;
     
 #ifdef PINOCCHIO_WITH_CPPAD_SUPPORT
     using CppAD::fabs;
     using CppAD::sqrt;
+    using CppAD::acos;
+    using CppAD::pow;
 #endif
   }
 }
diff --git a/src/math/quaternion.hpp b/src/math/quaternion.hpp
index d60bf2e0679c4fb0d5cdc346f8b1f3da3761b5a0..4a685b4832dfa6e7f06b2b94a105faa1ce68e0a7 100644
--- a/src/math/quaternion.hpp
+++ b/src/math/quaternion.hpp
@@ -96,7 +96,7 @@ namespace se3
     const Scalar alpha = ((Scalar)3 - N2) / Scalar(2);
     const_cast <Eigen::QuaternionBase<D> &> (q).coeffs() *= alpha;
 #ifndef NDEBUG
-    const Scalar M = Scalar(3) * std::pow(Scalar(1)-epsilon, ((Scalar)-Scalar(5))/Scalar(2)) / Scalar(4);
+    const Scalar M = Scalar(3) * math::pow(Scalar(1)-epsilon, ((Scalar)-Scalar(5))/Scalar(2)) / Scalar(4);
     assert(math::fabs(q.norm() - Scalar(1)) <=
         std::max(M * sqrt(N2) * (N2 - Scalar(1))*(N2 - Scalar(1)) / Scalar(2), Eigen::NumTraits<Scalar>::dummy_precision()));
 #endif