fsqr
fsqr(bb, b); // BB = B^2
fsqr(aa, a); // AA = A^2
fsqr(z2, dacb_m); // (DA - CB)^2
fsqr(x3, dacb_p); // x3 = (DA + CB)^2
fsqr(a0, i);
fsqr(t00, a0);
fsqr(o, inp, tmp);
fsqr(o, o, tmp);