libf2c/libF77/pow_zz.c

   1 #include "f2c.h"
   2
   3 #undef abs
   4 #include <math.h>
   5 extern double f__cabs (double, double);
   6 void
   7 pow_zz (doublecomplex * r, doublecomplex * a, doublecomplex * b)
   8 {
   9   double logr, logi, x, y;
  10
  11   logr = log (f__cabs (a->r, a->i));
  12   logi = atan2 (a->i, a->r);
  13
  14   x = exp (logr * b->r - logi * b->i);
  15   y = logr * b->i + logi * b->r;
  16
  17   r->r = x * cos (y);
  18   r->i = x * sin (y);
  19 }