← Back to branch summary

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

« back to all changes in this revision

Viewing changes to avidemux/libMpeg2Dec/motion_comp_altivec.c

  • Committer: Bazaar Package Importer
  • Author(s): Matvey Kozhev
  • Date: 2007-12-18 13:53:04 UTC
  • mfrom: (1.1.7 upstream)
  • Revision ID: james.westby@ubuntu.com-20071218135304-cdqec2lg2bglyz15
Tags: 1:2.4~preview3-0.0ubuntu1
https://launchpad.net/bugs/163287
https://launchpad.net/bugs/126572
https://launchpad.net/bugs/138854
* Upload to Ubuntu. (LP: #163287, LP: #126572)
* debian/changelog: re-added Ubuntu releases.
* debian/control:
  - Require debhelper >= 5.0.51 (for dh_icons) and imagemagick.
  - Build-depend on libsdl1.2-dev instead of libsdl-dev.
  - Build against newer libx264-dev. (LP: #138854)
  - Removed libamrnb-dev, not in Ubuntu yet.
* debian/rules:
  - Install all icon sizes, using convert (upstream installs none).
  - Added missing calls to dh_installmenu, dh_installman, dh_icons and
    dh_desktop.
* debian/menu, debian/avidemux-qt.menu:
  - Corrected package and executable names.
* debian/avidemux-common.install: Install icons.
* debian/avidemux.common.manpages: Install man/avidemux.1.
* debian/links, debian/avidemux-cli.links, debian/avidemux-gtk.links:
  - Link manpages to avidemux.1.gz.
* debian/install, debian/avidemux-qt.install, debian/avidemux-gtk.desktop,
  debian/avidemux-qt.desktop: Install desktop files.

Show diffs side-by-side

added added

removed removed

Lines of Context:
avidemux/libMpeg2Dec/motion_comp_altivec.c
1
 
/*
2
 
 * motion_comp_altivec.c
3
 
 * Copyright (C) 2000-2003 Michel Lespinasse <walken@zoy.org>
4
 
 * Copyright (C) 1999-2000 Aaron Holtzman <aholtzma@ess.engr.uvic.ca>
5
 
 *
6
 
 * This file is part of mpeg2dec, a free MPEG-2 video stream decoder.
7
 
 * See http://libmpeg2.sourceforge.net/ for updates.
8
 
 *
9
 
 * mpeg2dec is free software; you can redistribute it and/or modify
10
 
 * it under the terms of the GNU General Public License as published by
11
 
 * the Free Software Foundation; either version 2 of the License, or
12
 
 * (at your option) any later version.
13
 
 *
14
 
 * mpeg2dec is distributed in the hope that it will be useful,
15
 
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16
 
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
17
 
 * GNU General Public License for more details.
18
 
 *
19
 
 * You should have received a copy of the GNU General Public License
20
 
 * along with this program; if not, write to the Free Software
21
 
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
22
 
 */
23
 
 
24
 
#include "config.h"
25
 
 
26
 
#ifdef ARCH_PPC
27
 
 
28
 
#ifdef HAVE_ALTIVEC_H
29
 
#include <altivec.h>
30
 
#endif
31
 
#include <inttypes.h>
32
 
 
33
 
#include "mpeg2.h"
34
 
#include "attributes.h"
35
 
#include "mpeg2_internal.h"
36
 
 
37
 
typedef vector signed char vector_s8_t;
38
 
typedef vector unsigned char vector_u8_t;
39
 
typedef vector signed short vector_s16_t;
40
 
typedef vector unsigned short vector_u16_t;
41
 
typedef vector signed int vector_s32_t;
42
 
typedef vector unsigned int vector_u32_t;
43
 
 
44
 
#ifndef COFFEE_BREAK    /* Workarounds for gcc suckage */
45
 
 
46
 
static inline vector_u8_t my_vec_ld (int const A, const uint8_t * const B)
47
 
{
48
 
    return vec_ld (A, (uint8_t *)B);
49
 
}
50
 
#undef vec_ld
51
 
#define vec_ld my_vec_ld
52
 
 
53
 
static inline vector_u8_t my_vec_and (vector_u8_t const A, vector_u8_t const B)
54
 
{
55
 
    return vec_and (A, B);
56
 
}
57
 
#undef vec_and
58
 
#define vec_and my_vec_and
59
 
 
60
 
static inline vector_u8_t my_vec_avg (vector_u8_t const A, vector_u8_t const B)
61
 
{
62
 
    return vec_avg (A, B);
63
 
}
64
 
#undef vec_avg
65
 
#define vec_avg my_vec_avg
66
 
 
67
 
#endif
68
 
 
69
 
static void MC_put_o_16_altivec (uint8_t * dest, const uint8_t * ref,
70
 
                                 const int stride, int height)
71
 
{
72
 
    vector_u8_t perm, ref0, ref1, tmp;
73
 
 
74
 
    perm = vec_lvsl (0, ref);
75
 
 
76
 
    height = (height >> 1) - 1;
77
 
 
78
 
    ref0 = vec_ld (0, ref);
79
 
    ref1 = vec_ld (15, ref);
80
 
    ref += stride;
81
 
    tmp = vec_perm (ref0, ref1, perm);
82
 
 
83
 
    do {
84
 
        ref0 = vec_ld (0, ref);
85
 
        ref1 = vec_ld (15, ref);
86
 
        ref += stride;
87
 
        vec_st (tmp, 0, dest);
88
 
        tmp = vec_perm (ref0, ref1, perm);
89
 
 
90
 
        ref0 = vec_ld (0, ref);
91
 
        ref1 = vec_ld (15, ref);
92
 
        ref += stride;
93
 
        vec_st (tmp, stride, dest);
94
 
        dest += 2*stride;
95
 
        tmp = vec_perm (ref0, ref1, perm);
96
 
    } while (--height);
97
 
 
98
 
    ref0 = vec_ld (0, ref);
99
 
    ref1 = vec_ld (15, ref);
100
 
    vec_st (tmp, 0, dest);
101
 
    tmp = vec_perm (ref0, ref1, perm);
102
 
    vec_st (tmp, stride, dest);
103
 
}
104
 
 
105
 
static void MC_put_o_8_altivec (uint8_t * dest, const uint8_t * ref,
106
 
                                const int stride, int height)
107
 
{
108
 
    vector_u8_t perm0, perm1, tmp0, tmp1, ref0, ref1;
109
 
 
110
 
    tmp0 = vec_lvsl (0, ref);
111
 
    tmp0 = vec_mergeh (tmp0, tmp0);
112
 
    perm0 = vec_pack ((vector_u16_t)tmp0, (vector_u16_t)tmp0);
113
 
    tmp1 = vec_lvsl (stride, ref);
114
 
    tmp1 = vec_mergeh (tmp1, tmp1);
115
 
    perm1 = vec_pack ((vector_u16_t)tmp1, (vector_u16_t)tmp1);
116
 
 
117
 
    height = (height >> 1) - 1;
118
 
 
119
 
    ref0 = vec_ld (0, ref);
120
 
    ref1 = vec_ld (7, ref);
121
 
    ref += stride;
122
 
    tmp0 = vec_perm (ref0, ref1, perm0);
123
 
 
124
 
    do {
125
 
        ref0 = vec_ld (0, ref);
126
 
        ref1 = vec_ld (7, ref);
127
 
        ref += stride;
128
 
        vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
129
 
        vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
130
 
        dest += stride;
131
 
        tmp1 = vec_perm (ref0, ref1, perm1);
132
 
 
133
 
        ref0 = vec_ld (0, ref);
134
 
        ref1 = vec_ld (7, ref);
135
 
        ref += stride;
136
 
        vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
137
 
        vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
138
 
        dest += stride;
139
 
        tmp0 = vec_perm (ref0, ref1, perm0);
140
 
    } while (--height);
141
 
 
142
 
    ref0 = vec_ld (0, ref);
143
 
    ref1 = vec_ld (7, ref);
144
 
    vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
145
 
    vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
146
 
    dest += stride;
147
 
    tmp1 = vec_perm (ref0, ref1, perm1);
148
 
    vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
149
 
    vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
150
 
}
151
 
 
152
 
static void MC_put_x_16_altivec (uint8_t * dest, const uint8_t * ref,
153
 
                                 const int stride, int height)
154
 
{
155
 
    vector_u8_t permA, permB, ref0, ref1, tmp;
156
 
 
157
 
    permA = vec_lvsl (0, ref);
158
 
    permB = vec_add (permA, vec_splat_u8 (1));
159
 
 
160
 
    height = (height >> 1) - 1;
161
 
 
162
 
    ref0 = vec_ld (0, ref);
163
 
    ref1 = vec_ld (16, ref);
164
 
    ref += stride;
165
 
    tmp = vec_avg (vec_perm (ref0, ref1, permA),
166
 
                   vec_perm (ref0, ref1, permB));
167
 
 
168
 
    do {
169
 
        ref0 = vec_ld (0, ref);
170
 
        ref1 = vec_ld (16, ref);
171
 
        ref += stride;
172
 
        vec_st (tmp, 0, dest);
173
 
        tmp = vec_avg (vec_perm (ref0, ref1, permA),
174
 
                       vec_perm (ref0, ref1, permB));
175
 
 
176
 
        ref0 = vec_ld (0, ref);
177
 
        ref1 = vec_ld (16, ref);
178
 
        ref += stride;
179
 
        vec_st (tmp, stride, dest);
180
 
        dest += 2*stride;
181
 
        tmp = vec_avg (vec_perm (ref0, ref1, permA),
182
 
                       vec_perm (ref0, ref1, permB));
183
 
    } while (--height);
184
 
 
185
 
    ref0 = vec_ld (0, ref);
186
 
    ref1 = vec_ld (16, ref);
187
 
    vec_st (tmp, 0, dest);
188
 
    tmp = vec_avg (vec_perm (ref0, ref1, permA),
189
 
                   vec_perm (ref0, ref1, permB));
190
 
    vec_st (tmp, stride, dest);
191
 
}
192
 
 
193
 
static void MC_put_x_8_altivec (uint8_t * dest, const uint8_t * ref,
194
 
                                const int stride, int height)
195
 
{
196
 
    vector_u8_t perm0A, perm0B, perm1A, perm1B, ones, tmp0, tmp1, ref0, ref1;
197
 
 
198
 
    ones = vec_splat_u8 (1);
199
 
    tmp0 = vec_lvsl (0, ref);
200
 
    tmp0 = vec_mergeh (tmp0, tmp0);
201
 
    perm0A = vec_pack ((vector_u16_t)tmp0, (vector_u16_t)tmp0);
202
 
    perm0B = vec_add (perm0A, ones);
203
 
    tmp1 = vec_lvsl (stride, ref);
204
 
    tmp1 = vec_mergeh (tmp1, tmp1);
205
 
    perm1A = vec_pack ((vector_u16_t)tmp1, (vector_u16_t)tmp1);
206
 
    perm1B = vec_add (perm1A, ones);
207
 
 
208
 
    height = (height >> 1) - 1;
209
 
 
210
 
    ref0 = vec_ld (0, ref);
211
 
    ref1 = vec_ld (8, ref);
212
 
    ref += stride;
213
 
    tmp0 = vec_avg (vec_perm (ref0, ref1, perm0A),
214
 
                    vec_perm (ref0, ref1, perm0B));
215
 
 
216
 
    do {
217
 
        ref0 = vec_ld (0, ref);
218
 
        ref1 = vec_ld (8, ref);
219
 
        ref += stride;
220
 
        vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
221
 
        vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
222
 
        dest += stride;
223
 
        tmp1 = vec_avg (vec_perm (ref0, ref1, perm1A),
224
 
                        vec_perm (ref0, ref1, perm1B));
225
 
 
226
 
        ref0 = vec_ld (0, ref);
227
 
        ref1 = vec_ld (8, ref);
228
 
        ref += stride;
229
 
        vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
230
 
        vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
231
 
        dest += stride;
232
 
        tmp0 = vec_avg (vec_perm (ref0, ref1, perm0A),
233
 
                        vec_perm (ref0, ref1, perm0B));
234
 
    } while (--height);
235
 
 
236
 
    ref0 = vec_ld (0, ref);
237
 
    ref1 = vec_ld (8, ref);
238
 
    vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
239
 
    vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
240
 
    dest += stride;
241
 
    tmp1 = vec_avg (vec_perm (ref0, ref1, perm1A),
242
 
                    vec_perm (ref0, ref1, perm1B));
243
 
    vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
244
 
    vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
245
 
}
246
 
 
247
 
static void MC_put_y_16_altivec (uint8_t * dest, const uint8_t * ref,
248
 
                                 const int stride, int height)
249
 
{
250
 
    vector_u8_t perm, ref0, ref1, tmp0, tmp1, tmp;
251
 
 
252
 
    perm = vec_lvsl (0, ref);
253
 
 
254
 
    height = (height >> 1) - 1;
255
 
 
256
 
    ref0 = vec_ld (0, ref);
257
 
    ref1 = vec_ld (15, ref);
258
 
    ref += stride;
259
 
    tmp0 = vec_perm (ref0, ref1, perm);
260
 
    ref0 = vec_ld (0, ref);
261
 
    ref1 = vec_ld (15, ref);
262
 
    ref += stride;
263
 
    tmp1 = vec_perm (ref0, ref1, perm);
264
 
    tmp = vec_avg (tmp0, tmp1);
265
 
 
266
 
    do {
267
 
        ref0 = vec_ld (0, ref);
268
 
        ref1 = vec_ld (15, ref);
269
 
        ref += stride;
270
 
        vec_st (tmp, 0, dest);
271
 
        tmp0 = vec_perm (ref0, ref1, perm);
272
 
        tmp = vec_avg (tmp0, tmp1);
273
 
 
274
 
        ref0 = vec_ld (0, ref);
275
 
        ref1 = vec_ld (15, ref);
276
 
        ref += stride;
277
 
        vec_st (tmp, stride, dest);
278
 
        dest += 2*stride;
279
 
        tmp1 = vec_perm (ref0, ref1, perm);
280
 
        tmp = vec_avg (tmp0, tmp1);
281
 
    } while (--height);
282
 
 
283
 
    ref0 = vec_ld (0, ref);
284
 
    ref1 = vec_ld (15, ref);
285
 
    vec_st (tmp, 0, dest);
286
 
    tmp0 = vec_perm (ref0, ref1, perm);
287
 
    tmp = vec_avg (tmp0, tmp1);
288
 
    vec_st (tmp, stride, dest);
289
 
}
290
 
 
291
 
static void MC_put_y_8_altivec (uint8_t * dest, const uint8_t * ref,
292
 
                                const int stride, int height)
293
 
{
294
 
    vector_u8_t perm0, perm1, tmp0, tmp1, tmp, ref0, ref1;
295
 
 
296
 
    tmp0 = vec_lvsl (0, ref);
297
 
    tmp0 = vec_mergeh (tmp0, tmp0);
298
 
    perm0 = vec_pack ((vector_u16_t)tmp0, (vector_u16_t)tmp0);
299
 
    tmp1 = vec_lvsl (stride, ref);
300
 
    tmp1 = vec_mergeh (tmp1, tmp1);
301
 
    perm1 = vec_pack ((vector_u16_t)tmp1, (vector_u16_t)tmp1);
302
 
 
303
 
    height = (height >> 1) - 1;
304
 
 
305
 
    ref0 = vec_ld (0, ref);
306
 
    ref1 = vec_ld (7, ref);
307
 
    ref += stride;
308
 
    tmp0 = vec_perm (ref0, ref1, perm0);
309
 
    ref0 = vec_ld (0, ref);
310
 
    ref1 = vec_ld (7, ref);
311
 
    ref += stride;
312
 
    tmp1 = vec_perm (ref0, ref1, perm1);
313
 
    tmp = vec_avg (tmp0, tmp1);
314
 
 
315
 
    do {
316
 
        ref0 = vec_ld (0, ref);
317
 
        ref1 = vec_ld (7, ref);
318
 
        ref += stride;
319
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
320
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
321
 
        dest += stride;
322
 
        tmp0 = vec_perm (ref0, ref1, perm0);
323
 
        tmp = vec_avg (tmp0, tmp1);
324
 
 
325
 
        ref0 = vec_ld (0, ref);
326
 
        ref1 = vec_ld (7, ref);
327
 
        ref += stride;
328
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
329
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
330
 
        dest += stride;
331
 
        tmp1 = vec_perm (ref0, ref1, perm1);
332
 
        tmp = vec_avg (tmp0, tmp1);
333
 
    } while (--height);
334
 
 
335
 
    ref0 = vec_ld (0, ref);
336
 
    ref1 = vec_ld (7, ref);
337
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
338
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
339
 
    dest += stride;
340
 
    tmp0 = vec_perm (ref0, ref1, perm0);
341
 
    tmp = vec_avg (tmp0, tmp1);
342
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
343
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
344
 
}
345
 
 
346
 
static void MC_put_xy_16_altivec (uint8_t * dest, const uint8_t * ref,
347
 
                                  const int stride, int height)
348
 
{
349
 
    vector_u8_t permA, permB, ref0, ref1, A, B, avg0, avg1, xor0, xor1, tmp;
350
 
    vector_u8_t ones;
351
 
 
352
 
    ones = vec_splat_u8 (1);
353
 
    permA = vec_lvsl (0, ref);
354
 
    permB = vec_add (permA, ones);
355
 
 
356
 
    height = (height >> 1) - 1;
357
 
 
358
 
    ref0 = vec_ld (0, ref);
359
 
    ref1 = vec_ld (16, ref);
360
 
    ref += stride;
361
 
    A = vec_perm (ref0, ref1, permA);
362
 
    B = vec_perm (ref0, ref1, permB);
363
 
    avg0 = vec_avg (A, B);
364
 
    xor0 = vec_xor (A, B);
365
 
 
366
 
    ref0 = vec_ld (0, ref);
367
 
    ref1 = vec_ld (16, ref);
368
 
    ref += stride;
369
 
    A = vec_perm (ref0, ref1, permA);
370
 
    B = vec_perm (ref0, ref1, permB);
371
 
    avg1 = vec_avg (A, B);
372
 
    xor1 = vec_xor (A, B);
373
 
    tmp = vec_sub (vec_avg (avg0, avg1),
374
 
                   vec_and (vec_and (ones, vec_or (xor0, xor1)),
375
 
                            vec_xor (avg0, avg1)));
376
 
 
377
 
    do {
378
 
        ref0 = vec_ld (0, ref);
379
 
        ref1 = vec_ld (16, ref);
380
 
        ref += stride;
381
 
        vec_st (tmp, 0, dest);
382
 
        A = vec_perm (ref0, ref1, permA);
383
 
        B = vec_perm (ref0, ref1, permB);
384
 
        avg0 = vec_avg (A, B);
385
 
        xor0 = vec_xor (A, B);
386
 
        tmp = vec_sub (vec_avg (avg0, avg1),
387
 
                       vec_and (vec_and (ones, vec_or (xor0, xor1)),
388
 
                                vec_xor (avg0, avg1)));
389
 
 
390
 
        ref0 = vec_ld (0, ref);
391
 
        ref1 = vec_ld (16, ref);
392
 
        ref += stride;
393
 
        vec_st (tmp, stride, dest);
394
 
        dest += 2*stride;
395
 
        A = vec_perm (ref0, ref1, permA);
396
 
        B = vec_perm (ref0, ref1, permB);
397
 
        avg1 = vec_avg (A, B);
398
 
        xor1 = vec_xor (A, B);
399
 
        tmp = vec_sub (vec_avg (avg0, avg1),
400
 
                       vec_and (vec_and (ones, vec_or (xor0, xor1)),
401
 
                                vec_xor (avg0, avg1)));
402
 
    } while (--height);
403
 
 
404
 
    ref0 = vec_ld (0, ref);
405
 
    ref1 = vec_ld (16, ref);
406
 
    vec_st (tmp, 0, dest);
407
 
    A = vec_perm (ref0, ref1, permA);
408
 
    B = vec_perm (ref0, ref1, permB);
409
 
    avg0 = vec_avg (A, B);
410
 
    xor0 = vec_xor (A, B);
411
 
    tmp = vec_sub (vec_avg (avg0, avg1),
412
 
                   vec_and (vec_and (ones, vec_or (xor0, xor1)),
413
 
                            vec_xor (avg0, avg1)));
414
 
    vec_st (tmp, stride, dest);
415
 
}
416
 
 
417
 
static void MC_put_xy_8_altivec (uint8_t * dest, const uint8_t * ref,
418
 
                                 const int stride, int height)
419
 
{
420
 
    vector_u8_t perm0A, perm0B, perm1A, perm1B, ref0, ref1, A, B;
421
 
    vector_u8_t avg0, avg1, xor0, xor1, tmp, ones;
422
 
 
423
 
    ones = vec_splat_u8 (1);
424
 
    perm0A = vec_lvsl (0, ref);
425
 
    perm0A = vec_mergeh (perm0A, perm0A);
426
 
    perm0A = vec_pack ((vector_u16_t)perm0A, (vector_u16_t)perm0A);
427
 
    perm0B = vec_add (perm0A, ones);
428
 
    perm1A = vec_lvsl (stride, ref);
429
 
    perm1A = vec_mergeh (perm1A, perm1A);
430
 
    perm1A = vec_pack ((vector_u16_t)perm1A, (vector_u16_t)perm1A);
431
 
    perm1B = vec_add (perm1A, ones);
432
 
 
433
 
    height = (height >> 1) - 1;
434
 
 
435
 
    ref0 = vec_ld (0, ref);
436
 
    ref1 = vec_ld (8, ref);
437
 
    ref += stride;
438
 
    A = vec_perm (ref0, ref1, perm0A);
439
 
    B = vec_perm (ref0, ref1, perm0B);
440
 
    avg0 = vec_avg (A, B);
441
 
    xor0 = vec_xor (A, B);
442
 
 
443
 
    ref0 = vec_ld (0, ref);
444
 
    ref1 = vec_ld (8, ref);
445
 
    ref += stride;
446
 
    A = vec_perm (ref0, ref1, perm1A);
447
 
    B = vec_perm (ref0, ref1, perm1B);
448
 
    avg1 = vec_avg (A, B);
449
 
    xor1 = vec_xor (A, B);
450
 
    tmp = vec_sub (vec_avg (avg0, avg1),
451
 
                   vec_and (vec_and (ones, vec_or (xor0, xor1)),
452
 
                            vec_xor (avg0, avg1)));
453
 
 
454
 
    do {
455
 
        ref0 = vec_ld (0, ref);
456
 
        ref1 = vec_ld (8, ref);
457
 
        ref += stride;
458
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
459
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
460
 
        dest += stride;
461
 
        A = vec_perm (ref0, ref1, perm0A);
462
 
        B = vec_perm (ref0, ref1, perm0B);
463
 
        avg0 = vec_avg (A, B);
464
 
        xor0 = vec_xor (A, B);
465
 
        tmp = vec_sub (vec_avg (avg0, avg1),
466
 
                       vec_and (vec_and (ones, vec_or (xor0, xor1)),
467
 
                                vec_xor (avg0, avg1)));
468
 
 
469
 
        ref0 = vec_ld (0, ref);
470
 
        ref1 = vec_ld (8, ref);
471
 
        ref += stride;
472
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
473
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
474
 
        dest += stride;
475
 
        A = vec_perm (ref0, ref1, perm1A);
476
 
        B = vec_perm (ref0, ref1, perm1B);
477
 
        avg1 = vec_avg (A, B);
478
 
        xor1 = vec_xor (A, B);
479
 
        tmp = vec_sub (vec_avg (avg0, avg1),
480
 
                       vec_and (vec_and (ones, vec_or (xor0, xor1)),
481
 
                                vec_xor (avg0, avg1)));
482
 
    } while (--height);
483
 
 
484
 
    ref0 = vec_ld (0, ref);
485
 
    ref1 = vec_ld (8, ref);
486
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
487
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
488
 
    dest += stride;
489
 
    A = vec_perm (ref0, ref1, perm0A);
490
 
    B = vec_perm (ref0, ref1, perm0B);
491
 
    avg0 = vec_avg (A, B);
492
 
    xor0 = vec_xor (A, B);
493
 
    tmp = vec_sub (vec_avg (avg0, avg1),
494
 
                   vec_and (vec_and (ones, vec_or (xor0, xor1)),
495
 
                            vec_xor (avg0, avg1)));
496
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
497
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
498
 
}
499
 
 
500
 
#if 0
501
 
static void MC_put_xy_8_altivec (uint8_t * dest, const uint8_t * ref,
502
 
                                 const int stride, int height)
503
 
{
504
 
    vector_u8_t permA, permB, ref0, ref1, A, B, C, D, tmp, zero, ones;
505
 
    vector_u16_t splat2, temp;
506
 
 
507
 
    ones = vec_splat_u8 (1);
508
 
    permA = vec_lvsl (0, ref);
509
 
    permB = vec_add (permA, ones);
510
 
 
511
 
    zero = vec_splat_u8 (0);
512
 
    splat2 = vec_splat_u16 (2);
513
 
 
514
 
    do {
515
 
        ref0 = vec_ld (0, ref);
516
 
        ref1 = vec_ld (8, ref);
517
 
        ref += stride;
518
 
        A = vec_perm (ref0, ref1, permA);
519
 
        B = vec_perm (ref0, ref1, permB);
520
 
        ref0 = vec_ld (0, ref);
521
 
        ref1 = vec_ld (8, ref);
522
 
        C = vec_perm (ref0, ref1, permA);
523
 
        D = vec_perm (ref0, ref1, permB);
524
 
 
525
 
        temp = vec_add (vec_add ((vector_u16_t)vec_mergeh (zero, A),
526
 
                                (vector_u16_t)vec_mergeh (zero, B)),
527
 
                       vec_add ((vector_u16_t)vec_mergeh (zero, C),
528
 
                                (vector_u16_t)vec_mergeh (zero, D)));
529
 
        temp = vec_sr (vec_add (temp, splat2), splat2);
530
 
        tmp = vec_pack (temp, temp);
531
 
 
532
 
        vec_st (tmp, 0, dest);
533
 
        dest += stride;
534
 
        tmp = vec_avg (vec_perm (ref0, ref1, permA),
535
 
                       vec_perm (ref0, ref1, permB));
536
 
    } while (--height);
537
 
}
538
 
#endif
539
 
 
540
 
static void MC_avg_o_16_altivec (uint8_t * dest, const uint8_t * ref,
541
 
                                 const int stride, int height)
542
 
{
543
 
    vector_u8_t perm, ref0, ref1, tmp, prev;
544
 
 
545
 
    perm = vec_lvsl (0, ref);
546
 
 
547
 
    height = (height >> 1) - 1;
548
 
 
549
 
    ref0 = vec_ld (0, ref);
550
 
    ref1 = vec_ld (15, ref);
551
 
    ref += stride;
552
 
    prev = vec_ld (0, dest);
553
 
    tmp = vec_avg (prev, vec_perm (ref0, ref1, perm));
554
 
 
555
 
    do {
556
 
        ref0 = vec_ld (0, ref);
557
 
        ref1 = vec_ld (15, ref);
558
 
        ref += stride;
559
 
        prev = vec_ld (stride, dest);
560
 
        vec_st (tmp, 0, dest);
561
 
        tmp = vec_avg (prev, vec_perm (ref0, ref1, perm));
562
 
 
563
 
        ref0 = vec_ld (0, ref);
564
 
        ref1 = vec_ld (15, ref);
565
 
        ref += stride;
566
 
        prev = vec_ld (2*stride, dest);
567
 
        vec_st (tmp, stride, dest);
568
 
        dest += 2*stride;
569
 
        tmp = vec_avg (prev, vec_perm (ref0, ref1, perm));
570
 
    } while (--height);
571
 
 
572
 
    ref0 = vec_ld (0, ref);
573
 
    ref1 = vec_ld (15, ref);
574
 
    prev = vec_ld (stride, dest);
575
 
    vec_st (tmp, 0, dest);
576
 
    tmp = vec_avg (prev, vec_perm (ref0, ref1, perm));
577
 
    vec_st (tmp, stride, dest);
578
 
}
579
 
 
580
 
static void MC_avg_o_8_altivec (uint8_t * dest, const uint8_t * ref,
581
 
                                const int stride, int height)
582
 
{
583
 
    vector_u8_t perm0, perm1, tmp0, tmp1, ref0, ref1, prev;
584
 
 
585
 
    tmp0 = vec_lvsl (0, ref);
586
 
    tmp0 = vec_mergeh (tmp0, tmp0);
587
 
    perm0 = vec_pack ((vector_u16_t)tmp0, (vector_u16_t)tmp0);
588
 
    tmp1 = vec_lvsl (stride, ref);
589
 
    tmp1 = vec_mergeh (tmp1, tmp1);
590
 
    perm1 = vec_pack ((vector_u16_t)tmp1, (vector_u16_t)tmp1);
591
 
 
592
 
    height = (height >> 1) - 1;
593
 
 
594
 
    ref0 = vec_ld (0, ref);
595
 
    ref1 = vec_ld (7, ref);
596
 
    ref += stride;
597
 
    prev = vec_ld (0, dest);
598
 
    tmp0 = vec_avg (prev, vec_perm (ref0, ref1, perm0));
599
 
 
600
 
    do {
601
 
        ref0 = vec_ld (0, ref);
602
 
        ref1 = vec_ld (7, ref);
603
 
        ref += stride;
604
 
        prev = vec_ld (stride, dest);
605
 
        vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
606
 
        vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
607
 
        dest += stride;
608
 
        tmp1 = vec_avg (prev, vec_perm (ref0, ref1, perm1));
609
 
 
610
 
        ref0 = vec_ld (0, ref);
611
 
        ref1 = vec_ld (7, ref);
612
 
        ref += stride;
613
 
        prev = vec_ld (stride, dest);
614
 
        vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
615
 
        vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
616
 
        dest += stride;
617
 
        tmp0 = vec_avg (prev, vec_perm (ref0, ref1, perm0));
618
 
    } while (--height);
619
 
 
620
 
    ref0 = vec_ld (0, ref);
621
 
    ref1 = vec_ld (7, ref);
622
 
    prev = vec_ld (stride, dest);
623
 
    vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
624
 
    vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
625
 
    dest += stride;
626
 
    tmp1 = vec_avg (prev, vec_perm (ref0, ref1, perm1));
627
 
    vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
628
 
    vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
629
 
}
630
 
 
631
 
static void MC_avg_x_16_altivec (uint8_t * dest, const uint8_t * ref,
632
 
                                 const int stride, int height)
633
 
{
634
 
    vector_u8_t permA, permB, ref0, ref1, tmp, prev;
635
 
 
636
 
    permA = vec_lvsl (0, ref);
637
 
    permB = vec_add (permA, vec_splat_u8 (1));
638
 
 
639
 
    height = (height >> 1) - 1;
640
 
 
641
 
    ref0 = vec_ld (0, ref);
642
 
    ref1 = vec_ld (16, ref);
643
 
    prev = vec_ld (0, dest);
644
 
    ref += stride;
645
 
    tmp = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, permA),
646
 
                                  vec_perm (ref0, ref1, permB)));
647
 
 
648
 
    do {
649
 
        ref0 = vec_ld (0, ref);
650
 
        ref1 = vec_ld (16, ref);
651
 
        ref += stride;
652
 
        prev = vec_ld (stride, dest);
653
 
        vec_st (tmp, 0, dest);
654
 
        tmp = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, permA),
655
 
                                      vec_perm (ref0, ref1, permB)));
656
 
 
657
 
        ref0 = vec_ld (0, ref);
658
 
        ref1 = vec_ld (16, ref);
659
 
        ref += stride;
660
 
        prev = vec_ld (2*stride, dest);
661
 
        vec_st (tmp, stride, dest);
662
 
        dest += 2*stride;
663
 
        tmp = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, permA),
664
 
                                      vec_perm (ref0, ref1, permB)));
665
 
    } while (--height);
666
 
 
667
 
    ref0 = vec_ld (0, ref);
668
 
    ref1 = vec_ld (16, ref);
669
 
    prev = vec_ld (stride, dest);
670
 
    vec_st (tmp, 0, dest);
671
 
    tmp = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, permA),
672
 
                                  vec_perm (ref0, ref1, permB)));
673
 
    vec_st (tmp, stride, dest);
674
 
}
675
 
 
676
 
static void MC_avg_x_8_altivec (uint8_t * dest, const uint8_t * ref,
677
 
                                const int stride, int height)
678
 
{
679
 
    vector_u8_t perm0A, perm0B, perm1A, perm1B, ones, tmp0, tmp1, ref0, ref1;
680
 
    vector_u8_t prev;
681
 
 
682
 
    ones = vec_splat_u8 (1);
683
 
    tmp0 = vec_lvsl (0, ref);
684
 
    tmp0 = vec_mergeh (tmp0, tmp0);
685
 
    perm0A = vec_pack ((vector_u16_t)tmp0, (vector_u16_t)tmp0);
686
 
    perm0B = vec_add (perm0A, ones);
687
 
    tmp1 = vec_lvsl (stride, ref);
688
 
    tmp1 = vec_mergeh (tmp1, tmp1);
689
 
    perm1A = vec_pack ((vector_u16_t)tmp1, (vector_u16_t)tmp1);
690
 
    perm1B = vec_add (perm1A, ones);
691
 
 
692
 
    height = (height >> 1) - 1;
693
 
 
694
 
    ref0 = vec_ld (0, ref);
695
 
    ref1 = vec_ld (8, ref);
696
 
    prev = vec_ld (0, dest);
697
 
    ref += stride;
698
 
    tmp0 = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, perm0A),
699
 
                                   vec_perm (ref0, ref1, perm0B)));
700
 
 
701
 
    do {
702
 
        ref0 = vec_ld (0, ref);
703
 
        ref1 = vec_ld (8, ref);
704
 
        ref += stride;
705
 
        prev = vec_ld (stride, dest);
706
 
        vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
707
 
        vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
708
 
        dest += stride;
709
 
        tmp1 = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, perm1A),
710
 
                                       vec_perm (ref0, ref1, perm1B)));
711
 
 
712
 
        ref0 = vec_ld (0, ref);
713
 
        ref1 = vec_ld (8, ref);
714
 
        ref += stride;
715
 
        prev = vec_ld (stride, dest);
716
 
        vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
717
 
        vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
718
 
        dest += stride;
719
 
        tmp0 = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, perm0A),
720
 
                                       vec_perm (ref0, ref1, perm0B)));
721
 
    } while (--height);
722
 
 
723
 
    ref0 = vec_ld (0, ref);
724
 
    ref1 = vec_ld (8, ref);
725
 
    prev = vec_ld (stride, dest);
726
 
    vec_ste ((vector_u32_t)tmp0, 0, (unsigned int *)dest);
727
 
    vec_ste ((vector_u32_t)tmp0, 4, (unsigned int *)dest);
728
 
    dest += stride;
729
 
    tmp1 = vec_avg (prev, vec_avg (vec_perm (ref0, ref1, perm1A),
730
 
                                   vec_perm (ref0, ref1, perm1B)));
731
 
    vec_ste ((vector_u32_t)tmp1, 0, (unsigned int *)dest);
732
 
    vec_ste ((vector_u32_t)tmp1, 4, (unsigned int *)dest);
733
 
}
734
 
 
735
 
static void MC_avg_y_16_altivec (uint8_t * dest, const uint8_t * ref,
736
 
                                 const int stride, int height)
737
 
{
738
 
    vector_u8_t perm, ref0, ref1, tmp0, tmp1, tmp, prev;
739
 
 
740
 
    perm = vec_lvsl (0, ref);
741
 
 
742
 
    height = (height >> 1) - 1;
743
 
 
744
 
    ref0 = vec_ld (0, ref);
745
 
    ref1 = vec_ld (15, ref);
746
 
    ref += stride;
747
 
    tmp0 = vec_perm (ref0, ref1, perm);
748
 
    ref0 = vec_ld (0, ref);
749
 
    ref1 = vec_ld (15, ref);
750
 
    ref += stride;
751
 
    prev = vec_ld (0, dest);
752
 
    tmp1 = vec_perm (ref0, ref1, perm);
753
 
    tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
754
 
 
755
 
    do {
756
 
        ref0 = vec_ld (0, ref);
757
 
        ref1 = vec_ld (15, ref);
758
 
        ref += stride;
759
 
        prev = vec_ld (stride, dest);
760
 
        vec_st (tmp, 0, dest);
761
 
        tmp0 = vec_perm (ref0, ref1, perm);
762
 
        tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
763
 
 
764
 
        ref0 = vec_ld (0, ref);
765
 
        ref1 = vec_ld (15, ref);
766
 
        ref += stride;
767
 
        prev = vec_ld (2*stride, dest);
768
 
        vec_st (tmp, stride, dest);
769
 
        dest += 2*stride;
770
 
        tmp1 = vec_perm (ref0, ref1, perm);
771
 
        tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
772
 
    } while (--height);
773
 
 
774
 
    ref0 = vec_ld (0, ref);
775
 
    ref1 = vec_ld (15, ref);
776
 
    prev = vec_ld (stride, dest);
777
 
    vec_st (tmp, 0, dest);
778
 
    tmp0 = vec_perm (ref0, ref1, perm);
779
 
    tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
780
 
    vec_st (tmp, stride, dest);
781
 
}
782
 
 
783
 
static void MC_avg_y_8_altivec (uint8_t * dest, const uint8_t * ref,
784
 
                                const int stride, int height)
785
 
{
786
 
    vector_u8_t perm0, perm1, tmp0, tmp1, tmp, ref0, ref1, prev;
787
 
 
788
 
    tmp0 = vec_lvsl (0, ref);
789
 
    tmp0 = vec_mergeh (tmp0, tmp0);
790
 
    perm0 = vec_pack ((vector_u16_t)tmp0, (vector_u16_t)tmp0);
791
 
    tmp1 = vec_lvsl (stride, ref);
792
 
    tmp1 = vec_mergeh (tmp1, tmp1);
793
 
    perm1 = vec_pack ((vector_u16_t)tmp1, (vector_u16_t)tmp1);
794
 
 
795
 
    height = (height >> 1) - 1;
796
 
 
797
 
    ref0 = vec_ld (0, ref);
798
 
    ref1 = vec_ld (7, ref);
799
 
    ref += stride;
800
 
    tmp0 = vec_perm (ref0, ref1, perm0);
801
 
    ref0 = vec_ld (0, ref);
802
 
    ref1 = vec_ld (7, ref);
803
 
    ref += stride;
804
 
    prev = vec_ld (0, dest);
805
 
    tmp1 = vec_perm (ref0, ref1, perm1);
806
 
    tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
807
 
 
808
 
    do {
809
 
        ref0 = vec_ld (0, ref);
810
 
        ref1 = vec_ld (7, ref);
811
 
        ref += stride;
812
 
        prev = vec_ld (stride, dest);
813
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
814
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
815
 
        dest += stride;
816
 
        tmp0 = vec_perm (ref0, ref1, perm0);
817
 
        tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
818
 
 
819
 
        ref0 = vec_ld (0, ref);
820
 
        ref1 = vec_ld (7, ref);
821
 
        ref += stride;
822
 
        prev = vec_ld (stride, dest);
823
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
824
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
825
 
        dest += stride;
826
 
        tmp1 = vec_perm (ref0, ref1, perm1);
827
 
        tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
828
 
    } while (--height);
829
 
 
830
 
    ref0 = vec_ld (0, ref);
831
 
    ref1 = vec_ld (7, ref);
832
 
    prev = vec_ld (stride, dest);
833
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
834
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
835
 
    dest += stride;
836
 
    tmp0 = vec_perm (ref0, ref1, perm0);
837
 
    tmp = vec_avg (prev, vec_avg (tmp0, tmp1));
838
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
839
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
840
 
}
841
 
 
842
 
static void MC_avg_xy_16_altivec (uint8_t * dest, const uint8_t * ref,
843
 
                                  const int stride, int height)
844
 
{
845
 
    vector_u8_t permA, permB, ref0, ref1, A, B, avg0, avg1, xor0, xor1, tmp;
846
 
    vector_u8_t ones, prev;
847
 
 
848
 
    ones = vec_splat_u8 (1);
849
 
    permA = vec_lvsl (0, ref);
850
 
    permB = vec_add (permA, ones);
851
 
 
852
 
    height = (height >> 1) - 1;
853
 
 
854
 
    ref0 = vec_ld (0, ref);
855
 
    ref1 = vec_ld (16, ref);
856
 
    ref += stride;
857
 
    A = vec_perm (ref0, ref1, permA);
858
 
    B = vec_perm (ref0, ref1, permB);
859
 
    avg0 = vec_avg (A, B);
860
 
    xor0 = vec_xor (A, B);
861
 
 
862
 
    ref0 = vec_ld (0, ref);
863
 
    ref1 = vec_ld (16, ref);
864
 
    ref += stride;
865
 
    prev = vec_ld (0, dest);
866
 
    A = vec_perm (ref0, ref1, permA);
867
 
    B = vec_perm (ref0, ref1, permB);
868
 
    avg1 = vec_avg (A, B);
869
 
    xor1 = vec_xor (A, B);
870
 
    tmp = vec_avg (prev, vec_sub (vec_avg (avg0, avg1),
871
 
                                  vec_and (vec_and (ones, vec_or (xor0, xor1)),
872
 
                                           vec_xor (avg0, avg1))));
873
 
 
874
 
    do {
875
 
        ref0 = vec_ld (0, ref);
876
 
        ref1 = vec_ld (16, ref);
877
 
        ref += stride;
878
 
        prev = vec_ld (stride, dest);
879
 
        vec_st (tmp, 0, dest);
880
 
        A = vec_perm (ref0, ref1, permA);
881
 
        B = vec_perm (ref0, ref1, permB);
882
 
        avg0 = vec_avg (A, B);
883
 
        xor0 = vec_xor (A, B);
884
 
        tmp = vec_avg (prev,
885
 
                       vec_sub (vec_avg (avg0, avg1),
886
 
                                vec_and (vec_and (ones, vec_or (xor0, xor1)),
887
 
                                         vec_xor (avg0, avg1))));
888
 
 
889
 
        ref0 = vec_ld (0, ref);
890
 
        ref1 = vec_ld (16, ref);
891
 
        ref += stride;
892
 
        prev = vec_ld (2*stride, dest);
893
 
        vec_st (tmp, stride, dest);
894
 
        dest += 2*stride;
895
 
        A = vec_perm (ref0, ref1, permA);
896
 
        B = vec_perm (ref0, ref1, permB);
897
 
        avg1 = vec_avg (A, B);
898
 
        xor1 = vec_xor (A, B);
899
 
        tmp = vec_avg (prev,
900
 
                       vec_sub (vec_avg (avg0, avg1),
901
 
                                vec_and (vec_and (ones, vec_or (xor0, xor1)),
902
 
                                         vec_xor (avg0, avg1))));
903
 
    } while (--height);
904
 
 
905
 
    ref0 = vec_ld (0, ref);
906
 
    ref1 = vec_ld (16, ref);
907
 
    prev = vec_ld (stride, dest);
908
 
    vec_st (tmp, 0, dest);
909
 
    A = vec_perm (ref0, ref1, permA);
910
 
    B = vec_perm (ref0, ref1, permB);
911
 
    avg0 = vec_avg (A, B);
912
 
    xor0 = vec_xor (A, B);
913
 
    tmp = vec_avg (prev, vec_sub (vec_avg (avg0, avg1),
914
 
                                  vec_and (vec_and (ones, vec_or (xor0, xor1)),
915
 
                                           vec_xor (avg0, avg1))));
916
 
    vec_st (tmp, stride, dest);
917
 
}
918
 
 
919
 
static void MC_avg_xy_8_altivec (uint8_t * dest, const uint8_t * ref,
920
 
                                 const int stride, int height)
921
 
{
922
 
    vector_u8_t perm0A, perm0B, perm1A, perm1B, ref0, ref1, A, B;
923
 
    vector_u8_t avg0, avg1, xor0, xor1, tmp, ones, prev;
924
 
 
925
 
    ones = vec_splat_u8 (1);
926
 
    perm0A = vec_lvsl (0, ref);
927
 
    perm0A = vec_mergeh (perm0A, perm0A);
928
 
    perm0A = vec_pack ((vector_u16_t)perm0A, (vector_u16_t)perm0A);
929
 
    perm0B = vec_add (perm0A, ones);
930
 
    perm1A = vec_lvsl (stride, ref);
931
 
    perm1A = vec_mergeh (perm1A, perm1A);
932
 
    perm1A = vec_pack ((vector_u16_t)perm1A, (vector_u16_t)perm1A);
933
 
    perm1B = vec_add (perm1A, ones);
934
 
 
935
 
    height = (height >> 1) - 1;
936
 
 
937
 
    ref0 = vec_ld (0, ref);
938
 
    ref1 = vec_ld (8, ref);
939
 
    ref += stride;
940
 
    A = vec_perm (ref0, ref1, perm0A);
941
 
    B = vec_perm (ref0, ref1, perm0B);
942
 
    avg0 = vec_avg (A, B);
943
 
    xor0 = vec_xor (A, B);
944
 
 
945
 
    ref0 = vec_ld (0, ref);
946
 
    ref1 = vec_ld (8, ref);
947
 
    ref += stride;
948
 
    prev = vec_ld (0, dest);
949
 
    A = vec_perm (ref0, ref1, perm1A);
950
 
    B = vec_perm (ref0, ref1, perm1B);
951
 
    avg1 = vec_avg (A, B);
952
 
    xor1 = vec_xor (A, B);
953
 
    tmp = vec_avg (prev, vec_sub (vec_avg (avg0, avg1),
954
 
                                  vec_and (vec_and (ones, vec_or (xor0, xor1)),
955
 
                                           vec_xor (avg0, avg1))));
956
 
 
957
 
    do {
958
 
        ref0 = vec_ld (0, ref);
959
 
        ref1 = vec_ld (8, ref);
960
 
        ref += stride;
961
 
        prev = vec_ld (stride, dest);
962
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
963
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
964
 
        dest += stride;
965
 
        A = vec_perm (ref0, ref1, perm0A);
966
 
        B = vec_perm (ref0, ref1, perm0B);
967
 
        avg0 = vec_avg (A, B);
968
 
        xor0 = vec_xor (A, B);
969
 
        tmp = vec_avg (prev,
970
 
                       vec_sub (vec_avg (avg0, avg1),
971
 
                                vec_and (vec_and (ones, vec_or (xor0, xor1)),
972
 
                                         vec_xor (avg0, avg1))));
973
 
 
974
 
        ref0 = vec_ld (0, ref);
975
 
        ref1 = vec_ld (8, ref);
976
 
        ref += stride;
977
 
        prev = vec_ld (stride, dest);
978
 
        vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
979
 
        vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
980
 
        dest += stride;
981
 
        A = vec_perm (ref0, ref1, perm1A);
982
 
        B = vec_perm (ref0, ref1, perm1B);
983
 
        avg1 = vec_avg (A, B);
984
 
        xor1 = vec_xor (A, B);
985
 
        tmp = vec_avg (prev,
986
 
                       vec_sub (vec_avg (avg0, avg1),
987
 
                                vec_and (vec_and (ones, vec_or (xor0, xor1)),
988
 
                                         vec_xor (avg0, avg1))));
989
 
    } while (--height);
990
 
 
991
 
    ref0 = vec_ld (0, ref);
992
 
    ref1 = vec_ld (8, ref);
993
 
    prev = vec_ld (stride, dest);
994
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
995
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
996
 
    dest += stride;
997
 
    A = vec_perm (ref0, ref1, perm0A);
998
 
    B = vec_perm (ref0, ref1, perm0B);
999
 
    avg0 = vec_avg (A, B);
1000
 
    xor0 = vec_xor (A, B);
1001
 
    tmp = vec_avg (prev, vec_sub (vec_avg (avg0, avg1),
1002
 
                                  vec_and (vec_and (ones, vec_or (xor0, xor1)),
1003
 
                                           vec_xor (avg0, avg1))));
1004
 
    vec_ste ((vector_u32_t)tmp, 0, (unsigned int *)dest);
1005
 
    vec_ste ((vector_u32_t)tmp, 4, (unsigned int *)dest);
1006
 
}
1007
 
 
1008
 
MPEG2_MC_EXTERN (altivec)
1009
 
 
1010
 
#endif