| 236 |
|
diffY = 0; |
| 237 |
|
setcolor(cdiff, 0, 0, 0); |
| 238 |
|
} |
| 239 |
< |
/* need projected solid angles */ |
| 239 |
> |
/* need projected solid angle */ |
| 240 |
|
omega *= fabs(vsrc[2]); |
| 241 |
– |
ec = SDsizeBSDF(&tomega, ndp->vray, vsrc, SDqueryMin, ndp->sd); |
| 242 |
– |
if (ec) |
| 243 |
– |
goto baderror; |
| 241 |
|
/* check indirect over-counting */ |
| 242 |
|
if ((vsrc[2] > 0) ^ (ndp->vray[2] > 0) && bright(ndp->cthru) > FTINY) { |
| 243 |
< |
double dx = vsrc[0] + ndp->vray[0]; |
| 244 |
< |
double dy = vsrc[1] + ndp->vray[1]; |
| 245 |
< |
if (dx*dx + dy*dy <= (2.5*4./PI)*(omega + tomega + |
| 246 |
< |
2.*sqrt(omega*tomega))) |
| 243 |
> |
double dx = vsrc[0] + ndp->vray[0]; |
| 244 |
> |
double dy = vsrc[1] + ndp->vray[1]; |
| 245 |
> |
SDSpectralDF *dfp = (ndp->pr->rod > 0) ? |
| 246 |
> |
((ndp->sd->tf != NULL) ? ndp->sd->tf : ndp->sd->tb) : |
| 247 |
> |
((ndp->sd->tb != NULL) ? ndp->sd->tb : ndp->sd->tf) ; |
| 248 |
> |
|
| 249 |
> |
if (dx*dx + dy*dy <= (2.5*4./PI)*(omega + dfp->minProjSA + |
| 250 |
> |
2.*sqrt(omega*dfp->minProjSA))) |
| 251 |
|
return(0); |
| 252 |
|
} |
| 253 |
+ |
ec = SDsizeBSDF(&tomega, ndp->vray, vsrc, SDqueryMin, ndp->sd); |
| 254 |
+ |
if (ec) |
| 255 |
+ |
goto baderror; |
| 256 |
|
/* assign number of samples */ |
| 257 |
|
sf = specjitter * ndp->pr->rweight; |
| 258 |
|
if (tomega <= 0) |
