diff --git a/unsupported/Eigen/src/NonLinear/MathFunctions.h b/unsupported/Eigen/src/NonLinear/MathFunctions.h
index b2c696b3f..7002b4605 100644
--- a/unsupported/Eigen/src/NonLinear/MathFunctions.h
+++ b/unsupported/Eigen/src/NonLinear/MathFunctions.h
@@ -30,13 +30,13 @@
 template<typename Functor, typename Scalar>
 // TODO : fixe Scalar here
 int ei_hybrd1(
-        Eigen::Matrix< double, Eigen::Dynamic, 1 >  &x,
-        Eigen::Matrix< double, Eigen::Dynamic, 1 >  &fvec,
+        Eigen::Matrix< Scalar, Eigen::Dynamic, 1 >  &x,
+        Eigen::Matrix< Scalar, Eigen::Dynamic, 1 >  &fvec,
         Scalar tol = Eigen::ei_sqrt(Eigen::machine_epsilon<Scalar>())
         )
 {
     int lwa = (x.size()*(3*x.size()+13))/2;
-    Eigen::Matrix< double, Eigen::Dynamic, 1 > wa(lwa);
+    Eigen::Matrix< Scalar, Eigen::Dynamic, 1 > wa(lwa);
     fvec.resize(x.size());
     return hybrd1(Functor::f, 0, x.size(), x.data(), fvec.data(), tol, wa.data(), lwa);
 }