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.

341 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. #define PIXOP2(OPNAME, OP) \
  33. static inline void FUNC(OPNAME ## _no_rnd_pixels8_l2)(uint8_t *dst, \
  34. const uint8_t *src1, \
  35. const uint8_t *src2, \
  36. int dst_stride, \
  37. int src_stride1, \
  38. int src_stride2, \
  39. int h) \
  40. { \
  41. int i; \
  42. \
  43. for (i = 0; i < h; i++) { \
  44. pixel4 a, b; \
  45. a = AV_RN4P(&src1[i * src_stride1]); \
  46. b = AV_RN4P(&src2[i * src_stride2]); \
  47. OP(*((pixel4 *) &dst[i * dst_stride]), \
  48. no_rnd_avg_pixel4(a, b)); \
  49. a = AV_RN4P(&src1[i * src_stride1 + 4 * sizeof(pixel)]); \
  50. b = AV_RN4P(&src2[i * src_stride2 + 4 * sizeof(pixel)]); \
  51. OP(*((pixel4 *) &dst[i * dst_stride + 4 * sizeof(pixel)]), \
  52. no_rnd_avg_pixel4(a, b)); \
  53. } \
  54. } \
  55. \
  56. static inline void FUNCC(OPNAME ## _no_rnd_pixels8_x2)(uint8_t *block, \
  57. const uint8_t *pixels, \
  58. ptrdiff_t line_size, \
  59. int h) \
  60. { \
  61. FUNC(OPNAME ## _no_rnd_pixels8_l2)(block, pixels, pixels + sizeof(pixel), \
  62. line_size, line_size, line_size, h); \
  63. } \
  64. \
  65. static inline void FUNCC(OPNAME ## _pixels8_x2)(uint8_t *block, \
  66. const uint8_t *pixels, \
  67. ptrdiff_t line_size, \
  68. int h) \
  69. { \
  70. FUNC(OPNAME ## _pixels8_l2)(block, pixels, pixels + sizeof(pixel), \
  71. line_size, line_size, line_size, h); \
  72. } \
  73. \
  74. static inline void FUNCC(OPNAME ## _no_rnd_pixels8_y2)(uint8_t *block, \
  75. const uint8_t *pixels, \
  76. ptrdiff_t line_size, \
  77. int h) \
  78. { \
  79. FUNC(OPNAME ## _no_rnd_pixels8_l2)(block, pixels, pixels + line_size, \
  80. line_size, line_size, line_size, h); \
  81. } \
  82. \
  83. static inline void FUNCC(OPNAME ## _pixels8_y2)(uint8_t *block, \
  84. const uint8_t *pixels, \
  85. ptrdiff_t line_size, \
  86. int h) \
  87. { \
  88. FUNC(OPNAME ## _pixels8_l2)(block, pixels, pixels + line_size, \
  89. line_size, line_size, line_size, h); \
  90. } \
  91. \
  92. static inline void FUNCC(OPNAME ## _pixels4_x2)(uint8_t *block, \
  93. const uint8_t *pixels, \
  94. ptrdiff_t line_size, \
  95. int h) \
  96. { \
  97. FUNC(OPNAME ## _pixels4_l2)(block, pixels, pixels + sizeof(pixel), \
  98. line_size, line_size, line_size, h); \
  99. } \
  100. \
  101. static inline void FUNCC(OPNAME ## _pixels4_y2)(uint8_t *block, \
  102. const uint8_t *pixels, \
  103. ptrdiff_t line_size, \
  104. int h) \
  105. { \
  106. FUNC(OPNAME ## _pixels4_l2)(block, pixels, pixels + line_size, \
  107. line_size, line_size, line_size, h); \
  108. } \
  109. \
  110. static inline void FUNCC(OPNAME ## _pixels2_x2)(uint8_t *block, \
  111. const uint8_t *pixels, \
  112. ptrdiff_t line_size, \
  113. int h) \
  114. { \
  115. FUNC(OPNAME ## _pixels2_l2)(block, pixels, pixels + sizeof(pixel), \
  116. line_size, line_size, line_size, h); \
  117. } \
  118. \
  119. static inline void FUNCC(OPNAME ## _pixels2_y2)(uint8_t *block, \
  120. const uint8_t *pixels, \
  121. ptrdiff_t line_size, \
  122. int h) \
  123. { \
  124. FUNC(OPNAME ## _pixels2_l2)(block, pixels, pixels + line_size, \
  125. line_size, line_size, line_size, h); \
  126. } \
  127. \
  128. static inline void FUNCC(OPNAME ## _pixels2_xy2)(uint8_t *_block, \
  129. const uint8_t *_pixels, \
  130. ptrdiff_t line_size, \
  131. int h) \
  132. { \
  133. pixel *block = (pixel *) _block; \
  134. const pixel *pixels = (const pixel *) _pixels; \
  135. int i, a1, b1; \
  136. int a0 = pixels[0]; \
  137. int b0 = pixels[1] + 2; \
  138. \
  139. a0 += b0; \
  140. b0 += pixels[2]; \
  141. line_size /= sizeof(pixel); \
  142. pixels += line_size; \
  143. for (i = 0; i < h; i += 2) { \
  144. a1 = pixels[0]; \
  145. b1 = pixels[1]; \
  146. a1 += b1; \
  147. b1 += pixels[2]; \
  148. \
  149. block[0] = (a1 + a0) >> 2; /* FIXME non put */ \
  150. block[1] = (b1 + b0) >> 2; \
  151. \
  152. pixels += line_size; \
  153. block += line_size; \
  154. \
  155. a0 = pixels[0]; \
  156. b0 = pixels[1] + 2; \
  157. a0 += b0; \
  158. b0 += pixels[2]; \
  159. \
  160. block[0] = (a1 + a0) >> 2; \
  161. block[1] = (b1 + b0) >> 2; \
  162. pixels += line_size; \
  163. block += line_size; \
  164. } \
  165. } \
  166. \
  167. static inline void FUNCC(OPNAME ## _pixels4_xy2)(uint8_t *block, \
  168. const uint8_t *pixels, \
  169. ptrdiff_t line_size, \
  170. int h) \
  171. { \
  172. /* FIXME HIGH BIT DEPTH */ \
  173. int i; \
  174. const uint32_t a = AV_RN32(pixels); \
  175. const uint32_t b = AV_RN32(pixels + 1); \
  176. uint32_t l0 = (a & 0x03030303UL) + \
  177. (b & 0x03030303UL) + \
  178. 0x02020202UL; \
  179. uint32_t h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  180. ((b & 0xFCFCFCFCUL) >> 2); \
  181. uint32_t l1, h1; \
  182. \
  183. pixels += line_size; \
  184. for (i = 0; i < h; i += 2) { \
  185. uint32_t a = AV_RN32(pixels); \
  186. uint32_t b = AV_RN32(pixels + 1); \
  187. l1 = (a & 0x03030303UL) + \
  188. (b & 0x03030303UL); \
  189. h1 = ((a & 0xFCFCFCFCUL) >> 2) + \
  190. ((b & 0xFCFCFCFCUL) >> 2); \
  191. OP(*((uint32_t *) block), h0 + h1 + \
  192. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  193. pixels += line_size; \
  194. block += line_size; \
  195. a = AV_RN32(pixels); \
  196. b = AV_RN32(pixels + 1); \
  197. l0 = (a & 0x03030303UL) + \
  198. (b & 0x03030303UL) + \
  199. 0x02020202UL; \
  200. h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  201. ((b & 0xFCFCFCFCUL) >> 2); \
  202. OP(*((uint32_t *) block), h0 + h1 + \
  203. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  204. pixels += line_size; \
  205. block += line_size; \
  206. } \
  207. } \
  208. \
  209. static inline void FUNCC(OPNAME ## _pixels8_xy2)(uint8_t *block, \
  210. const uint8_t *pixels, \
  211. ptrdiff_t line_size, \
  212. int h) \
  213. { \
  214. /* FIXME HIGH BIT DEPTH */ \
  215. int j; \
  216. \
  217. for (j = 0; j < 2; j++) { \
  218. int i; \
  219. const uint32_t a = AV_RN32(pixels); \
  220. const uint32_t b = AV_RN32(pixels + 1); \
  221. uint32_t l0 = (a & 0x03030303UL) + \
  222. (b & 0x03030303UL) + \
  223. 0x02020202UL; \
  224. uint32_t h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  225. ((b & 0xFCFCFCFCUL) >> 2); \
  226. uint32_t l1, h1; \
  227. \
  228. pixels += line_size; \
  229. for (i = 0; i < h; i += 2) { \
  230. uint32_t a = AV_RN32(pixels); \
  231. uint32_t b = AV_RN32(pixels + 1); \
  232. l1 = (a & 0x03030303UL) + \
  233. (b & 0x03030303UL); \
  234. h1 = ((a & 0xFCFCFCFCUL) >> 2) + \
  235. ((b & 0xFCFCFCFCUL) >> 2); \
  236. OP(*((uint32_t *) block), h0 + h1 + \
  237. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  238. pixels += line_size; \
  239. block += line_size; \
  240. a = AV_RN32(pixels); \
  241. b = AV_RN32(pixels + 1); \
  242. l0 = (a & 0x03030303UL) + \
  243. (b & 0x03030303UL) + \
  244. 0x02020202UL; \
  245. h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  246. ((b & 0xFCFCFCFCUL) >> 2); \
  247. OP(*((uint32_t *) block), h0 + h1 + \
  248. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  249. pixels += line_size; \
  250. block += line_size; \
  251. } \
  252. pixels += 4 - line_size * (h + 1); \
  253. block += 4 - line_size * h; \
  254. } \
  255. } \
  256. \
  257. static inline void FUNCC(OPNAME ## _no_rnd_pixels8_xy2)(uint8_t *block, \
  258. const uint8_t *pixels, \
  259. ptrdiff_t line_size, \
  260. int h) \
  261. { \
  262. /* FIXME HIGH BIT DEPTH */ \
  263. int j; \
  264. \
  265. for (j = 0; j < 2; j++) { \
  266. int i; \
  267. const uint32_t a = AV_RN32(pixels); \
  268. const uint32_t b = AV_RN32(pixels + 1); \
  269. uint32_t l0 = (a & 0x03030303UL) + \
  270. (b & 0x03030303UL) + \
  271. 0x01010101UL; \
  272. uint32_t h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  273. ((b & 0xFCFCFCFCUL) >> 2); \
  274. uint32_t l1, h1; \
  275. \
  276. pixels += line_size; \
  277. for (i = 0; i < h; i += 2) { \
  278. uint32_t a = AV_RN32(pixels); \
  279. uint32_t b = AV_RN32(pixels + 1); \
  280. l1 = (a & 0x03030303UL) + \
  281. (b & 0x03030303UL); \
  282. h1 = ((a & 0xFCFCFCFCUL) >> 2) + \
  283. ((b & 0xFCFCFCFCUL) >> 2); \
  284. OP(*((uint32_t *) block), h0 + h1 + \
  285. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  286. pixels += line_size; \
  287. block += line_size; \
  288. a = AV_RN32(pixels); \
  289. b = AV_RN32(pixels + 1); \
  290. l0 = (a & 0x03030303UL) + \
  291. (b & 0x03030303UL) + \
  292. 0x01010101UL; \
  293. h0 = ((a & 0xFCFCFCFCUL) >> 2) + \
  294. ((b & 0xFCFCFCFCUL) >> 2); \
  295. OP(*((uint32_t *) block), h0 + h1 + \
  296. (((l0 + l1) >> 2) & 0x0F0F0F0FUL)); \
  297. pixels += line_size; \
  298. block += line_size; \
  299. } \
  300. pixels += 4 - line_size * (h + 1); \
  301. block += 4 - line_size * h; \
  302. } \
  303. } \
  304. \
  305. CALL_2X_PIXELS(FUNCC(OPNAME ## _pixels16_x2), \
  306. FUNCC(OPNAME ## _pixels8_x2), \
  307. 8 * sizeof(pixel)) \
  308. CALL_2X_PIXELS(FUNCC(OPNAME ## _pixels16_y2), \
  309. FUNCC(OPNAME ## _pixels8_y2), \
  310. 8 * sizeof(pixel)) \
  311. CALL_2X_PIXELS(FUNCC(OPNAME ## _pixels16_xy2), \
  312. FUNCC(OPNAME ## _pixels8_xy2), \
  313. 8 * sizeof(pixel)) \
  314. av_unused CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16), \
  315. FUNCC(OPNAME ## _pixels8), \
  316. 8 * sizeof(pixel)) \
  317. CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16_x2), \
  318. FUNCC(OPNAME ## _no_rnd_pixels8_x2), \
  319. 8 * sizeof(pixel)) \
  320. CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16_y2), \
  321. FUNCC(OPNAME ## _no_rnd_pixels8_y2), \
  322. 8 * sizeof(pixel)) \
  323. CALL_2X_PIXELS(FUNCC(OPNAME ## _no_rnd_pixels16_xy2), \
  324. FUNCC(OPNAME ## _no_rnd_pixels8_xy2), \
  325. 8 * sizeof(pixel)) \
  326. #define op_avg(a, b) a = rnd_avg_pixel4(a, b)
  327. #define op_put(a, b) a = b
  328. #if BIT_DEPTH == 8
  329. #define put_no_rnd_pixels8_8_c put_pixels8_8_c
  330. PIXOP2(avg, op_avg)
  331. PIXOP2(put, op_put)
  332. #endif
  333. #undef op_avg
  334. #undef op_put