It's slow compared to any efficient multivariate polynomial GCD algorithm. But it's a nice application of Groebner bases and it was already implemented so I kept it.