fix bug in TRDH

Author: MaxenceGollier

src/TRDH_alg.jl

@@ -308,6 +308,12 @@ function SolverCore.solve!(
   improper == true && return stats
 
   is_subsolver = h isa ShiftedProximableFunction # case TRDH is used as a subsolver
+
+  if is_subsolver
+    l_bound .= ψ.l
+    u_bound .= ψ.u
+  end
+
   if verbose > 0
     @info log_header(
       [:iter, :fx, :hx, :xi, :ρ, :Δ, :normx, :norms, :normD, :arrow],

src/TR_alg.jl

@@ -454,8 +454,6 @@ function SolverCore.solve!(
       if quasiNewtTest
         @. ∇fk⁻ = ∇fk - ∇fk⁻
         push!(nlp, s, ∇fk⁻) # update QN operator
-        #println(ψ)
-        #stats.iter == 1 && error("done")
       end
 
       solver.subpb.model.B = hess_op(nlp, xk)