603 |
|
/* record input range */ |
604 |
|
scale = myScaffold.wmin / (double)iwmax; |
605 |
|
for (i = myScaffold.nic; i--; ) { |
606 |
< |
cd->clim[i][0] = floor(pos[i]/scale) * scale; |
606 |
> |
cd->clim[i][0] = floor(pos[i+2*rev]/scale) * scale; |
607 |
|
cd->clim[i][1] = cd->clim[i][0] + scale; |
608 |
|
} |
609 |
|
if (cd->isodist) { /* avoid issue in SDqueryTreProjSA() */ |
728 |
|
} else { |
729 |
|
const SDTreCDst *cd = (const SDTreCDst *)SDgetTreCDist(v1, sdc); |
730 |
|
if (cd == NULL) |
731 |
< |
return SDEmemory; |
731 |
> |
cd = &SDemptyCD; |
732 |
|
myPSA[0] = M_PI * (cd->clim[0][1] - cd->clim[0][0]) * |
733 |
|
(cd->clim[1][1] - cd->clim[1][0]); |
734 |
|
myPSA[1] = cd->max_psa; |
1163 |
|
/* separate diffuse components */ |
1164 |
|
extract_diffuse(&sd->rLambFront, sd->rf); |
1165 |
|
extract_diffuse(&sd->rLambBack, sd->rb); |
1166 |
< |
extract_diffuse(&sd->tLamb, (sd->tf != NULL) ? sd->tf : sd->tb); |
1166 |
> |
if (sd->tf != NULL) |
1167 |
> |
extract_diffuse(&sd->tLamb, sd->tf); |
1168 |
> |
if (sd->tb != NULL) |
1169 |
> |
extract_diffuse(&sd->tLamb, sd->tb); |
1170 |
|
/* return success */ |
1171 |
|
return SDEnone; |
1172 |
|
} |