diff --git a/Registration_test/InertiaAlign_test.cpp b/Registration_test/InertiaAlign_test.cpp
index ea3053e..68d3dfc 100644
--- a/Registration_test/InertiaAlign_test.cpp
+++ b/Registration_test/InertiaAlign_test.cpp
@@ -68,7 +68,7 @@ using namespace modeling;
             //GenerateRigid (massSource, centerSource, &meshSource);
 
             type::Vec3d centerTarget ;
-            sofa::type::Vector3 translation;
+            sofa::type::Vec3 translation;
             defaulttype::Rigid3Mass massTarget;
 
             GenerateRigidMass::SPtr rigidSource = New<GenerateRigidMass>();
diff --git a/src/Registration/ClosestPointRegistrationForceField.inl b/src/Registration/ClosestPointRegistrationForceField.inl
index 1e1d920..a4c9304 100644
--- a/src/Registration/ClosestPointRegistrationForceField.inl
+++ b/src/Registration/ClosestPointRegistrationForceField.inl
@@ -385,12 +385,12 @@ void ClosestPointRegistrationForceField<DataTypes>::draw(const core::visual::Vis
     unsigned int nb = this->closestPos.size();
     if (vparams->displayFlags().getShowForceFields())
     {
-        std::vector< type::Vector3 > points;
+        std::vector< type::Vec3 > points;
         for (unsigned int i=0; i<nb; i++)
             if(!sourceIgnored[i])
             {
-                type::Vector3 point1 = DataTypes::getCPos(x[i]);
-                type::Vector3 point2 = DataTypes::getCPos(this->closestPos[i]);
+                type::Vec3 point1 = DataTypes::getCPos(x[i]);
+                type::Vec3 point2 = DataTypes::getCPos(this->closestPos[i]);
                 points.push_back(point1);
                 points.push_back(point2);
             }
@@ -409,8 +409,8 @@ void ClosestPointRegistrationForceField<DataTypes>::draw(const core::visual::Vis
         for (unsigned int i=0; i<nb; i++)
             if(!sourceIgnored[i])
             {
-                type::Vector3 point1 = DataTypes::getCPos(x[i]);
-                type::Vector3 point2 = DataTypes::getCPos(this->closestPos[i]);
+                type::Vec3 point1 = DataTypes::getCPos(x[i]);
+                type::Vec3 point2 = DataTypes::getCPos(this->closestPos[i]);
                 dists[i]=(point2-point1).norm();
             }
         Real max=0; for (unsigned int i=0; i<dists.size(); i++) if(max<dists[i]) max=dists[i];
diff --git a/src/Registration/InertiaAlign.cpp b/src/Registration/InertiaAlign.cpp
index 8a6223f..58a4d76 100644
--- a/src/Registration/InertiaAlign.cpp
+++ b/src/Registration/InertiaAlign.cpp
@@ -485,7 +485,7 @@ void InertiaAlign::init()
         positionDistSource = (*m_positions.beginEdit());
         for (size_t j = 0; j < waPositions.size(); j++)
         {
-            type::Vector4 pointS,pointT;
+            type::Vec4 pointS,pointT;
             pointS(0) = (*m_positions.beginEdit())[j][0];
             pointS(1) = (*m_positions.beginEdit())[j][1];
             pointS(2) = (*m_positions.beginEdit())[j][2];
@@ -577,7 +577,7 @@ void InertiaAlign::init()
         msg_info() << "The MTransformSourceMatrix is not a Transformation matrix";
     for (size_t i = 0; i < waPositions.size(); i++)
     {
-        type::Vector4 pointS,pointT;
+        type::Vec4 pointS,pointT;
         pointS(0) = (*m_positions.beginEdit())[i][0];
         pointS(1) = (*m_positions.beginEdit())[i][1];
         pointS(2) = (*m_positions.beginEdit())[i][2];
diff --git a/src/Registration/InertiaAlign.h b/src/Registration/InertiaAlign.h
index b4835fb..f9e68ba 100644
--- a/src/Registration/InertiaAlign.h
+++ b/src/Registration/InertiaAlign.h
@@ -52,8 +52,8 @@ class InertiaAlign: public sofa::core::objectmodel::BaseObject
       * Data Fields
       */
     /// input
-    Data <sofa::type::Vector3> targetC;
-    Data <sofa::type::Vector3> sourceC; ///< input: the gravity center of the source mesh
+    Data <sofa::type::Vec3> targetC;
+    Data <sofa::type::Vec3> sourceC; ///< input: the gravity center of the source mesh
 
     Data < Mat3x3 > targetInertiaMatrix; ///< input: the inertia matrix of the target mesh
     Data < Mat3x3 > sourceInertiaMatrix; ///< input: the inertia matrix of the source mesh
@@ -68,7 +68,7 @@ class InertiaAlign: public sofa::core::objectmodel::BaseObject
 
 protected:
 
-    typedef type::Vector3 Vector3;
+    typedef type::Vec3 Vector3;
     typedef type::Matrix4 Matrix4;
 
     SReal computeDistances(type::vector<sofa::type::Vec<3,SReal> >, type::vector<sofa::type::Vec<3,SReal> >);
diff --git a/src/Registration/RegistrationContactForceField.inl b/src/Registration/RegistrationContactForceField.inl
index ad29cf5..64f1f47 100644
--- a/src/Registration/RegistrationContactForceField.inl
+++ b/src/Registration/RegistrationContactForceField.inl
@@ -164,7 +164,7 @@ void RegistrationContactForceField<DataTypes>::draw(const core::visual::VisualPa
 
     glDisable(GL_LIGHTING);
 
-    std::vector< type::Vector3 > points[4];
+    std::vector< type::Vec3 > points[4];
 
     for (unsigned int i=0; i<cc.size(); i++)
 	{
@@ -199,7 +199,7 @@ void RegistrationContactForceField<DataTypes>::draw(const core::visual::VisualPa
     vparams->drawTool()->drawLines(points[2], 1, type::RGBAColor::red());
     vparams->drawTool()->drawLines(points[3], 1, type::RGBAColor::green());
 
-	std::vector< type::Vector3 > pointsN;
+    std::vector< type::Vec3 > pointsN;
 	if (vparams->displayFlags().getShowNormals())
 	{
         for (unsigned int i=0; i<cc.size(); i++)
diff --git a/src/Registration/RegistrationExporter.cpp b/src/Registration/RegistrationExporter.cpp
index 0fd617a..7055a71 100644
--- a/src/Registration/RegistrationExporter.cpp
+++ b/src/Registration/RegistrationExporter.cpp
@@ -73,12 +73,12 @@ void RegistrationExporter::init()
             if (this->f_printLog.getValue()) std::cout<<"RegistrationExporter: "<<this->inFileNames.back()<<"  ->  "<<this->outFileNames.back()<<std::endl;
 
             // get inverse transforms applied in loader
-                        type::Vector3 scale=loaders[l]->getScale();
+                        type::Vec3 scale=loaders[l]->getScale();
                         Mat4x4 m_scale; m_scale.fill(0);   for(unsigned int i=0;i<3;i++)	 m_scale[i][i]=1./scale[i]; m_scale[3][3]=1.;
                         type::Quat<SReal> q = type::Quat< SReal >::createQuaterFromEuler(type::Vec< 3, SReal >(loaders[l]->getRotation()) * M_PI / 180.0);
                         Mat4x4 m_rot;
                         q.inverse().toHomogeneousMatrix(m_rot);
-                        type::Vector3 translation=loaders[l]->getTranslation();
+                        type::Vec3 translation=loaders[l]->getTranslation();
                         Mat4x4 m_translation; m_translation.fill(0);   for(unsigned int i=0;i<3;i++)	 m_translation[i][3]=-translation[i];  for(unsigned int i=0;i<4;i++)	 m_translation[i][i]=1;
                         this->inverseTransforms.push_back(m_scale*m_rot*m_translation);
                         if (this->f_printLog.getValue()) std::cout<<"RegistrationExporter: transform = "<<this->inverseTransforms.back()<<std::endl;