X-Git-Url: https://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FMain.hs;h=9ad6199ec07edb2c0906235207ca45b7d59c1f56;hb=6fb9ab6b6068870323e996da931fc04c7710e3e4;hp=4467d174b94fa4d31d7273c47d9ef7afcdc5aa90;hpb=f97dbb52b73bd90a89940b9653ff274654aba9de;p=spline3.git diff --git a/src/Main.hs b/src/Main.hs index 4467d17..9ad6199 100644 --- 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)