X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FTetrahedron.hs;h=2342ba139d2b79cf64c2b5dbac60c3fa0de55b3a;hb=9220b061771e4ac6d2bd3adad3d926479b77f7c8;hp=e84b09103d7ad1bb850bed4b9380aeb47d6bb3b3;hpb=525f4e76394c74e5a6f6bcdc7fb50d0dc2b9ec2d;p=spline3.git

diff --git a/src/Tetrahedron.hs b/src/Tetrahedron.hs
index e84b091..2342ba1 100644
--- a/src/Tetrahedron.hs
+++ b/src/Tetrahedron.hs
@@ -281,16 +281,24 @@ c _ _ _ _ _ = error "coefficient index out of bounds"
 
 
 
+-- | Compute the determinant of the 4x4 matrix,
+--
+--   [1]
+--   [x]
+--   [y]
+--   [z]
+--
+--   where [1] = [1, 1, 1, 1],
+--         [x] = [x1,x2,x3,x4],
+--
+--   et cetera.
+--
 det :: Point -> Point -> Point -> Point -> Double
 det p0 p1 p2 p3 =
---  Both of these results are just copy/pasted from Sage. One of them
---  might be more numerically stable, faster, or both.
---
---  x1*y2*z4 - x1*y2*z3 + x1*y3*z2 - x1*y3*z4 - x1*y4*z2 + x1*y4*z3 +
---  x2*y1*z3 - x2*y1*z4 - x2*y3*z1 + x2*y3*z4 +
---  x2*y4*z1 - x2*y4*z3 - x3*y1*z2 + x3*y1*z4 + x3*y2*z1 - x3*y2*z4 - x3*y4*z1 +
---  x3*y4*z2 + x4*y1*z2 - x4*y1*z3 - x4*y2*z1 + x4*y2*z3 + x4*y3*z1 - x4*y3*z2
-  -((x2 - x3)*y1 - (x1 - x3)*y2 + (x1 - x2)*y3)*z4 + ((x2 - x4)*y1 - (x1 - x4)*y2 + (x1 - x2)*y4)*z3 + ((x3 - x4)*y2 - (x2 - x4)*y3 + (x2 - x3)*y4)*z1 - ((x3 - x4)*y1 - (x1 - x4)*y3 + (x1 - x3)*y4)*z2
+  x1*y2*z4 - x1*y2*z3 + x1*y3*z2 - x1*y3*z4 - x1*y4*z2 + x1*y4*z3 +
+  x2*y1*z3 - x2*y1*z4 - x2*y3*z1 + x2*y3*z4 + x2*y4*z1 + x3*y1*z4 +
+  x3*y2*z1 - x3*y2*z4 - x3*y4*z1 - x2*y4*z3 - x3*y1*z2 + x3*y4*z2 +
+  x4*y1*z2 - x4*y1*z3 - x4*y2*z1 + x4*y2*z3 + x4*y3*z1 - x4*y3*z2
   where
     (x1, y1, z1) = p0
     (x2, y2, z2) = p1