From 4c6e36a6e37ff5bc4da7655e125bc4e1c091dc2f Mon Sep 17 00:00:00 2001
From: Kevin Ryde <user42@zip.com.au>
Date: Fri, 29 Apr 2005 22:59:22 +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 08483d8f7..c60920b80 100644
--- a/libguile/numbers.c
+++ b/libguile/numbers.c
@@ -4622,7 +4622,7 @@ scm_i_divide (SCM x, SCM y, int inexact)
 	{
 	  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);
@@ -4694,7 +4694,7 @@ scm_i_divide (SCM x, SCM y, int inexact)
 	  {
 	    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);
@@ -4899,7 +4899,7 @@ scm_i_divide (SCM x, SCM y, int inexact)
 	{
 	  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);