| 1 |
+ |
#ifndef lint |
| 2 |
+ |
static const char RCSid[] = "$Id$"; |
| 3 |
+ |
#endif |
| 4 |
|
/* |
| 5 |
|
================================================================== |
| 6 |
|
Photon map support routines for emission from light sources |
| 10 |
|
(c) Lucerne University of Applied Sciences and Arts, |
| 11 |
|
supported by the Swiss National Science Foundation (SNSF, #147053) |
| 12 |
|
================================================================== |
| 13 |
< |
|
| 11 |
< |
$Id$ |
| 13 |
> |
|
| 14 |
|
*/ |
| 15 |
|
|
| 16 |
|
|
| 525 |
|
/* Find geometry declared as photon ports */ |
| 526 |
|
{ |
| 527 |
|
OBJECT i; |
| 528 |
< |
OBJREC* obj; |
| 528 |
> |
OBJREC *obj, *mat; |
| 529 |
|
|
| 530 |
|
/* Check for missing port modifiers */ |
| 531 |
|
if (!ambset [0]) |
| 532 |
< |
error(USER, "no photon ports found"); |
| 532 |
> |
error(USER, "no photon ports"); |
| 533 |
|
|
| 534 |
|
for (i = 0; i < nobjects; i++) { |
| 535 |
|
obj = objptr(i); |
| 536 |
+ |
mat = findmaterial(obj); |
| 537 |
|
|
| 538 |
< |
if (inset(ambset, obj -> omod)) { |
| 538 |
> |
/* Check if object is a surface and NOT a light source (duh) and |
| 539 |
> |
* resolve its material via any aliases, then check for inclusion in |
| 540 |
> |
* port modifier list */ |
| 541 |
> |
if (issurface(obj -> otype) && !islight(mat -> otype) && |
| 542 |
> |
inset(ambset, objndx(mat))) { |
| 543 |
|
/* Add photon port */ |
| 544 |
|
photonPorts = (SRCREC*)realloc(photonPorts, |
| 545 |
|
(numPhotonPorts + 1) * |
| 557 |
|
numPhotonPorts++; |
| 558 |
|
} |
| 559 |
|
} |
| 560 |
+ |
|
| 561 |
+ |
if (!numPhotonPorts) |
| 562 |
+ |
error(USER, "no valid photon ports found"); |
| 563 |
|
} |
| 564 |
|
|
| 565 |
|
|
| 602 |
|
unsigned i, t, p; |
| 603 |
|
double phi, cosTheta, sinTheta, du, dv, dOmega, thetaScale; |
| 604 |
|
EmissionSample* sample; |
| 605 |
< |
const OBJREC* mod = objptr(emap -> src -> so -> omod); |
| 605 |
> |
const OBJREC* mod = findmaterial(emap -> src -> so); |
| 606 |
|
static RAY r; |
| 607 |
|
#if 0 |
| 608 |
|
static double lastCosNorm = FHUGE; |
| 609 |
|
static SRCREC *lastSrc = NULL, *lastPort = NULL; |
| 610 |
|
#endif |
| 611 |
|
|
| 612 |
+ |
setcolor(emap -> partFlux, 0, 0, 0); |
| 613 |
+ |
|
| 614 |
|
photonOrigin [emap -> src -> so -> otype] (emap); |
| 615 |
|
cosTheta = DOT(emap -> ws, emap -> wh); |
| 616 |
|
|
| 633 |
|
#endif |
| 634 |
|
|
| 635 |
|
/* Need to recompute flux & PDF */ |
| 624 |
– |
setcolor(emap -> partFlux, 0, 0, 0); |
| 636 |
|
emap -> cdf = 0; |
| 637 |
|
emap -> numSamples = 0; |
| 638 |
|
|
| 681 |
|
|
| 682 |
|
VCOPY(r.rorg, emap -> photonOrg); |
| 683 |
|
VCOPY(r.rop, emap -> photonOrg); |
| 684 |
< |
r.rmax = FHUGE; |
| 684 |
> |
r.rmax = 0; |
| 685 |
|
|
| 686 |
|
for (t = 0; t < emap -> numTheta; t++) { |
| 687 |
|
for (p = 0; p < emap -> numPhi; p++) { |
| 751 |
|
unsigned long i, lo, hi; |
| 752 |
|
const EmissionSample* sample = emap -> samples; |
| 753 |
|
RREAL du, dv, cosTheta, cosThetaSqr, sinTheta, phi; |
| 754 |
< |
const OBJREC* mod = objptr(emap -> src -> so -> omod); |
| 754 |
> |
const OBJREC* mod = findmaterial(emap -> src -> so); |
| 755 |
|
|
| 756 |
|
/* Choose a new origin within current partition for every |
| 757 |
|
emitted photon to break up clustering artifacts */ |
| 758 |
|
photonOrigin [emap -> src -> so -> otype] ((EmissionMap*)emap); |
| 759 |
+ |
/* If we have a local glow source with a maximum radius, then |
| 760 |
+ |
restrict our photon to the specified distance (otherwise no limit) */ |
| 761 |
+ |
if (mod -> otype == MAT_GLOW && !(emap -> src -> sflags & SDISTANT) |
| 762 |
+ |
&& mod -> oargs.farg[3] > FTINY) |
| 763 |
+ |
ray -> rmax = mod -> oargs.farg[3]; |
| 764 |
+ |
else |
| 765 |
+ |
ray -> rmax = 0; |
| 766 |
|
rayorigin(ray, PRIMARY, NULL, NULL); |
| 749 |
– |
ray -> rmax = FHUGE; |
| 767 |
|
|
| 768 |
|
if (!emap -> numSamples) { |
| 769 |
|
/* Source is unmodified and has no port, and either local with |
