
5 Nov
2007
5 Nov
'07
4 p.m.
Update of /project/oct/cvsroot/oct In directory clnet:/tmp/cvs-serv7505 Modified Files: Tag: THREE-ARG-BRANCH qd.lisp Log Message: Cosmetic --- /project/oct/cvsroot/oct/qd.lisp 2007/11/04 16:31:45 1.60.2.4 +++ /project/oct/cvsroot/oct/qd.lisp 2007/11/05 16:00:39 1.60.2.5 @@ -870,7 +870,8 @@ (declaim (inline invert-qd)) -(defun invert-qd(v) ;; a quartic newton iteration for 1/v +(defun invert-qd (v) + ;; a quartic newton iteration for 1/v ;; to invert v, start with a good guess, x. ;; let h= 1-v*x ;; h is small ;; return x+ x*(h+h^2+h^3) . compute h3 in double-float