what does it mean to compute a galois group numerically? do you give a multiplication table? or do you identify it with a subgroup of S(n), or of a free group by generators and relations?
The algorithm output can be multiplication table or subgroup of [tex]S_n[/tex].
For given polynomial [tex]p[/tex], one can by automatic procedure find finite number of disjoint compacts [tex]K_1,\dots,K_s[/tex] (in [tex]\mathbb C[/tex]) such that any compact contains exactly any root of polynomial [tex]p[/tex] is contained in [tex]K_i[/tex] for some [tex]i[/tex]. Let we call the closed ball with rational radius and both real and imaginary part of center as fine ball. Compacts can be fine balls. Let we call this procedure as localization.
Any algebraic number [tex]a[/tex] can be expressed as root of an polynomial [tex]p\neq 0[/tex] over [tex]\mathbb Q[/tex] contained in certain fine ball which does not contain any other root of [tex]p[/tex]. Let we call this form of algebraic numbers as poly form. The algebraic number given in poly form can be computed numerical by automatic procedure.
We can calculate with algebraic numbers in poly form. Using previously explained procedure we can compute polynomial [tex]p\neq 0[/tex] over [tex]\mathbb Q[/tex] such that [tex]p(a+b)=0[/tex], where [tex]a,b[/tex] are algebraic numbers given in poly form. After localization of polynomial [tex]p[/tex] we can eliminate all compacts but one by numeric computing of [tex]a+b[/tex]. Then, poly form of [tex]a+b[/tex] is computed. The cases of multiplication, division and subtraction are analogous. Case of the n-th root is easer. these are left to reader.
Comparing algebraic numbers given in poly form is reduced to case of comparing with zero. If an algebraic number [tex]a[/tex] is given as root of polynomial [tex]p[/tex] contained in compact [tex]K[/tex], then [tex]a=0\Leftrightarrow(p(0)=0\land 0\in K)[/tex].
General case of computing Galois group of given polynomial can be easy reduced to the case of polynomial with no multiple roots.
Let [tex]a_1,\dots,a_n[/tex] are all different roots of polynomial [tex]p[/tex]. there is the automatic procedure for computing rationals [tex]r_1,\dos,r_n[/tex] such that [tex]b=r_1a_1+\cdots+r_na_n[/tex] be an primitive element. Then, one can express all conjugates of [tex]b[/tex] and all [tex]a_i[/tex] via [tex]b[/tex]. Automorphisms are determined by image of primitive element. Image of primitive element can be any of its conjugates. Image of [tex]a_i[/tex] can be determined by image of [tex]b[/tex] and compared with [tex]a_j[/tex]. Permutation of [tex]\{a_1,\dots,a_n\}[/tex] is computed.
This is a sketch. I hev no more time now.