From 9a3412eea8a2aa2daafa031f2e94f7049c1de202 Mon Sep 17 00:00:00 2001 From: Danny Smith Date: Sat, 25 Dec 2004 23:56:19 +0000 Subject: * mingwex/complex/(cabsf.c cacosf.c cacoshf.c cargf.c casinf.c casinhf.c catanf.c catanhf.c ccosf.c ccoshf.c cexpf.c cimagf.c clogf.c cpowf.c cprojf.c crealf.c csinf.c csinhf.c csqrtf.c ctanf.c ctanhf.c): New files. * mingwex/Makefile.in (COMPLEX_DISTFILES): Adjust. (COMPLEX_OBJS(: Adjust. * include/complex.h (cabsf, cacosf, cacoshf, cargf, casinf. casinhf, catanf, catanhf, ccosf, ccoshf, cexpf, cimagf, clogf, cpowf, cprojf, crealf, csinf, csinhf, csqrtf, ctanf, ctanhf): Declare. --- winsup/mingw/mingwex/complex/cpowf.c | 43 ++++++++++++++++++++++++++++++++++++ 1 file changed, 43 insertions(+) create mode 100755 winsup/mingw/mingwex/complex/cpowf.c (limited to 'winsup/mingw/mingwex/complex/cpowf.c') diff --git a/winsup/mingw/mingwex/complex/cpowf.c b/winsup/mingw/mingwex/complex/cpowf.c new file mode 100755 index 000000000..eeeed1b7f --- /dev/null +++ b/winsup/mingw/mingwex/complex/cpowf.c @@ -0,0 +1,43 @@ +/* cpowf.c */ +/* + Contributed by Danny Smith + 2004-12-24 +*/ + +/* cpow(X, Y) = cexp(X * clog(Y)) */ + +#include +#include + +float complex cpowf (float complex X, float complex Y) +{ + float complex Res; + float i; + float r = _hypot (__real__ X, __imag__ X); + if (r == 0.0f) + { + __real__ Res = __imag__ Res = 0.0; + } + else + { + float rho; + float theta; + i = cargf (X); + theta = i * __real__ Y; + + if (__imag__ Y == 0.0f) + /* This gives slightly more accurate results in these cases. */ + rho = powf (r, __real__ Y); + else + { + r = logf (r); + /* rearrangement of cexp(X * clog(Y)) */ + theta += r * __imag__ Y; + rho = expf (r * __real__ Y - i * __imag__ Y); + } + + __real__ Res = rho * cosf (theta); + __imag__ Res = rho * sinf (theta); + } + return Res; +} -- cgit v1.2.3