source 2005 Doche--Lange compute A = X1^2 compute B = A + Y1 Z1 compute C = X1 Z1 compute D = C^2 compute E = (B^2 + B C + a2 D) compute X3 = C E compute Y3 = (B + C) E + A^2 C compute Z3 = C D