Rename the xi function to domain_point, and make a new xi function which is an alias...

diff --git a/src/Tetrahedron.hs b/src/Tetrahedron.hs
index 032ce1eb9855b821d8964c44a824b911c1765a59..ad164d2426dce01e85cbf3d22cab0d0a19380410 100644 (file)
--- a/src/Tetrahedron.hs
+++ b/src/Tetrahedron.hs
@@ -54,10 +54,15 @@ polynomial t =
 
 
 -- | Returns the domain point of t with indices i,j,k,l.
+--   Simply an alias for the domain_point function.
 xi :: Tetrahedron -> Int -> Int -> Int -> Int -> Point
-xi t i j k l
+xi = domain_point
+
+-- | Returns the domain point of t with indices i,j,k,l.
+domain_point :: Tetrahedron -> Int -> Int -> Int -> Int -> Point
+domain_point t i j k l
    | i + j + k + l == 3 = weighted_sum `scale` (1/3)
-   | otherwise = error "xi index out of bounds"
+   | otherwise = error "domain point index out of bounds"
    where
      v0' = (v0 t) `scale` (fromIntegral i)
      v1' = (v1 t) `scale` (fromIntegral j)
@@ -66,7 +71,6 @@ xi t i j k l
      weighted_sum = v0' + v1' + v2' + v3'
 
 
-
 -- | The Bernstein polynomial on t with indices i,j,k,l. Denoted by a
 --   capital 'B' in the Sorokina/Zeilfelder paper.
 beta :: Tetrahedron -> Int -> Int -> Int -> Int -> (RealFunction Point)