source 2005 Doche--Lange assume Z2=1 compute O1 = Z1^2 compute B = X2 O1 compute D = Y2 O1 Z1 compute E = X1 + B compute F = Y1 + D compute Z3 = E Z1 compute H = F X2 + Z3 Y2 compute I = F + Z3 compute G = Z3^2 compute X3 = a2 G + F I + E E^2 compute Y3 = I X3 + G H