104 |
|
int |
105 |
|
use_symmetry(FVECT vec) |
106 |
|
{ |
107 |
< |
const double phi = get_phi360(vec); |
107 |
> |
double phi = get_phi360(vec); |
108 |
> |
/* because of -0. issue */ |
109 |
> |
while (phi >= 360.) phi -= 360.; |
110 |
> |
while (phi < 0.) phi += 360.; |
111 |
|
|
112 |
|
switch (inp_coverage) { |
113 |
|
case INP_QUAD1|INP_QUAD2|INP_QUAD3|INP_QUAD4: |
228 |
|
void |
229 |
|
ovec_from_pos(FVECT vec, int xpos, int ypos) |
230 |
|
{ |
231 |
< |
double uv[2]; |
231 |
> |
RREAL uv[2]; |
232 |
|
double r2; |
233 |
|
|
234 |
< |
SDsquare2disk(uv, (xpos+.5)/grid_res, (ypos+.5)/grid_res); |
234 |
> |
square2disk(uv, (xpos+.5)/grid_res, (ypos+.5)/grid_res); |
235 |
|
/* uniform hemispherical projection */ |
236 |
|
r2 = uv[0]*uv[0] + uv[1]*uv[1]; |
237 |
|
vec[0] = vec[1] = sqrt(2. - r2); |
244 |
|
void |
245 |
|
pos_from_vec(int pos[2], const FVECT vec) |
246 |
|
{ |
247 |
< |
double sq[2]; /* uniform hemispherical projection */ |
247 |
> |
RREAL sq[2]; /* uniform hemispherical projection */ |
248 |
|
double norm = 1./sqrt(1. + fabs(vec[2])); |
249 |
|
|
250 |
< |
SDdisk2square(sq, vec[0]*norm, vec[1]*norm); |
250 |
> |
disk2square(sq, vec[0]*norm, vec[1]*norm); |
251 |
|
|
252 |
|
pos[0] = (int)(sq[0]*grid_res); |
253 |
|
pos[1] = (int)(sq[1]*grid_res); |