| 1 |
< |
/* Copyright (c) 1991 Regents of the University of California */ |
| 1 |
> |
/* Copyright (c) 1992 Regents of the University of California */ |
| 2 |
|
|
| 3 |
|
#ifndef lint |
| 4 |
|
static char SCCSid[] = "$SunId$ LBL"; |
| 174 |
|
for (i = 0; i < 3; i++) |
| 175 |
|
src->ss[SU][i] = src->ss[SV][i] = src->ss[SW][i] = 0.0; |
| 176 |
|
for (i = 0; i < 3; i++) |
| 177 |
< |
src->ss[i][i] = .886227 * so->oargs.farg[3]; |
| 177 |
> |
src->ss[i][i] = .7236 * so->oargs.farg[3]; |
| 178 |
|
} |
| 179 |
|
|
| 180 |
|
|
| 212 |
|
co = getcone(so, 0); |
| 213 |
|
if (CO_R0(co) > .2*co->al) /* heuristic constraint */ |
| 214 |
|
objerror(so, WARNING, "source aspect too small"); |
| 215 |
+ |
src->sflags |= SCYL; |
| 216 |
|
for (i = 0; i < 3; i++) |
| 217 |
|
src->sloc[i] = .5 * (CO_P1(co)[i] + CO_P0(co)[i]); |
| 218 |
< |
src->srad = co->al; |
| 218 |
> |
src->srad = .5*co->al; |
| 219 |
|
src->ss2 = 2.*CO_R0(co)*co->al; |
| 220 |
|
/* set sampling vectors */ |
| 221 |
|
for (i = 0; i < 3; i++) |
| 228 |
|
fcross(src->ss[SW], src->ss[SV], co->ad); |
| 229 |
|
normalize(src->ss[SW]); |
| 230 |
|
for (i = 0; i < 3; i++) |
| 231 |
< |
src->ss[SW][i] *= .886227 * CO_R0(co); |
| 231 |
> |
src->ss[SW][i] *= .8559 * CO_R0(co); |
| 232 |
|
fcross(src->ss[SV], src->ss[SW], co->ad); |
| 233 |
|
} |
| 234 |
|
|
| 500 |
|
} |
| 501 |
|
|
| 502 |
|
|
| 503 |
< |
#define wrongsource(m, r) (r->rsrc>=0 && \ |
| 504 |
< |
source[r->rsrc].so!=r->ro && \ |
| 505 |
< |
(m->otype!=MAT_ILLUM || \ |
| 506 |
< |
objptr(source[r->rsrc].so->omod)->otype==MAT_ILLUM)) |
| 503 |
> |
/**************************************************************** |
| 504 |
> |
* The following macros were separated from the m_light() routine |
| 505 |
> |
* because they are very nasty and difficult to understand. |
| 506 |
> |
*/ |
| 507 |
|
|
| 508 |
+ |
/* wrongillum * |
| 509 |
+ |
* |
| 510 |
+ |
* We cannot allow an illum to pass to another illum, because that |
| 511 |
+ |
* would almost certainly constitute overcounting. |
| 512 |
+ |
* However, we do allow an illum to pass to another illum |
| 513 |
+ |
* that is actually going to relay to a virtual light source. |
| 514 |
+ |
*/ |
| 515 |
+ |
|
| 516 |
+ |
#define wrongillum(m, r) (!(source[r->rsrc].sflags&SVIRTUAL) && \ |
| 517 |
+ |
objptr(source[r->rsrc].so->omod)->otype==MAT_ILLUM) |
| 518 |
+ |
|
| 519 |
+ |
/* wrongsource * |
| 520 |
+ |
* |
| 521 |
+ |
* This source is the wrong source (ie. overcounted) if we are |
| 522 |
+ |
* aimed to a different source than the one we hit and the one |
| 523 |
+ |
* we hit is not an illum which should be passed. |
| 524 |
+ |
*/ |
| 525 |
+ |
|
| 526 |
+ |
#define wrongsource(m, r) (r->rsrc>=0 && source[r->rsrc].so!=r->ro && \ |
| 527 |
+ |
(m->otype!=MAT_ILLUM || wrongillum(m,r))) |
| 528 |
+ |
|
| 529 |
+ |
/* badspecular * |
| 530 |
+ |
* |
| 531 |
+ |
* Any undirected specular ray that hits a light source |
| 532 |
+ |
* should be discarded. This is because the source contribution to |
| 533 |
+ |
* specular components is calculated separately to reduce variance. |
| 534 |
+ |
*/ |
| 535 |
+ |
|
| 536 |
+ |
#define badspecular(m, r) (r->rsrc<0 && r->crtype&SPECULAR) |
| 537 |
+ |
|
| 538 |
+ |
/* distglow * |
| 539 |
+ |
* |
| 540 |
+ |
* A distant glow is an object that sometimes acts as a light source, |
| 541 |
+ |
* but is too far away from the test point to be one in this case. |
| 542 |
+ |
*/ |
| 543 |
+ |
|
| 544 |
|
#define distglow(m, r) (m->otype==MAT_GLOW && \ |
| 545 |
|
r->rot > m->oargs.farg[3]) |
| 546 |
|
|
| 547 |
+ |
/* badambient * |
| 548 |
+ |
* |
| 549 |
+ |
* We must avoid including counting light sources in the ambient calculation, |
| 550 |
+ |
* since the direct component is handled separately. Therefore, any |
| 551 |
+ |
* ambient ray which hits an active light source must be discarded. |
| 552 |
+ |
*/ |
| 553 |
+ |
|
| 554 |
|
#define badambient(m, r) ((r->crtype&(AMBIENT|SHADOW))==AMBIENT && \ |
| 555 |
|
!distglow(m, r)) |
| 556 |
|
|
| 557 |
+ |
/* overcount * |
| 558 |
+ |
* |
| 559 |
+ |
* All overcounting possibilities are contained here. |
| 560 |
+ |
*/ |
| 561 |
+ |
|
| 562 |
+ |
#define overcount(m, r) (badspecular(m,r) || wrongsource(m,r) || \ |
| 563 |
+ |
badambient(m,r)) |
| 564 |
+ |
|
| 565 |
+ |
/* passillum * |
| 566 |
+ |
* |
| 567 |
+ |
* An illum passes to another material type when we didn't hit it |
| 568 |
+ |
* on purpose (as part of a direct calculation), or it is relaying |
| 569 |
+ |
* a virtual light source. |
| 570 |
+ |
*/ |
| 571 |
+ |
|
| 572 |
|
#define passillum(m, r) (m->otype==MAT_ILLUM && \ |
| 573 |
< |
!(r->rsrc>=0&&source[r->rsrc].so==r->ro)) |
| 573 |
> |
(r->rsrc<0 || source[r->rsrc].so!=r->ro || \ |
| 574 |
> |
source[r->rsrc].sflags&SVIRTUAL)) |
| 575 |
|
|
| 576 |
+ |
/* srcignore * |
| 577 |
+ |
* |
| 578 |
+ |
* The -di flag renders light sources invisible, and here is the test. |
| 579 |
+ |
*/ |
| 580 |
+ |
|
| 581 |
|
#define srcignore(m, r) (directinvis && !(r->crtype&SHADOW) && \ |
| 582 |
|
!distglow(m, r)) |
| 583 |
|
|
| 587 |
|
register RAY *r; |
| 588 |
|
{ |
| 589 |
|
/* check for over-counting */ |
| 590 |
< |
if (wrongsource(m, r) || badambient(m, r)) |
| 590 |
> |
if (overcount(m, r)) |
| 591 |
|
return; |
| 592 |
|
/* check for passed illum */ |
| 593 |
|
if (passillum(m, r)) { |
