← Back to branch summary

~ubuntu-branches/ubuntu/hardy/avidemux/hardy

« back to all changes in this revision

Viewing changes to avidemux/ADM_lavcodec/i386/dsputil_mmx_avg.h

  • Committer: Bazaar Package Importer
  • Author(s): Daniel T Chen
  • Date: 2006-12-15 17:13:20 UTC
  • mfrom: (1.1.6 upstream)
  • Revision ID: james.westby@ubuntu.com-20061215171320-w79pvpehxx2fr217
Tags: 1:2.3.0-0.0ubuntu1
* Merge from debian-multimedia.org, remaining Ubuntu change:
  - desktop file,
  - no support for ccache and make -j.
* Closes Ubuntu: #69614.

Show diffs side-by-side

added added

removed removed

Lines of Context:
avidemux/ADM_lavcodec/i386/dsputil_mmx_avg.h
 
1
/*
 
2
 * DSP utils : average functions are compiled twice for 3dnow/mmx2
 
3
 * Copyright (c) 2000, 2001 Fabrice Bellard.
 
4
 * Copyright (c) 2002-2004 Michael Niedermayer
 
5
 *
 
6
 * This library is free software; you can redistribute it and/or
 
7
 * modify it under the terms of the GNU Lesser General Public
 
8
 * License as published by the Free Software Foundation; either
 
9
 * version 2 of the License, or (at your option) any later version.
 
10
 *
 
11
 * This library is distributed in the hope that it will be useful,
 
12
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 
14
 * Lesser General Public License for more details.
 
15
 *
 
16
 * You should have received a copy of the GNU Lesser General Public
 
17
 * License along with this library; if not, write to the Free Software
 
18
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 
19
 *
 
20
 * MMX optimization by Nick Kurshev <nickols_k@mail.ru>
 
21
 * mostly rewritten by Michael Niedermayer <michaelni@gmx.at>
 
22
 * and improved by Zdenek Kabelac <kabi@users.sf.net>
 
23
 */
 
24
 
 
25
/* XXX: we use explicit registers to avoid a gcc 2.95.2 register asm
 
26
   clobber bug - now it will work with 2.95.2 and also with -fPIC
 
27
 */
 
28
static void DEF(put_pixels8_x2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
29
{
 
30
    __asm __volatile(
 
31
        "lea (%3, %3), %%"REG_a"        \n\t"
 
32
        "1:                             \n\t"
 
33
        "movq (%1), %%mm0               \n\t"
 
34
        "movq (%1, %3), %%mm1           \n\t"
 
35
        PAVGB" 1(%1), %%mm0             \n\t"
 
36
        PAVGB" 1(%1, %3), %%mm1         \n\t"
 
37
        "movq %%mm0, (%2)               \n\t"
 
38
        "movq %%mm1, (%2, %3)           \n\t"
 
39
        "add %%"REG_a", %1              \n\t"
 
40
        "add %%"REG_a", %2              \n\t"
 
41
        "movq (%1), %%mm0               \n\t"
 
42
        "movq (%1, %3), %%mm1           \n\t"
 
43
        PAVGB" 1(%1), %%mm0             \n\t"
 
44
        PAVGB" 1(%1, %3), %%mm1         \n\t"
 
45
        "add %%"REG_a", %1              \n\t"
 
46
        "movq %%mm0, (%2)               \n\t"
 
47
        "movq %%mm1, (%2, %3)           \n\t"
 
48
        "add %%"REG_a", %2              \n\t"
 
49
        "subl $4, %0                    \n\t"
 
50
        "jnz 1b                         \n\t"
 
51
        :"+g"(h), "+S"(pixels), "+D"(block)
 
52
        :"r" ((long)line_size)
 
53
        :"%"REG_a, "memory");
 
54
}
 
55
 
 
56
static void DEF(put_pixels4_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
57
{
 
58
    __asm __volatile(
 
59
        "testl $1, %0                   \n\t"
 
60
            " jz 1f                     \n\t"
 
61
        "movd   (%1), %%mm0             \n\t"
 
62
        "movd   (%2), %%mm1             \n\t"
 
63
        "add    %4, %1                  \n\t"
 
64
        "add    $4, %2                  \n\t"
 
65
        PAVGB" %%mm1, %%mm0             \n\t"
 
66
        "movd   %%mm0, (%3)             \n\t"
 
67
        "add    %5, %3                  \n\t"
 
68
        "decl   %0                      \n\t"
 
69
        "1:                             \n\t"
 
70
        "movd   (%1), %%mm0             \n\t"
 
71
        "add    %4, %1                  \n\t"
 
72
        "movd   (%1), %%mm1             \n\t"
 
73
        "movd   (%2), %%mm2             \n\t"
 
74
        "movd   4(%2), %%mm3            \n\t"
 
75
        "add    %4, %1                  \n\t"
 
76
        PAVGB" %%mm2, %%mm0             \n\t"
 
77
        PAVGB" %%mm3, %%mm1             \n\t"
 
78
        "movd   %%mm0, (%3)             \n\t"
 
79
        "add    %5, %3                  \n\t"
 
80
        "movd   %%mm1, (%3)             \n\t"
 
81
        "add    %5, %3                  \n\t"
 
82
        "movd   (%1), %%mm0             \n\t"
 
83
        "add    %4, %1                  \n\t"
 
84
        "movd   (%1), %%mm1             \n\t"
 
85
        "movd   8(%2), %%mm2            \n\t"
 
86
        "movd   12(%2), %%mm3           \n\t"
 
87
        "add    %4, %1                  \n\t"
 
88
        PAVGB" %%mm2, %%mm0             \n\t"
 
89
        PAVGB" %%mm3, %%mm1             \n\t"
 
90
        "movd   %%mm0, (%3)             \n\t"
 
91
        "add    %5, %3                  \n\t"
 
92
        "movd   %%mm1, (%3)             \n\t"
 
93
        "add    %5, %3                  \n\t"
 
94
        "add    $16, %2                 \n\t"
 
95
        "subl   $4, %0                  \n\t"
 
96
        "jnz    1b                      \n\t"
 
97
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
98
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
99
#else
 
100
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
101
#endif
 
102
        :"S"((long)src1Stride), "D"((long)dstStride)
 
103
        :"memory");
 
104
}
 
105
 
 
106
 
 
107
static void DEF(put_pixels8_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
108
{
 
109
    __asm __volatile(
 
110
        "testl $1, %0                   \n\t"
 
111
            " jz 1f                     \n\t"
 
112
        "movq   (%1), %%mm0             \n\t"
 
113
        "movq   (%2), %%mm1             \n\t"
 
114
        "add    %4, %1                  \n\t"
 
115
        "add    $8, %2                  \n\t"
 
116
        PAVGB" %%mm1, %%mm0             \n\t"
 
117
        "movq   %%mm0, (%3)             \n\t"
 
118
        "add    %5, %3                  \n\t"
 
119
        "decl   %0                      \n\t"
 
120
        "1:                             \n\t"
 
121
        "movq   (%1), %%mm0             \n\t"
 
122
        "add    %4, %1                  \n\t"
 
123
        "movq   (%1), %%mm1             \n\t"
 
124
        "add    %4, %1                  \n\t"
 
125
        PAVGB" (%2), %%mm0              \n\t"
 
126
        PAVGB" 8(%2), %%mm1             \n\t"
 
127
        "movq   %%mm0, (%3)             \n\t"
 
128
        "add    %5, %3                  \n\t"
 
129
        "movq   %%mm1, (%3)             \n\t"
 
130
        "add    %5, %3                  \n\t"
 
131
        "movq   (%1), %%mm0             \n\t"
 
132
        "add    %4, %1                  \n\t"
 
133
        "movq   (%1), %%mm1             \n\t"
 
134
        "add    %4, %1                  \n\t"
 
135
        PAVGB" 16(%2), %%mm0            \n\t"
 
136
        PAVGB" 24(%2), %%mm1            \n\t"
 
137
        "movq   %%mm0, (%3)             \n\t"
 
138
        "add    %5, %3                  \n\t"
 
139
        "movq   %%mm1, (%3)             \n\t"
 
140
        "add    %5, %3                  \n\t"
 
141
        "add    $32, %2                 \n\t"
 
142
        "subl   $4, %0                  \n\t"
 
143
        "jnz    1b                      \n\t"
 
144
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
145
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
146
#else
 
147
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
148
#endif
 
149
        :"S"((long)src1Stride), "D"((long)dstStride)
 
150
        :"memory");
 
151
//the following should be used, though better not with gcc ...
 
152
/*        :"+g"(h), "+r"(src1), "+r"(src2), "+r"(dst)
 
153
        :"r"(src1Stride), "r"(dstStride)
 
154
        :"memory");*/
 
155
}
 
156
 
 
157
static void DEF(put_no_rnd_pixels8_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
158
{
 
159
    __asm __volatile(
 
160
        "pcmpeqb %%mm6, %%mm6           \n\t"
 
161
        "testl $1, %0                   \n\t"
 
162
            " jz 1f                     \n\t"
 
163
        "movq   (%1), %%mm0             \n\t"
 
164
        "movq   (%2), %%mm1             \n\t"
 
165
        "add    %4, %1                  \n\t"
 
166
        "add    $8, %2                  \n\t"
 
167
        "pxor %%mm6, %%mm0              \n\t"
 
168
        "pxor %%mm6, %%mm1              \n\t"
 
169
        PAVGB" %%mm1, %%mm0             \n\t"
 
170
        "pxor %%mm6, %%mm0              \n\t"
 
171
        "movq   %%mm0, (%3)             \n\t"
 
172
        "add    %5, %3                  \n\t"
 
173
        "decl   %0                      \n\t"
 
174
        "1:                             \n\t"
 
175
        "movq   (%1), %%mm0             \n\t"
 
176
        "add    %4, %1                  \n\t"
 
177
        "movq   (%1), %%mm1             \n\t"
 
178
        "add    %4, %1                  \n\t"
 
179
        "movq   (%2), %%mm2             \n\t"
 
180
        "movq   8(%2), %%mm3            \n\t"
 
181
        "pxor %%mm6, %%mm0              \n\t"
 
182
        "pxor %%mm6, %%mm1              \n\t"
 
183
        "pxor %%mm6, %%mm2              \n\t"
 
184
        "pxor %%mm6, %%mm3              \n\t"
 
185
        PAVGB" %%mm2, %%mm0             \n\t"
 
186
        PAVGB" %%mm3, %%mm1             \n\t"
 
187
        "pxor %%mm6, %%mm0              \n\t"
 
188
        "pxor %%mm6, %%mm1              \n\t"
 
189
        "movq   %%mm0, (%3)             \n\t"
 
190
        "add    %5, %3                  \n\t"
 
191
        "movq   %%mm1, (%3)             \n\t"
 
192
        "add    %5, %3                  \n\t"
 
193
        "movq   (%1), %%mm0             \n\t"
 
194
        "add    %4, %1                  \n\t"
 
195
        "movq   (%1), %%mm1             \n\t"
 
196
        "add    %4, %1                  \n\t"
 
197
        "movq   16(%2), %%mm2           \n\t"
 
198
        "movq   24(%2), %%mm3           \n\t"
 
199
        "pxor %%mm6, %%mm0              \n\t"
 
200
        "pxor %%mm6, %%mm1              \n\t"
 
201
        "pxor %%mm6, %%mm2              \n\t"
 
202
        "pxor %%mm6, %%mm3              \n\t"
 
203
        PAVGB" %%mm2, %%mm0             \n\t"
 
204
        PAVGB" %%mm3, %%mm1             \n\t"
 
205
        "pxor %%mm6, %%mm0              \n\t"
 
206
        "pxor %%mm6, %%mm1              \n\t"
 
207
        "movq   %%mm0, (%3)             \n\t"
 
208
        "add    %5, %3                  \n\t"
 
209
        "movq   %%mm1, (%3)             \n\t"
 
210
        "add    %5, %3                  \n\t"
 
211
        "add    $32, %2                 \n\t"
 
212
        "subl   $4, %0                  \n\t"
 
213
        "jnz    1b                      \n\t"
 
214
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
215
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
216
#else
 
217
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
218
#endif
 
219
        :"S"((long)src1Stride), "D"((long)dstStride)
 
220
        :"memory");
 
221
//the following should be used, though better not with gcc ...
 
222
/*        :"+g"(h), "+r"(src1), "+r"(src2), "+r"(dst)
 
223
        :"r"(src1Stride), "r"(dstStride)
 
224
        :"memory");*/
 
225
}
 
226
 
 
227
static void DEF(avg_pixels4_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
228
{
 
229
    __asm __volatile(
 
230
        "testl $1, %0                   \n\t"
 
231
            " jz 1f                     \n\t"
 
232
        "movd   (%1), %%mm0             \n\t"
 
233
        "movd   (%2), %%mm1             \n\t"
 
234
        "add    %4, %1                  \n\t"
 
235
        "add    $4, %2                  \n\t"
 
236
        PAVGB" %%mm1, %%mm0             \n\t"
 
237
        PAVGB" (%3), %%mm0              \n\t"
 
238
        "movd   %%mm0, (%3)             \n\t"
 
239
        "add    %5, %3                  \n\t"
 
240
        "decl   %0                      \n\t"
 
241
        "1:                             \n\t"
 
242
        "movd   (%1), %%mm0             \n\t"
 
243
        "add    %4, %1                  \n\t"
 
244
        "movd   (%1), %%mm1             \n\t"
 
245
        "add    %4, %1                  \n\t"
 
246
        PAVGB" (%2), %%mm0              \n\t"
 
247
        PAVGB" 4(%2), %%mm1             \n\t"
 
248
        PAVGB" (%3), %%mm0              \n\t"
 
249
        "movd   %%mm0, (%3)             \n\t"
 
250
        "add    %5, %3                  \n\t"
 
251
        PAVGB" (%3), %%mm1              \n\t"
 
252
        "movd   %%mm1, (%3)             \n\t"
 
253
        "add    %5, %3                  \n\t"
 
254
        "movd   (%1), %%mm0             \n\t"
 
255
        "add    %4, %1                  \n\t"
 
256
        "movd   (%1), %%mm1             \n\t"
 
257
        "add    %4, %1                  \n\t"
 
258
        PAVGB" 8(%2), %%mm0             \n\t"
 
259
        PAVGB" 12(%2), %%mm1            \n\t"
 
260
        PAVGB" (%3), %%mm0              \n\t"
 
261
        "movd   %%mm0, (%3)             \n\t"
 
262
        "add    %5, %3                  \n\t"
 
263
        PAVGB" (%3), %%mm1              \n\t"
 
264
        "movd   %%mm1, (%3)             \n\t"
 
265
        "add    %5, %3                  \n\t"
 
266
        "add    $16, %2                 \n\t"
 
267
        "subl   $4, %0                  \n\t"
 
268
        "jnz    1b                      \n\t"
 
269
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
270
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
271
#else
 
272
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
273
#endif
 
274
        :"S"((long)src1Stride), "D"((long)dstStride)
 
275
        :"memory");
 
276
}
 
277
 
 
278
 
 
279
static void DEF(avg_pixels8_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
280
{
 
281
    __asm __volatile(
 
282
        "testl $1, %0                   \n\t"
 
283
            " jz 1f                     \n\t"
 
284
        "movq   (%1), %%mm0             \n\t"
 
285
        "movq   (%2), %%mm1             \n\t"
 
286
        "add    %4, %1                  \n\t"
 
287
        "add    $8, %2                  \n\t"
 
288
        PAVGB" %%mm1, %%mm0             \n\t"
 
289
        PAVGB" (%3), %%mm0              \n\t"
 
290
        "movq   %%mm0, (%3)             \n\t"
 
291
        "add    %5, %3                  \n\t"
 
292
        "decl   %0                      \n\t"
 
293
        "1:                             \n\t"
 
294
        "movq   (%1), %%mm0             \n\t"
 
295
        "add    %4, %1                  \n\t"
 
296
        "movq   (%1), %%mm1             \n\t"
 
297
        "add    %4, %1                  \n\t"
 
298
        PAVGB" (%2), %%mm0              \n\t"
 
299
        PAVGB" 8(%2), %%mm1             \n\t"
 
300
        PAVGB" (%3), %%mm0              \n\t"
 
301
        "movq   %%mm0, (%3)             \n\t"
 
302
        "add    %5, %3                  \n\t"
 
303
        PAVGB" (%3), %%mm1              \n\t"
 
304
        "movq   %%mm1, (%3)             \n\t"
 
305
        "add    %5, %3                  \n\t"
 
306
        "movq   (%1), %%mm0             \n\t"
 
307
        "add    %4, %1                  \n\t"
 
308
        "movq   (%1), %%mm1             \n\t"
 
309
        "add    %4, %1                  \n\t"
 
310
        PAVGB" 16(%2), %%mm0            \n\t"
 
311
        PAVGB" 24(%2), %%mm1            \n\t"
 
312
        PAVGB" (%3), %%mm0              \n\t"
 
313
        "movq   %%mm0, (%3)             \n\t"
 
314
        "add    %5, %3                  \n\t"
 
315
        PAVGB" (%3), %%mm1              \n\t"
 
316
        "movq   %%mm1, (%3)             \n\t"
 
317
        "add    %5, %3                  \n\t"
 
318
        "add    $32, %2                 \n\t"
 
319
        "subl   $4, %0                  \n\t"
 
320
        "jnz    1b                      \n\t"
 
321
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
322
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
323
#else
 
324
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
325
#endif
 
326
        :"S"((long)src1Stride), "D"((long)dstStride)
 
327
        :"memory");
 
328
//the following should be used, though better not with gcc ...
 
329
/*        :"+g"(h), "+r"(src1), "+r"(src2), "+r"(dst)
 
330
        :"r"(src1Stride), "r"(dstStride)
 
331
        :"memory");*/
 
332
}
 
333
 
 
334
static void DEF(put_pixels16_x2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
335
{
 
336
    __asm __volatile(
 
337
        "lea (%3, %3), %%"REG_a"        \n\t"
 
338
        "1:                             \n\t"
 
339
        "movq (%1), %%mm0               \n\t"
 
340
        "movq (%1, %3), %%mm1           \n\t"
 
341
        "movq 8(%1), %%mm2              \n\t"
 
342
        "movq 8(%1, %3), %%mm3          \n\t"
 
343
        PAVGB" 1(%1), %%mm0             \n\t"
 
344
        PAVGB" 1(%1, %3), %%mm1         \n\t"
 
345
        PAVGB" 9(%1), %%mm2             \n\t"
 
346
        PAVGB" 9(%1, %3), %%mm3         \n\t"
 
347
        "movq %%mm0, (%2)               \n\t"
 
348
        "movq %%mm1, (%2, %3)           \n\t"
 
349
        "movq %%mm2, 8(%2)              \n\t"
 
350
        "movq %%mm3, 8(%2, %3)          \n\t"
 
351
        "add %%"REG_a", %1              \n\t"
 
352
        "add %%"REG_a", %2              \n\t"
 
353
        "movq (%1), %%mm0               \n\t"
 
354
        "movq (%1, %3), %%mm1           \n\t"
 
355
        "movq 8(%1), %%mm2              \n\t"
 
356
        "movq 8(%1, %3), %%mm3          \n\t"
 
357
        PAVGB" 1(%1), %%mm0             \n\t"
 
358
        PAVGB" 1(%1, %3), %%mm1         \n\t"
 
359
        PAVGB" 9(%1), %%mm2             \n\t"
 
360
        PAVGB" 9(%1, %3), %%mm3         \n\t"
 
361
        "add %%"REG_a", %1              \n\t"
 
362
        "movq %%mm0, (%2)               \n\t"
 
363
        "movq %%mm1, (%2, %3)           \n\t"
 
364
        "movq %%mm2, 8(%2)              \n\t"
 
365
        "movq %%mm3, 8(%2, %3)          \n\t"
 
366
        "add %%"REG_a", %2              \n\t"
 
367
        "subl $4, %0                    \n\t"
 
368
        "jnz 1b                         \n\t"
 
369
        :"+g"(h), "+S"(pixels), "+D"(block)
 
370
        :"r" ((long)line_size)
 
371
        :"%"REG_a, "memory");
 
372
}
 
373
 
 
374
static void DEF(put_pixels16_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
375
{
 
376
    __asm __volatile(
 
377
        "testl $1, %0                   \n\t"
 
378
            " jz 1f                     \n\t"
 
379
        "movq   (%1), %%mm0             \n\t"
 
380
        "movq   8(%1), %%mm1            \n\t"
 
381
        PAVGB" (%2), %%mm0              \n\t"
 
382
        PAVGB" 8(%2), %%mm1             \n\t"
 
383
        "add    %4, %1                  \n\t"
 
384
        "add    $16, %2                 \n\t"
 
385
        "movq   %%mm0, (%3)             \n\t"
 
386
        "movq   %%mm1, 8(%3)            \n\t"
 
387
        "add    %5, %3                  \n\t"
 
388
        "decl   %0                      \n\t"
 
389
        "1:                             \n\t"
 
390
        "movq   (%1), %%mm0             \n\t"
 
391
        "movq   8(%1), %%mm1            \n\t"
 
392
        "add    %4, %1                  \n\t"
 
393
        PAVGB" (%2), %%mm0              \n\t"
 
394
        PAVGB" 8(%2), %%mm1             \n\t"
 
395
        "movq   %%mm0, (%3)             \n\t"
 
396
        "movq   %%mm1, 8(%3)            \n\t"
 
397
        "add    %5, %3                  \n\t"
 
398
        "movq   (%1), %%mm0             \n\t"
 
399
        "movq   8(%1), %%mm1            \n\t"
 
400
        "add    %4, %1                  \n\t"
 
401
        PAVGB" 16(%2), %%mm0            \n\t"
 
402
        PAVGB" 24(%2), %%mm1            \n\t"
 
403
        "movq   %%mm0, (%3)             \n\t"
 
404
        "movq   %%mm1, 8(%3)            \n\t"
 
405
        "add    %5, %3                  \n\t"
 
406
        "add    $32, %2                 \n\t"
 
407
        "subl   $2, %0                  \n\t"
 
408
        "jnz    1b                      \n\t"
 
409
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
410
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
411
#else
 
412
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
413
#endif
 
414
        :"S"((long)src1Stride), "D"((long)dstStride)
 
415
        :"memory");
 
416
//the following should be used, though better not with gcc ...
 
417
/*        :"+g"(h), "+r"(src1), "+r"(src2), "+r"(dst)
 
418
        :"r"(src1Stride), "r"(dstStride)
 
419
        :"memory");*/
 
420
}
 
421
 
 
422
static void DEF(avg_pixels16_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
423
{
 
424
    __asm __volatile(
 
425
        "testl $1, %0                   \n\t"
 
426
            " jz 1f                     \n\t"
 
427
        "movq   (%1), %%mm0             \n\t"
 
428
        "movq   8(%1), %%mm1            \n\t"
 
429
        PAVGB" (%2), %%mm0              \n\t"
 
430
        PAVGB" 8(%2), %%mm1             \n\t"
 
431
        "add    %4, %1                  \n\t"
 
432
        "add    $16, %2                 \n\t"
 
433
        PAVGB" (%3), %%mm0              \n\t"
 
434
        PAVGB" 8(%3), %%mm1             \n\t"
 
435
        "movq   %%mm0, (%3)             \n\t"
 
436
        "movq   %%mm1, 8(%3)            \n\t"
 
437
        "add    %5, %3                  \n\t"
 
438
        "decl   %0                      \n\t"
 
439
        "1:                             \n\t"
 
440
        "movq   (%1), %%mm0             \n\t"
 
441
        "movq   8(%1), %%mm1            \n\t"
 
442
        "add    %4, %1                  \n\t"
 
443
        PAVGB" (%2), %%mm0              \n\t"
 
444
        PAVGB" 8(%2), %%mm1             \n\t"
 
445
        PAVGB" (%3), %%mm0              \n\t"
 
446
        PAVGB" 8(%3), %%mm1             \n\t"
 
447
        "movq   %%mm0, (%3)             \n\t"
 
448
        "movq   %%mm1, 8(%3)            \n\t"
 
449
        "add    %5, %3                  \n\t"
 
450
        "movq   (%1), %%mm0             \n\t"
 
451
        "movq   8(%1), %%mm1            \n\t"
 
452
        "add    %4, %1                  \n\t"
 
453
        PAVGB" 16(%2), %%mm0            \n\t"
 
454
        PAVGB" 24(%2), %%mm1            \n\t"
 
455
        PAVGB" (%3), %%mm0              \n\t"
 
456
        PAVGB" 8(%3), %%mm1             \n\t"
 
457
        "movq   %%mm0, (%3)             \n\t"
 
458
        "movq   %%mm1, 8(%3)            \n\t"
 
459
        "add    %5, %3                  \n\t"
 
460
        "add    $32, %2                 \n\t"
 
461
        "subl   $2, %0                  \n\t"
 
462
        "jnz    1b                      \n\t"
 
463
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
464
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
465
#else
 
466
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
467
#endif
 
468
        :"S"((long)src1Stride), "D"((long)dstStride)
 
469
        :"memory");
 
470
//the following should be used, though better not with gcc ...
 
471
/*        :"+g"(h), "+r"(src1), "+r"(src2), "+r"(dst)
 
472
        :"r"(src1Stride), "r"(dstStride)
 
473
        :"memory");*/
 
474
}
 
475
 
 
476
static void DEF(put_no_rnd_pixels16_l2)(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int src1Stride, int h)
 
477
{
 
478
    __asm __volatile(
 
479
        "pcmpeqb %%mm6, %%mm6           \n\t"
 
480
        "testl $1, %0                   \n\t"
 
481
            " jz 1f                     \n\t"
 
482
        "movq   (%1), %%mm0             \n\t"
 
483
        "movq   8(%1), %%mm1            \n\t"
 
484
        "movq   (%2), %%mm2             \n\t"
 
485
        "movq   8(%2), %%mm3            \n\t"
 
486
        "pxor %%mm6, %%mm0              \n\t"
 
487
        "pxor %%mm6, %%mm1              \n\t"
 
488
        "pxor %%mm6, %%mm2              \n\t"
 
489
        "pxor %%mm6, %%mm3              \n\t"
 
490
        PAVGB" %%mm2, %%mm0             \n\t"
 
491
        PAVGB" %%mm3, %%mm1             \n\t"
 
492
        "pxor %%mm6, %%mm0              \n\t"
 
493
        "pxor %%mm6, %%mm1              \n\t"
 
494
        "add    %4, %1                  \n\t"
 
495
        "add    $16, %2                 \n\t"
 
496
        "movq   %%mm0, (%3)             \n\t"
 
497
        "movq   %%mm1, 8(%3)            \n\t"
 
498
        "add    %5, %3                  \n\t"
 
499
        "decl   %0                      \n\t"
 
500
        "1:                             \n\t"
 
501
        "movq   (%1), %%mm0             \n\t"
 
502
        "movq   8(%1), %%mm1            \n\t"
 
503
        "add    %4, %1                  \n\t"
 
504
        "movq   (%2), %%mm2             \n\t"
 
505
        "movq   8(%2), %%mm3            \n\t"
 
506
        "pxor %%mm6, %%mm0              \n\t"
 
507
        "pxor %%mm6, %%mm1              \n\t"
 
508
        "pxor %%mm6, %%mm2              \n\t"
 
509
        "pxor %%mm6, %%mm3              \n\t"
 
510
        PAVGB" %%mm2, %%mm0             \n\t"
 
511
        PAVGB" %%mm3, %%mm1             \n\t"
 
512
        "pxor %%mm6, %%mm0              \n\t"
 
513
        "pxor %%mm6, %%mm1              \n\t"
 
514
        "movq   %%mm0, (%3)             \n\t"
 
515
        "movq   %%mm1, 8(%3)            \n\t"
 
516
        "add    %5, %3                  \n\t"
 
517
        "movq   (%1), %%mm0             \n\t"
 
518
        "movq   8(%1), %%mm1            \n\t"
 
519
        "add    %4, %1                  \n\t"
 
520
        "movq   16(%2), %%mm2           \n\t"
 
521
        "movq   24(%2), %%mm3           \n\t"
 
522
        "pxor %%mm6, %%mm0              \n\t"
 
523
        "pxor %%mm6, %%mm1              \n\t"
 
524
        "pxor %%mm6, %%mm2              \n\t"
 
525
        "pxor %%mm6, %%mm3              \n\t"
 
526
        PAVGB" %%mm2, %%mm0             \n\t"
 
527
        PAVGB" %%mm3, %%mm1             \n\t"
 
528
        "pxor %%mm6, %%mm0              \n\t"
 
529
        "pxor %%mm6, %%mm1              \n\t"
 
530
        "movq   %%mm0, (%3)             \n\t"
 
531
        "movq   %%mm1, 8(%3)            \n\t"
 
532
        "add    %5, %3                  \n\t"
 
533
        "add    $32, %2                 \n\t"
 
534
        "subl   $2, %0                  \n\t"
 
535
        "jnz    1b                      \n\t"
 
536
#ifdef PIC //Note "+bm" and "+mb" are buggy too (with gcc 3.2.2 at least) and cant be used
 
537
        :"+m"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
538
#else
 
539
        :"+b"(h), "+a"(src1), "+c"(src2), "+d"(dst)
 
540
#endif
 
541
        :"S"((long)src1Stride), "D"((long)dstStride)
 
542
        :"memory");
 
543
//the following should be used, though better not with gcc ...
 
544
/*        :"+g"(h), "+r"(src1), "+r"(src2), "+r"(dst)
 
545
        :"r"(src1Stride), "r"(dstStride)
 
546
        :"memory");*/
 
547
}
 
548
 
 
549
/* GL: this function does incorrect rounding if overflow */
 
550
static void DEF(put_no_rnd_pixels8_x2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
551
{
 
552
    MOVQ_BONE(mm6);
 
553
    __asm __volatile(
 
554
        "lea (%3, %3), %%"REG_a"        \n\t"
 
555
        "1:                             \n\t"
 
556
        "movq (%1), %%mm0               \n\t"
 
557
        "movq (%1, %3), %%mm2           \n\t"
 
558
        "movq 1(%1), %%mm1              \n\t"
 
559
        "movq 1(%1, %3), %%mm3          \n\t"
 
560
        "add %%"REG_a", %1              \n\t"
 
561
        "psubusb %%mm6, %%mm0           \n\t"
 
562
        "psubusb %%mm6, %%mm2           \n\t"
 
563
        PAVGB" %%mm1, %%mm0             \n\t"
 
564
        PAVGB" %%mm3, %%mm2             \n\t"
 
565
        "movq %%mm0, (%2)               \n\t"
 
566
        "movq %%mm2, (%2, %3)           \n\t"
 
567
        "movq (%1), %%mm0               \n\t"
 
568
        "movq 1(%1), %%mm1              \n\t"
 
569
        "movq (%1, %3), %%mm2           \n\t"
 
570
        "movq 1(%1, %3), %%mm3          \n\t"
 
571
        "add %%"REG_a", %2              \n\t"
 
572
        "add %%"REG_a", %1              \n\t"
 
573
        "psubusb %%mm6, %%mm0           \n\t"
 
574
        "psubusb %%mm6, %%mm2           \n\t"
 
575
        PAVGB" %%mm1, %%mm0             \n\t"
 
576
        PAVGB" %%mm3, %%mm2             \n\t"
 
577
        "movq %%mm0, (%2)               \n\t"
 
578
        "movq %%mm2, (%2, %3)           \n\t"
 
579
        "add %%"REG_a", %2              \n\t"
 
580
        "subl $4, %0                    \n\t"
 
581
        "jnz 1b                         \n\t"
 
582
        :"+g"(h), "+S"(pixels), "+D"(block)
 
583
        :"r" ((long)line_size)
 
584
        :"%"REG_a, "memory");
 
585
}
 
586
 
 
587
static void DEF(put_pixels8_y2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
588
{
 
589
    __asm __volatile(
 
590
        "lea (%3, %3), %%"REG_a"        \n\t"
 
591
        "movq (%1), %%mm0               \n\t"
 
592
        "sub %3, %2                     \n\t"
 
593
        "1:                             \n\t"
 
594
        "movq (%1, %3), %%mm1           \n\t"
 
595
        "movq (%1, %%"REG_a"), %%mm2    \n\t"
 
596
        "add %%"REG_a", %1              \n\t"
 
597
        PAVGB" %%mm1, %%mm0             \n\t"
 
598
        PAVGB" %%mm2, %%mm1             \n\t"
 
599
        "movq %%mm0, (%2, %3)           \n\t"
 
600
        "movq %%mm1, (%2, %%"REG_a")    \n\t"
 
601
        "movq (%1, %3), %%mm1           \n\t"
 
602
        "movq (%1, %%"REG_a"), %%mm0    \n\t"
 
603
        "add %%"REG_a", %2              \n\t"
 
604
        "add %%"REG_a", %1              \n\t"
 
605
        PAVGB" %%mm1, %%mm2             \n\t"
 
606
        PAVGB" %%mm0, %%mm1             \n\t"
 
607
        "movq %%mm2, (%2, %3)           \n\t"
 
608
        "movq %%mm1, (%2, %%"REG_a")    \n\t"
 
609
        "add %%"REG_a", %2              \n\t"
 
610
        "subl $4, %0                    \n\t"
 
611
        "jnz 1b                         \n\t"
 
612
        :"+g"(h), "+S"(pixels), "+D" (block)
 
613
        :"r" ((long)line_size)
 
614
        :"%"REG_a, "memory");
 
615
}
 
616
 
 
617
/* GL: this function does incorrect rounding if overflow */
 
618
static void DEF(put_no_rnd_pixels8_y2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
619
{
 
620
    MOVQ_BONE(mm6);
 
621
    __asm __volatile(
 
622
        "lea (%3, %3), %%"REG_a"        \n\t"
 
623
        "movq (%1), %%mm0               \n\t"
 
624
        "sub %3, %2                     \n\t"
 
625
        "1:                             \n\t"
 
626
        "movq (%1, %3), %%mm1           \n\t"
 
627
        "movq (%1, %%"REG_a"), %%mm2    \n\t"
 
628
        "add %%"REG_a", %1              \n\t"
 
629
        "psubusb %%mm6, %%mm1           \n\t"
 
630
        PAVGB" %%mm1, %%mm0             \n\t"
 
631
        PAVGB" %%mm2, %%mm1             \n\t"
 
632
        "movq %%mm0, (%2, %3)           \n\t"
 
633
        "movq %%mm1, (%2, %%"REG_a")    \n\t"
 
634
        "movq (%1, %3), %%mm1           \n\t"
 
635
        "movq (%1, %%"REG_a"), %%mm0    \n\t"
 
636
        "add %%"REG_a", %2              \n\t"
 
637
        "add %%"REG_a", %1              \n\t"
 
638
        "psubusb %%mm6, %%mm1           \n\t"
 
639
        PAVGB" %%mm1, %%mm2             \n\t"
 
640
        PAVGB" %%mm0, %%mm1             \n\t"
 
641
        "movq %%mm2, (%2, %3)           \n\t"
 
642
        "movq %%mm1, (%2, %%"REG_a")    \n\t"
 
643
        "add %%"REG_a", %2              \n\t"
 
644
        "subl $4, %0                    \n\t"
 
645
        "jnz 1b                         \n\t"
 
646
        :"+g"(h), "+S"(pixels), "+D" (block)
 
647
        :"r" ((long)line_size)
 
648
        :"%"REG_a, "memory");
 
649
}
 
650
 
 
651
static void DEF(avg_pixels8)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
652
{
 
653
    __asm __volatile(
 
654
        "lea (%3, %3), %%"REG_a"        \n\t"
 
655
        "1:                             \n\t"
 
656
        "movq (%2), %%mm0               \n\t"
 
657
        "movq (%2, %3), %%mm1           \n\t"
 
658
        PAVGB" (%1), %%mm0              \n\t"
 
659
        PAVGB" (%1, %3), %%mm1          \n\t"
 
660
        "movq %%mm0, (%2)               \n\t"
 
661
        "movq %%mm1, (%2, %3)           \n\t"
 
662
        "add %%"REG_a", %1              \n\t"
 
663
        "add %%"REG_a", %2              \n\t"
 
664
        "movq (%2), %%mm0               \n\t"
 
665
        "movq (%2, %3), %%mm1           \n\t"
 
666
        PAVGB" (%1), %%mm0              \n\t"
 
667
        PAVGB" (%1, %3), %%mm1          \n\t"
 
668
        "add %%"REG_a", %1              \n\t"
 
669
        "movq %%mm0, (%2)               \n\t"
 
670
        "movq %%mm1, (%2, %3)           \n\t"
 
671
        "add %%"REG_a", %2              \n\t"
 
672
        "subl $4, %0                    \n\t"
 
673
        "jnz 1b                         \n\t"
 
674
        :"+g"(h), "+S"(pixels), "+D"(block)
 
675
        :"r" ((long)line_size)
 
676
        :"%"REG_a, "memory");
 
677
}
 
678
 
 
679
static void DEF(avg_pixels8_x2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
680
{
 
681
    __asm __volatile(
 
682
        "lea (%3, %3), %%"REG_a"        \n\t"
 
683
        "1:                             \n\t"
 
684
        "movq (%1), %%mm0               \n\t"
 
685
        "movq (%1, %3), %%mm2           \n\t"
 
686
        PAVGB" 1(%1), %%mm0             \n\t"
 
687
        PAVGB" 1(%1, %3), %%mm2         \n\t"
 
688
        PAVGB" (%2), %%mm0              \n\t"
 
689
        PAVGB" (%2, %3), %%mm2          \n\t"
 
690
        "add %%"REG_a", %1              \n\t"
 
691
        "movq %%mm0, (%2)               \n\t"
 
692
        "movq %%mm2, (%2, %3)           \n\t"
 
693
        "movq (%1), %%mm0               \n\t"
 
694
        "movq (%1, %3), %%mm2           \n\t"
 
695
        PAVGB" 1(%1), %%mm0             \n\t"
 
696
        PAVGB" 1(%1, %3), %%mm2         \n\t"
 
697
        "add %%"REG_a", %2              \n\t"
 
698
        "add %%"REG_a", %1              \n\t"
 
699
        PAVGB" (%2), %%mm0              \n\t"
 
700
        PAVGB" (%2, %3), %%mm2          \n\t"
 
701
        "movq %%mm0, (%2)               \n\t"
 
702
        "movq %%mm2, (%2, %3)           \n\t"
 
703
        "add %%"REG_a", %2              \n\t"
 
704
        "subl $4, %0                    \n\t"
 
705
        "jnz 1b                         \n\t"
 
706
        :"+g"(h), "+S"(pixels), "+D"(block)
 
707
        :"r" ((long)line_size)
 
708
        :"%"REG_a, "memory");
 
709
}
 
710
 
 
711
static void DEF(avg_pixels8_y2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
712
{
 
713
    __asm __volatile(
 
714
        "lea (%3, %3), %%"REG_a"        \n\t"
 
715
        "movq (%1), %%mm0               \n\t"
 
716
        "sub %3, %2                     \n\t"
 
717
        "1:                             \n\t"
 
718
        "movq (%1, %3), %%mm1           \n\t"
 
719
        "movq (%1, %%"REG_a"), %%mm2    \n\t"
 
720
        "add %%"REG_a", %1              \n\t"
 
721
        PAVGB" %%mm1, %%mm0             \n\t"
 
722
        PAVGB" %%mm2, %%mm1             \n\t"
 
723
        "movq (%2, %3), %%mm3           \n\t"
 
724
        "movq (%2, %%"REG_a"), %%mm4    \n\t"
 
725
        PAVGB" %%mm3, %%mm0             \n\t"
 
726
        PAVGB" %%mm4, %%mm1             \n\t"
 
727
        "movq %%mm0, (%2, %3)           \n\t"
 
728
        "movq %%mm1, (%2, %%"REG_a")    \n\t"
 
729
        "movq (%1, %3), %%mm1           \n\t"
 
730
        "movq (%1, %%"REG_a"), %%mm0    \n\t"
 
731
        PAVGB" %%mm1, %%mm2             \n\t"
 
732
        PAVGB" %%mm0, %%mm1             \n\t"
 
733
        "add %%"REG_a", %2              \n\t"
 
734
        "add %%"REG_a", %1              \n\t"
 
735
        "movq (%2, %3), %%mm3           \n\t"
 
736
        "movq (%2, %%"REG_a"), %%mm4    \n\t"
 
737
        PAVGB" %%mm3, %%mm2             \n\t"
 
738
        PAVGB" %%mm4, %%mm1             \n\t"
 
739
        "movq %%mm2, (%2, %3)           \n\t"
 
740
        "movq %%mm1, (%2, %%"REG_a")    \n\t"
 
741
        "add %%"REG_a", %2              \n\t"
 
742
        "subl $4, %0                    \n\t"
 
743
        "jnz 1b                         \n\t"
 
744
        :"+g"(h), "+S"(pixels), "+D"(block)
 
745
        :"r" ((long)line_size)
 
746
        :"%"REG_a, "memory");
 
747
}
 
748
 
 
749
// Note this is not correctly rounded, but this function is only used for b frames so it doesnt matter
 
750
static void DEF(avg_pixels8_xy2)(uint8_t *block, const uint8_t *pixels, int line_size, int h)
 
751
{
 
752
    MOVQ_BONE(mm6);
 
753
    __asm __volatile(
 
754
        "lea (%3, %3), %%"REG_a"        \n\t"
 
755
        "movq (%1), %%mm0               \n\t"
 
756
        PAVGB" 1(%1), %%mm0             \n\t"
 
757
        ".balign 8                      \n\t"
 
758
        "1:                             \n\t"
 
759
        "movq (%1, %%"REG_a"), %%mm2    \n\t"
 
760
        "movq (%1, %3), %%mm1           \n\t"
 
761
        "psubusb %%mm6, %%mm2           \n\t"
 
762
        PAVGB" 1(%1, %3), %%mm1         \n\t"
 
763
        PAVGB" 1(%1, %%"REG_a"), %%mm2  \n\t"
 
764
        "add %%"REG_a", %1              \n\t"
 
765
        PAVGB" %%mm1, %%mm0             \n\t"
 
766
        PAVGB" %%mm2, %%mm1             \n\t"
 
767
        PAVGB" (%2), %%mm0              \n\t"
 
768
        PAVGB" (%2, %3), %%mm1          \n\t"
 
769
        "movq %%mm0, (%2)               \n\t"
 
770
        "movq %%mm1, (%2, %3)           \n\t"
 
771
        "movq (%1, %3), %%mm1           \n\t"
 
772
        "movq (%1, %%"REG_a"), %%mm0    \n\t"
 
773
        PAVGB" 1(%1, %3), %%mm1         \n\t"
 
774
        PAVGB" 1(%1, %%"REG_a"), %%mm0  \n\t"
 
775
        "add %%"REG_a", %2              \n\t"
 
776
        "add %%"REG_a", %1              \n\t"
 
777
        PAVGB" %%mm1, %%mm2             \n\t"
 
778
        PAVGB" %%mm0, %%mm1             \n\t"
 
779
        PAVGB" (%2), %%mm2              \n\t"
 
780
        PAVGB" (%2, %3), %%mm1          \n\t"
 
781
        "movq %%mm2, (%2)               \n\t"
 
782
        "movq %%mm1, (%2, %3)           \n\t"
 
783
        "add %%"REG_a", %2              \n\t"
 
784
        "subl $4, %0                    \n\t"
 
785
        "jnz 1b                         \n\t"
 
786
        :"+g"(h), "+S"(pixels), "+D"(block)
 
787
        :"r" ((long)line_size)
 
788
        :"%"REG_a,  "memory");
 
789
}
 
790
 
 
791
//FIXME the following could be optimized too ...
 
792
static void DEF(put_no_rnd_pixels16_x2)(uint8_t *block, const uint8_t *pixels, int line_size, int h){
 
793
    DEF(put_no_rnd_pixels8_x2)(block  , pixels  , line_size, h);
 
794
    DEF(put_no_rnd_pixels8_x2)(block+8, pixels+8, line_size, h);
 
795
}
 
796
static void DEF(put_pixels16_y2)(uint8_t *block, const uint8_t *pixels, int line_size, int h){
 
797
    DEF(put_pixels8_y2)(block  , pixels  , line_size, h);
 
798
    DEF(put_pixels8_y2)(block+8, pixels+8, line_size, h);
 
799
}
 
800
static void DEF(put_no_rnd_pixels16_y2)(uint8_t *block, const uint8_t *pixels, int line_size, int h){
 
801
    DEF(put_no_rnd_pixels8_y2)(block  , pixels  , line_size, h);
 
802
    DEF(put_no_rnd_pixels8_y2)(block+8, pixels+8, line_size, h);
 
803
}
 
804
static void DEF(avg_pixels16)(uint8_t *block, const uint8_t *pixels, int line_size, int h){
 
805
    DEF(avg_pixels8)(block  , pixels  , line_size, h);
 
806
    DEF(avg_pixels8)(block+8, pixels+8, line_size, h);
 
807
}
 
808
static void DEF(avg_pixels16_x2)(uint8_t *block, const uint8_t *pixels, int line_size, int h){
 
809
    DEF(avg_pixels8_x2)(block  , pixels  , line_size, h);
 
810
    DEF(avg_pixels8_x2)(block+8, pixels+8, line_size, h);
 
811
}
 
812
static void DEF(avg_pixels16_y2)(uint8_t *block, const uint8_t *pixels, int line_size, int h){
 
813
    DEF(avg_pixels8_y2)(block  , pixels  , line_size, h);
 
814
    DEF(avg_pixels8_y2)(block+8, pixels+8, line_size, h);
 
815
}
 
816
static void DEF(avg_pixels16_xy2)(uint8_t *block, const uint8_t *pixels, int line_size, int h){
 
817
    DEF(avg_pixels8_xy2)(block  , pixels  , line_size, h);
 
818
    DEF(avg_pixels8_xy2)(block+8, pixels+8, line_size, h);
 
819
}
 
820