1 files changed, 17 insertions, 0 deletions
diff --git a/lib/libF77/c_log.c b/lib/libF77/c_log.c
new file mode 100644
index 0000000..6ac990c
--- /dev/null
+++ b/lib/libF77/c_log.c
@@ -0,0 +1,17 @@
+#include "f2c.h"
+
+#ifdef KR_headers
+extern double log(), f__cabs(), atan2();
+VOID c_log(r, z) complex *r, *z;
+#else
+#undef abs
+#include "math.h"
+extern double f__cabs(double, double);
+
+void c_log(complex *r, complex *z)
+#endif
+{
+	double zi;
+	r->i = atan2(zi = z->i, z->r);
+	r->r = log( f__cabs(z->r, zi) );
+	}