X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FMain.hs;h=c4204c821f4f534f1df7c72784b155dbd7072f83;hb=f6d0c289ad3397cf392976c24f3afdb17da5d377;hp=8376ce954ff9f840d2b88f11b65a2393a283cd5f;hpb=89b8b6e94fcc944a1f4611811265f3c6217af850;p=spline3.git diff --git a/src/Main.hs b/src/Main.hs index 8376ce9..c4204c8 100644 --- a/src/Main.hs +++ b/src/Main.hs @@ -1,104 +1,27 @@ 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 ], - [ 1, 3, 5 ], - [ 1, 4, 7 ] ], - [ [ 1, 2, 3 ], - [ 1, 4, 7 ], - [ 1, 6, 11 ] ], - [ [ 1, 2, 3 ], - [ 1, 5, 9 ], - [ 1, 8, 15 ]]] -zeros :: [[[Double]]] -zeros = [ [ [ 0, 0, 0 ], - [ 0, 0, 0 ], - [ 0, 0, 0 ] ], - -- - [ [ 0, 0, 0 ], - [ 0, 0, 0 ], - [ 0, 0, 0 ] ], - -- - [ [ 0, 0, 0 ], - [ 0, 0, 0 ], - [ 0, 0, 0 ]]] -dummy :: [[[Double]]] -dummy = [ [ [ 0, 1, 2 ], - [ 3, 4, 5 ], - [ 6, 7, 8 ] ], - -- - [ [ 9, 10, 11 ], - [ 12, 13, 14 ], - [ 15, 16, 17 ] ], - -- - [ [ 18, 19, 20 ], - [ 21, 22, 23 ], - [ 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 = good_cubes !! 0 - good_tets = filter ((flip contains_point) p) (tetrahedrons target_cube) - target_tetrahedron = good_tets !! 0 - 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 - putStrLn $ show $ find_point_value (0,0,0) - putStrLn $ show $ find_point_value (1,0,0) - putStrLn $ show $ find_point_value (2,0,0) - putStrLn $ show $ find_point_value (0,1,0) - putStrLn $ show $ find_point_value (1,1,0) - putStrLn $ show $ find_point_value (2,1,0) - putStrLn $ show $ find_point_value (0,2,0) - putStrLn $ show $ find_point_value (1,2,0) - putStrLn $ show $ find_point_value (2,2,0) - putStrLn $ show $ find_point_value (0,0,1) - putStrLn $ show $ find_point_value (1,0,1) - putStrLn $ show $ find_point_value (2,0,1) - putStrLn $ show $ find_point_value (0,1,1) - putStrLn $ show $ find_point_value (1,1,1) - putStrLn $ show $ find_point_value (2,1,1) - putStrLn $ show $ find_point_value (0,2,1) - putStrLn $ show $ find_point_value (1,2,1) - putStrLn $ show $ find_point_value (2,2,1) - putStrLn $ show $ find_point_value (0,0,2) - putStrLn $ show $ find_point_value (1,0,2) - putStrLn $ show $ find_point_value (2,0,2) - putStrLn $ show $ find_point_value (0,1,2) - putStrLn $ show $ find_point_value (1,1,2) - putStrLn $ show $ find_point_value (2,1,2) - putStrLn $ show $ find_point_value (0,2,2) - putStrLn $ show $ find_point_value (1,2,2) - putStrLn $ show $ find_point_value (2,2,2) - -- let g0 = make_grid 1 trilinear - -- let the_cubes = flatten (cubes g0) - -- putStrLn $ show $ the_cubes - -- let p = (2, 0, 0) - -- let target_cubes = filter ((flip contains_point) p) the_cubes - -- putStrLn $ show $ target_cubes - -- let target_cube = (take 1 target_cubes) !! 0 - -- putStrLn $ show $ target_cube - -- let target_tetrahedra = filter ((flip contains_point) p) (tetrahedrons target_cube) - -- let target_tetrahedron = (take 1 target_tetrahedra) !! 0 - -- putStrLn $ show $ target_tetrahedron - -- let poly = polynomial target_tetrahedron - -- putStrLn $ show $ poly - -- putStrLn $ show $ poly p +main = putStrLn "Hello, World."