The function GIfacset computes the set of primary prime factors of its Gaussian integer argument n.

$\mathrm{with}\left(\mathrm{GaussInt}\right)\:$

$\mathrm{GIfacset}\left(72+100\mathrm{I}\right)$

$\left\{\mathrm{-3}-2\mathrm{I}\,\mathrm{-3}+8\mathrm{I}\,1+\mathrm{I}\right\}$

$\mathrm{GIfacset}\left(1574+368\mathrm{I}\right)$

$\left\{\mathrm{-7}+2\mathrm{I}\,1-4\mathrm{I}\,1+\mathrm{I}\,1+2\mathrm{I}\,5-2\mathrm{I}\right\}$