 	int Nconsecutive = 0;
 	int Nsteps = 0;
 	while (t + Roundoff < t_end) {
+		PRINT("Reject:%d\n", ANY(!Reject));
 		if(ANY(!Reject)) {
 			dydt (t, pr, y, F0);
 		}
+		PRINT("LU:%d\n", ANY(!SkipLU));
 		if(ANY(!SkipLU)) { 
 			//need to update Jac/LU
+			PRINT("Jac:%d\n", ANY(!SkipJac));
 			if(ANY(!SkipJac)) {
 				eval_jacob (t, pr, y, A);
 			}
 			RK_Decomp(H, E1, E2, A, ipiv1, ipiv2, &info);
+			PRINT("info:%d\n", ANY(!info));
 			if(ANY(info != 0)) {
 				Nconsecutive += 1;
 				if (Nconsecutive >= 5)


            NewtonDone = (NewtonRate * NewtonIncrement <= NewtonTol);
#ifndef NEWTON_UNROLL
            if (ALL (NewtonDone)) break;
#else //only break if it's at the end of the unroll
            if(ALL (NewtonDone && (NewtonIter + 1) % NewtonIter == 0)) break;

				return;
			}
		}
		PRINT("NotDone:%d\n", ANY(!NewtonDone));
		if(ANY (!NewtonDone)) {
			H = Fac * H;
			Reject = true;
