X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FFunctionValues.hs;h=9c789f3decde85e9628d8d00a89793e048e08fdf;hb=8c7f684c5b1a6941a552d2101aff70fe1f9a23ec;hp=b1fede7a414af4af42f9eb374e798786ac2e4366;hpb=3f7331f579118687cd73b977ce6aa7d401f88a09;p=spline3.git

diff --git a/src/FunctionValues.hs b/src/FunctionValues.hs
index b1fede7..9c789f3 100644
--- a/src/FunctionValues.hs
+++ b/src/FunctionValues.hs
@@ -1,3 +1,4 @@
+{-# LANGUAGE BangPatterns #-}
 -- | The FunctionValues module contains the 'FunctionValues' type and
 --   the functions used to manipulate it.
 module FunctionValues (
@@ -29,33 +30,33 @@ import Values (Values3D, dims, idx)
 --   cube. Each value of f can be accessed by the name of its
 --   direction.
 data FunctionValues =
-    FunctionValues { front  :: Double,
-                     back   :: Double,
-                     left   :: Double,
-                     right  :: Double,
-                     top    :: Double,
-                     down   :: Double,
-                     front_left :: Double,
-                     front_right :: Double,
-                     front_down :: Double,
-                     front_top :: Double,
-                     back_left :: Double,
-                     back_right :: Double,
-                     back_down :: Double,
-                     back_top :: Double,
-                     left_down :: Double,
-                     left_top :: Double,
-                     right_down :: Double,
-                     right_top :: Double,
-                     front_left_down :: Double,
-                     front_left_top :: Double,
-                     front_right_down :: Double,
-                     front_right_top :: Double,
-                     back_left_down :: Double,
-                     back_left_top :: Double,
-                     back_right_down :: Double,
-                     back_right_top :: Double,
-                     interior :: Double }
+    FunctionValues { front              :: !Double,
+                     back               :: !Double,
+                     left               :: !Double,
+                     right              :: !Double,
+                     top                :: !Double,
+                     down               :: !Double,
+                     front_left         :: !Double,
+                     front_right        :: !Double,
+                     front_down         :: !Double,
+                     front_top          :: !Double,
+                     back_left          :: !Double,
+                     back_right         :: !Double,
+                     back_down          :: !Double,
+                     back_top           :: !Double,
+                     left_down          :: !Double,
+                     left_top           :: !Double,
+                     right_down         :: !Double,
+                     right_top          :: !Double,
+                     front_left_down    :: !Double,
+                     front_left_top     :: !Double,
+                     front_right_down   :: !Double,
+                     front_right_top    :: !Double,
+                     back_left_down     :: !Double,
+                     back_left_top      :: !Double,
+                     back_right_down    :: !Double,
+                     back_right_top     :: !Double,
+                     interior           :: !Double }
       deriving (Eq, Show)
 
 
@@ -200,7 +201,7 @@ eval f (Quotient x y) = (eval f x) / (eval f y)
 --   5.0
 --
 value_at :: Values3D -> Int -> Int -> Int -> Double
-value_at v3d i j k
+value_at v3d !i !j !k
   -- Put the most common case first!
   | (valid_i i) && (valid_j j) && (valid_k k) =
       idx v3d i j k
@@ -249,14 +250,6 @@ value_at v3d i j k
           2*(value_at v3d i j 0) - (value_at v3d i j 1)
         else
           2*(value_at v3d i j (k-1)) - (value_at v3d i j (k-2))
-
-  | otherwise =
-      let istr = show i
-          jstr = show j
-          kstr = show k
-          coordstr = "(" ++ istr ++ "," ++ jstr ++ "," ++ kstr ++ ")"
-      in
-        error $ "value_at called outside of domain: " ++ coordstr
   where
     (dim_i, dim_j, dim_k) = dims v3d
 
@@ -275,7 +268,7 @@ value_at v3d i j k
 --   coordinates (i,j,k), constructs and returns the 'FunctionValues'
 --   object centered at (i,j,k)
 make_values :: Values3D -> Int -> Int -> Int -> FunctionValues
-make_values values i j k =
+make_values values !i !j !k =
     empty_values { front  = value_at values (i-1) j k,
                    back   = value_at values (i+1) j k,
                    left   = value_at values i (j-1) k,