| 13 |
|
|
| 14 |
|
#define _USE_MATH_DEFINES |
| 15 |
|
#include <math.h> |
| 16 |
+ |
#include "fvect.h" |
| 17 |
|
|
| 18 |
|
/* |
| 19 |
|
This transforms points on [0,1]^2 to points on unit disk centered at |
| 41 |
|
|
| 42 |
|
/* Map a [0,1]^2 square to a unit radius disk */ |
| 43 |
|
void |
| 44 |
< |
SDsquare2disk(double ds[2], double seedx, double seedy) |
| 44 |
> |
square2disk(RREAL ds[2], double seedx, double seedy) |
| 45 |
|
{ |
| 46 |
|
|
| 47 |
|
double phi, r; |
| 79 |
|
|
| 80 |
|
/* Map point on unit disk to a unit square in [0,1]^2 range */ |
| 81 |
|
void |
| 82 |
< |
SDdisk2square(double sq[2], double diskx, double disky) |
| 82 |
> |
disk2square(RREAL sq[2], double diskx, double disky) |
| 83 |
|
{ |
| 84 |
|
double r = sqrt( diskx*diskx + disky*disky ); |
| 85 |
|
double phi = atan2( disky, diskx ); |