From e21910acf39f7c7f173545939156154772d86cd6 Mon Sep 17 00:00:00 2001
From: Claus Fieker <fieker@mathematik.uni-kl.de>
Date: Wed, 31 Jan 2024 17:50:38 +0100
Subject: [PATCH] add support to adjoin many elements

---
 experimental/GModule/Misc.jl | 38 ++++++++++++++++++++++++++++++++++++
 1 file changed, 38 insertions(+)

diff --git a/experimental/GModule/Misc.jl b/experimental/GModule/Misc.jl
index 72bb48bc7f0f..25b22c05d99f 100644
--- a/experimental/GModule/Misc.jl
+++ b/experimental/GModule/Misc.jl
@@ -1,5 +1,39 @@
+module Misc
+using Oscar
+
 Hecke.minpoly(a::qqbar) = minpoly(Hecke.Globals.Qx, a)
 
+function primitive_element(a::Vector{qqbar})
+  pe = a[1]
+  f = minpoly(pe)
+  Qx = parent(f)
+  for i = 2:length(a)
+    g = minpoly(a[i])
+    f = minpoly(pe)
+    k, _ = number_field(f, check = false, cached = false)
+    lf = collect(keys(factor(k, g).fac))
+    for j = 1:length(lf)
+      h = map_coefficients(x->Qx(x)(pe), lf[j])
+      if is_zero(h(a[i]))
+        d = degree(f) * degree(h)
+        mu = 0
+        while degree(minpoly(pe+mu*a[i])) != d
+          mu += 1
+          if mu > 10
+            error("too bad")
+          end
+        end
+        pe += mu*a[i]
+      end
+    end
+  end
+  return pe
+end
+
+function Hecke.number_field(::QQField, a::Vector{qqbar}; cached::Bool = false)
+  return number_field(QQ, primitive_element(a))
+end
+
 function Hecke.number_field(::QQField, a::qqbar; cached::Bool = false)
   f = minpoly(a)
   k, b = number_field(f, check = false, cached = cached)
@@ -34,6 +68,7 @@ function Hecke.number_field(::QQField, a::qqbar; cached::Bool = false)
 end
 
 Base.getindex(::QQField, a::qqbar) = number_field(QQ, a)
+Base.getindex(::QQField, a::Vector{qqbar}) = number_field(QQ, a)
 
 function Hecke.numerator(f::QQPolyRingElem, parent::ZZPolyRing = Hecke.Globals.Zx)
   g = parent()
@@ -148,3 +183,6 @@ function cyclo_fixed_group_gens(A::AbstractArray{nf_elem})
   end
   return [(mR(sR(ms(x))), F) for x = gens(s)]
 end
+
+end # module
+using .Misc