| 26 |
|
|
| 27 |
|
#define CTRL(c) ((c)-'@') |
| 28 |
|
|
| 29 |
< |
#ifdef SMLFLT |
| 30 |
< |
#define sscanvec(s,v) (sscanf(s,"%f %f %f",v,v+1,v+2)==3) |
| 31 |
< |
#else |
| 32 |
< |
#define sscanvec(s,v) (sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3) |
| 33 |
< |
#endif |
| 29 |
> |
#define sscanvec(s,v) (sscanf(s,FVFORMAT,v,v+1,v+2)==3) |
| 30 |
|
|
| 31 |
|
extern char rifname[128]; /* rad input file name */ |
| 32 |
|
|
| 333 |
|
) |
| 334 |
|
{ |
| 335 |
|
VIEW nv = ourview; |
| 340 |
– |
FVECT v1; |
| 336 |
|
double angle, elev, zfact; |
| 337 |
|
|
| 338 |
|
elev = 0.0; zfact = 1.0; |
| 341 |
|
return; |
| 342 |
|
} |
| 343 |
|
spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.)); |
| 344 |
< |
if (elev != 0.0) { |
| 345 |
< |
fcross(v1, nv.vdir, ourview.vup); |
| 346 |
< |
normalize(v1); |
| 352 |
< |
spinvector(nv.vdir, nv.vdir, v1, elev*(PI/180.)); |
| 353 |
< |
} |
| 344 |
> |
if (elev != 0.0) |
| 345 |
> |
geodesic(nv.vdir, nv.vdir, nv.vup, elev*(PI/180.), GEOD_RAD); |
| 346 |
> |
|
| 347 |
|
zoomview(&nv, zfact); |
| 348 |
|
newview(&nv); |
| 349 |
|
} |
| 365 |
|
if (getinterest(sskip2(s,2), 0, vc, &mag) < 0) |
| 366 |
|
return; |
| 367 |
|
moveview(angle, elev, mag, vc); |
| 368 |
+ |
} |
| 369 |
+ |
|
| 370 |
+ |
|
| 371 |
+ |
void |
| 372 |
+ |
getorigin( /* origin viewpoint */ |
| 373 |
+ |
char *s |
| 374 |
+ |
) |
| 375 |
+ |
{ |
| 376 |
+ |
VIEW nv = ourview; |
| 377 |
+ |
double d; |
| 378 |
+ |
/* get new view origin */ |
| 379 |
+ |
if (sscanf(s, "%lf %lf", &d, &d) == 1) { |
| 380 |
+ |
/* just moving some distance */ |
| 381 |
+ |
VSUM(nv.vp, nv.vp, nv.vdir, d); |
| 382 |
+ |
} else if (!sscanvec(s, nv.vp)) { |
| 383 |
+ |
int x, y; /* need to pick origin */ |
| 384 |
+ |
RAY thisray; |
| 385 |
+ |
if (dev->getcur == NULL) |
| 386 |
+ |
return; |
| 387 |
+ |
(*dev->comout)("Pick point on surface for new origin\n"); |
| 388 |
+ |
if ((*dev->getcur)(&x, &y) == ABORT) |
| 389 |
+ |
return; |
| 390 |
+ |
if ((thisray.rmax = viewray(thisray.rorg, thisray.rdir, |
| 391 |
+ |
&ourview, (x+.5)/hresolu, (y+.5)/vresolu)) < -FTINY) { |
| 392 |
+ |
error(COMMAND, "not on image"); |
| 393 |
+ |
return; |
| 394 |
+ |
} |
| 395 |
+ |
rayorigin(&thisray, PRIMARY, NULL, NULL); |
| 396 |
+ |
if (!localhit(&thisray, &thescene)) { |
| 397 |
+ |
error(COMMAND, "not a local object"); |
| 398 |
+ |
return; |
| 399 |
+ |
} |
| 400 |
+ |
if (thisray.rod < 0.0) /* don't look through other side */ |
| 401 |
+ |
flipsurface(&thisray); |
| 402 |
+ |
VSUM(nv.vp, thisray.rop, thisray.ron, 20.0*FTINY); |
| 403 |
+ |
VCOPY(nv.vdir, thisray.ron); |
| 404 |
+ |
} else if (!sscanvec(sskip2(s,3), nv.vdir) || normalize(nv.vdir) == 0.0) |
| 405 |
+ |
VCOPY(nv.vdir, ourview.vdir); |
| 406 |
+ |
|
| 407 |
+ |
d = DOT(nv.vdir, nv.vup); /* need different up vector? */ |
| 408 |
+ |
if (d*d >= 1.-2.*FTINY) { |
| 409 |
+ |
int i; |
| 410 |
+ |
nv.vup[0] = nv.vup[1] = nv.vup[2] = 0.0; |
| 411 |
+ |
for (i = 3; i--; ) |
| 412 |
+ |
if (nv.vdir[i]*nv.vdir[i] < 0.34) |
| 413 |
+ |
break; |
| 414 |
+ |
nv.vup[i] = 1.; |
| 415 |
+ |
} |
| 416 |
+ |
newview(&nv); |
| 417 |
|
} |
| 418 |
|
|
| 419 |
|
|
