faster ``gf_edf_zassenhaus`` using the Frobenius map

For instance in this PR

``minpoly(sqrt(1 + root(2, 5)) + root(2, 5)*sqrt(1 + root(2, 10)))``

is 35% faster than in master.