1 |
+ |
#ifndef lint |
2 |
+ |
static const char RCSid[] = "$Id$"; |
3 |
+ |
#endif |
4 |
|
/* |
5 |
|
* Disk2Square.c |
6 |
|
* |
11 |
|
* Modified interface slightly (G. Ward) |
12 |
|
*/ |
13 |
|
|
14 |
+ |
#define _USE_MATH_DEFINES |
15 |
|
#include <math.h> |
16 |
|
|
17 |
|
/* |
70 |
|
phi = 0.; |
71 |
|
} |
72 |
|
} |
73 |
< |
|
73 |
> |
r *= 0.9999999999999; /* prophylactic against MS sin()/cos() impl. */ |
74 |
|
ds[0] = r * cos(phi); |
75 |
|
ds[1] = r * sin(phi); |
76 |
|
|
102 |
|
a = -(phi - 3*M_PI/2) * b / (M_PI/4); |
103 |
|
} |
104 |
|
|
105 |
< |
sq[0] = (a + 1) * 0.5; |
106 |
< |
sq[1] = (b + 1) * 0.5; |
105 |
> |
sq[0] = a*(0.5/0.9999999999999) + 0.5; |
106 |
> |
sq[1] = b*(0.5/0.9999999999999) + 0.5; |
107 |
|
} |