source GECC Algorithm 3.22 coords jacobian-3 assume Z2 = 1