The R Project SVN – R – Diff – /trunk/src/main/arithmetic.c

     R_NaReal = R_ValueOfNA();
     R_PosInf = 1.0/R_Zero_Hack;
     R_NegInf = -1.0/R_Zero_Hack;
+}
+/* Keep these two in step */
 static double myfmod(double x1, double x2)
+{
     double q = x1 / x2, tmp;
     if (x2 == 0.0) return R_NaN;
 	warning(_("probable complete loss of accuracy in modulus"));
     q = floor(tmp/x2);
     return tmp - q * x2;
+}
+static double myfloor(double x1, double x2)
+{
+    double q = x1 / x2, tmp;
+    if (x2 == 0.0) return q;
+    tmp = x1 - floor(q) * x2;
+    return floor(q) + floor(tmp/x2);
+}
 #ifdef HAVE_WORKING_LOG
 # define R_log	log
 #else
 double R_log(double x) { return(x > 0 ? log(x) : x < 0 ? R_NaN : R_NegInf); }
 	    REAL(ans)[i] = myfmod(REAL(s1)[i1], REAL(s2)[i2]);
+	}
 	break;
     case IDIVOP:
 	mod_iterate(n1, n2, i1, i2) {
-	    REAL(ans)[i] = floor(REAL(s1)[i1] / REAL(s2)[i2]);
+	    REAL(ans)[i] = myfloor(REAL(s1)[i1], REAL(s2)[i2]);
+	}
 	break;
+    }
     /* quick return if there are no attributes */

The R Project SVN R