 
 
 //  ~~~> New step size is bounded by FacMin <= Hnew/H <= FacMax
-            Fac  = fmin(FacMax,fmax(FacMin,FacSafe/pow(Err,ONE/ros_ELO)));
+            Fac  = fmin(FacMax,fmax(FacMin,FacSafe/pow(Err,ONE/4.0)));
             Hnew = H*Fac;
 
 //  ~~~>  Check the error magnitude and adjust step size

         *
         */
        ros_Integrator_ros4(var, fix, Tstart, Tend, Texit,
                //  Rosenbrock method coefficients
                ros_M, ros_E, ros_A, ros_C, 
                ros_Alpha, ros_Gamma, ros_ELO, ros_NewF, 
                //  Integration parameters
                autonomous, vectorTol, Max_no_steps, 
                roundoff, Hmin, Hmax, Hstart, Hexit, 

        Nsol = 0;
        Nsng = 0;

        /* Copy data from global memory to temporary array */
        /*
         * Optimization note: if we ever have enough constant
