summaryrefslogtreecommitdiff
path: root/libf2c/libF77/pow_zz.c
diff options
context:
space:
mode:
Diffstat (limited to 'libf2c/libF77/pow_zz.c')
-rw-r--r--libf2c/libF77/pow_zz.c25
1 files changed, 0 insertions, 25 deletions
diff --git a/libf2c/libF77/pow_zz.c b/libf2c/libF77/pow_zz.c
deleted file mode 100644
index d5cfbf3d931..00000000000
--- a/libf2c/libF77/pow_zz.c
+++ /dev/null
@@ -1,25 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-extern double f__cabs (double, double);
-void
-pow_zz (doublecomplex * r, doublecomplex * a, doublecomplex * b)
-{
- double logr, logi, x, y;
-
- if (a->r == 0.0 && a->i == 0.0)
- {
- /* Algorithm below doesn't cope. */
- r->r = r->i = 0.0;
- return;
- }
- logr = log (f__cabs (a->r, a->i));
- logi = atan2 (a->i, a->r);
-
- x = exp (logr * b->r - logi * b->i);
- y = logr * b->i + logi * b->r;
-
- r->r = x * cos (y);
- r->i = x * sin (y);
-}