362 |
|
nearestNeighbours(pmap, pos, norm, 1); |
363 |
|
} |
364 |
|
|
365 |
+ |
if (pmap -> maxDist < FTINY) { |
366 |
+ |
sprintf(errmsg, "itsy bitsy teeny weeny photon search radius %e", |
367 |
+ |
sqrt(pmap -> maxDist)); |
368 |
+ |
error(WARNING, errmsg); |
369 |
+ |
} |
370 |
+ |
|
371 |
|
if (pmap -> squeueEnd < pmap -> squeueSize * pmap -> gatherTolerance) { |
372 |
|
/* Short lookup; too few photons found */ |
373 |
|
if (pmap -> squeueEnd > PMAP_SHORT_LOOKUP_THRESH) { |
398 |
|
sprintf(errmsg, |
399 |
|
redo ? "restarting photon lookup with max radius %.1e" |
400 |
|
: "max photon lookup radius adjusted to %.1e", |
401 |
< |
pmap -> maxDist0); |
401 |
> |
sqrt(pmap -> maxDist0)); |
402 |
|
error(WARNING, errmsg); |
403 |
|
#endif |
404 |
|
} |
405 |
|
#ifdef PMAP_LOOKUP_REDO |
406 |
|
else { |
407 |
|
sprintf(errmsg, "max photon lookup radius clamped to %.1e", |
408 |
< |
pmap -> maxDist0); |
408 |
> |
sqrt(pmap -> maxDist0)); |
409 |
|
error(WARNING, errmsg); |
410 |
|
} |
411 |
|
#endif |