From: Michael Orlitzky Date: Fri, 8 Jul 2011 12:05:01 +0000 (-0400) Subject: Uncomment the (2.5) cijk1 identity and make a TODO item for it. X-Git-Tag: 0.0.1~268 X-Git-Url: https://gitweb.michael.orlitzky.com/?a=commitdiff_plain;h=c639d04341cb39f5b8e85a2d139c61017cbcaa66;p=spline3.git Uncomment the (2.5) cijk1 identity and make a TODO item for it. --- diff --git a/doc/TODO b/doc/TODO index a49de2b..2f8af09 100644 --- a/doc/TODO +++ b/doc/TODO @@ -6,3 +6,5 @@ * Add the section (2.9) tests. * Figure out why TC.prop_c0120_identity4 and friends are failing. + +* Figure out why the section (2.5) c_ijk1 identity is failing. diff --git a/src/Tests/Cube.hs b/src/Tests/Cube.hs index 631af5f..a772f4b 100644 --- a/src/Tests/Cube.hs +++ b/src/Tests/Cube.hs @@ -692,20 +692,19 @@ prop_c1011_identity cube = -- | 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 = tetrahedron0 cube --- t1 = tetrahedron1 cube - +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 = tetrahedron0 cube + t1 = tetrahedron1 cube -- | The function values at the interior should be the same for all tetrahedra.