From 60ff027e52faf66c111b3dc2779c3bf5c7dfb6b9 Mon Sep 17 00:00:00 2001
From: Michael Orlitzky <michael@orlitzky.com>
Date: Sun, 8 May 2011 16:56:09 -0400
Subject: [PATCH] Add two more function comments.

---
 src/Tetrahedron.hs | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/src/Tetrahedron.hs b/src/Tetrahedron.hs
index 9d90f71..bcf9a0b 100644
--- a/src/Tetrahedron.hs
+++ b/src/Tetrahedron.hs
@@ -76,6 +76,11 @@ beta t i j k l
     b3_term = (b3 t) `fexp` l
 
 
+-- | The coefficient function. c t i j k l returns the coefficient
+--   c_ijkl with respect to the tetrahedron t. The definition uses
+--   pattern matching to mimic the definitions given in Sorokina and
+--   Zeilfelder, pp. 84-86. If incorrect indices are supplied, the
+--   function will simply error.
 c :: Tetrahedron -> Int -> Int -> Int -> Int -> Double
 c t 0 0 3 0 = eval (fv t) $
               (1/8) * (I + F + L + T + LT + FL + FT + FLT)
@@ -197,6 +202,8 @@ c _ _ _ _ _ = error "coefficient index out of bounds"
 
 
 
+-- | The matrix used in the tetrahedron volume calculation as given in
+--   Lai & Schumaker, Definition 15.4, page 436.
 vol_matrix :: Tetrahedron -> Matrix Double
 vol_matrix t = (4><4)
                [1,  1,  1,  1,
-- 
2.43.2