diff --git a/Core/Geometry/src/BasicVector3D.cpp b/Core/Geometry/src/BasicVector3D.cpp
index e7eee932ddeedad0f401f25a6b14d6d2862a62d4..eb7b4141121b9b7fa50c7be258f43047d31718ce 100644
--- a/Core/Geometry/src/BasicVector3D.cpp
+++ b/Core/Geometry/src/BasicVector3D.cpp
@@ -109,14 +109,6 @@ double BasicVector3D<double>::cosTheta() const
     return std::abs(ma) == 0 ? 1 : z()/ma;
 }
 
-// //! @TODO eliminate this, it is plain wrong
-// template<>
-//complex_t BasicVector3D<complex_t>::cosTheta() const
-//{
-//   complex_t ma = mag();
-//   return std::abs(ma) == 0 ? 1 : z()/ma;
-//}
-
 // -----------------------------------------------------------------------------
 // Combine two vectors
 // -----------------------------------------------------------------------------
@@ -144,16 +136,6 @@ BasicVector3D<double> BasicVector3D<double>::cross(
                                  x()*v.y()-v.x()*y());
 }
 
-// //! @TODO check usage: unlikely to be correct
-// template<>
-// BasicVector3D<complex_t> BasicVector3D<complex_t>::cross(
-//    const BasicVector3D<complex_t>& v) const
-//{
-//    return BasicVector3D<complex_t>(y()*v.z()-v.y()*z(),
-//                                 z()*v.x()-v.z()*x(),
-//                                 x()*v.y()-v.x()*y());
-//}
-
 //! Returns square of transverse component with respect to given axis.
 template<>
 double BasicVector3D<double>::perp2(const BasicVector3D<double>& v) const