move the check

Author: zingale

integration/VODE/vode_dvnlsd.H

@@ -115,29 +115,6 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate)
             // sometime VODE goes way off tangent.  If these mass fractions
             // are really bad, then let's just bail now
 
-            if (integrator_rp::do_corrector_validation) {
-#ifdef SDC
-                const amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO];
-                const amrex::Real thresh = species_failure_tolerance * rho_current;
-#else
-                const amrex::Real thresh = species_failure_tolerance;
-#endif
-                bool fail{};
-                for (int i = 1; i <= NumSpec; ++i) {
-                    if (vstate.y(i) < -thresh) {
-                        fail = true;
-                        break;
-                    }
-                }
-
-                if (fail) {
-                    NFLAG = -1;
-                    vstate.ICF = 2;
-                    vstate.IPUP = 1;
-                    return ACNRM;
-                }
-            }
-
 #ifdef NEW_NETWORK_IMPLEMENTATION
             RHS::dgesl(vstate.jac, vstate.y);
 #else
@@ -171,6 +148,29 @@ amrex::Real dvnlsd (int& NFLAG, BurnT& state, DvodeT& vstate)
                 vstate.y(i) = vstate.yh(i,1) + vstate.acor(i);
             }
 
+            if (integrator_rp::do_corrector_validation) {
+#ifdef SDC
+                const amrex::Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO];
+                const amrex::Real thresh = species_failure_tolerance * rho_current;
+#else
+                const amrex::Real thresh = species_failure_tolerance;
+#endif
+                bool fail{};
+                for (int i = 1; i <= NumSpec; ++i) {
+                    if (vstate.y(i) < -thresh) {
+                        fail = true;
+                        break;
+                    }
+                }
+
+                if (fail) {
+                    NFLAG = -1;
+                    vstate.ICF = 2;
+                    vstate.IPUP = 1;
+                    return ACNRM;
+                }
+            }
+
             // Test for convergence. If M > 0, an estimate of the convergence
             // rate constant is stored in CRATE, and this is used in the test.