fix of divmod, divmodf

Alexander-Makaryev · Alexander-Makaryev · commit 9ab852e23efa · 2020-06-24T05:04:23.000-05:00
diff --git a/mkl_umath/src/loops_intel.c.src b/mkl_umath/src/loops_intel.c.src
@@ -122,6 +122,7 @@ static inline npy_float spacingf(npy_float x) {
 }
 
 static inline npy_double divmod(npy_double x1, npy_double x2, npy_double *mod) {
+    fesetround(FE_DOWNWARD);
     npy_double q = rint(x1 / x2);
 
     *mod = x1 - q * x2;
@@ -130,6 +131,7 @@ static inline npy_double divmod(npy_double x1, npy_double x2, npy_double *mod) {
 }
 
 static inline npy_float divmodf(npy_float x1, npy_float x2, npy_float *mod) {
+    fesetround(FE_DOWNWARD);
     npy_float q = rintf(x1 / x2);
 
     *mod = x1 - q * x2;

Original file line number	Diff line number	Diff line change
`@@ -122,6 +122,7 @@ static inline npy_float spacingf(npy_float x) {`
`122`	`122`	`}`
`123`	`123`
`124`	`124`	`static inline npy_double divmod(npy_double x1, npy_double x2, npy_double *mod) {`
	`125`	`+ fesetround(FE_DOWNWARD);`
`125`	`126`	`npy_double q = rint(x1 / x2);`
`126`	`127`
`127`	`128`	`mod = x1 - q x2;`
`@@ -130,6 +131,7 @@ static inline npy_double divmod(npy_double x1, npy_double x2, npy_double *mod) {`
`130`	`131`	`}`
`131`	`132`
`132`	`133`	`static inline npy_float divmodf(npy_float x1, npy_float x2, npy_float *mod) {`
	`134`	`+ fesetround(FE_DOWNWARD);`
`133`	`135`	`npy_float q = rintf(x1 / x2);`
`134`	`136`
`135`	`137`	`mod = x1 - q x2;`