diff --git a/edo/src/edoEstimatorNormalMulti.h b/edo/src/edoEstimatorNormalMulti.h
index 76d535da..f5666b8f 100644
--- a/edo/src/edoEstimatorNormalMulti.h
+++ b/edo/src/edoEstimatorNormalMulti.h
@@ -163,8 +163,10 @@ public:
     {
     public:
         typedef typename EOT::AtomType AtomType;
-        typedef Eigen::Matrix< AtomType, Eigen::Dynamic, 1> Vector;
-        typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic> Matrix;
+        // typedef typename edoNormalMulti<AtomType>::Vector Vector;
+        // typedef typename edoNormalMulti<AtomType>::Matrix Matrix;
+        typedef Eigen::Matrix< AtomType, 1, Eigen::Dynamic, Eigen::RowMajor> Vector;
+        typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> Matrix;
 
         CovMatrix( const eoPop< EOT >& pop )
         {
diff --git a/edo/src/edoNormalMulti.h b/edo/src/edoNormalMulti.h
index d5f76871..1f34bc83 100644
--- a/edo/src/edoNormalMulti.h
+++ b/edo/src/edoNormalMulti.h
@@ -83,8 +83,8 @@ class edoNormalMulti : public edoDistrib< EOT >
 {
 public:
     typedef typename EOT::AtomType AtomType;
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, 1> Vector; // Note: by default, Eigen is column-major
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic> Matrix;
+    typedef Eigen::Matrix< AtomType, 1, Eigen::Dynamic, Eigen::RowMajor> Vector;
+    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> Matrix;
 
     edoNormalMulti(
         const Vector & mean,
diff --git a/edo/src/edoNormalMultiCenter.h b/edo/src/edoNormalMultiCenter.h
index 8dac8461..07d5e3d3 100644
--- a/edo/src/edoNormalMultiCenter.h
+++ b/edo/src/edoNormalMultiCenter.h
@@ -57,11 +57,13 @@ class edoNormalMultiCenter : public edoModifierMass< edoNormalMulti< EOT > >
 {
 public:
     typedef typename EOT::AtomType AtomType;
+    // typedef typename edoNormalMulti<AtomType>::Vector Vector;
+    typedef Eigen::Matrix< AtomType, 1, Eigen::Dynamic, Eigen::RowMajor> Vector;
 
     void operator() ( edoNormalMulti< EOT >& distrib, EOT& mass )
     {
         assert( distrib.size() == mass.innerSize() );
-        Eigen::Matrix< AtomType, Eigen::Dynamic, 1 > mean( mass );
+        Vector mean( mass );
         distrib.mean() = mean;
     }
 };
diff --git a/edo/src/edoSamplerNormalMulti.h b/edo/src/edoSamplerNormalMulti.h
index b0a43d1a..702d05af 100644
--- a/edo/src/edoSamplerNormalMulti.h
+++ b/edo/src/edoSamplerNormalMulti.h
@@ -96,8 +96,10 @@ class edoSamplerNormalMulti : public edoSampler< EOD >
 {
 public:
     typedef typename EOT::AtomType AtomType;
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, 1> Vector;
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic> Matrix;
+    // typedef typename edoNormalMulti<AtomType>::Vector Vector;
+    // typedef typename edoNormalMulti<AtomType>::Matrix Matrix;
+    typedef Eigen::Matrix< AtomType, 1, Eigen::Dynamic, Eigen::RowMajor> Vector;
+    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> Matrix;
 
     edoSamplerNormalMulti( edoRepairer<EOT> & repairer ) 
         : edoSampler< EOD >( repairer)
diff --git a/edo/test/t-edoEstimatorNormalMulti.cpp b/edo/test/t-edoEstimatorNormalMulti.cpp
index c62f6dfa..e1fd8871 100644
--- a/edo/test/t-edoEstimatorNormalMulti.cpp
+++ b/edo/test/t-edoEstimatorNormalMulti.cpp
@@ -48,8 +48,10 @@ typedef EOT::AtomType AtomType;
 #else
 #ifdef WITH_EIGEN
 #include <Eigen/Dense>
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, 1> Vector;
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic> Matrix;
+    // typedef typename edoNormalMulti<AtomType>::Vector Vector;
+    // typedef typename edoNormalMulti<AtomType>::Matrix Matrix;
+    typedef Eigen::Matrix< AtomType, 1, Eigen::Dynamic, Eigen::RowMajor> Vector;
+    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> Matrix;
 #endif
 #endif
 
diff --git a/edo/test/t-mean-distance.cpp b/edo/test/t-mean-distance.cpp
index d119c1a9..45448e82 100644
--- a/edo/test/t-mean-distance.cpp
+++ b/edo/test/t-mean-distance.cpp
@@ -53,8 +53,10 @@ typedef EOT::AtomType AtomType;
 #else
 #ifdef WITH_EIGEN
 #include <Eigen/Dense>
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, 1> Vector;
-    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic> Matrix;
+    // typedef typename edoNormalMulti<AtomType>::Vector Vector;
+    // typedef typename edoNormalMulti<AtomType>::Matrix Matrix;
+    typedef Eigen::Matrix< AtomType, 1, Eigen::Dynamic, Eigen::RowMajor> Vector;
+    typedef Eigen::Matrix< AtomType, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> Matrix;
 #endif
 #endif