git: 9front

ref: c6fe3b0e8ac7d0ed478b1d48659f68cead87a215
dir: /sys/src/cmd/map/libmap/ccubrt.c/

View raw version 
#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);
}