annotate faandct.c @ 7246:c9dc58950620 libavcodec

10l, forgot to replace a / 256 for r14173. Fixes MS ADPCM regression test.
author reimar
date Sat, 12 Jul 2008 07:28:45 +0000
parents 322023e630a6
children 39b409447bfa
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
1 /*
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
2 * Floating point AAN DCT
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
3 * Copyright (c) 2003 Michael Niedermayer <michaelni@gmx.at>
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
4 *
5214
470601203f44 Group all copyright and author notices together.
diego
parents: 4283
diff changeset
5 * this implementation is based upon the IJG integer AAN DCT (see jfdctfst.c)
470601203f44 Group all copyright and author notices together.
diego
parents: 4283
diff changeset
6 *
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
7 * This file is part of FFmpeg.
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
8 *
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
9 * FFmpeg is free software; you can redistribute it and/or
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
10 * modify it under the terms of the GNU Lesser General Public
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
11 * License as published by the Free Software Foundation; either
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
12 * version 2.1 of the License, or (at your option) any later version.
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
13 *
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
14 * FFmpeg is distributed in the hope that it will be useful,
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
17 * Lesser General Public License for more details.
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
18 *
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
19 * You should have received a copy of the GNU Lesser General Public
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
20 * License along with FFmpeg; if not, write to the Free Software
3036
0b546eab515d Update licensing information: The FSF changed postal address.
diego
parents: 2979
diff changeset
21 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
6397
094f507b952b Allow use of my AAN DCT under BSD, AMD has agreed to write and contribute
michael
parents: 5214
diff changeset
22 *
094f507b952b Allow use of my AAN DCT under BSD, AMD has agreed to write and contribute
michael
parents: 5214
diff changeset
23 * The AAN DCT in this file except ff_faandct248() can also be used under the
6398
0a50abd54b7f Clarify which BSD license.
michael
parents: 6397
diff changeset
24 * new (3 clause) BSD license.
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
25 */
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
26
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
27 /**
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
28 * @file faandct.c
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
29 * @brief
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
30 * Floating point AAN DCT
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
31 * @author Michael Niedermayer <michaelni@gmx.at>
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
32 */
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
33
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
34 #include "dsputil.h"
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
35 #include "faandct.h"
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
36
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
37 #define FLOAT float
1562
bf452704100f optionally merge postscale into quantization table for the float aan dct
michael
parents: 1561
diff changeset
38 #ifdef FAAN_POSTSCALE
bf452704100f optionally merge postscale into quantization table for the float aan dct
michael
parents: 1561
diff changeset
39 # define SCALE(x) postscale[x]
bf452704100f optionally merge postscale into quantization table for the float aan dct
michael
parents: 1561
diff changeset
40 #else
bf452704100f optionally merge postscale into quantization table for the float aan dct
michael
parents: 1561
diff changeset
41 # define SCALE(x) 1
bf452704100f optionally merge postscale into quantization table for the float aan dct
michael
parents: 1561
diff changeset
42 #endif
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
43
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
44 //numbers generated by simple c code (not as accurate as they could be)
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
45 /*
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
46 for(i=0; i<8; i++){
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
47 printf("#define B%d %1.20llf\n", i, (long double)1.0/(cosl(i*acosl(-1.0)/(long double)16.0)*sqrtl(2)));
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
48 }
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
49 */
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
50 #define B0 1.00000000000000000000
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
51 #define B1 0.72095982200694791383 // (cos(pi*1/16)sqrt(2))^-1
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
52 #define B2 0.76536686473017954350 // (cos(pi*2/16)sqrt(2))^-1
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
53 #define B3 0.85043009476725644878 // (cos(pi*3/16)sqrt(2))^-1
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
54 #define B4 1.00000000000000000000 // (cos(pi*4/16)sqrt(2))^-1
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
55 #define B5 1.27275858057283393842 // (cos(pi*5/16)sqrt(2))^-1
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
56 #define B6 1.84775906502257351242 // (cos(pi*6/16)sqrt(2))^-1
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
57 #define B7 3.62450978541155137218 // (cos(pi*7/16)sqrt(2))^-1
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
58
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
59
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
60 #define A1 0.70710678118654752438 // cos(pi*4/16)
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
61 #define A2 0.54119610014619698435 // cos(pi*6/16)sqrt(2)
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
62 #define A5 0.38268343236508977170 // cos(pi*6/16)
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
63 #define A4 1.30656296487637652774 // cos(pi*2/16)sqrt(2)
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
64
7129
322023e630a6 mark read-only data as const
stefang
parents: 6471
diff changeset
65 static const FLOAT postscale[64]={
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
66 B0*B0, B0*B1, B0*B2, B0*B3, B0*B4, B0*B5, B0*B6, B0*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
67 B1*B0, B1*B1, B1*B2, B1*B3, B1*B4, B1*B5, B1*B6, B1*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
68 B2*B0, B2*B1, B2*B2, B2*B3, B2*B4, B2*B5, B2*B6, B2*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
69 B3*B0, B3*B1, B3*B2, B3*B3, B3*B4, B3*B5, B3*B6, B3*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
70 B4*B0, B4*B1, B4*B2, B4*B3, B4*B4, B4*B5, B4*B6, B4*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
71 B5*B0, B5*B1, B5*B2, B5*B3, B5*B4, B5*B5, B5*B6, B5*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
72 B6*B0, B6*B1, B6*B2, B6*B3, B6*B4, B6*B5, B6*B6, B6*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
73 B7*B0, B7*B1, B7*B2, B7*B3, B7*B4, B7*B5, B7*B6, B7*B7,
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
74 };
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
75
4283
d6f83e2f8804 rename always_inline to av_always_inline and move to common.h
mru
parents: 3947
diff changeset
76 static av_always_inline void row_fdct(FLOAT temp[64], DCTELEM * data)
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
77 {
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
78 FLOAT tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
79 FLOAT tmp10, tmp11, tmp12, tmp13;
6471
b72f6c4cee12 Fix unused variable warnings.
diego
parents: 6411
diff changeset
80 FLOAT z2, z4, z11, z13;
b72f6c4cee12 Fix unused variable warnings.
diego
parents: 6411
diff changeset
81 FLOAT av_unused z5;
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
82 int i;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
83
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
84 for (i=0; i<8*8; i+=8) {
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
85 tmp0= data[0 + i] + data[7 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
86 tmp7= data[0 + i] - data[7 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
87 tmp1= data[1 + i] + data[6 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
88 tmp6= data[1 + i] - data[6 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
89 tmp2= data[2 + i] + data[5 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
90 tmp5= data[2 + i] - data[5 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
91 tmp3= data[3 + i] + data[4 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
92 tmp4= data[3 + i] - data[4 + i];
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
93
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
94 tmp10= tmp0 + tmp3;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
95 tmp13= tmp0 - tmp3;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
96 tmp11= tmp1 + tmp2;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
97 tmp12= tmp1 - tmp2;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
98
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
99 temp[0 + i]= tmp10 + tmp11;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
100 temp[4 + i]= tmp10 - tmp11;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
101
6408
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
102 tmp12 += tmp13;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
103 tmp12 *= A1;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
104 temp[2 + i]= tmp13 + tmp12;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
105 temp[6 + i]= tmp13 - tmp12;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
106
6409
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
107 tmp4 += tmp5;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
108 tmp5 += tmp6;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
109 tmp6 += tmp7;
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
110
6411
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
111 #if 0
6409
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
112 z5= (tmp4 - tmp6) * A5;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
113 z2= tmp4*A2 + z5;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
114 z4= tmp6*A4 + z5;
6411
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
115 #else
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
116 z2= tmp4*(A2+A5) - tmp6*A5;
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
117 z4= tmp6*(A4-A5) + tmp4*A5;
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
118 #endif
6410
72e573046758 Get rid of another useless variable.
michael
parents: 6409
diff changeset
119 tmp5*=A1;
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
120
6410
72e573046758 Get rid of another useless variable.
michael
parents: 6409
diff changeset
121 z11= tmp7 + tmp5;
72e573046758 Get rid of another useless variable.
michael
parents: 6409
diff changeset
122 z13= tmp7 - tmp5;
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
123
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
124 temp[5 + i]= z13 + z2;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
125 temp[3 + i]= z13 - z2;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
126 temp[1 + i]= z11 + z4;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
127 temp[7 + i]= z11 - z4;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
128 }
1589
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
129 }
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
130
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
131 void ff_faandct(DCTELEM * data)
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
132 {
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
133 FLOAT tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
134 FLOAT tmp10, tmp11, tmp12, tmp13;
6471
b72f6c4cee12 Fix unused variable warnings.
diego
parents: 6411
diff changeset
135 FLOAT z2, z4, z11, z13;
b72f6c4cee12 Fix unused variable warnings.
diego
parents: 6411
diff changeset
136 FLOAT av_unused z5;
1589
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
137 FLOAT temp[64];
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
138 int i;
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
139
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
140 emms_c();
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
141
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
142 row_fdct(temp, data);
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
143
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
144 for (i=0; i<8; i++) {
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
145 tmp0= temp[8*0 + i] + temp[8*7 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
146 tmp7= temp[8*0 + i] - temp[8*7 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
147 tmp1= temp[8*1 + i] + temp[8*6 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
148 tmp6= temp[8*1 + i] - temp[8*6 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
149 tmp2= temp[8*2 + i] + temp[8*5 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
150 tmp5= temp[8*2 + i] - temp[8*5 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
151 tmp3= temp[8*3 + i] + temp[8*4 + i];
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
152 tmp4= temp[8*3 + i] - temp[8*4 + i];
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
153
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
154 tmp10= tmp0 + tmp3;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
155 tmp13= tmp0 - tmp3;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
156 tmp11= tmp1 + tmp2;
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
157 tmp12= tmp1 - tmp2;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
158
1563
820e06c6ca9b 4.9-RC FreeBSD doesnt like lrint() so change to lrintf(), not that bsd supports that but we emulate it ...
michael
parents: 1562
diff changeset
159 data[8*0 + i]= lrintf(SCALE(8*0 + i) * (tmp10 + tmp11));
820e06c6ca9b 4.9-RC FreeBSD doesnt like lrint() so change to lrintf(), not that bsd supports that but we emulate it ...
michael
parents: 1562
diff changeset
160 data[8*4 + i]= lrintf(SCALE(8*4 + i) * (tmp10 - tmp11));
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
161
6408
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
162 tmp12 += tmp13;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
163 tmp12 *= A1;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
164 data[8*2 + i]= lrintf(SCALE(8*2 + i) * (tmp13 + tmp12));
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
165 data[8*6 + i]= lrintf(SCALE(8*6 + i) * (tmp13 - tmp12));
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
166
6409
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
167 tmp4 += tmp5;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
168 tmp5 += tmp6;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
169 tmp6 += tmp7;
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
170
6411
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
171 #if 0
6409
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
172 z5= (tmp4 - tmp6) * A5;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
173 z2= tmp4*A2 + z5;
813ff53ed933 Avoid a=b+c but rather use a+=b (gcc is too stupid to do this itself),
michael
parents: 6408
diff changeset
174 z4= tmp6*A4 + z5;
6411
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
175 #else
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
176 z2= tmp4*(A2+A5) - tmp6*A5;
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
177 z4= tmp6*(A4-A5) + tmp4*A5;
5d50e8a7eb09 Alternative and faster (on gcc x86) plane rotation.
michael
parents: 6410
diff changeset
178 #endif
6410
72e573046758 Get rid of another useless variable.
michael
parents: 6409
diff changeset
179 tmp5*=A1;
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
180
6410
72e573046758 Get rid of another useless variable.
michael
parents: 6409
diff changeset
181 z11= tmp7 + tmp5;
72e573046758 Get rid of another useless variable.
michael
parents: 6409
diff changeset
182 z13= tmp7 - tmp5;
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
183
1563
820e06c6ca9b 4.9-RC FreeBSD doesnt like lrint() so change to lrintf(), not that bsd supports that but we emulate it ...
michael
parents: 1562
diff changeset
184 data[8*5 + i]= lrintf(SCALE(8*5 + i) * (z13 + z2));
820e06c6ca9b 4.9-RC FreeBSD doesnt like lrint() so change to lrintf(), not that bsd supports that but we emulate it ...
michael
parents: 1562
diff changeset
185 data[8*3 + i]= lrintf(SCALE(8*3 + i) * (z13 - z2));
820e06c6ca9b 4.9-RC FreeBSD doesnt like lrint() so change to lrintf(), not that bsd supports that but we emulate it ...
michael
parents: 1562
diff changeset
186 data[8*1 + i]= lrintf(SCALE(8*1 + i) * (z11 + z4));
820e06c6ca9b 4.9-RC FreeBSD doesnt like lrint() so change to lrintf(), not that bsd supports that but we emulate it ...
michael
parents: 1562
diff changeset
187 data[8*7 + i]= lrintf(SCALE(8*7 + i) * (z11 - z4));
1557
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
188 }
5d53c03186a1 floating point AAN DCT
michael
parents:
diff changeset
189 }
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
190
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
191 void ff_faandct248(DCTELEM * data)
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
192 {
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
193 FLOAT tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
194 FLOAT tmp10, tmp11, tmp12, tmp13;
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
195 FLOAT temp[64];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
196 int i;
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
197
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
198 emms_c();
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
199
1589
eb26d190cf5a move identical code into its own function
michael
parents: 1571
diff changeset
200 row_fdct(temp, data);
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
201
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
202 for (i=0; i<8; i++) {
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
203 tmp0 = temp[8*0 + i] + temp[8*1 + i];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
204 tmp1 = temp[8*2 + i] + temp[8*3 + i];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
205 tmp2 = temp[8*4 + i] + temp[8*5 + i];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
206 tmp3 = temp[8*6 + i] + temp[8*7 + i];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
207 tmp4 = temp[8*0 + i] - temp[8*1 + i];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
208 tmp5 = temp[8*2 + i] - temp[8*3 + i];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
209 tmp6 = temp[8*4 + i] - temp[8*5 + i];
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
210 tmp7 = temp[8*6 + i] - temp[8*7 + i];
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
211
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
212 tmp10 = tmp0 + tmp3;
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
213 tmp11 = tmp1 + tmp2;
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
214 tmp12 = tmp1 - tmp2;
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
215 tmp13 = tmp0 - tmp3;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
216
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
217 data[8*0 + i] = lrintf(SCALE(8*0 + i) * (tmp10 + tmp11));
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
218 data[8*4 + i] = lrintf(SCALE(8*4 + i) * (tmp10 - tmp11));
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
219
6408
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
220 tmp12 += tmp13;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
221 tmp12 *= A1;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
222 data[8*2 + i] = lrintf(SCALE(8*2 + i) * (tmp13 + tmp12));
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
223 data[8*6 + i] = lrintf(SCALE(8*6 + i) * (tmp13 - tmp12));
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2024
diff changeset
224
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
225 tmp10 = tmp4 + tmp7;
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
226 tmp11 = tmp5 + tmp6;
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
227 tmp12 = tmp5 - tmp6;
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
228 tmp13 = tmp4 - tmp7;
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
229
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
230 data[8*1 + i] = lrintf(SCALE(8*0 + i) * (tmp10 + tmp11));
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
231 data[8*5 + i] = lrintf(SCALE(8*4 + i) * (tmp10 - tmp11));
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
232
6408
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
233 tmp12 += tmp13;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
234 tmp12 *= A1;
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
235 data[8*3 + i] = lrintf(SCALE(8*2 + i) * (tmp13 + tmp12));
7af599600f2d One variable less. No real change to the actual code.
michael
parents: 6398
diff changeset
236 data[8*7 + i] = lrintf(SCALE(8*6 + i) * (tmp13 - tmp12));
1571
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
237 }
aa4dc16c0f18 * adding integer/floating point AAN implementations for DCT 2-4-8
romansh
parents: 1563
diff changeset
238 }