diff options
Diffstat (limited to 'Build/source/utils/pmx/pmx-src/libf2c/pow_zz.c')
-rw-r--r-- | Build/source/utils/pmx/pmx-src/libf2c/pow_zz.c | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/Build/source/utils/pmx/pmx-src/libf2c/pow_zz.c b/Build/source/utils/pmx/pmx-src/libf2c/pow_zz.c new file mode 100644 index 00000000000..b5ffd334835 --- /dev/null +++ b/Build/source/utils/pmx/pmx-src/libf2c/pow_zz.c @@ -0,0 +1,29 @@ +#include "f2c.h" + +#ifdef KR_headers +double log(), exp(), cos(), sin(), atan2(), f__cabs(); +VOID pow_zz(r,a,b) doublecomplex *r, *a, *b; +#else +#undef abs +#include "math.h" +#ifdef __cplusplus +extern "C" { +#endif +extern double f__cabs(double,double); +void pow_zz(doublecomplex *r, doublecomplex *a, doublecomplex *b) +#endif +{ +double logr, logi, x, y; + +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); +} +#ifdef __cplusplus +} +#endif |