1 |
|
|
/* |
2 |
|
|
* JPEG 2000 DSP functions |
3 |
|
|
* Copyright (c) 2007 Kamil Nowosad |
4 |
|
|
* Copyright (c) 2013 Nicolas Bertrand <nicoinattendu@gmail.com> |
5 |
|
|
* |
6 |
|
|
* This file is part of FFmpeg. |
7 |
|
|
* |
8 |
|
|
* FFmpeg is free software; you can redistribute it and/or |
9 |
|
|
* modify it under the terms of the GNU Lesser General Public |
10 |
|
|
* License as published by the Free Software Foundation; either |
11 |
|
|
* version 2.1 of the License, or (at your option) any later version. |
12 |
|
|
* |
13 |
|
|
* FFmpeg is distributed in the hope that it will be useful, |
14 |
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of |
15 |
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
16 |
|
|
* Lesser General Public License for more details. |
17 |
|
|
* |
18 |
|
|
* You should have received a copy of the GNU Lesser General Public |
19 |
|
|
* License along with FFmpeg; if not, write to the Free Software |
20 |
|
|
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
21 |
|
|
*/ |
22 |
|
|
|
23 |
|
|
#include "config.h" |
24 |
|
|
#include "libavutil/attributes.h" |
25 |
|
|
#include "jpeg2000dsp.h" |
26 |
|
|
|
27 |
|
|
/* Inverse ICT parameters in float and integer. |
28 |
|
|
* int value = (float value) * (1<<16) */ |
29 |
|
|
static const float f_ict_params[4] = { |
30 |
|
|
1.402f, |
31 |
|
|
0.34413f, |
32 |
|
|
0.71414f, |
33 |
|
|
1.772f |
34 |
|
|
}; |
35 |
|
|
|
36 |
|
|
static const int i_ict_params[4] = { |
37 |
|
|
91881, |
38 |
|
|
22553, |
39 |
|
|
46802, |
40 |
|
|
116130 |
41 |
|
|
}; |
42 |
|
|
|
43 |
|
3 |
static void ict_float(void *_src0, void *_src1, void *_src2, int csize) |
44 |
|
|
{ |
45 |
|
3 |
float *src0 = _src0, *src1 = _src1, *src2 = _src2; |
46 |
|
|
float i0f, i1f, i2f; |
47 |
|
|
int i; |
48 |
|
|
|
49 |
✓✓ |
1539 |
for (i = 0; i < csize; i++) { |
50 |
|
1536 |
i0f = *src0 + (f_ict_params[0] * *src2); |
51 |
|
1536 |
i1f = *src0 - (f_ict_params[1] * *src1) |
52 |
|
1536 |
- (f_ict_params[2] * *src2); |
53 |
|
1536 |
i2f = *src0 + (f_ict_params[3] * *src1); |
54 |
|
1536 |
*src0++ = i0f; |
55 |
|
1536 |
*src1++ = i1f; |
56 |
|
1536 |
*src2++ = i2f; |
57 |
|
|
} |
58 |
|
3 |
} |
59 |
|
|
|
60 |
|
3 |
static void ict_int(void *_src0, void *_src1, void *_src2, int csize) |
61 |
|
|
{ |
62 |
|
3 |
int32_t *src0 = _src0, *src1 = _src1, *src2 = _src2; |
63 |
|
|
int32_t i0, i1, i2; |
64 |
|
|
int i; |
65 |
|
|
|
66 |
✓✓ |
6220803 |
for (i = 0; i < csize; i++) { |
67 |
|
6220800 |
i0 = *src0 + *src2 + ((int)((26345U * *src2) + (1 << 15)) >> 16); |
68 |
|
6220800 |
i1 = *src0 - ((int)(((unsigned)i_ict_params[1] * *src1) + (1 << 15)) >> 16) |
69 |
|
6220800 |
- ((int)(((unsigned)i_ict_params[2] * *src2) + (1 << 15)) >> 16); |
70 |
|
6220800 |
i2 = *src0 + (2 * *src1) + ((int)((-14942U * *src1) + (1 << 15)) >> 16); |
71 |
|
6220800 |
*src0++ = i0; |
72 |
|
6220800 |
*src1++ = i1; |
73 |
|
6220800 |
*src2++ = i2; |
74 |
|
|
} |
75 |
|
3 |
} |
76 |
|
|
|
77 |
|
3 |
static void rct_int(void *_src0, void *_src1, void *_src2, int csize) |
78 |
|
|
{ |
79 |
|
3 |
int32_t *src0 = _src0, *src1 = _src1, *src2 = _src2; |
80 |
|
|
int32_t i0, i1, i2; |
81 |
|
|
int i; |
82 |
|
|
|
83 |
✓✓ |
1539 |
for (i = 0; i < csize; i++) { |
84 |
|
1536 |
i1 = *src0 - (*src2 + *src1 >> 2); |
85 |
|
1536 |
i0 = i1 + *src2; |
86 |
|
1536 |
i2 = i1 + *src1; |
87 |
|
1536 |
*src0++ = i0; |
88 |
|
1536 |
*src1++ = i1; |
89 |
|
1536 |
*src2++ = i2; |
90 |
|
|
} |
91 |
|
3 |
} |
92 |
|
|
|
93 |
|
32 |
av_cold void ff_jpeg2000dsp_init(Jpeg2000DSPContext *c) |
94 |
|
|
{ |
95 |
|
32 |
c->mct_decode[FF_DWT97] = ict_float; |
96 |
|
32 |
c->mct_decode[FF_DWT53] = rct_int; |
97 |
|
32 |
c->mct_decode[FF_DWT97_INT] = ict_int; |
98 |
|
|
|
99 |
|
|
if (ARCH_X86) |
100 |
|
32 |
ff_jpeg2000dsp_init_x86(c); |
101 |
|
32 |
} |