|
|
1.1 root 1: #include "map.h"
2: double cubrt();
3:
4: ccubrt(zr,zi,wr,wi)
5: double zr,zi;
6: double *wr,*wi;
7: {
8: double r, theta;
9: theta = atan2(zi,zr);
10: r = cubrt(hypot(zr,zi));
11: *wr = r*cos(theta/3);
12: *wi = r*sin(theta/3);
13: }
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.