1 files changed, 37 insertions, 0 deletions
diff --git a/lib/libF77/c_div.c b/lib/libF77/c_div.c
new file mode 100644
index 0000000..57139a0
--- /dev/null
+++ b/lib/libF77/c_div.c
@@ -0,0 +1,37 @@
+#include "f2c.h"
+
+#ifdef KR_headers
+extern VOID sig_die();
+VOID c_div(c, a, b)
+complex *a, *b, *c;
+#else
+extern void sig_die(char*,int);
+void c_div(complex *c, complex *a, complex *b)
+#endif
+{
+	double ratio, den;
+	double abr, abi;
+	double ai = a->i, ar = a->r, bi = b->i, br = b->r;
+
+	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)bi / br ;
+		den = br * (1 + ratio*ratio);
+		c->r = (ar + ai*ratio) / den;
+		c->i = (ai - ar*ratio) / den;
+		}
+	}