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