src/cmd/map/libmap/ccubrt.c


1
2
3
4
5
6
7
8
9
10
11
12
13

#include <u.h>
#include <libc.h>
#include "map.h"

void
ccubrt(double zr, double zi, double *wr, double *wi)
{
	double r, theta;
	theta = atan2(zi,zr);
	r = cubrt(hypot(zr,zi));
	*wr = r*cos(theta/3);
	*wi = r*sin(theta/3);
}