Index: distance_pythagoras.hpp
===================================================================
--- distance_pythagoras.hpp	(revision 74597)
+++ distance_pythagoras.hpp	(working copy)
@@ -43,8 +43,8 @@
 {
     static inline T apply(Point1 const& p1, Point2 const& p2)
     {
-        T const c1 = boost::numeric_cast<T>(get<I-1>(p2));
-        T const c2 = boost::numeric_cast<T>(get<I-1>(p1));
+        T const c1 = boost::numeric_cast<T>(get<I-1>(p1));
+        T const c2 = boost::numeric_cast<T>(get<I-1>(p2));
         T const d = c1 - c2;
         return d * d + compute_pythagoras<Point1, Point2, I-1, T>::apply(p1, p2);
     }
@@ -99,12 +99,14 @@
 
         assert_dimension_equal<Point1, Point2>();
 
-        return detail::compute_pythagoras
+        calculation_type const t = detail::compute_pythagoras
             <
                 Point1, Point2,
                 dimension<Point1>::value,
                 calculation_type
             >::apply(p1, p2);
+
+        return sqrt(t);
     }
 };