16 |
|
|
17 |
|
char *progname; /* global argv[0] */ |
18 |
|
|
19 |
< |
/* Data types for matrix loading */ |
19 |
> |
/* Data types for file loading */ |
20 |
|
enum {DTfromHeader, DTascii, DTfloat, DTdouble, DTrgbe, DTxyze}; |
21 |
|
|
22 |
|
/* A color coefficient matrix -- vectors have ncols==1 */ |
125 |
|
break; |
126 |
|
default: |
127 |
|
error(USER, "unexpected data type in cm_load()"); |
128 |
– |
return(NULL); |
128 |
|
} |
129 |
|
if (nrows <= 0) { /* don't know length? */ |
130 |
|
int guessrows = 147; /* usually big enough */ |
134 |
|
long endpos = ftell(fp); |
135 |
|
long elemsiz = 3*(dtype==DTfloat ? |
136 |
|
sizeof(float) : sizeof(double)); |
137 |
< |
guessrows = (endpos - startpos)/elemsiz; |
137 |
> |
|
138 |
> |
if ((endpos - startpos) % (ncols*elemsiz)) { |
139 |
> |
sprintf(errmsg, |
140 |
> |
"improper length for binary file '%s'", |
141 |
> |
fname); |
142 |
> |
error(USER, errmsg); |
143 |
> |
} |
144 |
> |
guessrows = (endpos - startpos)/(ncols*elemsiz); |
145 |
|
if (fseek(fp, startpos, SEEK_SET) < 0) { |
146 |
|
sprintf(errmsg, |
147 |
|
"fseek() error on file '%s'", |
148 |
|
fname); |
149 |
|
error(SYSTEM, errmsg); |
150 |
|
} |
151 |
< |
if ((endpos - startpos) % elemsiz == 0) |
146 |
< |
nrows = guessrows; /* confident */ |
151 |
> |
nrows = guessrows; /* we're confident */ |
152 |
|
} |
153 |
|
} |
154 |
|
cm = cm_alloc(guessrows, ncols); |
160 |
|
int maxrow = (nrows > 0 ? nrows : 32000); |
161 |
|
int r, c; |
162 |
|
for (r = 0; r < maxrow; r++) { |
163 |
< |
if (r >= cm->nrows) /* need more space? */ |
163 |
> |
if (r >= cm->nrows) /* need more space? */ |
164 |
|
cm = cm_resize(cm, 2*cm->nrows); |
165 |
|
for (c = 0; c < ncols; c++) { |
166 |
|
COLORV *cv = cm_lval(cm,r,c); |
167 |
|
if (fscanf(fp, COLSPEC, cv, cv+1, cv+2) != 3) |
168 |
< |
if ((nrows <= 0) & (r > 0) & (c == 0)) { |
168 |
> |
if ((nrows <= 0) & (r > 0) & !c) { |
169 |
|
cm = cm_resize(cm, maxrow=r); |
170 |
|
break; |
171 |
|
} else |
193 |
|
if (nread == cm->nrows*cm->ncols) |
194 |
|
/* need more space? */ |
195 |
|
cm = cm_resize(cm, 2*cm->nrows); |
196 |
< |
else if (nread % cm->ncols == 0) |
196 |
> |
else if (nread && !(nread % cm->ncols)) |
197 |
|
/* seem to be done */ |
198 |
|
cm = cm_resize(cm, nread/cm->ncols); |
199 |
|
else /* ended mid-row */ |
215 |
|
copycolor(cvp, dc); |
216 |
|
cvp += 3; |
217 |
|
} |
213 |
– |
|
218 |
|
} else /* dtype == DTfloat */ { |
219 |
|
float fc[3]; /* load from float */ |
220 |
|
COLORV *cvp = cm->cmem; |
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); |
302 |
> |
for (c = 0; c < cm->ncols; c++) { |
303 |
> |
float dom = getBSDF_incohm(bsdf,c); |
304 |
> |
FVECT v; |
305 |
> |
|
306 |
> |
if (dom <= .0) { |
307 |
> |
nbadohm++; |
308 |
> |
continue; |
309 |
> |
} |
310 |
> |
if (!getBSDF_incvec(v,bsdf,c) || v[2] > FTINY) |
311 |
> |
error(USER, "illegal incoming BTDF direction"); |
312 |
> |
dom *= -v[2]; |
313 |
> |
|
314 |
> |
for (r = 0; r < cm->nrows; r++) { |
315 |
> |
float f = BSDF_value(bsdf,c,r); |
316 |
> |
|
317 |
> |
if (f <= .0) { |
318 |
> |
nneg += (f < -FTINY); |
319 |
> |
f = .0f; |
320 |
> |
} |
321 |
> |
mp = cm_lval(cm,r,c); |
322 |
> |
mp[0] = mp[1] = mp[2] = f * dom; |
323 |
> |
} |
324 |
> |
} |
325 |
> |
if (nneg || nbadohm) { |
326 |
> |
sprintf(errmsg, |
327 |
> |
"BTDF has %d negatives and %d bad incoming solid angles", |
328 |
> |
nneg, nbadohm); |
329 |
> |
error(WARNING, errmsg); |
330 |
> |
} |
331 |
|
return(cm); |
332 |
|
} |
333 |
|
|
370 |
|
error(SYSTEM, "out of memory in sum_images()"); |
371 |
|
pmat = cm_alloc(myYR, myXR); |
372 |
|
memset(pmat->cmem, 0, sizeof(COLOR)*myXR*myYR); |
373 |
+ |
/* finish header */ |
374 |
+ |
fputformat(myDT==DTrgbe ? COLRFMT : CIEFMT, stdout); |
375 |
+ |
fputc('\n', stdout); |
376 |
+ |
fprtresolu(myXR, myYR, stdout); |
377 |
+ |
fflush(stdout); |
378 |
|
} else if ((dt != myDT) | (xr != myXR) | (yr != myYR)) { |
379 |
|
sprintf(errmsg, "picture '%s' format/size mismatch", |
380 |
|
fname); |
397 |
|
fclose(fp); /* done this picture */ |
398 |
|
} |
399 |
|
free(scanline); |
363 |
– |
/* finish header */ |
364 |
– |
fputformat(myDT==DTrgbe ? COLRFMT : CIEFMT, stdout); |
365 |
– |
fputc('\n', stdout); |
366 |
– |
fprtresolu(myXR, myYR, stdout); |
400 |
|
/* write scanlines */ |
401 |
|
for (y = 0; y < myYR; y++) |
402 |
|
if (fwritescan((COLOR *)cm_lval(pmat, y, 0), myXR, stdout) < 0) |
423 |
|
int |
424 |
|
main(int argc, char *argv[]) |
425 |
|
{ |
426 |
< |
CMATRIX *tvec, *Dmat, *Tmat, *vtmp, *cvec; |
426 |
> |
CMATRIX *tvec, *Dmat, *Tmat, *ivec, *cvec; |
427 |
|
struct BSDF_data *btdf; |
428 |
|
|
429 |
|
progname = argv[0]; |
433 |
|
progname); |
434 |
|
return(1); |
435 |
|
} |
436 |
+ |
/* load Tregenza vector */ |
437 |
|
tvec = cm_load(argv[4], 0, 1, DTascii); /* argv[4]==NULL iff argc==4 */ |
438 |
< |
Dmat = cm_load(argv[3], 0, tvec->nrows, DTfromHeader); |
438 |
> |
/* load BTDF */ |
439 |
|
btdf = load_BSDF(argv[2]); |
440 |
|
if (btdf == NULL) |
441 |
|
return(1); |
442 |
< |
if (btdf->ninc != Dmat->nrows) { |
443 |
< |
sprintf(errmsg, "Incoming BTDF dir (%d) mismatch to D (%d)", |
410 |
< |
btdf->ninc, Dmat->nrows); |
411 |
< |
error(USER, errmsg); |
412 |
< |
} |
442 |
> |
/* load Daylight matrix */ |
443 |
> |
Dmat = cm_load(argv[3], btdf->ninc, tvec->nrows, DTfromHeader); |
444 |
|
/* multiply vector through */ |
445 |
< |
vtmp = cm_multiply(Dmat, tvec); |
445 |
> |
ivec = cm_multiply(Dmat, tvec); |
446 |
|
cm_free(Dmat); cm_free(tvec); |
447 |
|
Tmat = cm_bsdf(btdf); /* convert BTDF to matrix */ |
448 |
|
free_BSDF(btdf); |
449 |
< |
cvec = cm_multiply(Tmat, vtmp); /* cvec = component vector */ |
450 |
< |
cm_free(Tmat); cm_free(vtmp); |
449 |
> |
cvec = cm_multiply(Tmat, ivec); /* cvec = component vector */ |
450 |
> |
cm_free(Tmat); cm_free(ivec); |
451 |
|
if (hasDecimalSpec(argv[1])) { /* generating image */ |
452 |
|
SET_FILE_BINARY(stdout); |
453 |
|
newheader("RADIANCE", stdout); |