implement squaring formula and reduce number of temporaries in fmpcb_mul