t0 = X1^2
t1 = Y1^2
t2 = Z1^2
t3 = X1 * Y1
t4 = Y1 * Z1
t3 = t3 + t3
Z = Z1 * X1
Z = Z * 2
Y = t2 * b
Y = Y - Z
Y = 3 * Y
X = t1 - Y
Y = Y + t1
Y = Y * X
X = X * t3
t2 = t2 * 3
Z = Z * b
Z = Z - t2
Z = Z - t0
Z = Z * 3
t0 = t0 * 3
t0 = t0 - t2
t0 = t0 * Z
Y3 = Y + t0
t4 = t4 + t4
Z = Z * t4
X3 = X - Z
Z = t4 * t1
Z3 = Z * 4