fgoldi_sub(&t0, Y, &t0); /* t0 = J = E - 2 H */
fgoldi_mul(X, X, &t0); /* X = c (B - E) J */
fgoldi_mul(Z, Y, &t0); /* Z = E J */
fgoldi_sub(&t0, Y, &t0); /* t0 = J = E - 2 H */
fgoldi_mul(X, X, &t0); /* X = c (B - E) J */
fgoldi_mul(Z, Y, &t0); /* Z = E J */