Upgrade to the 1997/02/26 version.

1 files changed, 23 insertions, 22 deletions

diff --git a/lib/libF77/c_div.c b/lib/libF77/c_div.c
index 0bb56b4..57139a0 100644
--- a/lib/libF77/c_div.c
+++ b/lib/libF77/c_div.c
@@ -9,28 +9,29 @@ extern void sig_die(char*,int);
 void c_div(complex *c, complex *a, complex *b)
 #endif
 {
-double ratio, den;
-double abr, abi;
+	double ratio, den;
+	double abr, abi;
+	double ai = a->i, ar = a->r, bi = b->i, br = b->r;
 
-if( (abr = b->r) < 0.)
-	abr = - abr;
-if( (abi = b->i) < 0.)
-	abi = - abi;
-if( abr <= abi )
-	{
-	if(abi == 0)
-		sig_die("complex division by zero", 1);
-	ratio = (double)b->r / b->i ;
-	den = b->i * (1 + ratio*ratio);
-	c->r = (a->r*ratio + a->i) / den;
-	c->i = (a->i*ratio - a->r) / den;
-	}
+	if( (abr = br) < 0.)
+		abr = - abr;
+	if( (abi = bi) < 0.)
+		abi = - abi;
+	if( abr <= abi )
+		{
+		if(abi == 0)
+			sig_die("complex division by zero", 1);
+		ratio = (double)br / bi ;
+		den = bi * (1 + ratio*ratio);
+		c->r = (ar*ratio + ai) / den;
+		c->i = (ai*ratio - ar) / den;
+		}
 
-else
-	{
-	ratio = (double)b->i / b->r ;
-	den = b->r * (1 + ratio*ratio);
-	c->r = (a->r + a->i*ratio) / den;
-	c->i = (a->i - a->r*ratio) / den;
+	else
+		{
+		ratio = (double)bi / br ;
+		den = br * (1 + ratio*ratio);
+		c->r = (ar + ai*ratio) / den;
+		c->i = (ai - ar*ratio) / den;
+		}
 	}
-}