source GECC Algorithm 3.21
parameter half
assume half = 1/2
coords jacobian-3
