| 24 |
|
printf(" i\t\t\t\t Report general information (metadata)\n"); |
| 25 |
|
printf(" c\t\t\t\t Report diffuse and specular components\n"); |
| 26 |
|
printf(" q theta_i phi_i theta_o phi_o\t Query BSDF for given path (CIE-XYZ)\n"); |
| 27 |
< |
printf(" s N theta phi\t\t\t Generate N ray directions at given incidence\n"); |
| 28 |
< |
printf(" h theta phi\t\t\t Report hemispherical total at given incidence\n"); |
| 29 |
< |
printf(" r theta phi\t\t\t Report hemispherical reflection at given incidence\n"); |
| 30 |
< |
printf(" t theta phi\t\t\t Report hemispherical transmission at given incidence\n"); |
| 27 |
> |
printf(" s[r|t][s|d] N theta phi\t Generate N ray directions & colors at given incidence\n"); |
| 28 |
> |
printf(" h[s|d] theta phi\t\t Report hemispherical scattering at given incidence\n"); |
| 29 |
> |
printf(" r[s|d] theta phi\t\t Report hemispherical reflection at given incidence\n"); |
| 30 |
> |
printf(" t[s|d] theta phi\t\t Report hemispherical transmission at given incidence\n"); |
| 31 |
|
printf(" a theta phi [t2 p2]\t\t Report resolution (in proj. steradians) for given direction(s)\n"); |
| 32 |
|
printf(" ^D\t\t\t\t Quit program\n"); |
| 33 |
|
} |
| 77 |
|
char *cp2; |
| 78 |
|
FVECT vin, vout; |
| 79 |
|
double proja[2]; |
| 80 |
< |
int n; |
| 80 |
> |
int n, i; |
| 81 |
|
SDValue val; |
| 82 |
|
|
| 83 |
|
while (isspace(*cp)) cp++; |
| 142 |
|
goto noBSDFerr; |
| 143 |
|
if (!*sskip2(cp,3)) |
| 144 |
|
break; |
| 145 |
< |
n = atoi(sskip2(cp,1)); |
| 145 |
> |
if (toupper(cp[1]) == 'R') { |
| 146 |
> |
sflags &= ~SDsampT; |
| 147 |
> |
++cp; |
| 148 |
> |
} else if (toupper(cp[1]) == 'T') { |
| 149 |
> |
sflags &= ~SDsampR; |
| 150 |
> |
++cp; |
| 151 |
> |
} |
| 152 |
> |
if (toupper(cp[1]) == 'S') |
| 153 |
> |
sflags &= ~SDsampDf; |
| 154 |
> |
else if (toupper(cp[1]) == 'D') |
| 155 |
> |
sflags &= ~SDsampSp; |
| 156 |
> |
i = n = atoi(sskip2(cp,1)); |
| 157 |
|
vec_from_deg(vin, atof(sskip2(cp,2)), atof(sskip2(cp,3))); |
| 158 |
< |
while (n-- > 0) { |
| 159 |
< |
if (SDreportError(SDsampBSDF(&val, vin, |
| 160 |
< |
rand()*(1./(RAND_MAX+.5)), |
| 158 |
> |
while (i-- > 0) { |
| 159 |
> |
VCOPY(vout, vin); |
| 160 |
> |
if (SDreportError(SDsampBSDF(&val, vout, |
| 161 |
> |
(i+rand()*(1./(RAND_MAX+.5)))/(double)n, |
| 162 |
|
sflags, bsdf), stderr)) |
| 163 |
|
break; |
| 164 |
< |
printf("%.8f %.8f %.8f\n", vin[0], vin[1], vin[2]); |
| 164 |
> |
printf("%.8f %.8f %.8f ", vout[0], vout[1], vout[2]); |
| 165 |
> |
printXYZ("", &val); |
| 166 |
|
} |
| 167 |
|
continue; |
| 168 |
< |
case 'H': /* hemispherical totals */ |
| 168 |
> |
case 'H': /* hemispherical values */ |
| 169 |
|
case 'R': |
| 170 |
|
case 'T': |
| 171 |
|
if (!bsdf) |
| 172 |
|
goto noBSDFerr; |
| 173 |
|
if (!*sskip2(cp,2)) |
| 174 |
|
break; |
| 175 |
< |
if (tolower(*cp) == 'r') |
| 175 |
> |
if (toupper(cp[0]) == 'R') |
| 176 |
|
sflags &= ~SDsampT; |
| 177 |
< |
else if (tolower(*cp) == 't') |
| 177 |
> |
else if (toupper(cp[0]) == 'T') |
| 178 |
|
sflags &= ~SDsampR; |
| 179 |
+ |
if (toupper(cp[1]) == 'S') |
| 180 |
+ |
sflags &= ~SDsampDf; |
| 181 |
+ |
else if (toupper(cp[1]) == 'D') |
| 182 |
+ |
sflags &= ~SDsampSp; |
| 183 |
|
vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2))); |
| 184 |
|
printf("%.4e\n", SDdirectHemi(vin, sflags, bsdf)); |
| 185 |
|
continue; |