-
Notifications
You must be signed in to change notification settings - Fork 27
/
trle.c
125 lines (112 loc) · 5.25 KB
/
trle.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
/**
Copyright (C) powturbo 2015-2023
SPDX-License-Identifier: GPL v2 License
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License along
with this program; if not, write to the Free Software Foundation, Inc.,
51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
- email : powturbo [AT] gmail.com
- github : https://github.com/powturbo
- homepage : https://sites.google.com/site/powturbo/
- twitter : https://twitter.com/powturbo
TurboRLE - "Most efficient and fastest Run Length Encoding"
**/
#include <string.h>
#include <stdio.h>
#include <stdlib.h>
#ifdef __APPLE__
#include <sys/malloc.h>
#else
#include <malloc.h>
#endif
#ifdef _MSC_VER
#include "vs/getopt.h"
#else
#include <getopt.h>
#endif
#include "include/trle.h"
#include "include_/time_.h"
#include "ext/mrle.c"
int memcheck(unsigned char *in, unsigned n, unsigned char *cpy) {
int i;
for(i = 0; i < n; i++)
if(in[i] != cpy[i]) {
printf("ERROR in[%d]=%x, dec[%d]=%x\n", i, in[i], i, cpy[i]);
return i+1;
}
return 0;
}
void libmemcpy(unsigned char *dst, unsigned char *src, int len) {
void *(*memcpy_ptr)(void *, const void *, size_t) = memcpy;
if (time(NULL) == 1)
memcpy_ptr = NULL;
memcpy_ptr(dst, src, len);
}
#define ID_MEMCPY 8
unsigned bench(unsigned char *in, unsigned n, unsigned char *out, unsigned char *cpy, int id) {
unsigned l;
memrcpy(cpy,in,n);
switch(id) {
case 1: TMBENCH("\ntrle ", l=trlec( in, n, out) ,n); printf("%10u %5.1f%%", l, (double)l*100.0/n); TMBENCH("",trled( out, l, cpy, n) ,n); break;
case 2: TMBENCH("\nsrle 0 ", l=srlec( in, n, out) ,n); printf("%10u %5.1f%%", l, (double)l*100.0/n); TMBENCH("",srled( out, l, cpy, n) ,n); break;
case 3: TMBENCH("\nsrle 8 ", l=srlec8( in, n, out, 0xda) ,n); printf("%10u %5.1f%%", l, (double)l*100.0/n); TMBENCH("",srled8( out, l, cpy, n, 0xda) ,n); break;
case 4: TMBENCH("\nsrle 16", l=srlec16(in, n, out, 0xdada) ,n); printf("%10u %5.1f%%", l, (double)l*100.0/n); TMBENCH("",srled16(out, l, cpy, n, 0xdada) ,n); break;
case 5: TMBENCH("\nsrle 32", l=srlec32(in, n, out, 0xdadadadau) ,n); printf("%10u %5.1f%%", l, (double)l*100.0/n); TMBENCH("",srled32(out, l, cpy, n, 0xdadadadau) ,n); break;
case 6: TMBENCH("\nsrle 64", l=srlec64(in, n, out, 0xdadadadadadadadaull) ,n); printf("%10u %5.1f%%", l, (double)l*100.0/n); TMBENCH("",srled64(out, l, cpy, n, 0xdadadadadadadadaull) ,n); break;
case 7: TMBENCH("\nmrle ", l=mrlec( in, n, out), n); printf("%10u %5.1f%%", l, (double)l*100.0/n); TMBENCH("",mrled( out, cpy, n), n); break;
case 8: TMBENCH("\nmemcpy ", libmemcpy(in,out,n) ,n); printf("%10u %5.1f%%", n, (double)100.0); return n;
default: return l;
}
memcheck(in,n,cpy);
return l;
}
int main(int argc, char* argv[]) {
unsigned trips = 4,cmp=1, b = 1 << 30, esize=4, lz=0, fno,id=0;
unsigned char *in,*out,*cpy;
int c, digit_optind = 0, this_option_optind = optind ? optind : 1, option_index = 0;
static struct option long_options[] = { {"blocsize", 0, 0, 'b'}, {0, 0, 0} };
for(;;) {
if((c = getopt_long(argc, argv, "B:ce:i:I:", long_options, &option_index)) == -1) break;
switch(c) {
case 0 : printf("Option %s", long_options[option_index].name); if(optarg) printf (" with arg %s", optarg); printf ("\n"); break;
case 'e': id = atoi(optarg); break;
case 'i': tm_rep = atoi(optarg); if(!tm_rep) tm_rep=1,trips=1; break;
case 'I': trips = atoi(optarg); if(!trips) trips=1; break;
case 'B': b = argtoi(optarg,1); break;
case 'c': cmp++; break;
}
}
if(argc - optind < 1) { fprintf(stderr, "File not specified\n"); exit(-1); }
for(fno = optind; fno < argc; fno++) {
char *inname = argv[fno];
int n,i;
long long flen;
FILE *fi = fopen(inname, "rb"); if(!fi ) { perror(inname); continue; }
fseek(fi, 0, SEEK_END);
flen = ftell(fi);
fseek(fi, 0, SEEK_SET);
if(flen > b) flen = b;
n = flen;
if(!(in = (unsigned char*)malloc(n+1024))) { fprintf(stderr, "malloc error\n"); exit(-1); } cpy = in;
if(!(out = (unsigned char*)malloc(flen*4/3+1024))) { fprintf(stderr, "malloc error\n"); exit(-1); }
if(cmp && !(cpy = (unsigned char*)malloc(n+1024))) { fprintf(stderr, "malloc error\n"); exit(-1); }
n = fread(in, 1, n, fi); printf("File='%s' Length=%u\n", inname, n);
fclose(fi);
if(n <= 0) exit(0);
if(fno == optind)
tm_init(trips, 1);
if(!id) {
printf("function E MB/s size ratio D MB/s");
for(i=1; i <= ID_MEMCPY; i++) bench(in,n,out,cpy,i);
} else
bench(in,n,out,cpy,id);
printf("\n");
}
}