| 8 |
|
*/ |
| 9 |
|
|
| 10 |
|
#define _USE_MATH_DEFINES |
| 11 |
– |
#include <stdio.h> |
| 11 |
|
#include <stdlib.h> |
| 13 |
– |
#include <string.h> |
| 12 |
|
#include <math.h> |
| 13 |
+ |
#include <ctype.h> |
| 14 |
|
#include "random.h" |
| 15 |
|
#include "platform.h" |
| 16 |
|
#include "paths.h" |
| 525 |
|
static int |
| 526 |
|
wrap_up(void) |
| 527 |
|
{ |
| 528 |
< |
char cmd[8192]; |
| 528 |
> |
char cmd[32700]; |
| 529 |
|
|
| 530 |
|
if (bsdf_manuf[0]) { |
| 531 |
|
add_wbsdf("-f", 1); |
| 576 |
|
} |
| 577 |
|
#endif |
| 578 |
|
|
| 579 |
+ |
#define HEAD_BUFLEN 10240 |
| 580 |
+ |
static char head_buf[HEAD_BUFLEN]; |
| 581 |
+ |
static int cur_headlen = 0; |
| 582 |
+ |
|
| 583 |
+ |
/* Record header line as comment associated with this SIR input */ |
| 584 |
+ |
static int |
| 585 |
+ |
record2header(char *s) |
| 586 |
+ |
{ |
| 587 |
+ |
int len = strlen(s); |
| 588 |
+ |
|
| 589 |
+ |
if (cur_headlen+len >= HEAD_BUFLEN-6) |
| 590 |
+ |
return(0); |
| 591 |
+ |
/* includes EOL */ |
| 592 |
+ |
strcpy(head_buf+cur_headlen, s); |
| 593 |
+ |
cur_headlen += len; |
| 594 |
+ |
|
| 595 |
+ |
#if defined(_WIN32) || defined(_WIN64) |
| 596 |
+ |
if (head_buf[cur_headlen-1] == '\n') |
| 597 |
+ |
head_buf[cur_headlen-1] = '\t'; |
| 598 |
+ |
#endif |
| 599 |
+ |
return(1); |
| 600 |
+ |
} |
| 601 |
+ |
|
| 602 |
+ |
/* Finish off header for this file */ |
| 603 |
+ |
static void |
| 604 |
+ |
done_header(void) |
| 605 |
+ |
{ |
| 606 |
+ |
while (cur_headlen > 0 && isspace(head_buf[cur_headlen-1])) |
| 607 |
+ |
--cur_headlen; |
| 608 |
+ |
head_buf[cur_headlen] = '\0'; |
| 609 |
+ |
if (!cur_headlen) |
| 610 |
+ |
return; |
| 611 |
+ |
add_wbsdf("-C", 1); |
| 612 |
+ |
add_wbsdf(head_buf, 0); |
| 613 |
+ |
head_buf[cur_headlen=0] = '\0'; |
| 614 |
+ |
} |
| 615 |
+ |
|
| 616 |
|
/* Read in BSDF and interpolate as Klems matrix representation */ |
| 617 |
|
int |
| 618 |
|
main(int argc, char *argv[]) |
| 619 |
|
{ |
| 620 |
|
int dofwd = 0, dobwd = 1; |
| 621 |
< |
char buf[2048]; |
| 621 |
> |
char buf[1024]; |
| 622 |
|
char *cp; |
| 623 |
|
int i, na; |
| 624 |
|
|
| 730 |
|
if (i < argc) { /* open input files if given */ |
| 731 |
|
int nbsdf = 0; |
| 732 |
|
for ( ; i < argc; i++) { /* interpolate each component */ |
| 697 |
– |
char pbuf[256]; |
| 733 |
|
FILE *fpin = fopen(argv[i], "rb"); |
| 734 |
|
if (fpin == NULL) { |
| 735 |
|
fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n", |
| 736 |
|
progname, argv[i]); |
| 737 |
|
return(1); |
| 738 |
|
} |
| 739 |
+ |
sprintf(buf, "%s:\n", argv[i]); |
| 740 |
+ |
record2header(buf); |
| 741 |
+ |
sir_headshare = &record2header; |
| 742 |
|
if (!load_bsdf_rep(fpin)) |
| 743 |
|
return(1); |
| 744 |
|
fclose(fpin); |
| 745 |
< |
sprintf(pbuf, "Interpolating component '%s'", argv[i]); |
| 746 |
< |
prog_start(pbuf); |
| 745 |
> |
done_header(); |
| 746 |
> |
sprintf(buf, "Interpolating component '%s'", argv[i]); |
| 747 |
> |
prog_start(buf); |
| 748 |
|
eval_rbf(); |
| 749 |
|
} |
| 750 |
|
return(wrap_up()); |
| 751 |
|
} |
| 752 |
|
SET_FILE_BINARY(stdin); /* load from stdin */ |
| 753 |
+ |
record2header("<stdin>:\n"); |
| 754 |
+ |
sir_headshare = &record2header; |
| 755 |
|
if (!load_bsdf_rep(stdin)) |
| 756 |
|
return(1); |
| 757 |
+ |
done_header(); |
| 758 |
|
prog_start("Interpolating from standard input"); |
| 759 |
|
eval_rbf(); /* resample dist. */ |
| 760 |
|
return(wrap_up()); |
