annotate ffmpeg/libavcodec/sh4/dsputil_align.c @ 13:844d341cf643 tip

Back up before ISMIR
author Yading Song <yading.song@eecs.qmul.ac.uk>
date Thu, 31 Oct 2013 13:17:06 +0000
parents 6840f77b83aa
children
rev   line source
yading@10 1 /*
yading@10 2 * aligned/packed access motion
yading@10 3 *
yading@10 4 * Copyright (c) 2001-2003 BERO <bero@geocities.co.jp>
yading@10 5 *
yading@10 6 * This file is part of FFmpeg.
yading@10 7 *
yading@10 8 * FFmpeg is free software; you can redistribute it and/or
yading@10 9 * modify it under the terms of the GNU Lesser General Public
yading@10 10 * License as published by the Free Software Foundation; either
yading@10 11 * version 2.1 of the License, or (at your option) any later version.
yading@10 12 *
yading@10 13 * FFmpeg is distributed in the hope that it will be useful,
yading@10 14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@10 15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
yading@10 16 * Lesser General Public License for more details.
yading@10 17 *
yading@10 18 * You should have received a copy of the GNU Lesser General Public
yading@10 19 * License along with FFmpeg; if not, write to the Free Software
yading@10 20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
yading@10 21 */
yading@10 22
yading@10 23 #include "libavutil/attributes.h"
yading@10 24 #include "libavcodec/avcodec.h"
yading@10 25 #include "libavcodec/dsputil.h"
yading@10 26 #include "libavcodec/rnd_avg.h"
yading@10 27 #include "dsputil_sh4.h"
yading@10 28
yading@10 29
yading@10 30 #define LP(p) *(uint32_t*)(p)
yading@10 31 #define LPC(p) *(const uint32_t*)(p)
yading@10 32
yading@10 33
yading@10 34 #define UNPACK(ph,pl,tt0,tt1) do { \
yading@10 35 uint32_t t0,t1; t0=tt0;t1=tt1; \
yading@10 36 ph = ( (t0 & ~BYTE_VEC32(0x03))>>2) + ( (t1 & ~BYTE_VEC32(0x03))>>2); \
yading@10 37 pl = (t0 & BYTE_VEC32(0x03)) + (t1 & BYTE_VEC32(0x03)); } while(0)
yading@10 38
yading@10 39 #define rnd_PACK(ph,pl,nph,npl) ph + nph + (((pl + npl + BYTE_VEC32(0x02))>>2) & BYTE_VEC32(0x03))
yading@10 40 #define no_rnd_PACK(ph,pl,nph,npl) ph + nph + (((pl + npl + BYTE_VEC32(0x01))>>2) & BYTE_VEC32(0x03))
yading@10 41
yading@10 42 /* little-endian */
yading@10 43 #define MERGE1(a,b,ofs) (ofs==0)?a:( ((a)>>(8*ofs))|((b)<<(32-8*ofs)) )
yading@10 44 #define MERGE2(a,b,ofs) (ofs==3)?b:( ((a)>>(8*(ofs+1)))|((b)<<(32-8*(ofs+1))) )
yading@10 45 /* big
yading@10 46 #define MERGE1(a,b,ofs) (ofs==0)?a:( ((a)<<(8*ofs))|((b)>>(32-8*ofs)) )
yading@10 47 #define MERGE2(a,b,ofs) (ofs==3)?b:( ((a)<<(8+8*ofs))|((b)>>(32-8-8*ofs)) )
yading@10 48 */
yading@10 49
yading@10 50
yading@10 51 #define put(d,s) d = s
yading@10 52 #define avg(d,s) d = rnd_avg32(s,d)
yading@10 53
yading@10 54 #define OP_C4(ofs) \
yading@10 55 ref-=ofs; \
yading@10 56 do { \
yading@10 57 OP(LP(dest),MERGE1(LPC(ref),LPC(ref+4),ofs)); \
yading@10 58 ref+=stride; \
yading@10 59 dest+=stride; \
yading@10 60 } while(--height)
yading@10 61
yading@10 62 #define OP_C40() \
yading@10 63 do { \
yading@10 64 OP(LP(dest),LPC(ref)); \
yading@10 65 ref+=stride; \
yading@10 66 dest+=stride; \
yading@10 67 } while(--height)
yading@10 68
yading@10 69 #define OP_C(ofs,sz,avg2) \
yading@10 70 { \
yading@10 71 ref-=ofs; \
yading@10 72 do { \
yading@10 73 uint32_t t0,t1; \
yading@10 74 t0 = LPC(ref+0); \
yading@10 75 t1 = LPC(ref+4); \
yading@10 76 OP(LP(dest+0), MERGE1(t0,t1,ofs)); \
yading@10 77 t0 = LPC(ref+8); \
yading@10 78 OP(LP(dest+4), MERGE1(t1,t0,ofs)); \
yading@10 79 if (sz==16) { \
yading@10 80 t1 = LPC(ref+12); \
yading@10 81 OP(LP(dest+8), MERGE1(t0,t1,ofs)); \
yading@10 82 t0 = LPC(ref+16); \
yading@10 83 OP(LP(dest+12), MERGE1(t1,t0,ofs)); \
yading@10 84 } \
yading@10 85 ref+=stride; \
yading@10 86 dest+= stride; \
yading@10 87 } while(--height); \
yading@10 88 }
yading@10 89
yading@10 90 /* aligned */
yading@10 91 #define OP_C0(sz,avg2) \
yading@10 92 { \
yading@10 93 do { \
yading@10 94 OP(LP(dest+0), LPC(ref+0)); \
yading@10 95 OP(LP(dest+4), LPC(ref+4)); \
yading@10 96 if (sz==16) { \
yading@10 97 OP(LP(dest+8), LPC(ref+8)); \
yading@10 98 OP(LP(dest+12), LPC(ref+12)); \
yading@10 99 } \
yading@10 100 ref+=stride; \
yading@10 101 dest+= stride; \
yading@10 102 } while(--height); \
yading@10 103 }
yading@10 104
yading@10 105 #define OP_X(ofs,sz,avg2) \
yading@10 106 { \
yading@10 107 ref-=ofs; \
yading@10 108 do { \
yading@10 109 uint32_t t0,t1; \
yading@10 110 t0 = LPC(ref+0); \
yading@10 111 t1 = LPC(ref+4); \
yading@10 112 OP(LP(dest+0), avg2(MERGE1(t0,t1,ofs),MERGE2(t0,t1,ofs))); \
yading@10 113 t0 = LPC(ref+8); \
yading@10 114 OP(LP(dest+4), avg2(MERGE1(t1,t0,ofs),MERGE2(t1,t0,ofs))); \
yading@10 115 if (sz==16) { \
yading@10 116 t1 = LPC(ref+12); \
yading@10 117 OP(LP(dest+8), avg2(MERGE1(t0,t1,ofs),MERGE2(t0,t1,ofs))); \
yading@10 118 t0 = LPC(ref+16); \
yading@10 119 OP(LP(dest+12), avg2(MERGE1(t1,t0,ofs),MERGE2(t1,t0,ofs))); \
yading@10 120 } \
yading@10 121 ref+=stride; \
yading@10 122 dest+= stride; \
yading@10 123 } while(--height); \
yading@10 124 }
yading@10 125
yading@10 126 /* aligned */
yading@10 127 #define OP_Y0(sz,avg2) \
yading@10 128 { \
yading@10 129 uint32_t t0,t1,t2,t3,t; \
yading@10 130 \
yading@10 131 t0 = LPC(ref+0); \
yading@10 132 t1 = LPC(ref+4); \
yading@10 133 if (sz==16) { \
yading@10 134 t2 = LPC(ref+8); \
yading@10 135 t3 = LPC(ref+12); \
yading@10 136 } \
yading@10 137 do { \
yading@10 138 ref += stride; \
yading@10 139 \
yading@10 140 t = LPC(ref+0); \
yading@10 141 OP(LP(dest+0), avg2(t0,t)); t0 = t; \
yading@10 142 t = LPC(ref+4); \
yading@10 143 OP(LP(dest+4), avg2(t1,t)); t1 = t; \
yading@10 144 if (sz==16) { \
yading@10 145 t = LPC(ref+8); \
yading@10 146 OP(LP(dest+8), avg2(t2,t)); t2 = t; \
yading@10 147 t = LPC(ref+12); \
yading@10 148 OP(LP(dest+12), avg2(t3,t)); t3 = t; \
yading@10 149 } \
yading@10 150 dest+= stride; \
yading@10 151 } while(--height); \
yading@10 152 }
yading@10 153
yading@10 154 #define OP_Y(ofs,sz,avg2) \
yading@10 155 { \
yading@10 156 uint32_t t0,t1,t2,t3,t,w0,w1; \
yading@10 157 \
yading@10 158 ref-=ofs; \
yading@10 159 w0 = LPC(ref+0); \
yading@10 160 w1 = LPC(ref+4); \
yading@10 161 t0 = MERGE1(w0,w1,ofs); \
yading@10 162 w0 = LPC(ref+8); \
yading@10 163 t1 = MERGE1(w1,w0,ofs); \
yading@10 164 if (sz==16) { \
yading@10 165 w1 = LPC(ref+12); \
yading@10 166 t2 = MERGE1(w0,w1,ofs); \
yading@10 167 w0 = LPC(ref+16); \
yading@10 168 t3 = MERGE1(w1,w0,ofs); \
yading@10 169 } \
yading@10 170 do { \
yading@10 171 ref += stride; \
yading@10 172 \
yading@10 173 w0 = LPC(ref+0); \
yading@10 174 w1 = LPC(ref+4); \
yading@10 175 t = MERGE1(w0,w1,ofs); \
yading@10 176 OP(LP(dest+0), avg2(t0,t)); t0 = t; \
yading@10 177 w0 = LPC(ref+8); \
yading@10 178 t = MERGE1(w1,w0,ofs); \
yading@10 179 OP(LP(dest+4), avg2(t1,t)); t1 = t; \
yading@10 180 if (sz==16) { \
yading@10 181 w1 = LPC(ref+12); \
yading@10 182 t = MERGE1(w0,w1,ofs); \
yading@10 183 OP(LP(dest+8), avg2(t2,t)); t2 = t; \
yading@10 184 w0 = LPC(ref+16); \
yading@10 185 t = MERGE1(w1,w0,ofs); \
yading@10 186 OP(LP(dest+12), avg2(t3,t)); t3 = t; \
yading@10 187 } \
yading@10 188 dest+=stride; \
yading@10 189 } while(--height); \
yading@10 190 }
yading@10 191
yading@10 192 #define OP_X0(sz,avg2) OP_X(0,sz,avg2)
yading@10 193 #define OP_XY0(sz,PACK) OP_XY(0,sz,PACK)
yading@10 194 #define OP_XY(ofs,sz,PACK) \
yading@10 195 { \
yading@10 196 uint32_t t2,t3,w0,w1; \
yading@10 197 uint32_t a0,a1,a2,a3,a4,a5,a6,a7; \
yading@10 198 \
yading@10 199 ref -= ofs; \
yading@10 200 w0 = LPC(ref+0); \
yading@10 201 w1 = LPC(ref+4); \
yading@10 202 UNPACK(a0,a1,MERGE1(w0,w1,ofs),MERGE2(w0,w1,ofs)); \
yading@10 203 w0 = LPC(ref+8); \
yading@10 204 UNPACK(a2,a3,MERGE1(w1,w0,ofs),MERGE2(w1,w0,ofs)); \
yading@10 205 if (sz==16) { \
yading@10 206 w1 = LPC(ref+12); \
yading@10 207 UNPACK(a4,a5,MERGE1(w0,w1,ofs),MERGE2(w0,w1,ofs)); \
yading@10 208 w0 = LPC(ref+16); \
yading@10 209 UNPACK(a6,a7,MERGE1(w1,w0,ofs),MERGE2(w1,w0,ofs)); \
yading@10 210 } \
yading@10 211 do { \
yading@10 212 ref+=stride; \
yading@10 213 w0 = LPC(ref+0); \
yading@10 214 w1 = LPC(ref+4); \
yading@10 215 UNPACK(t2,t3,MERGE1(w0,w1,ofs),MERGE2(w0,w1,ofs)); \
yading@10 216 OP(LP(dest+0),PACK(a0,a1,t2,t3)); \
yading@10 217 a0 = t2; a1 = t3; \
yading@10 218 w0 = LPC(ref+8); \
yading@10 219 UNPACK(t2,t3,MERGE1(w1,w0,ofs),MERGE2(w1,w0,ofs)); \
yading@10 220 OP(LP(dest+4),PACK(a2,a3,t2,t3)); \
yading@10 221 a2 = t2; a3 = t3; \
yading@10 222 if (sz==16) { \
yading@10 223 w1 = LPC(ref+12); \
yading@10 224 UNPACK(t2,t3,MERGE1(w0,w1,ofs),MERGE2(w0,w1,ofs)); \
yading@10 225 OP(LP(dest+8),PACK(a4,a5,t2,t3)); \
yading@10 226 a4 = t2; a5 = t3; \
yading@10 227 w0 = LPC(ref+16); \
yading@10 228 UNPACK(t2,t3,MERGE1(w1,w0,ofs),MERGE2(w1,w0,ofs)); \
yading@10 229 OP(LP(dest+12),PACK(a6,a7,t2,t3)); \
yading@10 230 a6 = t2; a7 = t3; \
yading@10 231 } \
yading@10 232 dest+=stride; \
yading@10 233 } while(--height); \
yading@10 234 }
yading@10 235
yading@10 236 #define put_pixels8_c ff_put_rnd_pixels8_o
yading@10 237 #define put_pixels16_c ff_put_rnd_pixels16_o
yading@10 238 #define avg_pixels8_c ff_avg_rnd_pixels8_o
yading@10 239 #define avg_pixels16_c ff_avg_rnd_pixels16_o
yading@10 240 #define put_no_rnd_pixels8_c ff_put_rnd_pixels8_o
yading@10 241 #define put_no_rnd_pixels16_c ff_put_rnd_pixels16_o
yading@10 242 #define avg_no_rnd_pixels16_c ff_avg_rnd_pixels16_o
yading@10 243
yading@10 244 #if CONFIG_HPELDSP
yading@10 245
yading@10 246 #include "qpel.c"
yading@10 247
yading@10 248 #endif
yading@10 249
yading@10 250 av_cold void ff_dsputil_init_align(DSPContext *c, AVCodecContext *avctx)
yading@10 251 {
yading@10 252 #if CONFIG_HPELDSP
yading@10 253
yading@10 254 #define dspfunc(PFX, IDX, NUM) \
yading@10 255 c->PFX ## _pixels_tab[IDX][ 0] = PFX ## NUM ## _mc00_sh4; \
yading@10 256 c->PFX ## _pixels_tab[IDX][ 1] = PFX ## NUM ## _mc10_sh4; \
yading@10 257 c->PFX ## _pixels_tab[IDX][ 2] = PFX ## NUM ## _mc20_sh4; \
yading@10 258 c->PFX ## _pixels_tab[IDX][ 3] = PFX ## NUM ## _mc30_sh4; \
yading@10 259 c->PFX ## _pixels_tab[IDX][ 4] = PFX ## NUM ## _mc01_sh4; \
yading@10 260 c->PFX ## _pixels_tab[IDX][ 5] = PFX ## NUM ## _mc11_sh4; \
yading@10 261 c->PFX ## _pixels_tab[IDX][ 6] = PFX ## NUM ## _mc21_sh4; \
yading@10 262 c->PFX ## _pixels_tab[IDX][ 7] = PFX ## NUM ## _mc31_sh4; \
yading@10 263 c->PFX ## _pixels_tab[IDX][ 8] = PFX ## NUM ## _mc02_sh4; \
yading@10 264 c->PFX ## _pixels_tab[IDX][ 9] = PFX ## NUM ## _mc12_sh4; \
yading@10 265 c->PFX ## _pixels_tab[IDX][10] = PFX ## NUM ## _mc22_sh4; \
yading@10 266 c->PFX ## _pixels_tab[IDX][11] = PFX ## NUM ## _mc32_sh4; \
yading@10 267 c->PFX ## _pixels_tab[IDX][12] = PFX ## NUM ## _mc03_sh4; \
yading@10 268 c->PFX ## _pixels_tab[IDX][13] = PFX ## NUM ## _mc13_sh4; \
yading@10 269 c->PFX ## _pixels_tab[IDX][14] = PFX ## NUM ## _mc23_sh4; \
yading@10 270 c->PFX ## _pixels_tab[IDX][15] = PFX ## NUM ## _mc33_sh4
yading@10 271
yading@10 272 dspfunc(put_qpel, 0, 16);
yading@10 273 dspfunc(put_no_rnd_qpel, 0, 16);
yading@10 274
yading@10 275 dspfunc(avg_qpel, 0, 16);
yading@10 276 /* dspfunc(avg_no_rnd_qpel, 0, 16); */
yading@10 277
yading@10 278 dspfunc(put_qpel, 1, 8);
yading@10 279 dspfunc(put_no_rnd_qpel, 1, 8);
yading@10 280
yading@10 281 dspfunc(avg_qpel, 1, 8);
yading@10 282 /* dspfunc(avg_no_rnd_qpel, 1, 8); */
yading@10 283
yading@10 284 #undef dspfunc
yading@10 285
yading@10 286 c->put_mspel_pixels_tab[0]= put_mspel8_mc00_sh4;
yading@10 287 c->put_mspel_pixels_tab[1]= put_mspel8_mc10_sh4;
yading@10 288 c->put_mspel_pixels_tab[2]= put_mspel8_mc20_sh4;
yading@10 289 c->put_mspel_pixels_tab[3]= put_mspel8_mc30_sh4;
yading@10 290 c->put_mspel_pixels_tab[4]= put_mspel8_mc02_sh4;
yading@10 291 c->put_mspel_pixels_tab[5]= put_mspel8_mc12_sh4;
yading@10 292 c->put_mspel_pixels_tab[6]= put_mspel8_mc22_sh4;
yading@10 293 c->put_mspel_pixels_tab[7]= put_mspel8_mc32_sh4;
yading@10 294
yading@10 295 c->gmc1 = gmc1_c;
yading@10 296
yading@10 297 #endif
yading@10 298 }