diff --git a/freegsnke/GSstaticsolver.py b/freegsnke/GSstaticsolver.py
index 61210cd..49522e4 100644
--- a/freegsnke/GSstaticsolver.py
+++ b/freegsnke/GSstaticsolver.py
@@ -938,6 +938,7 @@ def inverse_solve(
 
         check_equilibrium = False
         try:
+            rel_change_full=0
             GS_residual = self.F_function(
                 tokamak_psi=self.tokamak_psi.reshape(-1),
                 plasma_psi=eq.plasma_psi.reshape(-1),