add condition on B before resetting

MaxenceGollier · MaxenceGollier · commit 3c898ec6f756 · 2026-02-13T16:55:16.000-05:00
diff --git a/src/R2N.jl b/src/R2N.jl
@@ -98,7 +98,8 @@ end
 
 function SolverCore.reset!(solver::R2NSolver)
   _reset_power_method!(solver.v0)
-  LinearOperators.reset!(solver.subpb.model.B)
+  B = solver.subpb.model.B
+  isa(B, AbstractLinearOperator) && LinearOperators.reset!(B)
 end
 
 SolverCore.reset!(solver::R2NSolver, model) = SolverCore.reset!(solver)
diff --git a/src/TR_alg.jl b/src/TR_alg.jl
@@ -97,7 +97,8 @@ end
 
 function SolverCore.reset!(solver::TRSolver)
   _reset_power_method!(solver.v0)
-  LinearOperators.reset!(solver.subpb.model.B)
+  reset_data!(solver.subpb.model)
+  LinearOperators.reset!(solver.subpb.model)
 end
 
 SolverCore.reset!(solver::TRSolver, model) = SolverCore.reset!(solver)
