rem_pio2: actually return medium value for x ~<= 5pi/4

m1el · m1el · commit f42e57077266 · 2019-05-13T12:14:03.000+03:00
diff --git a/src/math/rem_pio2.rs b/src/math/rem_pio2.rs
@@ -85,7 +85,7 @@ pub(crate) fn rem_pio2(x: f64) -> (i32, f64, f64) {
         /* |x| ~<= 5pi/4 */
         if (ix & 0xfffff) == 0x921fb {
             /* |x| ~= pi/2 or 2pi/2 */
-            medium(x, ix); /* cancellation -- use medium case */
+            return medium(x, ix); /* cancellation -- use medium case */
         }
         if ix <= 0x4002d97c {
             /* |x| ~<= 3pi/4 */

Original file line number	Diff line number	Diff line change
`@@ -85,7 +85,7 @@ pub(crate) fn rem_pio2(x: f64) -> (i32, f64, f64) {`
`85`	`85`	`/* \|x\| ~<= 5pi/4 */`
`86`	`86`	`if (ix & 0xfffff) == 0x921fb {`
`87`	`87`	`/* \|x\| ~= pi/2 or 2pi/2 */`
`88`		`- medium(x, ix); /* cancellation -- use medium case */`
	`88`	`+ return medium(x, ix); /* cancellation -- use medium case */`
`89`	`89`	`}`
`90`	`90`	`if ix <= 0x4002d97c {`
`91`	`91`	`/* \|x\| ~<= 3pi/4 */`