Update of /project/oct/cvsroot/oct In directory clnet:/tmp/cvs-serv27836
Modified Files: qd-fun.lisp Log Message: Compare against pi/4 quad-double, not double-float.
--- /project/oct/cvsroot/oct/qd-fun.lisp 2007/10/17 03:44:12 1.89 +++ /project/oct/cvsroot/oct/qd-fun.lisp 2007/10/18 14:38:56 1.90 @@ -678,7 +678,7 @@ ;; We want the remainder part to be <= pi/4 because the trig ;; functions want that. So if the fraction is too big, adjust ;; it, and mod value. - (if (<= (abs (qd-0 f)) (/ pi 4)) + (if (qd-<= (abs-qd f) +qd-pi/4+) (values mod f) (values (mod (1+ mod) 4) (sub-qd f +qd-pi/2+))))))
@@ -723,7 +723,7 @@ +qd-pi/2+))) ;; If we did this right, (ash prod (- frac-bits)) should be 2 ;; bits long at most. - (if (<= (abs (qd-0 f)) (/ pi 4)) + (if (qd-<= (abs-qd f) +qd-pi/4+) (values mod f) (values (mod (1+ mod) 4) (sub-qd f +qd-pi/2+))))))