X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FCube.hs;h=7b139a4741aff80a5bae666395aa61c5d45aff5d;hb=2813bba7169c1457e630b3870858e96fd4bf61da;hp=8f22266049538cc697d3c0cbd83eb492654f28d3;hpb=a06c2cdbc189c6b16add4ca06856f0727c12a617;p=spline3.git diff --git a/src/Cube.hs b/src/Cube.hs index 8f22266..7b139a4 100644 --- a/src/Cube.hs +++ b/src/Cube.hs @@ -701,7 +701,7 @@ prop_v0_all_equal cube = (v0 t0) == (v0 t1) -- | Given in Sorokina and Zeilfelder, p. 79, (2.6). Note that the --- third and fourth indices of c-t1 have been switched. This is +-- third and fourth indices of c-t3 have been switched. This is -- because we store the triangles oriented such that their volume is -- positive. If T and T-tilde share \ and v3,v3-tilde point -- in opposite directions, one of them has to have negative volume! @@ -750,8 +750,8 @@ prop_c0120_identity5 cube = t4 = tetrahedron cube 4 t5 = tetrahedron cube 5 --- -- | Given in Sorokina and Zeilfelder, p. 79, (2.6). Repeats --- -- 'prop_c0120_identity1' with tetrahedrons 5 and 6. +-- | Given in Sorokina and Zeilfelder, p. 79, (2.6). Repeats +-- 'prop_c0120_identity1' with tetrahedrons 5 and 6. prop_c0120_identity6 :: Cube -> Bool prop_c0120_identity6 cube = c t6 0 1 2 0 ~= (c t6 0 0 2 1 + c t5 0 0 1 2) / 2 @@ -760,8 +760,8 @@ prop_c0120_identity6 cube = t6 = tetrahedron cube 6 --- -- | Given in Sorokina and Zeilfelder, p. 79, (2.6). Repeats --- -- 'prop_c0120_identity1' with tetrahedrons 6 and 7. +-- | Given in Sorokina and Zeilfelder, p. 79, (2.6). Repeats +-- 'prop_c0120_identity1' with tetrahedrons 6 and 7. prop_c0120_identity7 :: Cube -> Bool prop_c0120_identity7 cube = c t7 0 1 2 0 ~= (c t7 0 0 2 1 + c t6 0 0 1 2) / 2 @@ -955,23 +955,6 @@ prop_c1011_identity cube = t6 = tetrahedron cube 6 - --- | Given in Sorokina and Zeilfelder, p. 78. -prop_cijk1_identity :: Cube -> Bool -prop_cijk1_identity cube = - and [ c t0 i j k 1 ~= - (c t1 (i+1) j k 0) * ((b0 t0) (v3 t1)) + - (c t1 i (j+1) k 0) * ((b1 t0) (v3 t1)) + - (c t1 i j (k+1) 0) * ((b2 t0) (v3 t1)) + - (c t1 i j k 1) * ((b3 t0) (v3 t1)) | i <- [0..2], - j <- [0..2], - k <- [0..2], - i + j + k == 2] - where - t0 = tetrahedron cube 0 - t1 = tetrahedron cube 1 - - -- | The function values at the interior should be the same for all -- tetrahedra. prop_interior_values_all_identical :: Cube -> Bool @@ -1139,14 +1122,6 @@ prop_t7_shares_edge_with_t20 cube = t20 = tetrahedron cube 20 - - - -p78_25_properties :: Test.Framework.Test -p78_25_properties = - testGroup "p. 78, Section (2.5) Properties" [ - testProperty "c_ijk1 identity" prop_cijk1_identity ] - p79_26_properties :: Test.Framework.Test p79_26_properties = testGroup "p. 79, Section (2.6) Properties" [ @@ -1207,7 +1182,6 @@ edge_incidence_tests = cube_properties :: Test.Framework.Test cube_properties = testGroup "Cube Properties" [ - p78_25_properties, p79_26_properties, p79_27_properties, p79_28_properties,