11 import ThreeDimensional
13 trilinear :: [[[Double]]]
14 trilinear = [ [ [ 1, 2, 3 ],
25 zeros = [ [ [ 0, 0, 0 ],
38 dummy = [ [ [ 0, 1, 2 ],
51 find_point_value :: RealFunction Point
52 find_point_value p = poly p
54 g0 = make_grid 1 trilinear
55 the_cubes = flatten (cubes g0)
56 good_cubes = filter ((flip contains_point) p) the_cubes
57 target_cube = head good_cubes
58 good_tets = filter ((flip contains_point) p) (tetrahedrons target_cube)
59 target_tetrahedron = head good_tets
60 poly = polynomial target_tetrahedron
64 print $ find_point_value (0,0,0)
65 print $ find_point_value (1,0,0)
66 print $ find_point_value (2,0,0)
67 print $ find_point_value (0,1,0)
68 print $ find_point_value (1,1,0)
69 print $ find_point_value (2,1,0)
70 print $ find_point_value (0,2,0)
71 print $ find_point_value (1,2,0)
72 print $ find_point_value (2,2,0)
73 print $ find_point_value (0,0,1)
74 print $ find_point_value (1,0,1)
75 print $ find_point_value (2,0,1)
76 print $ find_point_value (0,1,1)
77 print $ find_point_value (1,1,1)
78 print $ find_point_value (2,1,1)
79 print $ find_point_value (0,2,1)
80 print $ find_point_value (1,2,1)
81 print $ find_point_value (2,2,1)
82 print $ find_point_value (0,0,2)
83 print $ find_point_value (1,0,2)
84 print $ find_point_value (2,0,2)
85 print $ find_point_value (0,1,2)
86 print $ find_point_value (1,1,2)
87 print $ find_point_value (2,1,2)
88 print $ find_point_value (0,2,2)
89 print $ find_point_value (1,2,2)
90 print $ find_point_value (2,2,2)