| 295 |
|
{ |
| 296 |
|
CMATRIX *cm = cm_alloc(bsdf->nout, bsdf->ninc); |
| 297 |
|
COLORV *mp = cm->cmem; |
| 298 |
+ |
int nbadohm = 0; |
| 299 |
+ |
int nneg = 0; |
| 300 |
|
int r, c; |
| 301 |
|
|
| 302 |
|
for (r = 0; r < cm->nrows; r++) |
| 303 |
< |
for (c = 0; c < cm->ncols; c++, mp += 3) |
| 304 |
< |
mp[0] = mp[1] = mp[2] = BSDF_value(bsdf,c,r); |
| 303 |
> |
for (c = 0; c < cm->ncols; c++, mp += 3) { |
| 304 |
> |
float f = BSDF_value(bsdf,c,r); |
| 305 |
> |
float dom = getBSDF_incohm(bsdf,c); |
| 306 |
> |
FVECT v; |
| 307 |
> |
|
| 308 |
> |
if (f <= .0) { |
| 309 |
> |
nneg += (f < -FTINY); |
| 310 |
> |
continue; |
| 311 |
> |
} |
| 312 |
> |
if (dom <= .0) { |
| 313 |
> |
nbadohm++; |
| 314 |
> |
continue; |
| 315 |
> |
} |
| 316 |
> |
if (!getBSDF_incvec(v,bsdf,c) || v[2] > FTINY) |
| 317 |
> |
error(USER, "illegal incoming BTDF direction"); |
| 318 |
> |
|
| 319 |
> |
mp[0] = mp[1] = mp[2] = f * dom * -v[2]; |
| 320 |
> |
} |
| 321 |
> |
if (nbadohm | nneg) { |
| 322 |
> |
sprintf(errmsg, |
| 323 |
> |
"BTDF has %d negatives and %d bad incoming solid angles"); |
| 324 |
> |
error(WARNING, errmsg); |
| 325 |
> |
} |
| 326 |
|
return(cm); |
| 327 |
|
} |
| 328 |
|
|