The R Project SVN – R – Diff – /branches/R-exp-R5/src/main/seq.c

 /*
    'along' has to be used on an unevaluated argument, and evalList
    tries to evaluate language objects.
  */
-#define FEPS 1e-7
+#define FEPS 1e-10
 /* to match seq.default */
 SEXP attribute_hidden do_seq(SEXP call, SEXP op, SEXP args, SEXP rho)
+{
     SEXP ans = R_NilValue /* -Wall */, ap, tmp, from, to, by, len, along;
     int i, nargs = length(args), lf, lout = NA_INTEGER;
+	    }
 	    if(n > (double) INT_MAX)
 		errorcall(call, _("'by' argument is much too small"));
 	    if(n < - FEPS)
 		errorcall(call, _("wrong sign in 'by' argument"));
-	    nn = (int)(n + FEPS);
 	    if(TYPEOF(from) == INTSXP &&
 	       TYPEOF(to) == INTSXP &&
 	       TYPEOF(by) == INTSXP) {
 		int *ia, ifrom = asInteger(from), iby = asInteger(by);
+		/* With the current limits on integers and FEPS
+		   reduced below 1/INT_MAX this is the same as the
+		   next, so this is future-proofing against longer integers.
+		*/
+		nn = (int)n;
-		/* seq.default gives integer result from from + (0:n)*by */
+		/* seq.default gives integer result from
+		   from + (0:n)*by
+		*/
 		ans = allocVector(INTSXP, nn+1);
 		ia = INTEGER(ans);
 		for(i = 0; i <= nn; i++)
 		    ia[i] = ifrom + i * iby;
 	    } else {
+		nn = (int)(n + FEPS);
 		ans = allocVector(REALSXP, nn+1);
 		ra = REAL(ans);
 		for(i = 0; i <= nn; i++)
 		    ra[i] = rfrom + i * rby;
 		/* Added in 2.9.0 */

The R Project SVN R