X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FPoint.hs;h=fd3ac58dec3ceef5dfd5c8aa1615e25a9bc0af53;hb=e725cfe579b9d05ac040efc08f1ad47e5060de38;hp=94bcd70f6e8c586f28a70e655c6bbd227ba1c801;hpb=58cf11569acb270995d2de924dda03ef526647e2;p=spline3.git

diff --git a/src/Point.hs b/src/Point.hs
index 94bcd70..fd3ac58 100644
--- a/src/Point.hs
+++ b/src/Point.hs
@@ -1,50 +1,18 @@
-{-# LANGUAGE TypeSynonymInstances #-}
+{-# LANGUAGE FlexibleInstances #-}
 
-module Point
+module Point (
+  Point,
+  dot,
+  scale
+  )
 where
 
-import Comparisons
-
-
 type Point = (Double, Double, Double)
 
-x_coord :: Point -> Double
-x_coord (x, _, _) = x
-
-y_coord :: Point -> Double
-y_coord (_, y, _) = y
-
-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)
@@ -55,19 +23,8 @@ 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)
-    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
-
-
--- | Returns 'True' if p1 'is_close' to p2, 'False' otherwise.
-is_close :: Point -> Point -> Bool
-is_close p1 p2 = (distance p1 p2) ~= 0
+-- | Returns the dot product of two points (taken as three-vectors).
+{-# INLINE dot #-}
+dot :: Point -> Point -> Double
+dot (x1, y1, z1) (x2, y2, z2) =
+    (x2 - x1)^(2::Int) + (y2 - y1)^(2::Int) + (z2 - z1)^(2::Int)