8 |
|
*/ |
9 |
|
|
10 |
|
#define _USE_MATH_DEFINES |
11 |
– |
#include <stdio.h> |
11 |
|
#include <stdlib.h> |
12 |
|
#include <math.h> |
13 |
|
#include <ctype.h> |
23 |
|
printf(" i\t\t\t\t Report general information (metadata)\n"); |
24 |
|
printf(" c\t\t\t\t Report diffuse and specular components\n"); |
25 |
|
printf(" q theta_i phi_i theta_o phi_o\t Query BSDF for given path (CIE-XYZ)\n"); |
26 |
< |
printf(" s[r|t][s|d] N theta phi\t\t\t Generate N ray directions & colors at given incidence\n"); |
27 |
< |
printf(" h[s|d] theta phi\t\t\t Report hemispherical scatteromg at given incidence\n"); |
28 |
< |
printf(" r[s|d] theta phi\t\t\t Report hemispherical reflection at given incidence\n"); |
29 |
< |
printf(" t[s|d] theta phi\t\t\t Report hemispherical transmission at given incidence\n"); |
26 |
> |
printf(" s[r|t][s|d] N theta phi\t Generate N ray directions & colors at given incidence\n"); |
27 |
> |
printf(" h[s|d] theta phi\t\t Report hemispherical scattering at given incidence\n"); |
28 |
> |
printf(" r[s|d] theta phi\t\t Report hemispherical reflection at given incidence\n"); |
29 |
> |
printf(" t[s|d] theta phi\t\t Report hemispherical transmission at given incidence\n"); |
30 |
|
printf(" a theta phi [t2 p2]\t\t Report resolution (in proj. steradians) for given direction(s)\n"); |
31 |
|
printf(" ^D\t\t\t\t Quit program\n"); |
32 |
|
} |
46 |
|
static void |
47 |
|
printXYZ(const char *intro, const SDValue *vp) |
48 |
|
{ |
49 |
+ |
if (vp->cieY <= 1e-9) { |
50 |
+ |
printf("%s0 0 0\n", intro); |
51 |
+ |
return; |
52 |
+ |
} |
53 |
|
printf("%s%.3e %.3e %.3e\n", intro, |
54 |
|
vp->spec.cx/vp->spec.cy*vp->cieY, |
55 |
|
vp->cieY, |
107 |
|
printf("Manufacturer: '%s'\n", bsdf->makr); |
108 |
|
printf("Width, Height, Thickness (m): %.4e, %.4e, %.4e\n", |
109 |
|
bsdf->dim[0], bsdf->dim[1], bsdf->dim[2]); |
110 |
< |
printf("Has geometry: %s\n", bsdf->mgf ? "yes" : "no"); |
111 |
< |
continue; |
110 |
> |
if (bsdf->mgf) |
111 |
> |
printf("Has geometry: %lu bytes\n", |
112 |
> |
(unsigned long)strlen(bsdf->mgf)); |
113 |
> |
else |
114 |
> |
printf("Has geometry: no\n"); |
115 |
> |
break; |
116 |
|
case 'C': /* report constant values */ |
117 |
|
if (!bsdf) |
118 |
|
goto noBSDFerr; |
133 |
|
printXYZ("Diffuse Front Reflectance: ", &bsdf->rLambFront); |
134 |
|
printXYZ("Diffuse Back Reflectance: ", &bsdf->rLambBack); |
135 |
|
printXYZ("Diffuse Transmittance: ", &bsdf->tLamb); |
136 |
< |
continue; |
136 |
> |
break; |
137 |
|
case 'Q': /* query BSDF value */ |
138 |
|
if (!bsdf) |
139 |
|
goto noBSDFerr; |
143 |
|
vec_from_deg(vout, atof(sskip2(cp,3)), atof(sskip2(cp,4))); |
144 |
|
if (!SDreportError(SDevalBSDF(&val, vout, vin, bsdf), stderr)) |
145 |
|
printXYZ("", &val); |
146 |
< |
continue; |
146 |
> |
break; |
147 |
|
case 'S': /* sample BSDF */ |
148 |
|
if (!bsdf) |
149 |
|
goto noBSDFerr; |
171 |
|
printf("%.8f %.8f %.8f ", vout[0], vout[1], vout[2]); |
172 |
|
printXYZ("", &val); |
173 |
|
} |
174 |
< |
continue; |
174 |
> |
break; |
175 |
|
case 'H': /* hemispherical values */ |
176 |
|
case 'R': |
177 |
|
case 'T': |
189 |
|
sflags &= ~SDsampSp; |
190 |
|
vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2))); |
191 |
|
printf("%.4e\n", SDdirectHemi(vin, sflags, bsdf)); |
192 |
< |
continue; |
192 |
> |
break; |
193 |
|
case 'A': /* resolution in proj. steradians */ |
194 |
|
if (!bsdf) |
195 |
|
goto noBSDFerr; |
198 |
|
vec_from_deg(vin, atof(sskip2(cp,1)), atof(sskip2(cp,2))); |
199 |
|
if (*sskip2(cp,4)) { |
200 |
|
vec_from_deg(vout, atof(sskip2(cp,3)), atof(sskip2(cp,4))); |
201 |
< |
if (SDreportError(SDsizeBSDF(proja, vin, vout, |
201 |
> |
if (SDreportError(SDsizeBSDF(proja, vout, vin, |
202 |
|
SDqueryMin+SDqueryMax, bsdf), stderr)) |
203 |
|
continue; |
204 |
|
} else if (SDreportError(SDsizeBSDF(proja, vin, NULL, |
205 |
|
SDqueryMin+SDqueryMax, bsdf), stderr)) |
206 |
|
continue; |
207 |
|
printf("%.4e %.4e\n", proja[0], proja[1]); |
208 |
< |
continue; |
208 |
> |
break; |
209 |
> |
default: |
210 |
> |
Usage(argv[0]); |
211 |
> |
break; |
212 |
|
} |
213 |
< |
Usage(argv[0]); |
213 |
> |
fflush(stdout); /* in case we're on remote */ |
214 |
|
continue; |
215 |
|
noBSDFerr: |
216 |
|
fprintf(stderr, "%s: First, use 'L' command to load BSDF\n", argv[0]); |