-- | Make sure we can reproduce a 9x9x9 trilinear from the 3x3x3 one.
+-- Use (t <- tetrahedra c0) for a much slower but comprehensive
+-- test.
test_trilinear9x9x9_reproduced :: Assertion
test_trilinear9x9x9_reproduced =
assertTrue "trilinear 9x9x9 is reproduced correctly" $
| i <- [0..8],
j <- [0..8],
k <- [0..8],
- t <- tetrahedra c0,
+ t <- [head $ tetrahedra c0],
let p = polynomial t,
let i' = (fromIntegral i) * 0.5,
let j' = (fromIntegral j) * 0.5,