Use isapprox for calculating convergence

SECQUOIA · Jan 9, 2025 · b56a6ee · b56a6ee
1 parent c2739fd
commit b56a6ee
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/src/solver.jl b/src/solver.jl
@@ -22,11 +22,10 @@ mutable struct ConvergenceObserver <: AbstractObserver
 end
 
 function ITensorMPS.checkdone!(o::ConvergenceObserver; energy, sweep, kwargs...)
-  abs_err = abs(energy - o.previous_energy)
-  rel_err = abs_err / abs(energy)
+  converged = isapprox(energy, o.previous_energy; atol = o.atol, rtol = o.rtol)
   o.previous_energy = energy
 
-  return abs_err < o.atol || rel_err < o.rtol
+  return converged
 end