395 |
|
static int |
396 |
|
filter_klems_matrix(FILE *fp) |
397 |
|
{ |
398 |
– |
#define MAX_COLUMNS 145 |
398 |
|
const char *bn = klems_basis_name[angle_basis]; |
400 |
– |
float col_corr[MAX_COLUMNS]; |
399 |
|
int i, j, n = nabases; |
400 |
|
/* get angle basis */ |
401 |
|
while (n-- > 0) |
403 |
|
break; |
404 |
|
if (n < 0) |
405 |
|
return 0; |
408 |
– |
if (abase_list[n].nangles > MAX_COLUMNS) { |
409 |
– |
fputs("Internal error - too many Klems columns!\n", stderr); |
410 |
– |
return 0; |
411 |
– |
} |
412 |
– |
/* get correction factors */ |
413 |
– |
for (j = abase_list[n].nangles; j--; ) |
414 |
– |
col_corr[j] = 1.f / io_getohm(j, &abase_list[n]); |
406 |
|
/* read/correct/write matrix */ |
407 |
|
for (i = 0; i < abase_list[n].nangles; i++) { |
408 |
+ |
const double corr = 1./io_getohm(i, &abase_list[n]); |
409 |
|
for (j = 0; j < abase_list[n].nangles; j++) { |
410 |
|
double d; |
411 |
|
if (fscanf(fp, "%lf", &d) != 1) |
414 |
|
fputs("Negative BSDF data!\n", stderr); |
415 |
|
return 0; |
416 |
|
} |
417 |
< |
printf(" %.3e", d*col_corr[j]*(d > 0)); |
417 |
> |
printf(" %.3e", d*corr*(d > 0)); |
418 |
|
} |
419 |
|
fputc('\n', stdout); |
420 |
|
} |
424 |
|
return 0; |
425 |
|
} |
426 |
|
return 1; /* all is good */ |
435 |
– |
#undef MAX_COLUMNS |
427 |
|
} |
428 |
|
|
429 |
|
/* Write out BSDF data block with surrounding tags */ |