Function fp2_mul 

pub fn fp2_mul(x: &Fp2, y: &Fp2, p: &BigInt) -> Fp2

(a+bi)(c+di) = (ac − bd) + (ad + bc)i.