  */
 __device__ void integrate (const double t_start, const double t_end, const double pr, double* y) {
 	double Hmin = 0;
-	double Hmax = fabs(t_end - t_start);
 	double Hold = 0;
 #ifdef Gustafsson
 	double Hacc = 0;

#ifndef NEWTON_UNROLL
            if (NewtonDone) break;
#else //only break if it's at the end of the unroll
            if(NewtonDone && (NewtonIter + 1) % NewtonIter == 0) break;
#endif
            if (NewtonIter >= NewtonMaxit)
            {

			}
			scale(y, y0, sc);
			safe_memcpy(y0, y);
			Hnew = fmin(fmax(Hnew, Hmin), t_end - t);
			if (Reject) {
				Hnew = fmin(Hnew, H);
			}
