aboutsummaryrefslogtreecommitdiffstatshomepage
path: root/libraries/openjpeg-libsl/libopenjpeg/t1_generate_luts.c
diff options
context:
space:
mode:
Diffstat (limited to 'libraries/openjpeg-libsl/libopenjpeg/t1_generate_luts.c')
-rw-r--r--libraries/openjpeg-libsl/libopenjpeg/t1_generate_luts.c295
1 files changed, 295 insertions, 0 deletions
diff --git a/libraries/openjpeg-libsl/libopenjpeg/t1_generate_luts.c b/libraries/openjpeg-libsl/libopenjpeg/t1_generate_luts.c
new file mode 100644
index 0000000..62cb158
--- /dev/null
+++ b/libraries/openjpeg-libsl/libopenjpeg/t1_generate_luts.c
@@ -0,0 +1,295 @@
1/*
2 * Copyright (c) 2002-2007, Communications and Remote Sensing Laboratory, Universite catholique de Louvain (UCL), Belgium
3 * Copyright (c) 2002-2007, Professor Benoit Macq
4 * Copyright (c) 2001-2003, David Janssens
5 * Copyright (c) 2002-2003, Yannick Verschueren
6 * Copyright (c) 2003-2007, Francois-Olivier Devaux and Antonin Descampe
7 * Copyright (c) 2005, Herve Drolon, FreeImage Team
8 * Copyright (c) 2007, Callum Lerwick <seg@haxxed.com>
9 * All rights reserved.
10 *
11 * Redistribution and use in source and binary forms, with or without
12 * modification, are permitted provided that the following conditions
13 * are met:
14 * 1. Redistributions of source code must retain the above copyright
15 * notice, this list of conditions and the following disclaimer.
16 * 2. Redistributions in binary form must reproduce the above copyright
17 * notice, this list of conditions and the following disclaimer in the
18 * documentation and/or other materials provided with the distribution.
19 *
20 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
21 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
22 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
23 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
24 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
25 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
26 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
27 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
28 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
29 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
30 * POSSIBILITY OF SUCH DAMAGE.
31 */
32
33#include "opj_includes.h"
34#include <math.h>
35
36static int t1_init_ctxno_zc(int f, int orient) {
37 int h, v, d, n, t, hv;
38 n = 0;
39 h = ((f & T1_SIG_W) != 0) + ((f & T1_SIG_E) != 0);
40 v = ((f & T1_SIG_N) != 0) + ((f & T1_SIG_S) != 0);
41 d = ((f & T1_SIG_NW) != 0) + ((f & T1_SIG_NE) != 0) + ((f & T1_SIG_SE) != 0) + ((f & T1_SIG_SW) != 0);
42
43 switch (orient) {
44 case 2:
45 t = h;
46 h = v;
47 v = t;
48 case 0:
49 case 1:
50 if (!h) {
51 if (!v) {
52 if (!d)
53 n = 0;
54 else if (d == 1)
55 n = 1;
56 else
57 n = 2;
58 } else if (v == 1) {
59 n = 3;
60 } else {
61 n = 4;
62 }
63 } else if (h == 1) {
64 if (!v) {
65 if (!d)
66 n = 5;
67 else
68 n = 6;
69 } else {
70 n = 7;
71 }
72 } else
73 n = 8;
74 break;
75 case 3:
76 hv = h + v;
77 if (!d) {
78 if (!hv) {
79 n = 0;
80 } else if (hv == 1) {
81 n = 1;
82 } else {
83 n = 2;
84 }
85 } else if (d == 1) {
86 if (!hv) {
87 n = 3;
88 } else if (hv == 1) {
89 n = 4;
90 } else {
91 n = 5;
92 }
93 } else if (d == 2) {
94 if (!hv) {
95 n = 6;
96 } else {
97 n = 7;
98 }
99 } else {
100 n = 8;
101 }
102 break;
103 }
104
105 return (T1_CTXNO_ZC + n);
106}
107
108static int t1_init_ctxno_sc(int f) {
109 int hc, vc, n;
110 n = 0;
111
112 hc = int_min(((f & (T1_SIG_E | T1_SGN_E)) ==
113 T1_SIG_E) + ((f & (T1_SIG_W | T1_SGN_W)) == T1_SIG_W),
114 1) - int_min(((f & (T1_SIG_E | T1_SGN_E)) ==
115 (T1_SIG_E | T1_SGN_E)) +
116 ((f & (T1_SIG_W | T1_SGN_W)) ==
117 (T1_SIG_W | T1_SGN_W)), 1);
118
119 vc = int_min(((f & (T1_SIG_N | T1_SGN_N)) ==
120 T1_SIG_N) + ((f & (T1_SIG_S | T1_SGN_S)) == T1_SIG_S),
121 1) - int_min(((f & (T1_SIG_N | T1_SGN_N)) ==
122 (T1_SIG_N | T1_SGN_N)) +
123 ((f & (T1_SIG_S | T1_SGN_S)) ==
124 (T1_SIG_S | T1_SGN_S)), 1);
125
126 if (hc < 0) {
127 hc = -hc;
128 vc = -vc;
129 }
130 if (!hc) {
131 if (vc == -1)
132 n = 1;
133 else if (!vc)
134 n = 0;
135 else
136 n = 1;
137 } else if (hc == 1) {
138 if (vc == -1)
139 n = 2;
140 else if (!vc)
141 n = 3;
142 else
143 n = 4;
144 }
145
146 return (T1_CTXNO_SC + n);
147}
148
149static int t1_init_ctxno_mag(int f) {
150 int n;
151 if (!(f & T1_REFINE))
152 n = (f & (T1_SIG_OTH)) ? 1 : 0;
153 else
154 n = 2;
155
156 return (T1_CTXNO_MAG + n);
157}
158
159static int t1_init_spb(int f) {
160 int hc, vc, n;
161
162 hc = int_min(((f & (T1_SIG_E | T1_SGN_E)) ==
163 T1_SIG_E) + ((f & (T1_SIG_W | T1_SGN_W)) == T1_SIG_W),
164 1) - int_min(((f & (T1_SIG_E | T1_SGN_E)) ==
165 (T1_SIG_E | T1_SGN_E)) +
166 ((f & (T1_SIG_W | T1_SGN_W)) ==
167 (T1_SIG_W | T1_SGN_W)), 1);
168
169 vc = int_min(((f & (T1_SIG_N | T1_SGN_N)) ==
170 T1_SIG_N) + ((f & (T1_SIG_S | T1_SGN_S)) == T1_SIG_S),
171 1) - int_min(((f & (T1_SIG_N | T1_SGN_N)) ==
172 (T1_SIG_N | T1_SGN_N)) +
173 ((f & (T1_SIG_S | T1_SGN_S)) ==
174 (T1_SIG_S | T1_SGN_S)), 1);
175
176 if (!hc && !vc)
177 n = 0;
178 else
179 n = (!(hc > 0 || (!hc && vc > 0)));
180
181 return n;
182}
183
184void dump_array16(int array[],int size){
185 int i;
186 --size;
187 for (i = 0; i < size; ++i) {
188 printf("0x%04x, ", array[i]);
189 if(!((i+1)&0x7))
190 printf("\n ");
191 }
192 printf("0x%04x\n};\n\n", array[size]);
193}
194
195int main(){
196 int i, j;
197 double u, v, t;
198
199 int lut_ctxno_zc[1024];
200 int lut_ctxno_mag[4096];
201 int lut_nmsedec_sig[1 << T1_NMSEDEC_BITS];
202 int lut_nmsedec_sig0[1 << T1_NMSEDEC_BITS];
203 int lut_nmsedec_ref[1 << T1_NMSEDEC_BITS];
204 int lut_nmsedec_ref0[1 << T1_NMSEDEC_BITS];
205
206 printf("/* This file was automatically generated by t1_generate_luts.c */\n\n");
207
208 // lut_ctxno_zc
209 for (j = 0; j < 4; ++j) {
210 for (i = 0; i < 256; ++i) {
211 lut_ctxno_zc[(j << 8) | i] = t1_init_ctxno_zc(i, j);
212 }
213 }
214
215 printf("static int8_t lut_ctxno_zc[1024] = {\n ");
216 for (i = 0; i < 1023; ++i) {
217 printf("%i, ", lut_ctxno_zc[i]);
218 if(!((i+1)&0x1f))
219 printf("\n ");
220 }
221 printf("%i\n};\n\n", lut_ctxno_zc[1023]);
222
223 // lut_ctxno_sc
224 printf("static int8_t lut_ctxno_sc[256] = {\n ");
225 for (i = 0; i < 255; ++i) {
226 printf("0x%x, ", t1_init_ctxno_sc(i << 4));
227 if(!((i+1)&0xf))
228 printf("\n ");
229 }
230 printf("0x%x\n};\n\n", t1_init_ctxno_sc(255 << 4));
231
232 // lut_ctxno_mag
233 for (j = 0; j < 2; ++j) {
234 for (i = 0; i < 2048; ++i) {
235 lut_ctxno_mag[(j << 11) + i] = t1_init_ctxno_mag((j ? T1_REFINE : 0) | i);
236 }
237 }
238
239 printf("static int8_t lut_ctxno_mag[4096] = {\n ");
240 for (i = 0; i < 4095; ++i) {
241 printf("%i, ", lut_ctxno_mag[i]);
242 if(!((i+1)&0xf))
243 printf("\n ");
244 }
245 printf("%i\n};\n\n", lut_ctxno_mag[4095]);
246
247 // lut_spb
248 printf("static int8_t lut_spb[256] = {\n ");
249 for (i = 0; i < 255; ++i) {
250 printf("%i, ", t1_init_spb(i << 4));
251 if(!((i+1)&0x1f))
252 printf("\n ");
253 }
254 printf("%i\n};\n\n", t1_init_spb(255 << 4));
255
256 /* FIXME FIXME FIXME */
257 /* fprintf(stdout,"nmsedec luts:\n"); */
258 for (i = 0; i < (1 << T1_NMSEDEC_BITS); ++i) {
259 t = i / pow(2, T1_NMSEDEC_FRACBITS);
260 u = t;
261 v = t - 1.5;
262 lut_nmsedec_sig[i] =
263 int_max(0,
264 (int) (floor((u * u - v * v) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2, T1_NMSEDEC_FRACBITS) * 8192.0));
265 lut_nmsedec_sig0[i] =
266 int_max(0,
267 (int) (floor((u * u) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2, T1_NMSEDEC_FRACBITS) * 8192.0));
268 u = t - 1.0;
269 if (i & (1 << (T1_NMSEDEC_BITS - 1))) {
270 v = t - 1.5;
271 } else {
272 v = t - 0.5;
273 }
274 lut_nmsedec_ref[i] =
275 int_max(0,
276 (int) (floor((u * u - v * v) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2, T1_NMSEDEC_FRACBITS) * 8192.0));
277 lut_nmsedec_ref0[i] =
278 int_max(0,
279 (int) (floor((u * u) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2, T1_NMSEDEC_FRACBITS) * 8192.0));
280 }
281
282 printf("static int16_t lut_nmsedec_sig[1 << T1_NMSEDEC_BITS] = {\n ");
283 dump_array16(&lut_nmsedec_sig, 1 << T1_NMSEDEC_BITS);
284
285 printf("static int16_t lut_nmsedec_sig0[1 << T1_NMSEDEC_BITS] = {\n ");
286 dump_array16(&lut_nmsedec_sig0, 1 << T1_NMSEDEC_BITS);
287
288 printf("static int16_t lut_nmsedec_ref[1 << T1_NMSEDEC_BITS] = {\n ");
289 dump_array16(&lut_nmsedec_ref, 1 << T1_NMSEDEC_BITS);
290
291 printf("static int16_t lut_nmsedec_ref0[1 << T1_NMSEDEC_BITS] = {\n ");
292 dump_array16(&lut_nmsedec_ref0, 1 << T1_NMSEDEC_BITS);
293
294 return 0;
295}