You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

342 lines
23KB

  1. /*
  2. * Half-pel DSP functions
  3. *
  4. * Copyright (c) 2000, 2001 Fabrice Bellard
  5. * Copyright (c) 2002-2004 Michael Niedermayer <michaelni@gmx.at>
  6. *
  7. * gmc & q-pel & 32/64 bit based MC by Michael Niedermayer <michaelni@gmx.at>
  8. *
  9. * This file is part of Libav.
  10. *
  11. * Libav is free software; you can redistribute it and/or
  12. * modify it under the terms of the GNU Lesser General Public
  13. * License as published by the Free Software Foundation; either
  14. * version 2.1 of the License, or (at your option) any later version.
  15. *
  16. * Libav is distributed in the hope that it will be useful,
  17. * but WITHOUT ANY WARRANTY; without even the implied warranty of
  18. * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  19. * Lesser General Public License for more details.
  20. *
  21. * You should have received a copy of the GNU Lesser General Public
  22. * License along with Libav; if not, write to the Free Software
  23. * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  24. */
  25. /**
  26. * @file
  27. * Half-pel DSP functions
  28. */
  29. #include "pixels.h"
  30. #include "bit_depth_template.c"
  31. #include "hpel_template.c"
  32. #include "tpel_template.c"
  33. #define PIXOP2(OPNAME, OP) \
  34. static inline void FUNC(OPNAME ## _no_rnd_pixels8_l2)(uint8_t *dst, \
  35. const uint8_t *src1, \
  36. const uint8_t *src2, \
  37. int dst_stride, \
  38. int src_stride1, \
  39. int src_stride2, \
  40. int h) \
  41. { \
  42. int i; \
  43. \
  44. for (i = 0; i < h; i++) { \
  45. pixel4 a, b; \
  46. a = AV_RN4P(&src1[i * src_stride1]); \
  47. b = AV_RN4P(&src2[i * src_stride2]); \
  48. OP(*((pixel4 *) &dst[i * dst_stride]), \
  49. no_rnd_avg_pixel4(a, b)); \
  50. a = AV_RN4P(&src1[i * src_stride1 + 4 * sizeof(pixel)]); \
  51. b = AV_RN4P(&src2[i * src_stride2 + 4 * sizeof(pixel)]); \
  52. OP(*((pixel4 *) &dst[i * dst_stride + 4 * sizeof(pixel)]), \
  53. no_rnd_avg_pixel4(a, b)); \
  54. } \
  55. } \
  56. \
  57. static inline void FUNCC(OPNAME ## _no_rnd_pixels8_x2)(uint8_t *block, \
  58. const uint8_t *pixels, \
  59. ptrdiff_t line_size, \
  60. int h) \
  61. { \
  62. FUNC(OPNAME ## _no_rnd_pixels8_l2)(block, pixels, pixels + sizeof(pixel), \
  63. line_size, line_size, line_size, h); \
  64. } \
  65. \
  66. static inline void FUNCC(OPNAME ## _pixels8_x2)(uint8_t *block, \
  67. const uint8_t *pixels, \
  68. ptrdiff_t line_size, \
  69. int h) \
  70. { \
  71. FUNC(OPNAME ## _pixels8_l2)(block, pixels, pixels + sizeof(pixel), \
  72. line_size, line_size, line_size, h); \
  73. } \
  74. \
  75. static inline void FUNCC(OPNAME ## _no_rnd_pixels8_y2)(uint8_t *block, \
  76. const uint8_t *pixels, \
  77. ptrdiff_t line_size, \
  78. int h) \
  79. { \
  80. FUNC(OPNAME ## _no_rnd_pixels8_l2)(block, pixels, pixels + line_size, \
  81. line_size, line_size, line_size, h); \
  82. } \
  83. \
  84. static inline void FUNCC(OPNAME ## _pixels8_y2)(uint8_t *block, \
  85. const uint8_t *pixels, \
  86. ptrdiff_t line_size, \
  87. int h) \
  88. { \
  89. FUNC(OPNAME ## _pixels8_l2)(block, pixels, pixels + line_size, \
  90. line_size, line_size, line_size, h); \
  91. } \
  92. \
  93. static inline void FUNCC(OPNAME ## _pixels4_x2)(uint8_t *block, \
  94. const uint8_t *pixels, \
  95. ptrdiff_t line_size, \
  96. int h) \
  97. { \
  98. FUNC(OPNAME ## _pixels4_l2)(block, pixels, pixels + sizeof(pixel), \
  99. line_size, line_size, line_size, h); \
  100. } \
  101. \
  102. static inline void FUNCC(OPNAME ## _pixels4_y2)(uint8_t *block, \
  103. const uint8_t *pixels, \
  104. ptrdiff_t line_size, \
  105. int h) \
  106. { \
  107. FUNC(OPNAME ## _pixels4_l2)(block, pixels, pixels + line_size, \
  108. line_size, line_size, line_size, h); \
  109. } \
  110. \
  111. static inline void FUNCC(OPNAME ## _pixels2_x2)(uint8_t *block, \
  112. const uint8_t *pixels, \
  113. ptrdiff_t line_size, \
  114. int h) \
  115. { \
  116. FUNC(OPNAME ## _pixels2_l2)(block, pixels, pixels + sizeof(pixel), \
  117. line_size, line_size, line_size, h); \
  118. } \
  119. \
  120. static inline void FUNCC(OPNAME ## _pixels2_y2)(uint8_t *block, \
  121. const uint8_t *pixels, \
  122. ptrdiff_t line_size, \
  123. int h) \
  124. { \
  125. FUNC(OPNAME ## _pixels2_l2)(block, pixels, pixels + line_size, \
  126. line_size, line_size, line_size, h); \
  127. } \
  128. \
  129. static inline void FUNCC(OPNAME ## _pixels2_xy2)(uint8_t *_block, \
  130. const uint8_t *_pixels, \
  131. ptrdiff_t line_size, \
  132. int h) \
  133. { \
  134. pixel *block = (pixel *) _block; \
  135. const pixel *pixels = (const pixel *) _pixels; \
  136. int i, a1, b1; \
  137. int a0 = pixels[0]; \
  138. int b0 = pixels[1] + 2; \
  139. \
  140. a0 += b0; \
  141. b0 += pixels[2]; \
  142. line_size /= sizeof(pixel); \
  143. pixels += line_size; \
  144. for (i = 0; i < h; i += 2) { \
  145. a1 = pixels[0]; \
  146. b1 = pixels[1]; \
  147. a1 += b1; \
  148. b1 += pixels[2]; \
  149. \
  150. block[0] = (a1 + a0) >> 2; /* FIXME non put */ \
  151. block[1] = (b1 + b0) >> 2; \
  152. \
  153. pixels += line_size; \
  154. block += line_size; \
  155. \
  156. a0 = pixels[0]; \
  157. b0 = pixels[1] + 2; \
  158. a0 += b0; \
  159. b0 += pixels[2]; \
  160. \
  161. block[0] = (a1 + a0) >> 2; \
  162. block[1] = (b1 + b0) >> 2; \
  163. pixels += line_size; \
  164. block += line_size; \
  165. } \
  166. } \
  167. \
  168. static inline void FUNCC(OPNAME ## _pixels4_xy2)(uint8_t *block, \
  169. const uint8_t *pixels, \
  170. ptrdiff_t line_size, \
  171. int h) \
  172. { \
  173. /* FIXME HIGH BIT DEPTH */ \
  174. int i; \
  175. const uint32_t a = AV_RN32(pixels); \
  176. const uint32_t b = AV_RN32(pixels + 1); \
  177. uint32_t l0 = (a & 0x03030303UL) + \
  178. (b & 0x03030303UL) + \
  179. 0x02020202UL; \
  180. uint32_t h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  181. ((b & 0xFCFCFCFCUL) >> 2); \
  182. uint32_t l1, h1; \
  183. \
  184. pixels += line_size; \
  185. for (i = 0; i < h; i += 2) { \
  186. uint32_t a = AV_RN32(pixels); \
  187. uint32_t b = AV_RN32(pixels + 1); \
  188. l1 = (a & 0x03030303UL) + \
  189. (b & 0x03030303UL); \
  190. h1 = ((a & 0xFCFCFCFCUL) >> 2) + \
  191. ((b & 0xFCFCFCFCUL) >> 2); \
  192. OP(*((uint32_t *) block), h0 + h1 + \
  193. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  194. pixels += line_size; \
  195. block += line_size; \
  196. a = AV_RN32(pixels); \
  197. b = AV_RN32(pixels + 1); \
  198. l0 = (a & 0x03030303UL) + \
  199. (b & 0x03030303UL) + \
  200. 0x02020202UL; \
  201. h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  202. ((b & 0xFCFCFCFCUL) >> 2); \
  203. OP(*((uint32_t *) block), h0 + h1 + \
  204. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  205. pixels += line_size; \
  206. block += line_size; \
  207. } \
  208. } \
  209. \
  210. static inline void FUNCC(OPNAME ## _pixels8_xy2)(uint8_t *block, \
  211. const uint8_t *pixels, \
  212. ptrdiff_t line_size, \
  213. int h) \
  214. { \
  215. /* FIXME HIGH BIT DEPTH */ \
  216. int j; \
  217. \
  218. for (j = 0; j < 2; j++) { \
  219. int i; \
  220. const uint32_t a = AV_RN32(pixels); \
  221. const uint32_t b = AV_RN32(pixels + 1); \
  222. uint32_t l0 = (a & 0x03030303UL) + \
  223. (b & 0x03030303UL) + \
  224. 0x02020202UL; \
  225. uint32_t h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  226. ((b & 0xFCFCFCFCUL) >> 2); \
  227. uint32_t l1, h1; \
  228. \
  229. pixels += line_size; \
  230. for (i = 0; i < h; i += 2) { \
  231. uint32_t a = AV_RN32(pixels); \
  232. uint32_t b = AV_RN32(pixels + 1); \
  233. l1 = (a & 0x03030303UL) + \
  234. (b & 0x03030303UL); \
  235. h1 = ((a & 0xFCFCFCFCUL) >> 2) + \
  236. ((b & 0xFCFCFCFCUL) >> 2); \
  237. OP(*((uint32_t *) block), h0 + h1 + \
  238. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  239. pixels += line_size; \
  240. block += line_size; \
  241. a = AV_RN32(pixels); \
  242. b = AV_RN32(pixels + 1); \
  243. l0 = (a & 0x03030303UL) + \
  244. (b & 0x03030303UL) + \
  245. 0x02020202UL; \
  246. h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  247. ((b & 0xFCFCFCFCUL) >> 2); \
  248. OP(*((uint32_t *) block), h0 + h1 + \
  249. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  250. pixels += line_size; \
  251. block += line_size; \
  252. } \
  253. pixels += 4 - line_size * (h + 1); \
  254. block += 4 - line_size * h; \
  255. } \
  256. } \
  257. \
  258. static inline void FUNCC(OPNAME ## _no_rnd_pixels8_xy2)(uint8_t *block, \
  259. const uint8_t *pixels, \
  260. ptrdiff_t line_size, \
  261. int h) \
  262. { \
  263. /* FIXME HIGH BIT DEPTH */ \
  264. int j; \
  265. \
  266. for (j = 0; j < 2; j++) { \
  267. int i; \
  268. const uint32_t a = AV_RN32(pixels); \
  269. const uint32_t b = AV_RN32(pixels + 1); \
  270. uint32_t l0 = (a & 0x03030303UL) + \
  271. (b & 0x03030303UL) + \
  272. 0x01010101UL; \
  273. uint32_t h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  274. ((b & 0xFCFCFCFCUL) >> 2); \
  275. uint32_t l1, h1; \
  276. \
  277. pixels += line_size; \
  278. for (i = 0; i < h; i += 2) { \
  279. uint32_t a = AV_RN32(pixels); \
  280. uint32_t b = AV_RN32(pixels + 1); \
  281. l1 = (a & 0x03030303UL) + \
  282. (b & 0x03030303UL); \
  283. h1 = ((a & 0xFCFCFCFCUL) >> 2) + \
  284. ((b & 0xFCFCFCFCUL) >> 2); \
  285. OP(*((uint32_t *) block), h0 + h1 + \
  286. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  287. pixels += line_size; \
  288. block += line_size; \
  289. a = AV_RN32(pixels); \
  290. b = AV_RN32(pixels + 1); \
  291. l0 = (a & 0x03030303UL) + \
  292. (b & 0x03030303UL) + \
  293. 0x01010101UL; \
  294. h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  295. ((b & 0xFCFCFCFCUL) >> 2); \
  296. OP(*((uint32_t *) block), h0 + h1 + \
  297. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  298. pixels += line_size; \
  299. block += line_size; \
  300. } \
  301. pixels += 4 - line_size * (h + 1); \
  302. block += 4 - line_size * h; \
  303. } \
  304. } \
  305. \
  306. CALL_2X_PIXELS(FUNCC(OPNAME ## _pixels16_x2), \
  307. FUNCC(OPNAME ## _pixels8_x2), \
  308. 8 * sizeof(pixel)) \
  309. CALL_2X_PIXELS(FUNCC(OPNAME ## _pixels16_y2), \
  310. FUNCC(OPNAME ## _pixels8_y2), \
  311. 8 * sizeof(pixel)) \
  312. CALL_2X_PIXELS(FUNCC(OPNAME ## _pixels16_xy2), \
  313. FUNCC(OPNAME ## _pixels8_xy2), \
  314. 8 * sizeof(pixel)) \
  315. av_unused CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16), \
  316. FUNCC(OPNAME ## _pixels8), \
  317. 8 * sizeof(pixel)) \
  318. CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16_x2), \
  319. FUNCC(OPNAME ## _no_rnd_pixels8_x2), \
  320. 8 * sizeof(pixel)) \
  321. CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16_y2), \
  322. FUNCC(OPNAME ## _no_rnd_pixels8_y2), \
  323. 8 * sizeof(pixel)) \
  324. CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16_xy2), \
  325. FUNCC(OPNAME ## _no_rnd_pixels8_xy2), \
  326. 8 * sizeof(pixel)) \
  327. #define op_avg(a, b) a = rnd_avg_pixel4(a, b)
  328. #define op_put(a, b) a = b
  329. #if BIT_DEPTH == 8
  330. #define put_no_rnd_pixels8_8_c put_pixels8_8_c
  331. PIXOP2(avg, op_avg)
  332. PIXOP2(put, op_put)
  333. #endif
  334. #undef op_avg
  335. #undef op_put