def METHOD_NAME():  # FIXME
    x = SpatialCoordinate(tetrahedron)
    xv = (5, 7, 9)
    a, b, c = 6.5, -4, 3
    d, e, f = 2, 3, 4
    xx = as_matrix(((x[0], x[1], x[2]),
                    (a, b, c),
                    (d, e, f)))
    s = det(2 * xx)
    e = s(xv)
    v = 2 ** 3 * \
        (xv[0] * (b * f - e * c) - xv[1] *
         (a * f - c * d) + xv[2] * (a * e - b * d))
    assert e == v