Raymond Toy pushed to branch issue-456-more-accurate-complex-div at cmucl / cmucl
Commits:
-
6f6a8758
by Raymond Toy at 2025-12-12T09:13:22-08:00
1 changed file:
Changes:
| ... | ... | @@ -632,14 +632,8 @@ |
| 632 | 632 | ;; = (a + b*r)/(c + d*r).
|
| 633 | 633 | ;;
|
| 634 | 634 | ;; Thus tt = (c + d*r).
|
| 635 | - #+nil
|
|
| 636 | - (progn
|
|
| 637 | - (format t "a,b,c,d = ~A ~A ~A ~A~%" a b c d)
|
|
| 638 | - (format t " r, tt = ~A ~A~%" r tt))
|
|
| 639 | 635 | (cond ((>= (abs r) +rmin+)
|
| 640 | 636 | (let ((br (* b r)))
|
| 641 | - #+nil
|
|
| 642 | - (format t "br = ~A~%" br)
|
|
| 643 | 637 | (if (/= br 0)
|
| 644 | 638 | (/ (+ a br) tt)
|
| 645 | 639 | ;; b*r underflows. Instead, compute
|
| ... | ... | @@ -654,11 +648,6 @@ |
| 654 | 648 | (t
|
| 655 | 649 | ;; r = 0 so d is very tiny compared to c.
|
| 656 | 650 | ;;
|
| 657 | - ;; (a + b*r)/tt = (a + b*(d/c))/tt
|
|
| 658 | - #+nil
|
|
| 659 | - (progn
|
|
| 660 | - (format t "r = 0~%")
|
|
| 661 | - (format t "a*tt = ~A~%" (* a tt)))
|
|
| 662 | 651 | (/ (+ a (* d (/ b c)))
|
| 663 | 652 | tt))))
|
| 664 | 653 | (robust-subinternal (a b c d)
|