lib/libm/cpowf.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

#include "__complex.h"

float complex cpowf(float complex a, float complex z)
{
	float x, y, r, theta, absa, arga;
	x = crealf(z);
	y = cimagf(z);
	absa = cabsf(a);
	if (absa == 0.0f) {
		return (0.0f + 0.0f * I);
	}
	arga = cargf(a);
	r = powf(absa, x);
	theta = x * arga;
	if (y != 0.0f) {
		r = r * expf(-y * arga);
		theta = theta + y * logf(absa);
	}
	return r * cosf(theta) + (r * sinf(theta)) * I;
}