11 |
|
#include "platform.h" |
12 |
|
#include "rtio.h" |
13 |
|
#include "resolu.h" |
14 |
< |
#ifdef _WIN32 |
15 |
< |
#undef ftello |
16 |
< |
#define ftello ftell |
17 |
< |
#undef ssize_t |
18 |
< |
#define ssize_t size_t |
14 |
> |
#if defined(_WIN32) || defined(_WIN64) |
15 |
> |
#undef ftello |
16 |
> |
#define ftello ftell |
17 |
> |
#undef ssize_t |
18 |
> |
#define ssize_t size_t |
19 |
|
#else |
20 |
< |
#include <sys/mman.h> |
20 |
> |
#include <sys/mman.h> |
21 |
|
#endif |
22 |
|
|
23 |
– |
#ifdef getc_unlocked /* avoid horrendous overhead of flockfile */ |
24 |
– |
#undef getc |
25 |
– |
#undef putc |
26 |
– |
#define getc getc_unlocked |
27 |
– |
#define putc putc_unlocked |
28 |
– |
#endif |
29 |
– |
|
23 |
|
typedef struct { |
24 |
|
void *base; /* pointer to base memory */ |
25 |
|
size_t len; /* allocated memory length */ |
50 |
|
mp->len = 0; |
51 |
|
} |
52 |
|
|
53 |
+ |
/* load memory from an input stream, starting from current position */ |
54 |
+ |
static int |
55 |
+ |
load_stream(MEMLOAD *mp, FILE *fp) |
56 |
+ |
{ |
57 |
+ |
size_t alloced = 0; |
58 |
+ |
char buf[8192]; |
59 |
+ |
size_t nr; |
60 |
+ |
|
61 |
+ |
if (mp == NULL) |
62 |
+ |
return(-1); |
63 |
+ |
mp->base = NULL; |
64 |
+ |
mp->len = 0; |
65 |
+ |
mp->mapped = 0; |
66 |
+ |
if (fp == NULL) |
67 |
+ |
return(-1); |
68 |
+ |
while ((nr = fread(buf, 1, sizeof(buf), fp)) > 0) { |
69 |
+ |
if (!alloced) |
70 |
+ |
mp->base = malloc(alloced = nr); |
71 |
+ |
else if (mp->len+nr > alloced) |
72 |
+ |
mp->base = realloc(mp->base, |
73 |
+ |
alloced = alloced*(2+(nr==sizeof(buf)))/2+nr); |
74 |
+ |
if (mp->base == NULL) |
75 |
+ |
return(-1); |
76 |
+ |
memcpy((char *)mp->base + mp->len, buf, nr); |
77 |
+ |
mp->len += nr; |
78 |
+ |
} |
79 |
+ |
if (ferror(fp)) { |
80 |
+ |
free_load(mp); |
81 |
+ |
return(-1); |
82 |
+ |
} |
83 |
+ |
if (alloced > mp->len*5/4) /* don't waste too much space */ |
84 |
+ |
mp->base = realloc(mp->base, mp->len); |
85 |
+ |
return(mp->len > 0); |
86 |
+ |
} |
87 |
+ |
|
88 |
|
/* load a file into memory */ |
89 |
|
static int |
90 |
|
load_file(MEMLOAD *mp, FILE *fp) |
92 |
|
int fd; |
93 |
|
off_t skip, flen; |
94 |
|
|
95 |
+ |
#if defined(_WIN32) || defined(_WIN64) |
96 |
+ |
/* too difficult to fix this */ |
97 |
+ |
return load_stream(mp, fp); |
98 |
+ |
#endif |
99 |
|
if (mp == NULL) |
100 |
|
return(-1); |
101 |
|
mp->base = NULL; |
131 |
|
return(1); |
132 |
|
} |
133 |
|
|
102 |
– |
/* load memory from an input stream, starting from current position */ |
103 |
– |
static int |
104 |
– |
load_stream(MEMLOAD *mp, FILE *fp) |
105 |
– |
{ |
106 |
– |
size_t alloced = 0; |
107 |
– |
char buf[8192]; |
108 |
– |
size_t nr; |
109 |
– |
|
110 |
– |
if (mp == NULL) |
111 |
– |
return(-1); |
112 |
– |
mp->base = NULL; |
113 |
– |
mp->len = 0; |
114 |
– |
mp->mapped = 0; |
115 |
– |
if (fp == NULL) |
116 |
– |
return(-1); |
117 |
– |
while ((nr = fread(buf, 1, sizeof(buf), fp)) > 0) { |
118 |
– |
if (!alloced) |
119 |
– |
mp->base = malloc(nr); |
120 |
– |
else if (mp->len+nr > alloced) |
121 |
– |
mp->base = realloc(mp->base, |
122 |
– |
alloced = alloced*(2+(nr==sizeof(buf)))/2+nr); |
123 |
– |
if (mp->base == NULL) |
124 |
– |
return(-1); |
125 |
– |
memcpy((char *)mp->base + mp->len, buf, nr); |
126 |
– |
mp->len += nr; |
127 |
– |
} |
128 |
– |
if (ferror(fp)) { |
129 |
– |
free_load(mp); |
130 |
– |
return(-1); |
131 |
– |
} |
132 |
– |
if (alloced > mp->len*5/4) /* don't waste too much space */ |
133 |
– |
mp->base = realloc(mp->base, mp->len); |
134 |
– |
return(mp->len > 0); |
135 |
– |
} |
136 |
– |
|
134 |
|
/* free a record index */ |
135 |
|
#define free_records(rp) free(rp) |
136 |
|
|
234 |
|
|
235 |
|
if (fp == NULL) |
236 |
|
return(0); |
237 |
< |
fflush(stdout); /* assumes nothing in input buffer */ |
238 |
< |
while ((n = read(fileno(fp), buf, sizeof(buf))) > 0) |
237 |
> |
fflush(stdout); |
238 |
> |
while ((n = fread(buf, 1, sizeof(buf), fp)) > 0) |
239 |
|
if (write(fileno(stdout), buf, n) != n) |
240 |
|
return(0); |
241 |
< |
return(n >= 0); |
241 |
> |
return(!ferror(fp)); |
242 |
|
} |
243 |
|
|
244 |
|
/* get next word from stream, leaving stream on EOL or start of next word */ |
271 |
|
int ni_rows = 0; /* number of input rows */ |
272 |
|
int no_columns = 0; /* number of output columns */ |
273 |
|
int no_rows = 0; /* number of output rows */ |
274 |
+ |
int transpose = 0; /* transpose rows & cols? */ |
275 |
+ |
int i_header = 1; /* input header? */ |
276 |
+ |
int o_header = 1; /* output header? */ |
277 |
|
|
278 |
|
/* check settings and assign defaults */ |
279 |
|
static int |
288 |
|
comp_size = sizeof(double); |
289 |
|
else if (!strcmp(fmtid, "byte")) |
290 |
|
comp_size = 1; |
291 |
< |
else { |
291 |
> |
else if (strcmp(fmtid, "ascii")) { |
292 |
|
fprintf(stderr, "Unsupported format: %s\n", fmtid); |
293 |
|
return(0); |
294 |
|
} |
295 |
|
} |
296 |
+ |
if (transpose && (no_rows <= 0) & (no_columns <= 0)) { |
297 |
+ |
if (ni_rows > 0) no_columns = ni_rows; |
298 |
+ |
if (ni_columns > 0) no_rows = ni_columns; |
299 |
+ |
} else if ((no_rows <= 0) & (no_columns > 0) && |
300 |
+ |
!((ni_rows*ni_columns) % no_columns)) |
301 |
+ |
no_rows = ni_rows*ni_columns/no_columns; |
302 |
|
if (n_comp <= 0) |
303 |
|
n_comp = 3; |
304 |
|
return(1); |
378 |
|
int columns2go = no_columns; |
379 |
|
char word[256]; |
380 |
|
/* sanity checks */ |
381 |
< |
if (comp_size) |
382 |
< |
return(output_stream(fp)); /* binary data -- just copy */ |
381 |
> |
if (comp_size || (no_columns == ni_columns) & (no_rows == ni_rows)) |
382 |
> |
return(output_stream(fp)); /* no-op -- just copy */ |
383 |
|
if (no_columns <= 0) { |
384 |
|
fprintf(stderr, "Missing -oc specification\n"); |
385 |
|
return(0); |
466 |
|
if (!strncmp(s, "NCOMP=", 6)) { |
467 |
|
n = atoi(s+6); |
468 |
|
if ((n_comp > 0) & (n != n_comp)) { |
469 |
< |
fputs("Incorrect number of components", stderr); |
469 |
> |
fputs("Incorrect number of components\n", stderr); |
470 |
|
return(-1); |
471 |
|
} |
472 |
|
n_comp = n; |
473 |
|
return(0); |
474 |
|
} |
475 |
< |
fputs(s, stdout); /* copy header info. */ |
475 |
> |
if (o_header) |
476 |
> |
fputs(s, stdout); /* copy header info. */ |
477 |
|
return(0); |
478 |
|
} |
479 |
|
|
481 |
|
int |
482 |
|
main(int argc, char *argv[]) |
483 |
|
{ |
484 |
< |
int i_header = 1; /* input header? */ |
478 |
< |
int o_header = 1; /* output header? */ |
479 |
< |
int transpose = 0; /* transpose rows & cols? */ |
480 |
< |
int i; |
484 |
> |
int a; |
485 |
|
|
486 |
< |
for (i = 1; i < argc && argv[i][0] == '-'; i++) |
487 |
< |
switch (argv[i][1]) { |
486 |
> |
for (a = 1; a < argc && argv[a][0] == '-'; a++) |
487 |
> |
switch (argv[a][1]) { |
488 |
|
case 'i': /* input */ |
489 |
< |
if (argv[i][2] == 'c') /* columns */ |
490 |
< |
ni_columns = atoi(argv[++i]); |
491 |
< |
else if (argv[i][2] == 'r') |
492 |
< |
ni_rows = atoi(argv[++i]); |
489 |
> |
if (argv[a][2] == 'c') /* columns */ |
490 |
> |
ni_columns = atoi(argv[++a]); |
491 |
> |
else if (argv[a][2] == 'r') |
492 |
> |
ni_rows = atoi(argv[++a]); |
493 |
|
else |
494 |
|
goto userr; |
495 |
|
break; |
496 |
|
case 'o': /* output */ |
497 |
< |
if (argv[i][2] == 'c') /* columns */ |
498 |
< |
no_columns = atoi(argv[++i]); |
499 |
< |
else if (argv[i][2] == 'r') |
500 |
< |
no_rows = atoi(argv[++i]); |
497 |
> |
if (argv[a][2] == 'c') /* columns */ |
498 |
> |
no_columns = atoi(argv[++a]); |
499 |
> |
else if (argv[a][2] == 'r') |
500 |
> |
no_rows = atoi(argv[++a]); |
501 |
|
else |
502 |
|
goto userr; |
503 |
|
break; |
504 |
|
case 'h': /* turn off header */ |
505 |
< |
switch (argv[i][2]) { |
505 |
> |
switch (argv[a][2]) { |
506 |
|
case 'i': |
507 |
|
i_header = 0; |
508 |
|
break; |
520 |
|
transpose = !transpose; |
521 |
|
break; |
522 |
|
case 'f': /* format */ |
523 |
< |
switch (argv[i][2]) { |
523 |
> |
switch (argv[a][2]) { |
524 |
|
case 'a': /* ASCII */ |
525 |
|
case 'A': |
526 |
|
fmtid = "ascii"; |
544 |
|
default: |
545 |
|
goto userr; |
546 |
|
} |
547 |
< |
if (argv[i][3]) { |
548 |
< |
if (!isdigit(argv[i][3])) |
547 |
> |
if (argv[a][3]) { |
548 |
> |
if (!isdigit(argv[a][3])) |
549 |
|
goto userr; |
550 |
< |
n_comp = atoi(argv[i]+3); |
551 |
< |
} |
550 |
> |
n_comp = atoi(argv[a]+3); |
551 |
> |
} else |
552 |
> |
n_comp = 1; |
553 |
|
break; |
554 |
|
case 'w': /* warnings on/off */ |
555 |
|
warnings = !warnings; |
557 |
|
default: |
558 |
|
goto userr; |
559 |
|
} |
560 |
< |
if (i < argc-1) /* arg count OK? */ |
560 |
> |
if (a < argc-1) /* arg count OK? */ |
561 |
|
goto userr; |
562 |
|
/* open input file? */ |
563 |
< |
if (i == argc-1 && freopen(argv[i], "r", stdin) == NULL) { |
564 |
< |
fprintf(stderr, "%s: cannot open for reading\n", argv[i]); |
563 |
> |
if (a == argc-1 && freopen(argv[a], "r", stdin) == NULL) { |
564 |
> |
fprintf(stderr, "%s: cannot open for reading\n", argv[a]); |
565 |
|
return(1); |
566 |
|
} |
567 |
|
if (comp_size) { |
569 |
|
SET_FILE_BINARY(stdout); |
570 |
|
} |
571 |
|
/* check for no-op */ |
572 |
< |
if (!transpose & (i_header == o_header) && (comp_size || |
573 |
< |
(no_columns == ni_columns) & (no_rows == ni_rows))) { |
572 |
> |
if (!transpose & (i_header == o_header) && |
573 |
> |
(no_columns == ni_columns) & (no_rows == ni_rows)) { |
574 |
|
if (warnings) |
575 |
|
fprintf(stderr, "%s: no-op -- copying input verbatim\n", |
576 |
|
argv[0]); |
579 |
|
return(0); |
580 |
|
} |
581 |
|
if (i_header) { /* read header */ |
582 |
< |
if (getheader(stdin, &headline, NULL) < 0) |
582 |
> |
if (getheader(stdin, headline, NULL) < 0) |
583 |
|
return(1); |
584 |
|
if (!check_sizes()) |
585 |
|
return(1); |
590 |
|
} else if (!check_sizes()) |
591 |
|
return(1); |
592 |
|
if (o_header) { /* write header */ |
593 |
< |
printargs(argc, argv, stdout); |
589 |
< |
if (transpose && (no_rows <= 0) & (no_columns <= 0)) { |
590 |
< |
if (ni_rows > 0) no_columns = ni_rows; |
591 |
< |
if (ni_columns > 0) no_rows = ni_columns; |
592 |
< |
} |
593 |
> |
printargs(a, argv, stdout); |
594 |
|
if (no_rows > 0) |
595 |
|
printf("NROWS=%d\n", no_rows); |
596 |
|
if (no_columns > 0) |
600 |
|
fputc('\n', stdout); /* finish new header */ |
601 |
|
} |
602 |
|
if (transpose) { /* transposing rows & columns? */ |
603 |
< |
MEMLOAD myMem; /* need to load into memory */ |
604 |
< |
if (i == argc-1) { |
603 |
> |
MEMLOAD myMem; /* need to map into memory */ |
604 |
> |
if (a == argc-1) { |
605 |
|
if (load_file(&myMem, stdin) <= 0) { |
606 |
|
fprintf(stderr, "%s: error loading file into memory\n", |
607 |
< |
argv[i]); |
607 |
> |
argv[a]); |
608 |
|
return(1); |
609 |
|
} |
610 |
|
} else if (load_stream(&myMem, stdin) <= 0) { |
614 |
|
} |
615 |
|
if (!do_transpose(&myMem)) |
616 |
|
return(1); |
617 |
< |
/* free_load(&myMem); */ |
618 |
< |
} else if (!do_resize(stdin)) /* just reshaping input */ |
617 |
> |
/* free_load(&myMem); about to exit, so don't bother */ |
618 |
> |
} else if (!do_resize(stdin)) /* reshaping input */ |
619 |
|
return(1); |
620 |
|
return(0); |
621 |
|
userr: |