summaryrefslogtreecommitdiffstats
path: root/gnu/lib/libf2c/libF77/pow_zi.c
diff options
context:
space:
mode:
Diffstat (limited to 'gnu/lib/libf2c/libF77/pow_zi.c')
-rw-r--r--gnu/lib/libf2c/libF77/pow_zi.c50
1 files changed, 0 insertions, 50 deletions
diff --git a/gnu/lib/libf2c/libF77/pow_zi.c b/gnu/lib/libf2c/libF77/pow_zi.c
deleted file mode 100644
index 214db3d7a0f..00000000000
--- a/gnu/lib/libf2c/libF77/pow_zi.c
+++ /dev/null
@@ -1,50 +0,0 @@
-#include "f2c.h"
-
-extern void z_div (doublecomplex *, doublecomplex *, doublecomplex *);
-void
-pow_zi (doublecomplex * p, doublecomplex * a, integer * b) /* p = a**b */
-{
- integer n;
- unsigned long u;
- double t;
- doublecomplex q, x;
- static doublecomplex one = { 1.0, 0.0 };
-
- n = *b;
- q.r = 1;
- q.i = 0;
-
- if (n == 0)
- goto done;
- if (n < 0)
- {
- n = -n;
- z_div (&x, &one, a);
- }
- else
- {
- x.r = a->r;
- x.i = a->i;
- }
-
- for (u = n;;)
- {
- if (u & 01)
- {
- t = q.r * x.r - q.i * x.i;
- q.i = q.r * x.i + q.i * x.r;
- q.r = t;
- }
- if (u >>= 1)
- {
- t = x.r * x.r - x.i * x.i;
- x.i = 2 * x.r * x.i;
- x.r = t;
- }
- else
- break;
- }
-done:
- p->i = q.i;
- p->r = q.r;
-}