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
42static 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
119static 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
163static 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
190static 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
205int 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