| 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 |
|
} |