 			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
+				if(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)) {
+				if(NewtonPredictedErr >= NewtonTol) {
 					//Non-convergence of Newton: predicted error too large
 					double Qnewton = fmin(10.0, NewtonPredictedErr / NewtonTol);
                     Fac = 0.8 * pow(Qnewton, -ONE/((double)(NewtonMaxit-NewtonIter)));

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

#ifdef DEBUG
#include <stdio.h>
#define PRINT(X, Y)  \
	if (threadIdx.x == 0) {\
		printf("Block:%d\t", blockIdx.x); \
		printf(X, Y); \
	} \
	else \
		(void)0
#else
#define PRINT(X, Y)
#endif

#ifdef NEWTON_UNROLL
	#define NEWTON_UNROLLER	#pragma unroll NEWTON_UNROLL

				return;
			}
		}
		if(!NewtonDone) {
			H = Fac * H;
			Reject = true;
			SkipJac = true;
