The function GIfacpoly computes the factorization of a multivariate polynomial with Gaussian integer coefficients.

The result is returned in the form  where ,  is an irreducible polynomial over Gaussian ring with the gcd of its coefficients being 1,  is its exponent (multiplicity) and u is a Gaussian integer.