429 |
|
progname); |
430 |
|
return(1); |
431 |
|
} |
432 |
+ |
/* load Tregenza vector */ |
433 |
|
tvec = cm_load(argv[4], 0, 1, DTascii); /* argv[4]==NULL iff argc==4 */ |
434 |
< |
Dmat = cm_load(argv[3], 0, tvec->nrows, DTfromHeader); |
434 |
> |
/* load BTDF */ |
435 |
|
btdf = load_BSDF(argv[2]); |
436 |
|
if (btdf == NULL) |
437 |
|
return(1); |
438 |
< |
if (btdf->ninc != Dmat->nrows) { |
439 |
< |
sprintf(errmsg, "Incoming BTDF dir (%d) mismatch to D (%d)", |
439 |
< |
btdf->ninc, Dmat->nrows); |
440 |
< |
error(USER, errmsg); |
441 |
< |
} |
438 |
> |
/* load Daylight matrix */ |
439 |
> |
Dmat = cm_load(argv[3], btdf->ninc, tvec->nrows, DTfromHeader); |
440 |
|
/* multiply vector through */ |
441 |
|
ivec = cm_multiply(Dmat, tvec); |
442 |
|
cm_free(Dmat); cm_free(tvec); |