Begin overhauling the program to handle other tetrahedra. Main is

[spline3.git] / src / Main.hs

diff --git a/src/Main.hs b/src/Main.hs
index 4467d174b94fa4d31d7273c47d9ef7afcdc5aa90..9ad6199ec07edb2c0906235207ca45b7d59c1f56 100644 (file)
--- a/src/Main.hs
+++ b/src/Main.hs
@@ -1,14 +1,14 @@
 module Main
 where
 
-import Cube
-import Face
-import Grid
-import Misc (flatten)
-import Point
-import RealFunction
-import Tetrahedron
-import ThreeDimensional
+--import Cube
+--import Face
+--import Grid
+--import Misc (flatten)
+--import Point
+--import RealFunction
+--import Tetrahedron
+--import ThreeDimensional
 
 trilinear :: [[[Double]]]
 trilinear = [ [ [ 1, 2, 3 ],
@@ -48,43 +48,44 @@ dummy = [ [ [ 0, 1, 2 ],
             [ 24, 25, 26 ]]]
 
 
-find_point_value :: RealFunction Point
-find_point_value p = poly p
-    where
-      g0 = make_grid 1 trilinear
-      the_cubes = flatten (cubes g0)
-      good_cubes = filter ((flip contains_point) p) the_cubes
-      target_cube = head good_cubes
-      good_tets = filter ((flip contains_point) p) (tetrahedrons target_cube)
-      target_tetrahedron = head good_tets
-      poly = polynomial target_tetrahedron
+--find_point_value :: RealFunction Point
+--find_point_value p = poly p
+--    where
+--      g0 = make_grid 1 trilinear
+--      the_cubes = flatten (cubes g0)
+--      good_cubes = filter ((flip contains_point) p) the_cubes
+--       target_cube = head good_cubes
+--       good_tets = filter ((flip contains_point) p) (tetrahedrons target_cube)
+--       target_tetrahedron = head good_tets
+--       poly = polynomial target_tetrahedron
 
 main :: IO ()
 main = do
-  print $ find_point_value (0,0,0)
-  print $ find_point_value (1,0,0)
-  print $ find_point_value (2,0,0)
-  print $ find_point_value (0,1,0)
-  print $ find_point_value (1,1,0)
-  print $ find_point_value (2,1,0)
-  print $ find_point_value (0,2,0)
-  print $ find_point_value (1,2,0)
-  print $ find_point_value (2,2,0)
-  print $ find_point_value (0,0,1)
-  print $ find_point_value (1,0,1)
-  print $ find_point_value (2,0,1)
-  print $ find_point_value (0,1,1)
-  print $ find_point_value (1,1,1)
-  print $ find_point_value (2,1,1)
-  print $ find_point_value (0,2,1)
-  print $ find_point_value (1,2,1)
-  print $ find_point_value (2,2,1)
-  print $ find_point_value (0,0,2)
-  print $ find_point_value (1,0,2)
-  print $ find_point_value (2,0,2)
-  print $ find_point_value (0,1,2)
-  print $ find_point_value (1,1,2)
-  print $ find_point_value (2,1,2)
-  print $ find_point_value (0,2,2)
-  print $ find_point_value (1,2,2)
-  print $ find_point_value (2,2,2)
+  putStrLn "Hello, World."
+  -- print $ find_point_value (0,0,0)
+  -- print $ find_point_value (1,0,0)
+  -- print $ find_point_value (2,0,0)
+  -- print $ find_point_value (0,1,0)
+  -- print $ find_point_value (1,1,0)
+  -- print $ find_point_value (2,1,0)
+  -- print $ find_point_value (0,2,0)
+  -- print $ find_point_value (1,2,0)
+  -- print $ find_point_value (2,2,0)
+  -- print $ find_point_value (0,0,1)
+  -- print $ find_point_value (1,0,1)
+  -- print $ find_point_value (2,0,1)
+  -- print $ find_point_value (0,1,1)
+  -- print $ find_point_value (1,1,1)
+  -- print $ find_point_value (2,1,1)
+  -- print $ find_point_value (0,2,1)
+  -- print $ find_point_value (1,2,1)
+  -- print $ find_point_value (2,2,1)
+  -- print $ find_point_value (0,0,2)
+  -- print $ find_point_value (1,0,2)
+  -- print $ find_point_value (2,0,2)
+  -- print $ find_point_value (0,1,2)
+  -- print $ find_point_value (1,1,2)
+  -- print $ find_point_value (2,1,2)
+  -- print $ find_point_value (0,2,2)
+  -- print $ find_point_value (1,2,2)
+  -- print $ find_point_value (2,2,2)