1 /*
2 * The copyright in this software is being made available under the 2-clauses
3 * BSD License, included below. This software may be subject to other third
4 * party and contributor rights, including patent rights, and no such rights
5 * are granted under this license.
6 *
7 * Copyright (c) 2002-2014, Universite catholique de Louvain (UCL), Belgium
8 * Copyright (c) 2002-2014, Professor Benoit Macq
9 * Copyright (c) 2001-2003, David Janssens
10 * Copyright (c) 2002-2003, Yannick Verschueren
11 * Copyright (c) 2003-2007, Francois-Olivier Devaux
12 * Copyright (c) 2003-2014, Antonin Descampe
13 * Copyright (c) 2005, Herve Drolon, FreeImage Team
14 * Copyright (c) 2007, Callum Lerwick <seg@haxxed.com>
15 * Copyright (c) 2012, Carl Hetherington
16 * All rights reserved.
17 *
18 * Redistribution and use in source and binary forms, with or without
19 * modification, are permitted provided that the following conditions
20 * are met:
21 * 1. Redistributions of source code must retain the above copyright
22 * notice, this list of conditions and the following disclaimer.
23 * 2. Redistributions in binary form must reproduce the above copyright
24 * notice, this list of conditions and the following disclaimer in the
25 * documentation and/or other materials provided with the distribution.
26 *
27 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
28 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
29 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
30 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
31 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
32 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
33 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
34 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
35 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
36 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
37 * POSSIBILITY OF SUCH DAMAGE.
38 */
39
40 #include "opj_includes.h"
41
t1_init_ctxno_zc(OPJ_UINT32 f,OPJ_UINT32 orient)42 static int t1_init_ctxno_zc(OPJ_UINT32 f, OPJ_UINT32 orient)
43 {
44 int h, v, d, n, t, hv;
45 n = 0;
46 h = ((f & T1_SIGMA_3) != 0) + ((f & T1_SIGMA_5) != 0);
47 v = ((f & T1_SIGMA_1) != 0) + ((f & T1_SIGMA_7) != 0);
48 d = ((f & T1_SIGMA_0) != 0) + ((f & T1_SIGMA_2) != 0) + ((
49 f & T1_SIGMA_8) != 0) + ((f & T1_SIGMA_6) != 0);
50
51 switch (orient) {
52 case 2:
53 t = h;
54 h = v;
55 v = t;
56 case 0:
57 case 1:
58 if (!h) {
59 if (!v) {
60 if (!d) {
61 n = 0;
62 } else if (d == 1) {
63 n = 1;
64 } else {
65 n = 2;
66 }
67 } else if (v == 1) {
68 n = 3;
69 } else {
70 n = 4;
71 }
72 } else if (h == 1) {
73 if (!v) {
74 if (!d) {
75 n = 5;
76 } else {
77 n = 6;
78 }
79 } else {
80 n = 7;
81 }
82 } else {
83 n = 8;
84 }
85 break;
86 case 3:
87 hv = h + v;
88 if (!d) {
89 if (!hv) {
90 n = 0;
91 } else if (hv == 1) {
92 n = 1;
93 } else {
94 n = 2;
95 }
96 } else if (d == 1) {
97 if (!hv) {
98 n = 3;
99 } else if (hv == 1) {
100 n = 4;
101 } else {
102 n = 5;
103 }
104 } else if (d == 2) {
105 if (!hv) {
106 n = 6;
107 } else {
108 n = 7;
109 }
110 } else {
111 n = 8;
112 }
113 break;
114 }
115
116 return (T1_CTXNO_ZC + n);
117 }
118
t1_init_ctxno_sc(OPJ_UINT32 f)119 static int t1_init_ctxno_sc(OPJ_UINT32 f)
120 {
121 int hc, vc, n;
122 n = 0;
123
124 hc = opj_int_min(((f & (T1_LUT_SIG_E | T1_LUT_SGN_E)) ==
125 T1_LUT_SIG_E) + ((f & (T1_LUT_SIG_W | T1_LUT_SGN_W)) == T1_LUT_SIG_W),
126 1) - opj_int_min(((f & (T1_LUT_SIG_E | T1_LUT_SGN_E)) ==
127 (T1_LUT_SIG_E | T1_LUT_SGN_E)) +
128 ((f & (T1_LUT_SIG_W | T1_LUT_SGN_W)) ==
129 (T1_LUT_SIG_W | T1_LUT_SGN_W)), 1);
130
131 vc = opj_int_min(((f & (T1_LUT_SIG_N | T1_LUT_SGN_N)) ==
132 T1_LUT_SIG_N) + ((f & (T1_LUT_SIG_S | T1_LUT_SGN_S)) == T1_LUT_SIG_S),
133 1) - opj_int_min(((f & (T1_LUT_SIG_N | T1_LUT_SGN_N)) ==
134 (T1_LUT_SIG_N | T1_LUT_SGN_N)) +
135 ((f & (T1_LUT_SIG_S | T1_LUT_SGN_S)) ==
136 (T1_LUT_SIG_S | T1_LUT_SGN_S)), 1);
137
138 if (hc < 0) {
139 hc = -hc;
140 vc = -vc;
141 }
142 if (!hc) {
143 if (vc == -1) {
144 n = 1;
145 } else if (!vc) {
146 n = 0;
147 } else {
148 n = 1;
149 }
150 } else if (hc == 1) {
151 if (vc == -1) {
152 n = 2;
153 } else if (!vc) {
154 n = 3;
155 } else {
156 n = 4;
157 }
158 }
159
160 return (T1_CTXNO_SC + n);
161 }
162
t1_init_spb(OPJ_UINT32 f)163 static int t1_init_spb(OPJ_UINT32 f)
164 {
165 int hc, vc, n;
166
167 hc = opj_int_min(((f & (T1_LUT_SIG_E | T1_LUT_SGN_E)) ==
168 T1_LUT_SIG_E) + ((f & (T1_LUT_SIG_W | T1_LUT_SGN_W)) == T1_LUT_SIG_W),
169 1) - opj_int_min(((f & (T1_LUT_SIG_E | T1_LUT_SGN_E)) ==
170 (T1_LUT_SIG_E | T1_LUT_SGN_E)) +
171 ((f & (T1_LUT_SIG_W | T1_LUT_SGN_W)) ==
172 (T1_LUT_SIG_W | T1_LUT_SGN_W)), 1);
173
174 vc = opj_int_min(((f & (T1_LUT_SIG_N | T1_LUT_SGN_N)) ==
175 T1_LUT_SIG_N) + ((f & (T1_LUT_SIG_S | T1_LUT_SGN_S)) == T1_LUT_SIG_S),
176 1) - opj_int_min(((f & (T1_LUT_SIG_N | T1_LUT_SGN_N)) ==
177 (T1_LUT_SIG_N | T1_LUT_SGN_N)) +
178 ((f & (T1_LUT_SIG_S | T1_LUT_SGN_S)) ==
179 (T1_LUT_SIG_S | T1_LUT_SGN_S)), 1);
180
181 if (!hc && !vc) {
182 n = 0;
183 } else {
184 n = (!(hc > 0 || (!hc && vc > 0)));
185 }
186
187 return n;
188 }
189
dump_array16(int array[],int size)190 static void dump_array16(int array[], int size)
191 {
192 int i;
193 --size;
194 for (i = 0; i < size; ++i) {
195 printf("0x%04x,", array[i]);
196 if (!((i + 1) & 0x7)) {
197 printf("\n ");
198 } else {
199 printf(" ");
200 }
201 }
202 printf("0x%04x\n};\n\n", array[size]);
203 }
204
main(int argc,char ** argv)205 int main(int argc, char **argv)
206 {
207 unsigned int i, j;
208 double u, v, t;
209
210 int lut_ctxno_zc[2048];
211 int lut_nmsedec_sig[1 << T1_NMSEDEC_BITS];
212 int lut_nmsedec_sig0[1 << T1_NMSEDEC_BITS];
213 int lut_nmsedec_ref[1 << T1_NMSEDEC_BITS];
214 int lut_nmsedec_ref0[1 << T1_NMSEDEC_BITS];
215 (void)argc;
216 (void)argv;
217
218 printf("/* This file was automatically generated by t1_generate_luts.c */\n\n");
219
220 /* lut_ctxno_zc */
221 for (j = 0; j < 4; ++j) {
222 for (i = 0; i < 512; ++i) {
223 OPJ_UINT32 orient = j;
224 if (orient == 2) {
225 orient = 1;
226 } else if (orient == 1) {
227 orient = 2;
228 }
229 lut_ctxno_zc[(orient << 9) | i] = t1_init_ctxno_zc(i, j);
230 }
231 }
232
233 printf("static const OPJ_BYTE lut_ctxno_zc[2048] = {\n ");
234 for (i = 0; i < 2047; ++i) {
235 printf("%i,", lut_ctxno_zc[i]);
236 if (!((i + 1) & 0x1f)) {
237 printf("\n ");
238 } else {
239 printf(" ");
240 }
241 }
242 printf("%i\n};\n\n", lut_ctxno_zc[2047]);
243
244 /* lut_ctxno_sc */
245 printf("static const OPJ_BYTE lut_ctxno_sc[256] = {\n ");
246 for (i = 0; i < 255; ++i) {
247 printf("0x%x,", t1_init_ctxno_sc(i));
248 if (!((i + 1) & 0xf)) {
249 printf("\n ");
250 } else {
251 printf(" ");
252 }
253 }
254 printf("0x%x\n};\n\n", t1_init_ctxno_sc(255));
255
256 /* lut_spb */
257 printf("static const OPJ_BYTE lut_spb[256] = {\n ");
258 for (i = 0; i < 255; ++i) {
259 printf("%i,", t1_init_spb(i));
260 if (!((i + 1) & 0x1f)) {
261 printf("\n ");
262 } else {
263 printf(" ");
264 }
265 }
266 printf("%i\n};\n\n", t1_init_spb(255));
267
268 /* FIXME FIXME FIXME */
269 /* fprintf(stdout,"nmsedec luts:\n"); */
270 for (i = 0U; i < (1U << T1_NMSEDEC_BITS); ++i) {
271 t = i / pow(2, T1_NMSEDEC_FRACBITS);
272 u = t;
273 v = t - 1.5;
274 lut_nmsedec_sig[i] =
275 opj_int_max(0,
276 (int)(floor((u * u - v * v) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2,
277 T1_NMSEDEC_FRACBITS) * 8192.0));
278 lut_nmsedec_sig0[i] =
279 opj_int_max(0,
280 (int)(floor((u * u) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2,
281 T1_NMSEDEC_FRACBITS) * 8192.0));
282 u = t - 1.0;
283 if (i & (1 << (T1_NMSEDEC_BITS - 1))) {
284 v = t - 1.5;
285 } else {
286 v = t - 0.5;
287 }
288 lut_nmsedec_ref[i] =
289 opj_int_max(0,
290 (int)(floor((u * u - v * v) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2,
291 T1_NMSEDEC_FRACBITS) * 8192.0));
292 lut_nmsedec_ref0[i] =
293 opj_int_max(0,
294 (int)(floor((u * u) * pow(2, T1_NMSEDEC_FRACBITS) + 0.5) / pow(2,
295 T1_NMSEDEC_FRACBITS) * 8192.0));
296 }
297
298 printf("static const OPJ_INT16 lut_nmsedec_sig[1U << T1_NMSEDEC_BITS] = {\n ");
299 dump_array16(lut_nmsedec_sig, 1U << T1_NMSEDEC_BITS);
300
301 printf("static const OPJ_INT16 lut_nmsedec_sig0[1U << T1_NMSEDEC_BITS] = {\n ");
302 dump_array16(lut_nmsedec_sig0, 1U << T1_NMSEDEC_BITS);
303
304 printf("static const OPJ_INT16 lut_nmsedec_ref[1U << T1_NMSEDEC_BITS] = {\n ");
305 dump_array16(lut_nmsedec_ref, 1U << T1_NMSEDEC_BITS);
306
307 printf("static const OPJ_INT16 lut_nmsedec_ref0[1U << T1_NMSEDEC_BITS] = {\n ");
308 dump_array16(lut_nmsedec_ref0, 1U << T1_NMSEDEC_BITS);
309
310 return 0;
311 }
312