Fix bug when number of (solve steps % options.save_every) == 0 (#83)

* Fix bug when number of solve steps % save_every == 0

tdgl/solver/runner.py

@@ -449,6 +449,6 @@ def save_step(step):
                         self.logger.warning(msg.format("Cancelling"))
                         cancelled = True
                         break
-            if save:
+            if save and (i % self.options.save_every):
                 save_step(i)
             return not cancelled

tdgl/solver/solver.py

@@ -632,7 +632,7 @@ def update(
                 (current_A_applied - prev_A_applied) / dt,
                 self.normalized_directions,
             )
-            if xp.any(xp.absolute(dA_dt) > 0):
+            if not xp.allclose(current_A_applied, self.current_A_applied):
                 # Update the link exponents only if the applied vector potential
                 # has actually changed.
                 operators.set_link_exponents(current_A_applied)
@@ -642,10 +642,10 @@ def update(
         self.current_A_applied = current_A_applied
 
         # Update the value of epsilon
-        epsilon = self.epsilon
         if self.dynamic_epsilon:
-            epsilon = self.epsilon = self.update_epsilon(time)
+            self.epsilon = self.update_epsilon(time)
 
+        epsilon = self.epsilon
         old_sq_psi = xp.absolute(psi) ** 2
         screening_error = np.inf
         A_induced_vals = [A_induced]

tdgl/version.py

@@ -1,4 +1,4 @@
-__version_info__ = (0, 8, 2)
+__version_info__ = (0, 8, 3)
 __version__ = ".".join(map(str, __version_info__))