From 1c9c30b6b2d49d2e4816d90bb37685e188c566b4 Mon Sep 17 00:00:00 2001
From: Kevin Ryde <user42@zip.com.au>
Date: Fri, 29 Apr 2005 22:54:35 +0000
Subject: [PATCH] (scm_divide): Correction to 1/complex and <any>/complex, need
 to test abs(re)<abs(im) for choice of cases, otherwise divide by zero when
 re==0 and im<0.  Reported by Jean Crepeau.

---
 libguile/numbers.c | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/libguile/numbers.c b/libguile/numbers.c
index 471be67a9..fad41074c 100644
--- a/libguile/numbers.c
+++ b/libguile/numbers.c
@@ -3738,7 +3738,7 @@ scm_divide (SCM x, SCM y)
     } else if (SCM_COMPLEXP (x)) {
       double r = SCM_COMPLEX_REAL (x);
       double i = SCM_COMPLEX_IMAG (x);
-      if (r <= i) {
+      if (fabs(r) <= fabs(i)) {
 	double t = r / i;
 	double d = i * (1.0 + t * t);
 	return scm_make_complex (t / d, -1.0 / d);
@@ -3782,7 +3782,7 @@ scm_divide (SCM x, SCM y)
       {
 	double r = SCM_COMPLEX_REAL (y);
 	double i = SCM_COMPLEX_IMAG (y);
-	if (r <= i) {
+        if (fabs(r) <= fabs(i)) {
 	  double t = r / i;
 	  double d = i * (1.0 + t * t);
 	  return scm_make_complex ((a * t) / d,  -a / d);
@@ -3873,7 +3873,7 @@ scm_divide (SCM x, SCM y)
     } else if (SCM_COMPLEXP (y)) {
       double ry = SCM_COMPLEX_REAL (y);
       double iy = SCM_COMPLEX_IMAG (y);
-      if (ry <= iy) {
+      if (fabs(ry) <= fabs(iy)) {
 	double t = ry / iy;
 	double d = iy * (1.0 + t * t);
 	return scm_make_complex ((rx * t + ix) / d, (ix * t - rx) / d);