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 S_n.
For given polynomial p, one can by automatic procedure find finite number of disjoint compacts K_1,\dots,K_s (in \mathbb C) such that any compact contains exactly any root of polynomial p is contained in K_i for some i. 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 a can be expressed as root of an polynomial p\neq 0 over \mathbb Q contained in certain fine ball which does not contain any other root of p. 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 p\neq 0 over \mathbb Q such that p(a+b)=0, where a,b are algebraic numbers given in poly form. After localization of polynomial p we can eliminate all compacts but one by numeric computing of a+b. Then, poly form of a+b 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 a is given as root of polynomial p contained in compact K, then a=0\Leftrightarrow(p(0)=0\land 0\in K).
General case of computing Galois group of given polynomial can be easy reduced to the case of polynomial with no multiple roots.
Let a_1,\dots,a_n are all different roots of polynomial p. there is the automatic procedure for computing rationals r_1,\dos,r_n such that b=r_1a_1+\cdots+r_na_n be an primitive element. Then, one can express all conjugates of b and all a_i via b. Automorphisms are determined by image of primitive element. Image of primitive element can be any of its conjugates. Image of a_i can be determined by image of b and compared with a_j. Permutation of \{a_1,\dots,a_n\} is computed.
This is a sketch. I hev no more time now.