always impose the algebraic condition

SciML · Feb 28, 2024 · cb01bf3 · cb01bf3
1 parent c79951f
commit cb01bf3
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/perform_step/sdirk_perform_step.jl b/src/perform_step/sdirk_perform_step.jl
@@ -103,7 +103,7 @@ end
 
     integrator.fsallast = f(u, p, t + dt)
 
-    if integrator.opts.adaptive && integrator.success_iter > 0  && integrator.f.mass_matrix !== I
+    if integrator.opts.adaptive && integrator.f.mass_matrix !== I
         atmp = @. ifelse(!integrator.differential_vars, integrator.fsallast, false) ./
                     integrator.opts.abstol
         integrator.EEst += integrator.opts.internalnorm(atmp, t)
@@ -160,7 +160,7 @@ end
     integrator.stats.nf += 1
     f(integrator.fsallast, u, p, t + dt)
 
-    if integrator.opts.adaptive && integrator.success_iter > 0 && integrator.f.mass_matrix !== I
+    if integrator.opts.adaptive && integrator.f.mass_matrix !== I
         @.. broadcast=false atmp=ifelse(cache.algebraic_vars, integrator.fsallast, false) /
                                  integrator.opts.abstol
         integrator.EEst += integrator.opts.internalnorm(atmp, t)