Begin updating everything to use Repa arrays (Values3D).

[spline3.git] / src / Point.hs

diff --git a/src/Point.hs b/src/Point.hs
index 5d7954c810634f64fa65599afd352502fce40eca..a92041e89845147f31af704f53613383c4244ee6 100644 (file)
--- a/src/Point.hs
+++ b/src/Point.hs
@@ -18,42 +18,20 @@ z_coord :: Point -> Double
 z_coord (_, _, z) = z
 
 instance Num Point where
-    p1 + p2 = (x1+x2, y1+y2, z1+z2)
-        where
-          x1 = x_coord p1
-          x2 = x_coord p2
-          y1 = y_coord p1
-          y2 = y_coord p2
-          z1 = z_coord p1
-          z2 = z_coord p2
-
-    p1 - p2 = (x1-x2, y1-y2, z1-z2)
-        where
-          x1 = x_coord p1
-          x2 = x_coord p2
-          y1 = y_coord p1
-          y2 = y_coord p2
-          z1 = z_coord p1
-          z2 = z_coord p2
-
-    p1 * p2 = (x1*x2, y1*y2, z1*z2)
-        where
-          x1 = x_coord p1
-          x2 = x_coord p2
-          y1 = y_coord p1
-          y2 = y_coord p2
-          z1 = z_coord p1
-          z2 = z_coord p2
-
+    (x1,y1,z1) + (x2,y2,z2) = (x1+x2, y1+y2, z1+z2)
+    (x1,y1,z1) - (x2,y2,z2) = (x1-x2, y1-y2, z1-z2)
+    (x1,y1,z1) * (x2,y2,z2) = (x1*x2, y1*y2, z1*z2)
     abs (x, y, z) = (abs x, abs y, abs z)
     signum (x, y, z) = (signum x, signum y, signum z)
     fromInteger n = (fromInteger n, fromInteger n, fromInteger n)
 
 
+-- | Scale a point by a constant.
 scale :: Point -> Double -> Point
 scale (x, y, z) d = (x*d, y*d, z*d)
 
 
+-- | Returns the distance between p1 and p2.
 distance :: Point -> Point -> Double
 distance p1 p2 =
     sqrt $ (x2 - x1)^(2::Int) + (y2 - y1)^(2::Int) + (z2 - z1)^(2::Int)
@@ -66,5 +44,7 @@ distance p1 p2 =
       z2 = z_coord p2
 
 
+-- | Returns 'True' if p1 is close to (within 'epsilon' of) p2,
+--   'False' otherwise.
 is_close :: Point -> Point -> Bool
 is_close p1 p2 = (distance p1 p2) ~= 0