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 |
+ |
#ifdef _WIN32 /* too difficult to fix this */ |
96 |
+ |
return load_stream(mp, fp); |
97 |
+ |
#endif |
98 |
|
if (mp == NULL) |
99 |
|
return(-1); |
100 |
|
mp->base = NULL; |
130 |
|
return(1); |
131 |
|
} |
132 |
|
|
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 |
– |
|
133 |
|
/* free a record index */ |
134 |
|
#define free_records(rp) free(rp) |
135 |
|
|
270 |
|
int ni_rows = 0; /* number of input rows */ |
271 |
|
int no_columns = 0; /* number of output columns */ |
272 |
|
int no_rows = 0; /* number of output rows */ |
273 |
+ |
int transpose = 0; /* transpose rows & cols? */ |
274 |
+ |
int i_header = 1; /* input header? */ |
275 |
+ |
int o_header = 1; /* output header? */ |
276 |
|
|
277 |
|
/* check settings and assign defaults */ |
278 |
|
static int |
292 |
|
return(0); |
293 |
|
} |
294 |
|
} |
295 |
+ |
if (transpose && (no_rows <= 0) & (no_columns <= 0)) { |
296 |
+ |
if (ni_rows > 0) no_columns = ni_rows; |
297 |
+ |
if (ni_columns > 0) no_rows = ni_columns; |
298 |
+ |
} else if ((no_rows <= 0) & (no_columns > 0) && |
299 |
+ |
!((ni_rows*ni_columns) % no_columns)) |
300 |
+ |
no_rows = ni_rows*ni_columns/no_columns; |
301 |
|
if (n_comp <= 0) |
302 |
|
n_comp = 3; |
303 |
|
return(1); |
377 |
|
int columns2go = no_columns; |
378 |
|
char word[256]; |
379 |
|
/* sanity checks */ |
380 |
< |
if (comp_size) |
381 |
< |
return(output_stream(fp)); /* binary data -- just copy */ |
380 |
> |
if (comp_size || (no_columns == ni_columns) & (no_rows == ni_rows)) |
381 |
> |
return(output_stream(fp)); /* no-op -- just copy */ |
382 |
|
if (no_columns <= 0) { |
383 |
|
fprintf(stderr, "Missing -oc specification\n"); |
384 |
|
return(0); |
471 |
|
n_comp = n; |
472 |
|
return(0); |
473 |
|
} |
474 |
< |
fputs(s, stdout); /* copy header info. */ |
474 |
> |
if (o_header) |
475 |
> |
fputs(s, stdout); /* copy header info. */ |
476 |
|
return(0); |
477 |
|
} |
478 |
|
|
480 |
|
int |
481 |
|
main(int argc, char *argv[]) |
482 |
|
{ |
477 |
– |
int i_header = 1; /* input header? */ |
478 |
– |
int o_header = 1; /* output header? */ |
479 |
– |
int transpose = 0; /* transpose rows & cols? */ |
483 |
|
int a; |
484 |
|
|
485 |
|
for (a = 1; a < argc && argv[a][0] == '-'; a++) |
568 |
|
SET_FILE_BINARY(stdout); |
569 |
|
} |
570 |
|
/* check for no-op */ |
571 |
< |
if (!transpose & (i_header == o_header) && (comp_size || |
572 |
< |
(no_columns == ni_columns) & (no_rows == ni_rows))) { |
571 |
> |
if (!transpose & (i_header == o_header) && |
572 |
> |
(no_columns == ni_columns) & (no_rows == ni_rows)) { |
573 |
|
if (warnings) |
574 |
|
fprintf(stderr, "%s: no-op -- copying input verbatim\n", |
575 |
|
argv[0]); |
578 |
|
return(0); |
579 |
|
} |
580 |
|
if (i_header) { /* read header */ |
581 |
< |
if (getheader(stdin, &headline, NULL) < 0) |
581 |
> |
if (getheader(stdin, headline, NULL) < 0) |
582 |
|
return(1); |
583 |
|
if (!check_sizes()) |
584 |
|
return(1); |
590 |
|
return(1); |
591 |
|
if (o_header) { /* write header */ |
592 |
|
printargs(a, argv, stdout); |
590 |
– |
if (transpose && (no_rows <= 0) & (no_columns <= 0)) { |
591 |
– |
if (ni_rows > 0) no_columns = ni_rows; |
592 |
– |
if (ni_columns > 0) no_rows = ni_columns; |
593 |
– |
} |
593 |
|
if (no_rows > 0) |
594 |
|
printf("NROWS=%d\n", no_rows); |
595 |
|
if (no_columns > 0) |
599 |
|
fputc('\n', stdout); /* finish new header */ |
600 |
|
} |
601 |
|
if (transpose) { /* transposing rows & columns? */ |
602 |
< |
MEMLOAD myMem; /* need to load into memory */ |
602 |
> |
MEMLOAD myMem; /* need to map into memory */ |
603 |
|
if (a == argc-1) { |
604 |
|
if (load_file(&myMem, stdin) <= 0) { |
605 |
|
fprintf(stderr, "%s: error loading file into memory\n", |
613 |
|
} |
614 |
|
if (!do_transpose(&myMem)) |
615 |
|
return(1); |
616 |
< |
/* free_load(&myMem); */ |
617 |
< |
} else if (!do_resize(stdin)) /* just reshaping input */ |
616 |
> |
/* free_load(&myMem); about to exit, so don't bother */ |
617 |
> |
} else if (!do_resize(stdin)) /* reshaping input */ |
618 |
|
return(1); |
619 |
|
return(0); |
620 |
|
userr: |