date and time created 83/01/21 11:08:07 by dlw

SCCS-vsn: usr.bin/f77/libF77/c_div.c 1.1

diff --git a/usr/src/usr.bin/f77/libF77/c_div.c b/usr/src/usr.bin/f77/libF77/c_div.c
new file mode 100644 (file)
index 0000000..7147514
--- /dev/null
+++ b/usr/src/usr.bin/f77/libF77/c_div.c
@@ -0,0 +1,35 @@
+/*
+ *     "@(#)c_div.c    1.1"
+ */
+
+#include "complex"
+
+c_div(c, a, b)
+complex *a, *b, *c;
+{
+double ratio, den;
+double abr, abi;
+
+if( (abr = b->real) < 0.)
+       abr = - abr;
+if( (abi = b->imag) < 0.)
+       abi = - abi;
+if( abr <= abi )
+       {
+       if(abi == 0)
+               abort(); /* fatal("complex division by zero"); */
+       ratio = b->real / b->imag ;
+       den = b->imag * (1 + ratio*ratio);
+       c->real = (a->real*ratio + a->imag) / den;
+       c->imag = (a->imag*ratio - a->real) / den;
+       }
+
+else
+       {
+       ratio = b->imag / b->real ;
+       den = b->real * (1 + ratio*ratio);
+       c->real = (a->real + a->imag*ratio) / den;
+       c->imag = (a->imag - a->real*ratio) / den;
+       }
+
+}