| 16 |
|
#include "rtprocess.h" /* win_popen() */ |
| 17 |
|
#include "paths.h" |
| 18 |
|
#include "ray.h" |
| 19 |
– |
#include "source.h" |
| 19 |
|
#include "ambient.h" |
| 20 |
|
#include "otypes.h" |
| 21 |
+ |
#include "otspecial.h" |
| 22 |
|
#include "rpaint.h" |
| 23 |
|
|
| 24 |
|
extern int psample; /* pixel sample size */ |
| 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 |
|
|
| 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 |
+ |
|
| 420 |
+ |
void |
| 421 |
|
getexposure( /* get new exposure */ |
| 422 |
|
char *s |
| 423 |
|
) |
| 483 |
|
} |
| 484 |
|
|
| 485 |
|
typedef union {int i; double d; COLOR C;} *MyUptr; |
| 441 |
– |
#define FEQ(x,y) (fabs((x)-(y)) <= FTINY) |
| 486 |
|
|
| 487 |
|
int |
| 488 |
|
getparam( /* get variable from user */ |
| 520 |
|
if (sscanf(buf, "%lf", &d0) != 1) |
| 521 |
|
return(0); |
| 522 |
|
} |
| 523 |
< |
if (FEQ(ptr->d, d0)) |
| 523 |
> |
if (FABSEQ(ptr->d, d0)) |
| 524 |
|
return(0); |
| 525 |
|
ptr->d = d0; |
| 526 |
|
break; |
| 552 |
|
if (sscanf(buf, "%lf %lf %lf", &d0, &d1, &d2) != 3) |
| 553 |
|
return(0); |
| 554 |
|
} |
| 555 |
< |
if (FEQ(colval(ptr->C,RED), d0) && |
| 556 |
< |
FEQ(colval(ptr->C,GRN), d1) && |
| 557 |
< |
FEQ(colval(ptr->C,BLU), d2)) |
| 555 |
> |
if (FABSEQ(colval(ptr->C,RED), d0) && |
| 556 |
> |
FABSEQ(colval(ptr->C,GRN), d1) && |
| 557 |
> |
FABSEQ(colval(ptr->C,BLU), d2)) |
| 558 |
|
return(0); |
| 559 |
|
setcolor(ptr->C, d0, d1, d2); |
| 560 |
|
break; |
| 792 |
|
matspec[0] = '\0'; |
| 793 |
|
if (thisray.ro->omod != OVOID) { |
| 794 |
|
mod = objptr(thisray.ro->omod); |
| 795 |
< |
mat = findmaterial(mod); |
| 795 |
> |
mat = findmaterial(thisray.ro); |
| 796 |
|
} |
| 797 |
|
if (thisray.rod < 0.0) |
| 798 |
|
strcpy(matspec, "back of "); |
| 811 |
|
ofun[ino->otype].funame, ino->oname); |
| 812 |
|
(*dev->comout)(buf); |
| 813 |
|
(*dev->comin)(buf, NULL); |
| 814 |
< |
if (thisray.rot >= FHUGE) |
| 814 |
> |
if (thisray.rot >= FHUGE*.99) |
| 815 |
|
(*dev->comout)("at infinity"); |
| 816 |
|
else { |
| 817 |
|
sprintf(buf, "at (%.6g %.6g %.6g) (%.6g)", |
| 818 |
|
thisray.rop[0], thisray.rop[1], |
| 819 |
< |
thisray.rop[2], thisray.rt); |
| 819 |
> |
thisray.rop[2], raydistance(&thisray)); |
| 820 |
|
(*dev->comout)(buf); |
| 821 |
|
} |
| 822 |
|
(*dev->comin)(buf, NULL); |
