update s and t in Hensel validation

functions.jl

@@ -364,6 +364,7 @@ function hensel_validate(l::Int, f::PolyElem, g::PolyElem, h::PolyElem; ϵ=1e-10
     dg = degree(g)
     dh = degree(h)
     df = degree(f)
+    
     @show dg
     @show dh
 
@@ -374,6 +375,8 @@ function hensel_validate(l::Int, f::PolyElem, g::PolyElem, h::PolyElem; ϵ=1e-10
     #compute Hensel with Float64
     (g,h,s,t) = henseltruncate(l,fr,gr,hr)
     #convert resulting polynomials to thin intervals
+    ds = degree(s)
+    dt = degree(t)
     g = convert_float_to_balls(g)
     h = convert_float_to_balls(h)
     s = convert_float_to_balls(s)
@@ -426,6 +429,18 @@ function hensel_validate(l::Int, f::PolyElem, g::PolyElem, h::PolyElem; ϵ=1e-10
         end
         g = gg
         h = hh        
+
+        #update the cofactors
+        
+        b = s*g +t*h -1
+        BivPolMod!(b,l)
+        (c,d) = fastdivrem(s*b,h,l,validation = true)
+        s = s-d
+        s = truncate(s,ds+1)
+        BivPolMod!(s,l)
+        t = t-t*b-c*g
+        t = truncate(t,dt+1)
+        BivPolMod!(t,l)
     end
     print("no convergence\n")
     (g,h)