Use vectorize for initializing vector in dummy implementations

closes #59

src/OpenMesh/Core/Geometry/EigenVectorT.hh

@@ -89,9 +89,7 @@ namespace Eigen {
 
   template <typename Derived>
   MatrixBase<Derived> normalize(MatrixBase<Derived> &x) {
-      std::cerr << "Norm : " << x.norm() << std::endl;
       x /= x.norm();
-      std::cerr << "Norm : " << x << std::endl;
       return x;
   }
 
@@ -99,6 +97,7 @@ namespace Eigen {
   MatrixBase<Derived> &vectorize(MatrixBase<Derived> &x,
                                  typename Derived::Scalar const &val) {
       x.fill(val);
+      std::cerr << "Vec" << std::endl;
       return x;
   }
 

src/OpenMesh/Core/Mesh/PolyMeshT.cc

@@ -154,7 +154,15 @@ typename PolyMeshT<Kernel>::Normal
 PolyMeshT<Kernel>::calc_face_normal_impl(FaceHandle, PointIsNot3DTag) const
 {
   // Dummy fallback implementation
-  return Normal(typename Normal::value_type(0));
+  // Returns just an initialized all 0 normal
+  // This function is only used if we don't hate a matching implementation
+  // for normal computation with the current vector type defined in the mesh traits
+
+  assert(false);
+
+  Normal normal;
+  vectorize(normal,0);
+  return normal;
 }
 
 //-----------------------------------------------------------------------------
@@ -212,7 +220,17 @@ template <class Kernel>
 typename PolyMeshT<Kernel>::Normal
 PolyMeshT<Kernel>::calc_face_normal_impl(const Point&, const Point&, const Point&, PointIsNot3DTag) const
 {
-  return Normal(typename Normal::value_type(0));
+
+  // Dummy fallback implementation
+  // Returns just an initialized all 0 normal
+  // This function is only used if we don't hate a matching implementation
+  // for normal computation with the current vector type defined in the mesh traits
+
+  assert(false);
+
+  Normal normal;
+  vectorize(normal,0);
+  return normal;
 }
 
 //-----------------------------------------------------------------------------