 			if (NewtonIter > 0) 
 			{
 				Theta = NewtonIncrement / NewtonIncrementOld;
+				PRINT("Theta:%d\n", ANY(Theta >= 0.99));
 				if(ANY(Theta >= 0.99)) //! Non-convergence of Newton: Theta too large
 					break;
 				else
 					NewtonRate = Theta / (ONE - Theta);
 				//Predict error at the end of Newton process 
 				double NewtonPredictedErr = (NewtonIncrement * pow(Theta, (NewtonMaxit - NewtonIter - 1))) / (ONE - Theta);
+				PRINT("Pred:%d\n", ANY(NewtonPredictedErr >= NewtonTol));
 				if(ANY (NewtonPredictedErr >= NewtonTol)) {
 					//Non-convergence of Newton: predicted error too large
 					double Qnewton = fmin(10.0, NewtonPredictedErr / NewtonTol);

		sc[i] = 1.0 / (ATOL + fabs(y0[i]) * RTOL);
	}
}

#ifdef NEWTON_UNROLL
	#define NEWTON_UNROLLER	#pragma unroll NEWTON_UNROLL

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