diff options
Diffstat (limited to 'libf2c/libF77/pow_zz.c')
-rw-r--r-- | libf2c/libF77/pow_zz.c | 25 |
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); -} |