25 |
|
* from his, the first two interfaces and the style of some comments are very |
26 |
|
* much derived from his work. */ |
27 |
|
|
28 |
< |
|
28 |
> |
#include <stdlib.h> |
29 |
|
#include "hilbert.h" |
30 |
|
|
31 |
|
/* implementation of the hilbert functions */ |
1256 |
|
printf("How many points? "); |
1257 |
|
scanf("%d", &n); |
1258 |
|
a = (double*) malloc(2*n*sizeof(double)); |
1259 |
< |
for (i = 0; i < n; ++i) |
1260 |
< |
a[2*i] = drand48()-0.5, a[2*i+1] = drand48()-0.5; |
1261 |
< |
|
1259 |
> |
for (i = 0; i < n; ++i) { |
1260 |
> |
a[2*i] = rand()*(1./RAND_MAX) - 0.5; |
1261 |
> |
a[2*i+1] = rand()*(1./RAND_MAX) - 0.5; |
1262 |
> |
} |
1263 |
|
qsort(a, n, 2*sizeof(double), cmp); |
1264 |
|
|
1265 |
|
for (i = 0; i < n; ++i) |