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.

3538 lines
104KB

  1. /*
  2. * Copyright (C) 2003 David S. Miller <davem@redhat.com>
  3. *
  4. * This file is part of FFmpeg.
  5. *
  6. * FFmpeg 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.1 of the License, or (at your option) any later version.
  10. *
  11. * FFmpeg 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 FFmpeg; if not, write to the Free Software
  18. * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  19. */
  20. /* The *no_round* functions have been added by James A. Morrison, 2003,2004.
  21. The vis code from libmpeg2 was adapted for libavcodec by James A. Morrison.
  22. */
  23. #include "config.h"
  24. #include <inttypes.h>
  25. #include "libavcodec/hpeldsp.h"
  26. #include "libavutil/mem.h"
  27. #include "vis.h"
  28. /* The trick used in some of this file is the formula from the MMX
  29. * motion comp code, which is:
  30. *
  31. * (x+y+1)>>1 == (x|y)-((x^y)>>1)
  32. *
  33. * This allows us to average 8 bytes at a time in a 64-bit FPU reg.
  34. * We avoid overflows by masking before we do the shift, and we
  35. * implement the shift by multiplying by 1/2 using mul8x16. So in
  36. * VIS this is (assume 'x' is in f0, 'y' is in f2, a repeating mask
  37. * of '0xfe' is in f4, a repeating mask of '0x7f' is in f6, and
  38. * the value 0x80808080 is in f8):
  39. *
  40. * fxor f0, f2, f10
  41. * fand f10, f4, f10
  42. * fmul8x16 f8, f10, f10
  43. * fand f10, f6, f10
  44. * for f0, f2, f12
  45. * fpsub16 f12, f10, f10
  46. */
  47. #define DUP4(x) {x, x, x, x}
  48. #define DUP8(x) {x, x, x, x, x, x, x, x}
  49. DECLARE_ALIGNED(8, static const int16_t, constants1)[] = DUP4 (1);
  50. DECLARE_ALIGNED(8, static const int16_t, constants2)[] = DUP4 (2);
  51. DECLARE_ALIGNED(8, static const int16_t, constants3)[] = DUP4 (3);
  52. DECLARE_ALIGNED(8, static const int16_t, constants6)[] = DUP4 (6);
  53. DECLARE_ALIGNED(8, static const int8_t, constants_fe)[] = DUP8 (0xfe);
  54. DECLARE_ALIGNED(8, static const int8_t, constants_7f)[] = DUP8 (0x7f);
  55. DECLARE_ALIGNED(8, static const int8_t, constants128)[] = DUP8 (128);
  56. DECLARE_ALIGNED(8, static const int16_t, constants256_512)[] =
  57. {256, 512, 256, 512};
  58. DECLARE_ALIGNED(8, static const int16_t, constants256_1024)[] =
  59. {256, 1024, 256, 1024};
  60. #define REF_0 0
  61. #define REF_0_1 1
  62. #define REF_2 2
  63. #define REF_2_1 3
  64. #define REF_4 4
  65. #define REF_4_1 5
  66. #define REF_6 6
  67. #define REF_6_1 7
  68. #define REF_S0 8
  69. #define REF_S0_1 9
  70. #define REF_S2 10
  71. #define REF_S2_1 11
  72. #define REF_S4 12
  73. #define REF_S4_1 13
  74. #define REF_S6 14
  75. #define REF_S6_1 15
  76. #define DST_0 16
  77. #define DST_1 17
  78. #define DST_2 18
  79. #define DST_3 19
  80. #define CONST_1 20
  81. #define CONST_2 20
  82. #define CONST_3 20
  83. #define CONST_6 20
  84. #define MASK_fe 20
  85. #define CONST_128 22
  86. #define CONST_256 22
  87. #define CONST_512 22
  88. #define CONST_1024 22
  89. #define TMP0 24
  90. #define TMP1 25
  91. #define TMP2 26
  92. #define TMP3 27
  93. #define TMP4 28
  94. #define TMP5 29
  95. #define ZERO 30
  96. #define MASK_7f 30
  97. #define TMP6 32
  98. #define TMP8 34
  99. #define TMP10 36
  100. #define TMP12 38
  101. #define TMP14 40
  102. #define TMP16 42
  103. #define TMP18 44
  104. #define TMP20 46
  105. #define TMP22 48
  106. #define TMP24 50
  107. #define TMP26 52
  108. #define TMP28 54
  109. #define TMP30 56
  110. #define TMP32 58
  111. static void MC_put_o_16_vis (uint8_t * dest, const uint8_t * ref,
  112. const ptrdiff_t stride, int height)
  113. {
  114. ref = vis_alignaddr(ref);
  115. do { /* 5 cycles */
  116. vis_ld64(ref[0], TMP0);
  117. vis_ld64_2(ref, 8, TMP2);
  118. vis_ld64_2(ref, 16, TMP4);
  119. ref += stride;
  120. vis_faligndata(TMP0, TMP2, REF_0);
  121. vis_st64(REF_0, dest[0]);
  122. vis_faligndata(TMP2, TMP4, REF_2);
  123. vis_st64_2(REF_2, dest, 8);
  124. dest += stride;
  125. } while (--height);
  126. }
  127. static void MC_put_o_8_vis (uint8_t * dest, const uint8_t * ref,
  128. const ptrdiff_t stride, int height)
  129. {
  130. ref = vis_alignaddr(ref);
  131. do { /* 4 cycles */
  132. vis_ld64(ref[0], TMP0);
  133. vis_ld64(ref[8], TMP2);
  134. ref += stride;
  135. /* stall */
  136. vis_faligndata(TMP0, TMP2, REF_0);
  137. vis_st64(REF_0, dest[0]);
  138. dest += stride;
  139. } while (--height);
  140. }
  141. static void MC_avg_o_16_vis (uint8_t * dest, const uint8_t * ref,
  142. const ptrdiff_t stride, int height)
  143. {
  144. int stride_8 = stride + 8;
  145. ref = vis_alignaddr(ref);
  146. vis_ld64(ref[0], TMP0);
  147. vis_ld64(ref[8], TMP2);
  148. vis_ld64(ref[16], TMP4);
  149. vis_ld64(dest[0], DST_0);
  150. vis_ld64(dest[8], DST_2);
  151. vis_ld64(constants_fe[0], MASK_fe);
  152. vis_faligndata(TMP0, TMP2, REF_0);
  153. vis_ld64(constants_7f[0], MASK_7f);
  154. vis_faligndata(TMP2, TMP4, REF_2);
  155. vis_ld64(constants128[0], CONST_128);
  156. ref += stride;
  157. height = (height >> 1) - 1;
  158. do { /* 24 cycles */
  159. vis_ld64(ref[0], TMP0);
  160. vis_xor(DST_0, REF_0, TMP6);
  161. vis_ld64_2(ref, 8, TMP2);
  162. vis_and(TMP6, MASK_fe, TMP6);
  163. vis_ld64_2(ref, 16, TMP4);
  164. ref += stride;
  165. vis_mul8x16(CONST_128, TMP6, TMP6);
  166. vis_xor(DST_2, REF_2, TMP8);
  167. vis_and(TMP8, MASK_fe, TMP8);
  168. vis_or(DST_0, REF_0, TMP10);
  169. vis_ld64_2(dest, stride, DST_0);
  170. vis_mul8x16(CONST_128, TMP8, TMP8);
  171. vis_or(DST_2, REF_2, TMP12);
  172. vis_ld64_2(dest, stride_8, DST_2);
  173. vis_ld64(ref[0], TMP14);
  174. vis_and(TMP6, MASK_7f, TMP6);
  175. vis_and(TMP8, MASK_7f, TMP8);
  176. vis_psub16(TMP10, TMP6, TMP6);
  177. vis_st64(TMP6, dest[0]);
  178. vis_psub16(TMP12, TMP8, TMP8);
  179. vis_st64_2(TMP8, dest, 8);
  180. dest += stride;
  181. vis_ld64_2(ref, 8, TMP16);
  182. vis_faligndata(TMP0, TMP2, REF_0);
  183. vis_ld64_2(ref, 16, TMP18);
  184. vis_faligndata(TMP2, TMP4, REF_2);
  185. ref += stride;
  186. vis_xor(DST_0, REF_0, TMP20);
  187. vis_and(TMP20, MASK_fe, TMP20);
  188. vis_xor(DST_2, REF_2, TMP22);
  189. vis_mul8x16(CONST_128, TMP20, TMP20);
  190. vis_and(TMP22, MASK_fe, TMP22);
  191. vis_or(DST_0, REF_0, TMP24);
  192. vis_mul8x16(CONST_128, TMP22, TMP22);
  193. vis_or(DST_2, REF_2, TMP26);
  194. vis_ld64_2(dest, stride, DST_0);
  195. vis_faligndata(TMP14, TMP16, REF_0);
  196. vis_ld64_2(dest, stride_8, DST_2);
  197. vis_faligndata(TMP16, TMP18, REF_2);
  198. vis_and(TMP20, MASK_7f, TMP20);
  199. vis_and(TMP22, MASK_7f, TMP22);
  200. vis_psub16(TMP24, TMP20, TMP20);
  201. vis_st64(TMP20, dest[0]);
  202. vis_psub16(TMP26, TMP22, TMP22);
  203. vis_st64_2(TMP22, dest, 8);
  204. dest += stride;
  205. } while (--height);
  206. vis_ld64(ref[0], TMP0);
  207. vis_xor(DST_0, REF_0, TMP6);
  208. vis_ld64_2(ref, 8, TMP2);
  209. vis_and(TMP6, MASK_fe, TMP6);
  210. vis_ld64_2(ref, 16, TMP4);
  211. vis_mul8x16(CONST_128, TMP6, TMP6);
  212. vis_xor(DST_2, REF_2, TMP8);
  213. vis_and(TMP8, MASK_fe, TMP8);
  214. vis_or(DST_0, REF_0, TMP10);
  215. vis_ld64_2(dest, stride, DST_0);
  216. vis_mul8x16(CONST_128, TMP8, TMP8);
  217. vis_or(DST_2, REF_2, TMP12);
  218. vis_ld64_2(dest, stride_8, DST_2);
  219. vis_ld64(ref[0], TMP14);
  220. vis_and(TMP6, MASK_7f, TMP6);
  221. vis_and(TMP8, MASK_7f, TMP8);
  222. vis_psub16(TMP10, TMP6, TMP6);
  223. vis_st64(TMP6, dest[0]);
  224. vis_psub16(TMP12, TMP8, TMP8);
  225. vis_st64_2(TMP8, dest, 8);
  226. dest += stride;
  227. vis_faligndata(TMP0, TMP2, REF_0);
  228. vis_faligndata(TMP2, TMP4, REF_2);
  229. vis_xor(DST_0, REF_0, TMP20);
  230. vis_and(TMP20, MASK_fe, TMP20);
  231. vis_xor(DST_2, REF_2, TMP22);
  232. vis_mul8x16(CONST_128, TMP20, TMP20);
  233. vis_and(TMP22, MASK_fe, TMP22);
  234. vis_or(DST_0, REF_0, TMP24);
  235. vis_mul8x16(CONST_128, TMP22, TMP22);
  236. vis_or(DST_2, REF_2, TMP26);
  237. vis_and(TMP20, MASK_7f, TMP20);
  238. vis_and(TMP22, MASK_7f, TMP22);
  239. vis_psub16(TMP24, TMP20, TMP20);
  240. vis_st64(TMP20, dest[0]);
  241. vis_psub16(TMP26, TMP22, TMP22);
  242. vis_st64_2(TMP22, dest, 8);
  243. }
  244. static void MC_avg_o_8_vis (uint8_t * dest, const uint8_t * ref,
  245. const ptrdiff_t stride, int height)
  246. {
  247. ref = vis_alignaddr(ref);
  248. vis_ld64(ref[0], TMP0);
  249. vis_ld64(ref[8], TMP2);
  250. vis_ld64(dest[0], DST_0);
  251. vis_ld64(constants_fe[0], MASK_fe);
  252. vis_ld64(constants_7f[0], MASK_7f);
  253. vis_faligndata(TMP0, TMP2, REF_0);
  254. vis_ld64(constants128[0], CONST_128);
  255. ref += stride;
  256. height = (height >> 1) - 1;
  257. do { /* 12 cycles */
  258. vis_ld64(ref[0], TMP0);
  259. vis_xor(DST_0, REF_0, TMP4);
  260. vis_ld64(ref[8], TMP2);
  261. vis_and(TMP4, MASK_fe, TMP4);
  262. vis_or(DST_0, REF_0, TMP6);
  263. vis_ld64_2(dest, stride, DST_0);
  264. ref += stride;
  265. vis_mul8x16(CONST_128, TMP4, TMP4);
  266. vis_ld64(ref[0], TMP12);
  267. vis_faligndata(TMP0, TMP2, REF_0);
  268. vis_ld64(ref[8], TMP2);
  269. vis_xor(DST_0, REF_0, TMP0);
  270. ref += stride;
  271. vis_and(TMP0, MASK_fe, TMP0);
  272. vis_and(TMP4, MASK_7f, TMP4);
  273. vis_psub16(TMP6, TMP4, TMP4);
  274. vis_st64(TMP4, dest[0]);
  275. dest += stride;
  276. vis_mul8x16(CONST_128, TMP0, TMP0);
  277. vis_or(DST_0, REF_0, TMP6);
  278. vis_ld64_2(dest, stride, DST_0);
  279. vis_faligndata(TMP12, TMP2, REF_0);
  280. vis_and(TMP0, MASK_7f, TMP0);
  281. vis_psub16(TMP6, TMP0, TMP4);
  282. vis_st64(TMP4, dest[0]);
  283. dest += stride;
  284. } while (--height);
  285. vis_ld64(ref[0], TMP0);
  286. vis_xor(DST_0, REF_0, TMP4);
  287. vis_ld64(ref[8], TMP2);
  288. vis_and(TMP4, MASK_fe, TMP4);
  289. vis_or(DST_0, REF_0, TMP6);
  290. vis_ld64_2(dest, stride, DST_0);
  291. vis_mul8x16(CONST_128, TMP4, TMP4);
  292. vis_faligndata(TMP0, TMP2, REF_0);
  293. vis_xor(DST_0, REF_0, TMP0);
  294. vis_and(TMP0, MASK_fe, TMP0);
  295. vis_and(TMP4, MASK_7f, TMP4);
  296. vis_psub16(TMP6, TMP4, TMP4);
  297. vis_st64(TMP4, dest[0]);
  298. dest += stride;
  299. vis_mul8x16(CONST_128, TMP0, TMP0);
  300. vis_or(DST_0, REF_0, TMP6);
  301. vis_and(TMP0, MASK_7f, TMP0);
  302. vis_psub16(TMP6, TMP0, TMP4);
  303. vis_st64(TMP4, dest[0]);
  304. }
  305. static void MC_put_x_16_vis (uint8_t * dest, const uint8_t * ref,
  306. const ptrdiff_t stride, int height)
  307. {
  308. unsigned long off = (unsigned long) ref & 0x7;
  309. unsigned long off_plus_1 = off + 1;
  310. ref = vis_alignaddr(ref);
  311. vis_ld64(ref[0], TMP0);
  312. vis_ld64_2(ref, 8, TMP2);
  313. vis_ld64_2(ref, 16, TMP4);
  314. vis_ld64(constants_fe[0], MASK_fe);
  315. vis_ld64(constants_7f[0], MASK_7f);
  316. vis_faligndata(TMP0, TMP2, REF_0);
  317. vis_ld64(constants128[0], CONST_128);
  318. vis_faligndata(TMP2, TMP4, REF_4);
  319. if (off != 0x7) {
  320. vis_alignaddr_g0((void *)off_plus_1);
  321. vis_faligndata(TMP0, TMP2, REF_2);
  322. vis_faligndata(TMP2, TMP4, REF_6);
  323. } else {
  324. vis_src1(TMP2, REF_2);
  325. vis_src1(TMP4, REF_6);
  326. }
  327. ref += stride;
  328. height = (height >> 1) - 1;
  329. do { /* 34 cycles */
  330. vis_ld64(ref[0], TMP0);
  331. vis_xor(REF_0, REF_2, TMP6);
  332. vis_ld64_2(ref, 8, TMP2);
  333. vis_xor(REF_4, REF_6, TMP8);
  334. vis_ld64_2(ref, 16, TMP4);
  335. vis_and(TMP6, MASK_fe, TMP6);
  336. ref += stride;
  337. vis_ld64(ref[0], TMP14);
  338. vis_mul8x16(CONST_128, TMP6, TMP6);
  339. vis_and(TMP8, MASK_fe, TMP8);
  340. vis_ld64_2(ref, 8, TMP16);
  341. vis_mul8x16(CONST_128, TMP8, TMP8);
  342. vis_or(REF_0, REF_2, TMP10);
  343. vis_ld64_2(ref, 16, TMP18);
  344. ref += stride;
  345. vis_or(REF_4, REF_6, TMP12);
  346. vis_alignaddr_g0((void *)off);
  347. vis_faligndata(TMP0, TMP2, REF_0);
  348. vis_faligndata(TMP2, TMP4, REF_4);
  349. if (off != 0x7) {
  350. vis_alignaddr_g0((void *)off_plus_1);
  351. vis_faligndata(TMP0, TMP2, REF_2);
  352. vis_faligndata(TMP2, TMP4, REF_6);
  353. } else {
  354. vis_src1(TMP2, REF_2);
  355. vis_src1(TMP4, REF_6);
  356. }
  357. vis_and(TMP6, MASK_7f, TMP6);
  358. vis_and(TMP8, MASK_7f, TMP8);
  359. vis_psub16(TMP10, TMP6, TMP6);
  360. vis_st64(TMP6, dest[0]);
  361. vis_psub16(TMP12, TMP8, TMP8);
  362. vis_st64_2(TMP8, dest, 8);
  363. dest += stride;
  364. vis_xor(REF_0, REF_2, TMP6);
  365. vis_xor(REF_4, REF_6, TMP8);
  366. vis_and(TMP6, MASK_fe, TMP6);
  367. vis_mul8x16(CONST_128, TMP6, TMP6);
  368. vis_and(TMP8, MASK_fe, TMP8);
  369. vis_mul8x16(CONST_128, TMP8, TMP8);
  370. vis_or(REF_0, REF_2, TMP10);
  371. vis_or(REF_4, REF_6, TMP12);
  372. vis_alignaddr_g0((void *)off);
  373. vis_faligndata(TMP14, TMP16, REF_0);
  374. vis_faligndata(TMP16, TMP18, REF_4);
  375. if (off != 0x7) {
  376. vis_alignaddr_g0((void *)off_plus_1);
  377. vis_faligndata(TMP14, TMP16, REF_2);
  378. vis_faligndata(TMP16, TMP18, REF_6);
  379. } else {
  380. vis_src1(TMP16, REF_2);
  381. vis_src1(TMP18, REF_6);
  382. }
  383. vis_and(TMP6, MASK_7f, TMP6);
  384. vis_and(TMP8, MASK_7f, TMP8);
  385. vis_psub16(TMP10, TMP6, TMP6);
  386. vis_st64(TMP6, dest[0]);
  387. vis_psub16(TMP12, TMP8, TMP8);
  388. vis_st64_2(TMP8, dest, 8);
  389. dest += stride;
  390. } while (--height);
  391. vis_ld64(ref[0], TMP0);
  392. vis_xor(REF_0, REF_2, TMP6);
  393. vis_ld64_2(ref, 8, TMP2);
  394. vis_xor(REF_4, REF_6, TMP8);
  395. vis_ld64_2(ref, 16, TMP4);
  396. vis_and(TMP6, MASK_fe, TMP6);
  397. vis_mul8x16(CONST_128, TMP6, TMP6);
  398. vis_and(TMP8, MASK_fe, TMP8);
  399. vis_mul8x16(CONST_128, TMP8, TMP8);
  400. vis_or(REF_0, REF_2, TMP10);
  401. vis_or(REF_4, REF_6, TMP12);
  402. vis_alignaddr_g0((void *)off);
  403. vis_faligndata(TMP0, TMP2, REF_0);
  404. vis_faligndata(TMP2, TMP4, REF_4);
  405. if (off != 0x7) {
  406. vis_alignaddr_g0((void *)off_plus_1);
  407. vis_faligndata(TMP0, TMP2, REF_2);
  408. vis_faligndata(TMP2, TMP4, REF_6);
  409. } else {
  410. vis_src1(TMP2, REF_2);
  411. vis_src1(TMP4, REF_6);
  412. }
  413. vis_and(TMP6, MASK_7f, TMP6);
  414. vis_and(TMP8, MASK_7f, TMP8);
  415. vis_psub16(TMP10, TMP6, TMP6);
  416. vis_st64(TMP6, dest[0]);
  417. vis_psub16(TMP12, TMP8, TMP8);
  418. vis_st64_2(TMP8, dest, 8);
  419. dest += stride;
  420. vis_xor(REF_0, REF_2, TMP6);
  421. vis_xor(REF_4, REF_6, TMP8);
  422. vis_and(TMP6, MASK_fe, TMP6);
  423. vis_mul8x16(CONST_128, TMP6, TMP6);
  424. vis_and(TMP8, MASK_fe, TMP8);
  425. vis_mul8x16(CONST_128, TMP8, TMP8);
  426. vis_or(REF_0, REF_2, TMP10);
  427. vis_or(REF_4, REF_6, TMP12);
  428. vis_and(TMP6, MASK_7f, TMP6);
  429. vis_and(TMP8, MASK_7f, TMP8);
  430. vis_psub16(TMP10, TMP6, TMP6);
  431. vis_st64(TMP6, dest[0]);
  432. vis_psub16(TMP12, TMP8, TMP8);
  433. vis_st64_2(TMP8, dest, 8);
  434. }
  435. static void MC_put_x_8_vis (uint8_t * dest, const uint8_t * ref,
  436. const ptrdiff_t stride, int height)
  437. {
  438. unsigned long off = (unsigned long) ref & 0x7;
  439. unsigned long off_plus_1 = off + 1;
  440. ref = vis_alignaddr(ref);
  441. vis_ld64(ref[0], TMP0);
  442. vis_ld64(ref[8], TMP2);
  443. vis_ld64(constants_fe[0], MASK_fe);
  444. vis_ld64(constants_7f[0], MASK_7f);
  445. vis_ld64(constants128[0], CONST_128);
  446. vis_faligndata(TMP0, TMP2, REF_0);
  447. if (off != 0x7) {
  448. vis_alignaddr_g0((void *)off_plus_1);
  449. vis_faligndata(TMP0, TMP2, REF_2);
  450. } else {
  451. vis_src1(TMP2, REF_2);
  452. }
  453. ref += stride;
  454. height = (height >> 1) - 1;
  455. do { /* 20 cycles */
  456. vis_ld64(ref[0], TMP0);
  457. vis_xor(REF_0, REF_2, TMP4);
  458. vis_ld64_2(ref, 8, TMP2);
  459. vis_and(TMP4, MASK_fe, TMP4);
  460. ref += stride;
  461. vis_ld64(ref[0], TMP8);
  462. vis_or(REF_0, REF_2, TMP6);
  463. vis_mul8x16(CONST_128, TMP4, TMP4);
  464. vis_alignaddr_g0((void *)off);
  465. vis_ld64_2(ref, 8, TMP10);
  466. ref += stride;
  467. vis_faligndata(TMP0, TMP2, REF_0);
  468. if (off != 0x7) {
  469. vis_alignaddr_g0((void *)off_plus_1);
  470. vis_faligndata(TMP0, TMP2, REF_2);
  471. } else {
  472. vis_src1(TMP2, REF_2);
  473. }
  474. vis_and(TMP4, MASK_7f, TMP4);
  475. vis_psub16(TMP6, TMP4, DST_0);
  476. vis_st64(DST_0, dest[0]);
  477. dest += stride;
  478. vis_xor(REF_0, REF_2, TMP12);
  479. vis_and(TMP12, MASK_fe, TMP12);
  480. vis_or(REF_0, REF_2, TMP14);
  481. vis_mul8x16(CONST_128, TMP12, TMP12);
  482. vis_alignaddr_g0((void *)off);
  483. vis_faligndata(TMP8, TMP10, REF_0);
  484. if (off != 0x7) {
  485. vis_alignaddr_g0((void *)off_plus_1);
  486. vis_faligndata(TMP8, TMP10, REF_2);
  487. } else {
  488. vis_src1(TMP10, REF_2);
  489. }
  490. vis_and(TMP12, MASK_7f, TMP12);
  491. vis_psub16(TMP14, TMP12, DST_0);
  492. vis_st64(DST_0, dest[0]);
  493. dest += stride;
  494. } while (--height);
  495. vis_ld64(ref[0], TMP0);
  496. vis_xor(REF_0, REF_2, TMP4);
  497. vis_ld64_2(ref, 8, TMP2);
  498. vis_and(TMP4, MASK_fe, TMP4);
  499. vis_or(REF_0, REF_2, TMP6);
  500. vis_mul8x16(CONST_128, TMP4, TMP4);
  501. vis_alignaddr_g0((void *)off);
  502. vis_faligndata(TMP0, TMP2, REF_0);
  503. if (off != 0x7) {
  504. vis_alignaddr_g0((void *)off_plus_1);
  505. vis_faligndata(TMP0, TMP2, REF_2);
  506. } else {
  507. vis_src1(TMP2, REF_2);
  508. }
  509. vis_and(TMP4, MASK_7f, TMP4);
  510. vis_psub16(TMP6, TMP4, DST_0);
  511. vis_st64(DST_0, dest[0]);
  512. dest += stride;
  513. vis_xor(REF_0, REF_2, TMP12);
  514. vis_and(TMP12, MASK_fe, TMP12);
  515. vis_or(REF_0, REF_2, TMP14);
  516. vis_mul8x16(CONST_128, TMP12, TMP12);
  517. vis_and(TMP12, MASK_7f, TMP12);
  518. vis_psub16(TMP14, TMP12, DST_0);
  519. vis_st64(DST_0, dest[0]);
  520. dest += stride;
  521. }
  522. static void MC_avg_x_16_vis (uint8_t * dest, const uint8_t * ref,
  523. const ptrdiff_t stride, int height)
  524. {
  525. unsigned long off = (unsigned long) ref & 0x7;
  526. unsigned long off_plus_1 = off + 1;
  527. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  528. vis_ld64(constants3[0], CONST_3);
  529. vis_fzero(ZERO);
  530. vis_ld64(constants256_512[0], CONST_256);
  531. ref = vis_alignaddr(ref);
  532. do { /* 26 cycles */
  533. vis_ld64(ref[0], TMP0);
  534. vis_ld64(ref[8], TMP2);
  535. vis_alignaddr_g0((void *)off);
  536. vis_ld64(ref[16], TMP4);
  537. vis_ld64(dest[0], DST_0);
  538. vis_faligndata(TMP0, TMP2, REF_0);
  539. vis_ld64(dest[8], DST_2);
  540. vis_faligndata(TMP2, TMP4, REF_4);
  541. if (off != 0x7) {
  542. vis_alignaddr_g0((void *)off_plus_1);
  543. vis_faligndata(TMP0, TMP2, REF_2);
  544. vis_faligndata(TMP2, TMP4, REF_6);
  545. } else {
  546. vis_src1(TMP2, REF_2);
  547. vis_src1(TMP4, REF_6);
  548. }
  549. vis_mul8x16au(REF_0, CONST_256, TMP0);
  550. vis_pmerge(ZERO, REF_2, TMP4);
  551. vis_mul8x16au(REF_0_1, CONST_256, TMP2);
  552. vis_pmerge(ZERO, REF_2_1, TMP6);
  553. vis_padd16(TMP0, TMP4, TMP0);
  554. vis_mul8x16al(DST_0, CONST_512, TMP4);
  555. vis_padd16(TMP2, TMP6, TMP2);
  556. vis_mul8x16al(DST_1, CONST_512, TMP6);
  557. vis_mul8x16au(REF_6, CONST_256, TMP12);
  558. vis_padd16(TMP0, TMP4, TMP0);
  559. vis_mul8x16au(REF_6_1, CONST_256, TMP14);
  560. vis_padd16(TMP2, TMP6, TMP2);
  561. vis_mul8x16au(REF_4, CONST_256, TMP16);
  562. vis_padd16(TMP0, CONST_3, TMP8);
  563. vis_mul8x16au(REF_4_1, CONST_256, TMP18);
  564. vis_padd16(TMP2, CONST_3, TMP10);
  565. vis_pack16(TMP8, DST_0);
  566. vis_pack16(TMP10, DST_1);
  567. vis_padd16(TMP16, TMP12, TMP0);
  568. vis_st64(DST_0, dest[0]);
  569. vis_mul8x16al(DST_2, CONST_512, TMP4);
  570. vis_padd16(TMP18, TMP14, TMP2);
  571. vis_mul8x16al(DST_3, CONST_512, TMP6);
  572. vis_padd16(TMP0, CONST_3, TMP0);
  573. vis_padd16(TMP2, CONST_3, TMP2);
  574. vis_padd16(TMP0, TMP4, TMP0);
  575. vis_padd16(TMP2, TMP6, TMP2);
  576. vis_pack16(TMP0, DST_2);
  577. vis_pack16(TMP2, DST_3);
  578. vis_st64(DST_2, dest[8]);
  579. ref += stride;
  580. dest += stride;
  581. } while (--height);
  582. }
  583. static void MC_avg_x_8_vis (uint8_t * dest, const uint8_t * ref,
  584. const ptrdiff_t stride, int height)
  585. {
  586. unsigned long off = (unsigned long) ref & 0x7;
  587. unsigned long off_plus_1 = off + 1;
  588. int stride_times_2 = stride << 1;
  589. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  590. vis_ld64(constants3[0], CONST_3);
  591. vis_fzero(ZERO);
  592. vis_ld64(constants256_512[0], CONST_256);
  593. ref = vis_alignaddr(ref);
  594. height >>= 2;
  595. do { /* 47 cycles */
  596. vis_ld64(ref[0], TMP0);
  597. vis_ld64_2(ref, 8, TMP2);
  598. ref += stride;
  599. vis_alignaddr_g0((void *)off);
  600. vis_ld64(ref[0], TMP4);
  601. vis_faligndata(TMP0, TMP2, REF_0);
  602. vis_ld64_2(ref, 8, TMP6);
  603. ref += stride;
  604. vis_ld64(ref[0], TMP8);
  605. vis_ld64_2(ref, 8, TMP10);
  606. ref += stride;
  607. vis_faligndata(TMP4, TMP6, REF_4);
  608. vis_ld64(ref[0], TMP12);
  609. vis_ld64_2(ref, 8, TMP14);
  610. ref += stride;
  611. vis_faligndata(TMP8, TMP10, REF_S0);
  612. vis_faligndata(TMP12, TMP14, REF_S4);
  613. if (off != 0x7) {
  614. vis_alignaddr_g0((void *)off_plus_1);
  615. vis_ld64(dest[0], DST_0);
  616. vis_faligndata(TMP0, TMP2, REF_2);
  617. vis_ld64_2(dest, stride, DST_2);
  618. vis_faligndata(TMP4, TMP6, REF_6);
  619. vis_faligndata(TMP8, TMP10, REF_S2);
  620. vis_faligndata(TMP12, TMP14, REF_S6);
  621. } else {
  622. vis_ld64(dest[0], DST_0);
  623. vis_src1(TMP2, REF_2);
  624. vis_ld64_2(dest, stride, DST_2);
  625. vis_src1(TMP6, REF_6);
  626. vis_src1(TMP10, REF_S2);
  627. vis_src1(TMP14, REF_S6);
  628. }
  629. vis_pmerge(ZERO, REF_0, TMP0);
  630. vis_mul8x16au(REF_0_1, CONST_256, TMP2);
  631. vis_pmerge(ZERO, REF_2, TMP4);
  632. vis_mul8x16au(REF_2_1, CONST_256, TMP6);
  633. vis_padd16(TMP0, CONST_3, TMP0);
  634. vis_mul8x16al(DST_0, CONST_512, TMP16);
  635. vis_padd16(TMP2, CONST_3, TMP2);
  636. vis_mul8x16al(DST_1, CONST_512, TMP18);
  637. vis_padd16(TMP0, TMP4, TMP0);
  638. vis_mul8x16au(REF_4, CONST_256, TMP8);
  639. vis_padd16(TMP2, TMP6, TMP2);
  640. vis_mul8x16au(REF_4_1, CONST_256, TMP10);
  641. vis_padd16(TMP0, TMP16, TMP0);
  642. vis_mul8x16au(REF_6, CONST_256, TMP12);
  643. vis_padd16(TMP2, TMP18, TMP2);
  644. vis_mul8x16au(REF_6_1, CONST_256, TMP14);
  645. vis_padd16(TMP8, CONST_3, TMP8);
  646. vis_mul8x16al(DST_2, CONST_512, TMP16);
  647. vis_padd16(TMP8, TMP12, TMP8);
  648. vis_mul8x16al(DST_3, CONST_512, TMP18);
  649. vis_padd16(TMP10, TMP14, TMP10);
  650. vis_pack16(TMP0, DST_0);
  651. vis_pack16(TMP2, DST_1);
  652. vis_st64(DST_0, dest[0]);
  653. dest += stride;
  654. vis_padd16(TMP10, CONST_3, TMP10);
  655. vis_ld64_2(dest, stride, DST_0);
  656. vis_padd16(TMP8, TMP16, TMP8);
  657. vis_ld64_2(dest, stride_times_2, TMP4/*DST_2*/);
  658. vis_padd16(TMP10, TMP18, TMP10);
  659. vis_pack16(TMP8, DST_2);
  660. vis_pack16(TMP10, DST_3);
  661. vis_st64(DST_2, dest[0]);
  662. dest += stride;
  663. vis_mul8x16au(REF_S0_1, CONST_256, TMP2);
  664. vis_pmerge(ZERO, REF_S0, TMP0);
  665. vis_pmerge(ZERO, REF_S2, TMP24);
  666. vis_mul8x16au(REF_S2_1, CONST_256, TMP6);
  667. vis_padd16(TMP0, CONST_3, TMP0);
  668. vis_mul8x16au(REF_S4, CONST_256, TMP8);
  669. vis_padd16(TMP2, CONST_3, TMP2);
  670. vis_mul8x16au(REF_S4_1, CONST_256, TMP10);
  671. vis_padd16(TMP0, TMP24, TMP0);
  672. vis_mul8x16au(REF_S6, CONST_256, TMP12);
  673. vis_padd16(TMP2, TMP6, TMP2);
  674. vis_mul8x16au(REF_S6_1, CONST_256, TMP14);
  675. vis_padd16(TMP8, CONST_3, TMP8);
  676. vis_mul8x16al(DST_0, CONST_512, TMP16);
  677. vis_padd16(TMP10, CONST_3, TMP10);
  678. vis_mul8x16al(DST_1, CONST_512, TMP18);
  679. vis_padd16(TMP8, TMP12, TMP8);
  680. vis_mul8x16al(TMP4/*DST_2*/, CONST_512, TMP20);
  681. vis_mul8x16al(TMP5/*DST_3*/, CONST_512, TMP22);
  682. vis_padd16(TMP0, TMP16, TMP0);
  683. vis_padd16(TMP2, TMP18, TMP2);
  684. vis_pack16(TMP0, DST_0);
  685. vis_padd16(TMP10, TMP14, TMP10);
  686. vis_pack16(TMP2, DST_1);
  687. vis_st64(DST_0, dest[0]);
  688. dest += stride;
  689. vis_padd16(TMP8, TMP20, TMP8);
  690. vis_padd16(TMP10, TMP22, TMP10);
  691. vis_pack16(TMP8, DST_2);
  692. vis_pack16(TMP10, DST_3);
  693. vis_st64(DST_2, dest[0]);
  694. dest += stride;
  695. } while (--height);
  696. }
  697. static void MC_put_y_16_vis (uint8_t * dest, const uint8_t * ref,
  698. const ptrdiff_t stride, int height)
  699. {
  700. ref = vis_alignaddr(ref);
  701. vis_ld64(ref[0], TMP0);
  702. vis_ld64_2(ref, 8, TMP2);
  703. vis_ld64_2(ref, 16, TMP4);
  704. ref += stride;
  705. vis_ld64(ref[0], TMP6);
  706. vis_faligndata(TMP0, TMP2, REF_0);
  707. vis_ld64_2(ref, 8, TMP8);
  708. vis_faligndata(TMP2, TMP4, REF_4);
  709. vis_ld64_2(ref, 16, TMP10);
  710. ref += stride;
  711. vis_ld64(constants_fe[0], MASK_fe);
  712. vis_faligndata(TMP6, TMP8, REF_2);
  713. vis_ld64(constants_7f[0], MASK_7f);
  714. vis_faligndata(TMP8, TMP10, REF_6);
  715. vis_ld64(constants128[0], CONST_128);
  716. height = (height >> 1) - 1;
  717. do { /* 24 cycles */
  718. vis_ld64(ref[0], TMP0);
  719. vis_xor(REF_0, REF_2, TMP12);
  720. vis_ld64_2(ref, 8, TMP2);
  721. vis_xor(REF_4, REF_6, TMP16);
  722. vis_ld64_2(ref, 16, TMP4);
  723. ref += stride;
  724. vis_or(REF_0, REF_2, TMP14);
  725. vis_ld64(ref[0], TMP6);
  726. vis_or(REF_4, REF_6, TMP18);
  727. vis_ld64_2(ref, 8, TMP8);
  728. vis_faligndata(TMP0, TMP2, REF_0);
  729. vis_ld64_2(ref, 16, TMP10);
  730. ref += stride;
  731. vis_faligndata(TMP2, TMP4, REF_4);
  732. vis_and(TMP12, MASK_fe, TMP12);
  733. vis_and(TMP16, MASK_fe, TMP16);
  734. vis_mul8x16(CONST_128, TMP12, TMP12);
  735. vis_mul8x16(CONST_128, TMP16, TMP16);
  736. vis_xor(REF_0, REF_2, TMP0);
  737. vis_xor(REF_4, REF_6, TMP2);
  738. vis_or(REF_0, REF_2, TMP20);
  739. vis_and(TMP12, MASK_7f, TMP12);
  740. vis_and(TMP16, MASK_7f, TMP16);
  741. vis_psub16(TMP14, TMP12, TMP12);
  742. vis_st64(TMP12, dest[0]);
  743. vis_psub16(TMP18, TMP16, TMP16);
  744. vis_st64_2(TMP16, dest, 8);
  745. dest += stride;
  746. vis_or(REF_4, REF_6, TMP18);
  747. vis_and(TMP0, MASK_fe, TMP0);
  748. vis_and(TMP2, MASK_fe, TMP2);
  749. vis_mul8x16(CONST_128, TMP0, TMP0);
  750. vis_faligndata(TMP6, TMP8, REF_2);
  751. vis_mul8x16(CONST_128, TMP2, TMP2);
  752. vis_faligndata(TMP8, TMP10, REF_6);
  753. vis_and(TMP0, MASK_7f, TMP0);
  754. vis_and(TMP2, MASK_7f, TMP2);
  755. vis_psub16(TMP20, TMP0, TMP0);
  756. vis_st64(TMP0, dest[0]);
  757. vis_psub16(TMP18, TMP2, TMP2);
  758. vis_st64_2(TMP2, dest, 8);
  759. dest += stride;
  760. } while (--height);
  761. vis_ld64(ref[0], TMP0);
  762. vis_xor(REF_0, REF_2, TMP12);
  763. vis_ld64_2(ref, 8, TMP2);
  764. vis_xor(REF_4, REF_6, TMP16);
  765. vis_ld64_2(ref, 16, TMP4);
  766. vis_or(REF_0, REF_2, TMP14);
  767. vis_or(REF_4, REF_6, TMP18);
  768. vis_faligndata(TMP0, TMP2, REF_0);
  769. vis_faligndata(TMP2, TMP4, REF_4);
  770. vis_and(TMP12, MASK_fe, TMP12);
  771. vis_and(TMP16, MASK_fe, TMP16);
  772. vis_mul8x16(CONST_128, TMP12, TMP12);
  773. vis_mul8x16(CONST_128, TMP16, TMP16);
  774. vis_xor(REF_0, REF_2, TMP0);
  775. vis_xor(REF_4, REF_6, TMP2);
  776. vis_or(REF_0, REF_2, TMP20);
  777. vis_and(TMP12, MASK_7f, TMP12);
  778. vis_and(TMP16, MASK_7f, TMP16);
  779. vis_psub16(TMP14, TMP12, TMP12);
  780. vis_st64(TMP12, dest[0]);
  781. vis_psub16(TMP18, TMP16, TMP16);
  782. vis_st64_2(TMP16, dest, 8);
  783. dest += stride;
  784. vis_or(REF_4, REF_6, TMP18);
  785. vis_and(TMP0, MASK_fe, TMP0);
  786. vis_and(TMP2, MASK_fe, TMP2);
  787. vis_mul8x16(CONST_128, TMP0, TMP0);
  788. vis_mul8x16(CONST_128, TMP2, TMP2);
  789. vis_and(TMP0, MASK_7f, TMP0);
  790. vis_and(TMP2, MASK_7f, TMP2);
  791. vis_psub16(TMP20, TMP0, TMP0);
  792. vis_st64(TMP0, dest[0]);
  793. vis_psub16(TMP18, TMP2, TMP2);
  794. vis_st64_2(TMP2, dest, 8);
  795. }
  796. static void MC_put_y_8_vis (uint8_t * dest, const uint8_t * ref,
  797. const ptrdiff_t stride, int height)
  798. {
  799. ref = vis_alignaddr(ref);
  800. vis_ld64(ref[0], TMP0);
  801. vis_ld64_2(ref, 8, TMP2);
  802. ref += stride;
  803. vis_ld64(ref[0], TMP4);
  804. vis_ld64_2(ref, 8, TMP6);
  805. ref += stride;
  806. vis_ld64(constants_fe[0], MASK_fe);
  807. vis_faligndata(TMP0, TMP2, REF_0);
  808. vis_ld64(constants_7f[0], MASK_7f);
  809. vis_faligndata(TMP4, TMP6, REF_2);
  810. vis_ld64(constants128[0], CONST_128);
  811. height = (height >> 1) - 1;
  812. do { /* 12 cycles */
  813. vis_ld64(ref[0], TMP0);
  814. vis_xor(REF_0, REF_2, TMP4);
  815. vis_ld64_2(ref, 8, TMP2);
  816. ref += stride;
  817. vis_and(TMP4, MASK_fe, TMP4);
  818. vis_or(REF_0, REF_2, TMP6);
  819. vis_mul8x16(CONST_128, TMP4, TMP4);
  820. vis_faligndata(TMP0, TMP2, REF_0);
  821. vis_ld64(ref[0], TMP0);
  822. vis_ld64_2(ref, 8, TMP2);
  823. ref += stride;
  824. vis_xor(REF_0, REF_2, TMP12);
  825. vis_and(TMP4, MASK_7f, TMP4);
  826. vis_and(TMP12, MASK_fe, TMP12);
  827. vis_mul8x16(CONST_128, TMP12, TMP12);
  828. vis_or(REF_0, REF_2, TMP14);
  829. vis_psub16(TMP6, TMP4, DST_0);
  830. vis_st64(DST_0, dest[0]);
  831. dest += stride;
  832. vis_faligndata(TMP0, TMP2, REF_2);
  833. vis_and(TMP12, MASK_7f, TMP12);
  834. vis_psub16(TMP14, TMP12, DST_0);
  835. vis_st64(DST_0, dest[0]);
  836. dest += stride;
  837. } while (--height);
  838. vis_ld64(ref[0], TMP0);
  839. vis_xor(REF_0, REF_2, TMP4);
  840. vis_ld64_2(ref, 8, TMP2);
  841. vis_and(TMP4, MASK_fe, TMP4);
  842. vis_or(REF_0, REF_2, TMP6);
  843. vis_mul8x16(CONST_128, TMP4, TMP4);
  844. vis_faligndata(TMP0, TMP2, REF_0);
  845. vis_xor(REF_0, REF_2, TMP12);
  846. vis_and(TMP4, MASK_7f, TMP4);
  847. vis_and(TMP12, MASK_fe, TMP12);
  848. vis_mul8x16(CONST_128, TMP12, TMP12);
  849. vis_or(REF_0, REF_2, TMP14);
  850. vis_psub16(TMP6, TMP4, DST_0);
  851. vis_st64(DST_0, dest[0]);
  852. dest += stride;
  853. vis_and(TMP12, MASK_7f, TMP12);
  854. vis_psub16(TMP14, TMP12, DST_0);
  855. vis_st64(DST_0, dest[0]);
  856. }
  857. static void MC_avg_y_16_vis (uint8_t * dest, const uint8_t * ref,
  858. const ptrdiff_t stride, int height)
  859. {
  860. int stride_8 = stride + 8;
  861. int stride_16 = stride + 16;
  862. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  863. ref = vis_alignaddr(ref);
  864. vis_ld64(ref[ 0], TMP0);
  865. vis_fzero(ZERO);
  866. vis_ld64(ref[ 8], TMP2);
  867. vis_ld64(ref[16], TMP4);
  868. vis_ld64(constants3[0], CONST_3);
  869. vis_faligndata(TMP0, TMP2, REF_2);
  870. vis_ld64(constants256_512[0], CONST_256);
  871. vis_faligndata(TMP2, TMP4, REF_6);
  872. height >>= 1;
  873. do { /* 31 cycles */
  874. vis_ld64_2(ref, stride, TMP0);
  875. vis_pmerge(ZERO, REF_2, TMP12);
  876. vis_mul8x16au(REF_2_1, CONST_256, TMP14);
  877. vis_ld64_2(ref, stride_8, TMP2);
  878. vis_pmerge(ZERO, REF_6, TMP16);
  879. vis_mul8x16au(REF_6_1, CONST_256, TMP18);
  880. vis_ld64_2(ref, stride_16, TMP4);
  881. ref += stride;
  882. vis_ld64(dest[0], DST_0);
  883. vis_faligndata(TMP0, TMP2, REF_0);
  884. vis_ld64_2(dest, 8, DST_2);
  885. vis_faligndata(TMP2, TMP4, REF_4);
  886. vis_ld64_2(ref, stride, TMP6);
  887. vis_pmerge(ZERO, REF_0, TMP0);
  888. vis_mul8x16au(REF_0_1, CONST_256, TMP2);
  889. vis_ld64_2(ref, stride_8, TMP8);
  890. vis_pmerge(ZERO, REF_4, TMP4);
  891. vis_ld64_2(ref, stride_16, TMP10);
  892. ref += stride;
  893. vis_ld64_2(dest, stride, REF_S0/*DST_4*/);
  894. vis_faligndata(TMP6, TMP8, REF_2);
  895. vis_mul8x16au(REF_4_1, CONST_256, TMP6);
  896. vis_ld64_2(dest, stride_8, REF_S2/*DST_6*/);
  897. vis_faligndata(TMP8, TMP10, REF_6);
  898. vis_mul8x16al(DST_0, CONST_512, TMP20);
  899. vis_padd16(TMP0, CONST_3, TMP0);
  900. vis_mul8x16al(DST_1, CONST_512, TMP22);
  901. vis_padd16(TMP2, CONST_3, TMP2);
  902. vis_mul8x16al(DST_2, CONST_512, TMP24);
  903. vis_padd16(TMP4, CONST_3, TMP4);
  904. vis_mul8x16al(DST_3, CONST_512, TMP26);
  905. vis_padd16(TMP6, CONST_3, TMP6);
  906. vis_padd16(TMP12, TMP20, TMP12);
  907. vis_mul8x16al(REF_S0, CONST_512, TMP20);
  908. vis_padd16(TMP14, TMP22, TMP14);
  909. vis_mul8x16al(REF_S0_1, CONST_512, TMP22);
  910. vis_padd16(TMP16, TMP24, TMP16);
  911. vis_mul8x16al(REF_S2, CONST_512, TMP24);
  912. vis_padd16(TMP18, TMP26, TMP18);
  913. vis_mul8x16al(REF_S2_1, CONST_512, TMP26);
  914. vis_padd16(TMP12, TMP0, TMP12);
  915. vis_mul8x16au(REF_2, CONST_256, TMP28);
  916. vis_padd16(TMP14, TMP2, TMP14);
  917. vis_mul8x16au(REF_2_1, CONST_256, TMP30);
  918. vis_padd16(TMP16, TMP4, TMP16);
  919. vis_mul8x16au(REF_6, CONST_256, REF_S4);
  920. vis_padd16(TMP18, TMP6, TMP18);
  921. vis_mul8x16au(REF_6_1, CONST_256, REF_S6);
  922. vis_pack16(TMP12, DST_0);
  923. vis_padd16(TMP28, TMP0, TMP12);
  924. vis_pack16(TMP14, DST_1);
  925. vis_st64(DST_0, dest[0]);
  926. vis_padd16(TMP30, TMP2, TMP14);
  927. vis_pack16(TMP16, DST_2);
  928. vis_padd16(REF_S4, TMP4, TMP16);
  929. vis_pack16(TMP18, DST_3);
  930. vis_st64_2(DST_2, dest, 8);
  931. dest += stride;
  932. vis_padd16(REF_S6, TMP6, TMP18);
  933. vis_padd16(TMP12, TMP20, TMP12);
  934. vis_padd16(TMP14, TMP22, TMP14);
  935. vis_pack16(TMP12, DST_0);
  936. vis_padd16(TMP16, TMP24, TMP16);
  937. vis_pack16(TMP14, DST_1);
  938. vis_st64(DST_0, dest[0]);
  939. vis_padd16(TMP18, TMP26, TMP18);
  940. vis_pack16(TMP16, DST_2);
  941. vis_pack16(TMP18, DST_3);
  942. vis_st64_2(DST_2, dest, 8);
  943. dest += stride;
  944. } while (--height);
  945. }
  946. static void MC_avg_y_8_vis (uint8_t * dest, const uint8_t * ref,
  947. const ptrdiff_t stride, int height)
  948. {
  949. int stride_8 = stride + 8;
  950. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  951. ref = vis_alignaddr(ref);
  952. vis_ld64(ref[ 0], TMP0);
  953. vis_fzero(ZERO);
  954. vis_ld64(ref[ 8], TMP2);
  955. vis_ld64(constants3[0], CONST_3);
  956. vis_faligndata(TMP0, TMP2, REF_2);
  957. vis_ld64(constants256_512[0], CONST_256);
  958. height >>= 1;
  959. do { /* 20 cycles */
  960. vis_ld64_2(ref, stride, TMP0);
  961. vis_pmerge(ZERO, REF_2, TMP8);
  962. vis_mul8x16au(REF_2_1, CONST_256, TMP10);
  963. vis_ld64_2(ref, stride_8, TMP2);
  964. ref += stride;
  965. vis_ld64(dest[0], DST_0);
  966. vis_ld64_2(dest, stride, DST_2);
  967. vis_faligndata(TMP0, TMP2, REF_0);
  968. vis_ld64_2(ref, stride, TMP4);
  969. vis_mul8x16al(DST_0, CONST_512, TMP16);
  970. vis_pmerge(ZERO, REF_0, TMP12);
  971. vis_ld64_2(ref, stride_8, TMP6);
  972. ref += stride;
  973. vis_mul8x16al(DST_1, CONST_512, TMP18);
  974. vis_pmerge(ZERO, REF_0_1, TMP14);
  975. vis_padd16(TMP12, CONST_3, TMP12);
  976. vis_mul8x16al(DST_2, CONST_512, TMP24);
  977. vis_padd16(TMP14, CONST_3, TMP14);
  978. vis_mul8x16al(DST_3, CONST_512, TMP26);
  979. vis_faligndata(TMP4, TMP6, REF_2);
  980. vis_padd16(TMP8, TMP12, TMP8);
  981. vis_padd16(TMP10, TMP14, TMP10);
  982. vis_mul8x16au(REF_2, CONST_256, TMP20);
  983. vis_padd16(TMP8, TMP16, TMP0);
  984. vis_mul8x16au(REF_2_1, CONST_256, TMP22);
  985. vis_padd16(TMP10, TMP18, TMP2);
  986. vis_pack16(TMP0, DST_0);
  987. vis_pack16(TMP2, DST_1);
  988. vis_st64(DST_0, dest[0]);
  989. dest += stride;
  990. vis_padd16(TMP12, TMP20, TMP12);
  991. vis_padd16(TMP14, TMP22, TMP14);
  992. vis_padd16(TMP12, TMP24, TMP0);
  993. vis_padd16(TMP14, TMP26, TMP2);
  994. vis_pack16(TMP0, DST_2);
  995. vis_pack16(TMP2, DST_3);
  996. vis_st64(DST_2, dest[0]);
  997. dest += stride;
  998. } while (--height);
  999. }
  1000. static void MC_put_xy_16_vis (uint8_t * dest, const uint8_t * ref,
  1001. const ptrdiff_t stride, int height)
  1002. {
  1003. unsigned long off = (unsigned long) ref & 0x7;
  1004. unsigned long off_plus_1 = off + 1;
  1005. int stride_8 = stride + 8;
  1006. int stride_16 = stride + 16;
  1007. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  1008. ref = vis_alignaddr(ref);
  1009. vis_ld64(ref[ 0], TMP0);
  1010. vis_fzero(ZERO);
  1011. vis_ld64(ref[ 8], TMP2);
  1012. vis_ld64(ref[16], TMP4);
  1013. vis_ld64(constants2[0], CONST_2);
  1014. vis_faligndata(TMP0, TMP2, REF_S0);
  1015. vis_ld64(constants256_512[0], CONST_256);
  1016. vis_faligndata(TMP2, TMP4, REF_S4);
  1017. if (off != 0x7) {
  1018. vis_alignaddr_g0((void *)off_plus_1);
  1019. vis_faligndata(TMP0, TMP2, REF_S2);
  1020. vis_faligndata(TMP2, TMP4, REF_S6);
  1021. } else {
  1022. vis_src1(TMP2, REF_S2);
  1023. vis_src1(TMP4, REF_S6);
  1024. }
  1025. height >>= 1;
  1026. do {
  1027. vis_ld64_2(ref, stride, TMP0);
  1028. vis_mul8x16au(REF_S0, CONST_256, TMP12);
  1029. vis_pmerge(ZERO, REF_S0_1, TMP14);
  1030. vis_alignaddr_g0((void *)off);
  1031. vis_ld64_2(ref, stride_8, TMP2);
  1032. vis_mul8x16au(REF_S2, CONST_256, TMP16);
  1033. vis_pmerge(ZERO, REF_S2_1, TMP18);
  1034. vis_ld64_2(ref, stride_16, TMP4);
  1035. ref += stride;
  1036. vis_mul8x16au(REF_S4, CONST_256, TMP20);
  1037. vis_pmerge(ZERO, REF_S4_1, TMP22);
  1038. vis_ld64_2(ref, stride, TMP6);
  1039. vis_mul8x16au(REF_S6, CONST_256, TMP24);
  1040. vis_pmerge(ZERO, REF_S6_1, TMP26);
  1041. vis_ld64_2(ref, stride_8, TMP8);
  1042. vis_faligndata(TMP0, TMP2, REF_0);
  1043. vis_ld64_2(ref, stride_16, TMP10);
  1044. ref += stride;
  1045. vis_faligndata(TMP2, TMP4, REF_4);
  1046. vis_faligndata(TMP6, TMP8, REF_S0);
  1047. vis_faligndata(TMP8, TMP10, REF_S4);
  1048. if (off != 0x7) {
  1049. vis_alignaddr_g0((void *)off_plus_1);
  1050. vis_faligndata(TMP0, TMP2, REF_2);
  1051. vis_faligndata(TMP2, TMP4, REF_6);
  1052. vis_faligndata(TMP6, TMP8, REF_S2);
  1053. vis_faligndata(TMP8, TMP10, REF_S6);
  1054. } else {
  1055. vis_src1(TMP2, REF_2);
  1056. vis_src1(TMP4, REF_6);
  1057. vis_src1(TMP8, REF_S2);
  1058. vis_src1(TMP10, REF_S6);
  1059. }
  1060. vis_mul8x16au(REF_0, CONST_256, TMP0);
  1061. vis_pmerge(ZERO, REF_0_1, TMP2);
  1062. vis_mul8x16au(REF_2, CONST_256, TMP4);
  1063. vis_pmerge(ZERO, REF_2_1, TMP6);
  1064. vis_padd16(TMP0, CONST_2, TMP8);
  1065. vis_mul8x16au(REF_4, CONST_256, TMP0);
  1066. vis_padd16(TMP2, CONST_2, TMP10);
  1067. vis_mul8x16au(REF_4_1, CONST_256, TMP2);
  1068. vis_padd16(TMP8, TMP4, TMP8);
  1069. vis_mul8x16au(REF_6, CONST_256, TMP4);
  1070. vis_padd16(TMP10, TMP6, TMP10);
  1071. vis_mul8x16au(REF_6_1, CONST_256, TMP6);
  1072. vis_padd16(TMP12, TMP8, TMP12);
  1073. vis_padd16(TMP14, TMP10, TMP14);
  1074. vis_padd16(TMP12, TMP16, TMP12);
  1075. vis_padd16(TMP14, TMP18, TMP14);
  1076. vis_pack16(TMP12, DST_0);
  1077. vis_pack16(TMP14, DST_1);
  1078. vis_st64(DST_0, dest[0]);
  1079. vis_padd16(TMP0, CONST_2, TMP12);
  1080. vis_mul8x16au(REF_S0, CONST_256, TMP0);
  1081. vis_padd16(TMP2, CONST_2, TMP14);
  1082. vis_mul8x16au(REF_S0_1, CONST_256, TMP2);
  1083. vis_padd16(TMP12, TMP4, TMP12);
  1084. vis_mul8x16au(REF_S2, CONST_256, TMP4);
  1085. vis_padd16(TMP14, TMP6, TMP14);
  1086. vis_mul8x16au(REF_S2_1, CONST_256, TMP6);
  1087. vis_padd16(TMP20, TMP12, TMP20);
  1088. vis_padd16(TMP22, TMP14, TMP22);
  1089. vis_padd16(TMP20, TMP24, TMP20);
  1090. vis_padd16(TMP22, TMP26, TMP22);
  1091. vis_pack16(TMP20, DST_2);
  1092. vis_pack16(TMP22, DST_3);
  1093. vis_st64_2(DST_2, dest, 8);
  1094. dest += stride;
  1095. vis_padd16(TMP0, TMP4, TMP24);
  1096. vis_mul8x16au(REF_S4, CONST_256, TMP0);
  1097. vis_padd16(TMP2, TMP6, TMP26);
  1098. vis_mul8x16au(REF_S4_1, CONST_256, TMP2);
  1099. vis_padd16(TMP24, TMP8, TMP24);
  1100. vis_padd16(TMP26, TMP10, TMP26);
  1101. vis_pack16(TMP24, DST_0);
  1102. vis_pack16(TMP26, DST_1);
  1103. vis_st64(DST_0, dest[0]);
  1104. vis_pmerge(ZERO, REF_S6, TMP4);
  1105. vis_pmerge(ZERO, REF_S6_1, TMP6);
  1106. vis_padd16(TMP0, TMP4, TMP0);
  1107. vis_padd16(TMP2, TMP6, TMP2);
  1108. vis_padd16(TMP0, TMP12, TMP0);
  1109. vis_padd16(TMP2, TMP14, TMP2);
  1110. vis_pack16(TMP0, DST_2);
  1111. vis_pack16(TMP2, DST_3);
  1112. vis_st64_2(DST_2, dest, 8);
  1113. dest += stride;
  1114. } while (--height);
  1115. }
  1116. static void MC_put_xy_8_vis (uint8_t * dest, const uint8_t * ref,
  1117. const ptrdiff_t stride, int height)
  1118. {
  1119. unsigned long off = (unsigned long) ref & 0x7;
  1120. unsigned long off_plus_1 = off + 1;
  1121. int stride_8 = stride + 8;
  1122. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  1123. ref = vis_alignaddr(ref);
  1124. vis_ld64(ref[ 0], TMP0);
  1125. vis_fzero(ZERO);
  1126. vis_ld64(ref[ 8], TMP2);
  1127. vis_ld64(constants2[0], CONST_2);
  1128. vis_ld64(constants256_512[0], CONST_256);
  1129. vis_faligndata(TMP0, TMP2, REF_S0);
  1130. if (off != 0x7) {
  1131. vis_alignaddr_g0((void *)off_plus_1);
  1132. vis_faligndata(TMP0, TMP2, REF_S2);
  1133. } else {
  1134. vis_src1(TMP2, REF_S2);
  1135. }
  1136. height >>= 1;
  1137. do { /* 26 cycles */
  1138. vis_ld64_2(ref, stride, TMP0);
  1139. vis_mul8x16au(REF_S0, CONST_256, TMP8);
  1140. vis_pmerge(ZERO, REF_S2, TMP12);
  1141. vis_alignaddr_g0((void *)off);
  1142. vis_ld64_2(ref, stride_8, TMP2);
  1143. ref += stride;
  1144. vis_mul8x16au(REF_S0_1, CONST_256, TMP10);
  1145. vis_pmerge(ZERO, REF_S2_1, TMP14);
  1146. vis_ld64_2(ref, stride, TMP4);
  1147. vis_ld64_2(ref, stride_8, TMP6);
  1148. ref += stride;
  1149. vis_faligndata(TMP0, TMP2, REF_S4);
  1150. vis_pmerge(ZERO, REF_S4, TMP18);
  1151. vis_pmerge(ZERO, REF_S4_1, TMP20);
  1152. vis_faligndata(TMP4, TMP6, REF_S0);
  1153. if (off != 0x7) {
  1154. vis_alignaddr_g0((void *)off_plus_1);
  1155. vis_faligndata(TMP0, TMP2, REF_S6);
  1156. vis_faligndata(TMP4, TMP6, REF_S2);
  1157. } else {
  1158. vis_src1(TMP2, REF_S6);
  1159. vis_src1(TMP6, REF_S2);
  1160. }
  1161. vis_padd16(TMP18, CONST_2, TMP18);
  1162. vis_mul8x16au(REF_S6, CONST_256, TMP22);
  1163. vis_padd16(TMP20, CONST_2, TMP20);
  1164. vis_mul8x16au(REF_S6_1, CONST_256, TMP24);
  1165. vis_mul8x16au(REF_S0, CONST_256, TMP26);
  1166. vis_pmerge(ZERO, REF_S0_1, TMP28);
  1167. vis_mul8x16au(REF_S2, CONST_256, TMP30);
  1168. vis_padd16(TMP18, TMP22, TMP18);
  1169. vis_mul8x16au(REF_S2_1, CONST_256, TMP32);
  1170. vis_padd16(TMP20, TMP24, TMP20);
  1171. vis_padd16(TMP8, TMP18, TMP8);
  1172. vis_padd16(TMP10, TMP20, TMP10);
  1173. vis_padd16(TMP8, TMP12, TMP8);
  1174. vis_padd16(TMP10, TMP14, TMP10);
  1175. vis_pack16(TMP8, DST_0);
  1176. vis_pack16(TMP10, DST_1);
  1177. vis_st64(DST_0, dest[0]);
  1178. dest += stride;
  1179. vis_padd16(TMP18, TMP26, TMP18);
  1180. vis_padd16(TMP20, TMP28, TMP20);
  1181. vis_padd16(TMP18, TMP30, TMP18);
  1182. vis_padd16(TMP20, TMP32, TMP20);
  1183. vis_pack16(TMP18, DST_2);
  1184. vis_pack16(TMP20, DST_3);
  1185. vis_st64(DST_2, dest[0]);
  1186. dest += stride;
  1187. } while (--height);
  1188. }
  1189. static void MC_avg_xy_16_vis (uint8_t * dest, const uint8_t * ref,
  1190. const ptrdiff_t stride, int height)
  1191. {
  1192. unsigned long off = (unsigned long) ref & 0x7;
  1193. unsigned long off_plus_1 = off + 1;
  1194. int stride_8 = stride + 8;
  1195. int stride_16 = stride + 16;
  1196. vis_set_gsr(4 << VIS_GSR_SCALEFACT_SHIFT);
  1197. ref = vis_alignaddr(ref);
  1198. vis_ld64(ref[ 0], TMP0);
  1199. vis_fzero(ZERO);
  1200. vis_ld64(ref[ 8], TMP2);
  1201. vis_ld64(ref[16], TMP4);
  1202. vis_ld64(constants6[0], CONST_6);
  1203. vis_faligndata(TMP0, TMP2, REF_S0);
  1204. vis_ld64(constants256_1024[0], CONST_256);
  1205. vis_faligndata(TMP2, TMP4, REF_S4);
  1206. if (off != 0x7) {
  1207. vis_alignaddr_g0((void *)off_plus_1);
  1208. vis_faligndata(TMP0, TMP2, REF_S2);
  1209. vis_faligndata(TMP2, TMP4, REF_S6);
  1210. } else {
  1211. vis_src1(TMP2, REF_S2);
  1212. vis_src1(TMP4, REF_S6);
  1213. }
  1214. height >>= 1;
  1215. do { /* 55 cycles */
  1216. vis_ld64_2(ref, stride, TMP0);
  1217. vis_mul8x16au(REF_S0, CONST_256, TMP12);
  1218. vis_pmerge(ZERO, REF_S0_1, TMP14);
  1219. vis_alignaddr_g0((void *)off);
  1220. vis_ld64_2(ref, stride_8, TMP2);
  1221. vis_mul8x16au(REF_S2, CONST_256, TMP16);
  1222. vis_pmerge(ZERO, REF_S2_1, TMP18);
  1223. vis_ld64_2(ref, stride_16, TMP4);
  1224. ref += stride;
  1225. vis_mul8x16au(REF_S4, CONST_256, TMP20);
  1226. vis_pmerge(ZERO, REF_S4_1, TMP22);
  1227. vis_ld64_2(ref, stride, TMP6);
  1228. vis_mul8x16au(REF_S6, CONST_256, TMP24);
  1229. vis_pmerge(ZERO, REF_S6_1, TMP26);
  1230. vis_ld64_2(ref, stride_8, TMP8);
  1231. vis_faligndata(TMP0, TMP2, REF_0);
  1232. vis_ld64_2(ref, stride_16, TMP10);
  1233. ref += stride;
  1234. vis_faligndata(TMP2, TMP4, REF_4);
  1235. vis_ld64(dest[0], DST_0);
  1236. vis_faligndata(TMP6, TMP8, REF_S0);
  1237. vis_ld64_2(dest, 8, DST_2);
  1238. vis_faligndata(TMP8, TMP10, REF_S4);
  1239. if (off != 0x7) {
  1240. vis_alignaddr_g0((void *)off_plus_1);
  1241. vis_faligndata(TMP0, TMP2, REF_2);
  1242. vis_faligndata(TMP2, TMP4, REF_6);
  1243. vis_faligndata(TMP6, TMP8, REF_S2);
  1244. vis_faligndata(TMP8, TMP10, REF_S6);
  1245. } else {
  1246. vis_src1(TMP2, REF_2);
  1247. vis_src1(TMP4, REF_6);
  1248. vis_src1(TMP8, REF_S2);
  1249. vis_src1(TMP10, REF_S6);
  1250. }
  1251. vis_mul8x16al(DST_0, CONST_1024, TMP30);
  1252. vis_pmerge(ZERO, REF_0, TMP0);
  1253. vis_mul8x16al(DST_1, CONST_1024, TMP32);
  1254. vis_pmerge(ZERO, REF_0_1, TMP2);
  1255. vis_mul8x16au(REF_2, CONST_256, TMP4);
  1256. vis_pmerge(ZERO, REF_2_1, TMP6);
  1257. vis_mul8x16al(DST_2, CONST_1024, REF_0);
  1258. vis_padd16(TMP0, CONST_6, TMP0);
  1259. vis_mul8x16al(DST_3, CONST_1024, REF_2);
  1260. vis_padd16(TMP2, CONST_6, TMP2);
  1261. vis_padd16(TMP0, TMP4, TMP0);
  1262. vis_mul8x16au(REF_4, CONST_256, TMP4);
  1263. vis_padd16(TMP2, TMP6, TMP2);
  1264. vis_mul8x16au(REF_4_1, CONST_256, TMP6);
  1265. vis_padd16(TMP12, TMP0, TMP12);
  1266. vis_mul8x16au(REF_6, CONST_256, TMP8);
  1267. vis_padd16(TMP14, TMP2, TMP14);
  1268. vis_mul8x16au(REF_6_1, CONST_256, TMP10);
  1269. vis_padd16(TMP12, TMP16, TMP12);
  1270. vis_mul8x16au(REF_S0, CONST_256, REF_4);
  1271. vis_padd16(TMP14, TMP18, TMP14);
  1272. vis_mul8x16au(REF_S0_1, CONST_256, REF_6);
  1273. vis_padd16(TMP12, TMP30, TMP12);
  1274. vis_padd16(TMP14, TMP32, TMP14);
  1275. vis_pack16(TMP12, DST_0);
  1276. vis_pack16(TMP14, DST_1);
  1277. vis_st64(DST_0, dest[0]);
  1278. vis_padd16(TMP4, CONST_6, TMP4);
  1279. vis_ld64_2(dest, stride, DST_0);
  1280. vis_padd16(TMP6, CONST_6, TMP6);
  1281. vis_mul8x16au(REF_S2, CONST_256, TMP12);
  1282. vis_padd16(TMP4, TMP8, TMP4);
  1283. vis_mul8x16au(REF_S2_1, CONST_256, TMP14);
  1284. vis_padd16(TMP6, TMP10, TMP6);
  1285. vis_padd16(TMP20, TMP4, TMP20);
  1286. vis_padd16(TMP22, TMP6, TMP22);
  1287. vis_padd16(TMP20, TMP24, TMP20);
  1288. vis_padd16(TMP22, TMP26, TMP22);
  1289. vis_padd16(TMP20, REF_0, TMP20);
  1290. vis_mul8x16au(REF_S4, CONST_256, REF_0);
  1291. vis_padd16(TMP22, REF_2, TMP22);
  1292. vis_pack16(TMP20, DST_2);
  1293. vis_pack16(TMP22, DST_3);
  1294. vis_st64_2(DST_2, dest, 8);
  1295. dest += stride;
  1296. vis_ld64_2(dest, 8, DST_2);
  1297. vis_mul8x16al(DST_0, CONST_1024, TMP30);
  1298. vis_pmerge(ZERO, REF_S4_1, REF_2);
  1299. vis_mul8x16al(DST_1, CONST_1024, TMP32);
  1300. vis_padd16(REF_4, TMP0, TMP8);
  1301. vis_mul8x16au(REF_S6, CONST_256, REF_4);
  1302. vis_padd16(REF_6, TMP2, TMP10);
  1303. vis_mul8x16au(REF_S6_1, CONST_256, REF_6);
  1304. vis_padd16(TMP8, TMP12, TMP8);
  1305. vis_padd16(TMP10, TMP14, TMP10);
  1306. vis_padd16(TMP8, TMP30, TMP8);
  1307. vis_padd16(TMP10, TMP32, TMP10);
  1308. vis_pack16(TMP8, DST_0);
  1309. vis_pack16(TMP10, DST_1);
  1310. vis_st64(DST_0, dest[0]);
  1311. vis_padd16(REF_0, TMP4, REF_0);
  1312. vis_mul8x16al(DST_2, CONST_1024, TMP30);
  1313. vis_padd16(REF_2, TMP6, REF_2);
  1314. vis_mul8x16al(DST_3, CONST_1024, TMP32);
  1315. vis_padd16(REF_0, REF_4, REF_0);
  1316. vis_padd16(REF_2, REF_6, REF_2);
  1317. vis_padd16(REF_0, TMP30, REF_0);
  1318. /* stall */
  1319. vis_padd16(REF_2, TMP32, REF_2);
  1320. vis_pack16(REF_0, DST_2);
  1321. vis_pack16(REF_2, DST_3);
  1322. vis_st64_2(DST_2, dest, 8);
  1323. dest += stride;
  1324. } while (--height);
  1325. }
  1326. static void MC_avg_xy_8_vis (uint8_t * dest, const uint8_t * ref,
  1327. const ptrdiff_t stride, int height)
  1328. {
  1329. unsigned long off = (unsigned long) ref & 0x7;
  1330. unsigned long off_plus_1 = off + 1;
  1331. int stride_8 = stride + 8;
  1332. vis_set_gsr(4 << VIS_GSR_SCALEFACT_SHIFT);
  1333. ref = vis_alignaddr(ref);
  1334. vis_ld64(ref[0], TMP0);
  1335. vis_fzero(ZERO);
  1336. vis_ld64_2(ref, 8, TMP2);
  1337. vis_ld64(constants6[0], CONST_6);
  1338. vis_ld64(constants256_1024[0], CONST_256);
  1339. vis_faligndata(TMP0, TMP2, REF_S0);
  1340. if (off != 0x7) {
  1341. vis_alignaddr_g0((void *)off_plus_1);
  1342. vis_faligndata(TMP0, TMP2, REF_S2);
  1343. } else {
  1344. vis_src1(TMP2, REF_S2);
  1345. }
  1346. height >>= 1;
  1347. do { /* 31 cycles */
  1348. vis_ld64_2(ref, stride, TMP0);
  1349. vis_mul8x16au(REF_S0, CONST_256, TMP8);
  1350. vis_pmerge(ZERO, REF_S0_1, TMP10);
  1351. vis_ld64_2(ref, stride_8, TMP2);
  1352. ref += stride;
  1353. vis_mul8x16au(REF_S2, CONST_256, TMP12);
  1354. vis_pmerge(ZERO, REF_S2_1, TMP14);
  1355. vis_alignaddr_g0((void *)off);
  1356. vis_ld64_2(ref, stride, TMP4);
  1357. vis_faligndata(TMP0, TMP2, REF_S4);
  1358. vis_ld64_2(ref, stride_8, TMP6);
  1359. ref += stride;
  1360. vis_ld64(dest[0], DST_0);
  1361. vis_faligndata(TMP4, TMP6, REF_S0);
  1362. vis_ld64_2(dest, stride, DST_2);
  1363. if (off != 0x7) {
  1364. vis_alignaddr_g0((void *)off_plus_1);
  1365. vis_faligndata(TMP0, TMP2, REF_S6);
  1366. vis_faligndata(TMP4, TMP6, REF_S2);
  1367. } else {
  1368. vis_src1(TMP2, REF_S6);
  1369. vis_src1(TMP6, REF_S2);
  1370. }
  1371. vis_mul8x16al(DST_0, CONST_1024, TMP30);
  1372. vis_pmerge(ZERO, REF_S4, TMP22);
  1373. vis_mul8x16al(DST_1, CONST_1024, TMP32);
  1374. vis_pmerge(ZERO, REF_S4_1, TMP24);
  1375. vis_mul8x16au(REF_S6, CONST_256, TMP26);
  1376. vis_pmerge(ZERO, REF_S6_1, TMP28);
  1377. vis_mul8x16au(REF_S0, CONST_256, REF_S4);
  1378. vis_padd16(TMP22, CONST_6, TMP22);
  1379. vis_mul8x16au(REF_S0_1, CONST_256, REF_S6);
  1380. vis_padd16(TMP24, CONST_6, TMP24);
  1381. vis_mul8x16al(DST_2, CONST_1024, REF_0);
  1382. vis_padd16(TMP22, TMP26, TMP22);
  1383. vis_mul8x16al(DST_3, CONST_1024, REF_2);
  1384. vis_padd16(TMP24, TMP28, TMP24);
  1385. vis_mul8x16au(REF_S2, CONST_256, TMP26);
  1386. vis_padd16(TMP8, TMP22, TMP8);
  1387. vis_mul8x16au(REF_S2_1, CONST_256, TMP28);
  1388. vis_padd16(TMP10, TMP24, TMP10);
  1389. vis_padd16(TMP8, TMP12, TMP8);
  1390. vis_padd16(TMP10, TMP14, TMP10);
  1391. vis_padd16(TMP8, TMP30, TMP8);
  1392. vis_padd16(TMP10, TMP32, TMP10);
  1393. vis_pack16(TMP8, DST_0);
  1394. vis_pack16(TMP10, DST_1);
  1395. vis_st64(DST_0, dest[0]);
  1396. dest += stride;
  1397. vis_padd16(REF_S4, TMP22, TMP12);
  1398. vis_padd16(REF_S6, TMP24, TMP14);
  1399. vis_padd16(TMP12, TMP26, TMP12);
  1400. vis_padd16(TMP14, TMP28, TMP14);
  1401. vis_padd16(TMP12, REF_0, TMP12);
  1402. vis_padd16(TMP14, REF_2, TMP14);
  1403. vis_pack16(TMP12, DST_2);
  1404. vis_pack16(TMP14, DST_3);
  1405. vis_st64(DST_2, dest[0]);
  1406. dest += stride;
  1407. } while (--height);
  1408. }
  1409. /* End of rounding code */
  1410. /* Start of no rounding code */
  1411. /* The trick used in some of this file is the formula from the MMX
  1412. * motion comp code, which is:
  1413. *
  1414. * (x+y)>>1 == (x&y)+((x^y)>>1)
  1415. *
  1416. * This allows us to average 8 bytes at a time in a 64-bit FPU reg.
  1417. * We avoid overflows by masking before we do the shift, and we
  1418. * implement the shift by multiplying by 1/2 using mul8x16. So in
  1419. * VIS this is (assume 'x' is in f0, 'y' is in f2, a repeating mask
  1420. * of '0xfe' is in f4, a repeating mask of '0x7f' is in f6, and
  1421. * the value 0x80808080 is in f8):
  1422. *
  1423. * fxor f0, f2, f10
  1424. * fand f10, f4, f10
  1425. * fmul8x16 f8, f10, f10
  1426. * fand f10, f6, f10
  1427. * fand f0, f2, f12
  1428. * fpadd16 f12, f10, f10
  1429. */
  1430. static void MC_put_no_round_o_16_vis (uint8_t * dest, const uint8_t * ref,
  1431. const ptrdiff_t stride, int height)
  1432. {
  1433. ref = vis_alignaddr(ref);
  1434. do { /* 5 cycles */
  1435. vis_ld64(ref[0], TMP0);
  1436. vis_ld64_2(ref, 8, TMP2);
  1437. vis_ld64_2(ref, 16, TMP4);
  1438. ref += stride;
  1439. vis_faligndata(TMP0, TMP2, REF_0);
  1440. vis_st64(REF_0, dest[0]);
  1441. vis_faligndata(TMP2, TMP4, REF_2);
  1442. vis_st64_2(REF_2, dest, 8);
  1443. dest += stride;
  1444. } while (--height);
  1445. }
  1446. static void MC_put_no_round_o_8_vis (uint8_t * dest, const uint8_t * ref,
  1447. const ptrdiff_t stride, int height)
  1448. {
  1449. ref = vis_alignaddr(ref);
  1450. do { /* 4 cycles */
  1451. vis_ld64(ref[0], TMP0);
  1452. vis_ld64(ref[8], TMP2);
  1453. ref += stride;
  1454. /* stall */
  1455. vis_faligndata(TMP0, TMP2, REF_0);
  1456. vis_st64(REF_0, dest[0]);
  1457. dest += stride;
  1458. } while (--height);
  1459. }
  1460. static void MC_avg_no_round_o_16_vis (uint8_t * dest, const uint8_t * ref,
  1461. const ptrdiff_t stride, int height)
  1462. {
  1463. int stride_8 = stride + 8;
  1464. ref = vis_alignaddr(ref);
  1465. vis_ld64(ref[0], TMP0);
  1466. vis_ld64(ref[8], TMP2);
  1467. vis_ld64(ref[16], TMP4);
  1468. vis_ld64(dest[0], DST_0);
  1469. vis_ld64(dest[8], DST_2);
  1470. vis_ld64(constants_fe[0], MASK_fe);
  1471. vis_faligndata(TMP0, TMP2, REF_0);
  1472. vis_ld64(constants_7f[0], MASK_7f);
  1473. vis_faligndata(TMP2, TMP4, REF_2);
  1474. vis_ld64(constants128[0], CONST_128);
  1475. ref += stride;
  1476. height = (height >> 1) - 1;
  1477. do { /* 24 cycles */
  1478. vis_ld64(ref[0], TMP0);
  1479. vis_xor(DST_0, REF_0, TMP6);
  1480. vis_ld64_2(ref, 8, TMP2);
  1481. vis_and(TMP6, MASK_fe, TMP6);
  1482. vis_ld64_2(ref, 16, TMP4);
  1483. ref += stride;
  1484. vis_mul8x16(CONST_128, TMP6, TMP6);
  1485. vis_xor(DST_2, REF_2, TMP8);
  1486. vis_and(TMP8, MASK_fe, TMP8);
  1487. vis_and(DST_0, REF_0, TMP10);
  1488. vis_ld64_2(dest, stride, DST_0);
  1489. vis_mul8x16(CONST_128, TMP8, TMP8);
  1490. vis_and(DST_2, REF_2, TMP12);
  1491. vis_ld64_2(dest, stride_8, DST_2);
  1492. vis_ld64(ref[0], TMP14);
  1493. vis_and(TMP6, MASK_7f, TMP6);
  1494. vis_and(TMP8, MASK_7f, TMP8);
  1495. vis_padd16(TMP10, TMP6, TMP6);
  1496. vis_st64(TMP6, dest[0]);
  1497. vis_padd16(TMP12, TMP8, TMP8);
  1498. vis_st64_2(TMP8, dest, 8);
  1499. dest += stride;
  1500. vis_ld64_2(ref, 8, TMP16);
  1501. vis_faligndata(TMP0, TMP2, REF_0);
  1502. vis_ld64_2(ref, 16, TMP18);
  1503. vis_faligndata(TMP2, TMP4, REF_2);
  1504. ref += stride;
  1505. vis_xor(DST_0, REF_0, TMP20);
  1506. vis_and(TMP20, MASK_fe, TMP20);
  1507. vis_xor(DST_2, REF_2, TMP22);
  1508. vis_mul8x16(CONST_128, TMP20, TMP20);
  1509. vis_and(TMP22, MASK_fe, TMP22);
  1510. vis_and(DST_0, REF_0, TMP24);
  1511. vis_mul8x16(CONST_128, TMP22, TMP22);
  1512. vis_and(DST_2, REF_2, TMP26);
  1513. vis_ld64_2(dest, stride, DST_0);
  1514. vis_faligndata(TMP14, TMP16, REF_0);
  1515. vis_ld64_2(dest, stride_8, DST_2);
  1516. vis_faligndata(TMP16, TMP18, REF_2);
  1517. vis_and(TMP20, MASK_7f, TMP20);
  1518. vis_and(TMP22, MASK_7f, TMP22);
  1519. vis_padd16(TMP24, TMP20, TMP20);
  1520. vis_st64(TMP20, dest[0]);
  1521. vis_padd16(TMP26, TMP22, TMP22);
  1522. vis_st64_2(TMP22, dest, 8);
  1523. dest += stride;
  1524. } while (--height);
  1525. vis_ld64(ref[0], TMP0);
  1526. vis_xor(DST_0, REF_0, TMP6);
  1527. vis_ld64_2(ref, 8, TMP2);
  1528. vis_and(TMP6, MASK_fe, TMP6);
  1529. vis_ld64_2(ref, 16, TMP4);
  1530. vis_mul8x16(CONST_128, TMP6, TMP6);
  1531. vis_xor(DST_2, REF_2, TMP8);
  1532. vis_and(TMP8, MASK_fe, TMP8);
  1533. vis_and(DST_0, REF_0, TMP10);
  1534. vis_ld64_2(dest, stride, DST_0);
  1535. vis_mul8x16(CONST_128, TMP8, TMP8);
  1536. vis_and(DST_2, REF_2, TMP12);
  1537. vis_ld64_2(dest, stride_8, DST_2);
  1538. vis_ld64(ref[0], TMP14);
  1539. vis_and(TMP6, MASK_7f, TMP6);
  1540. vis_and(TMP8, MASK_7f, TMP8);
  1541. vis_padd16(TMP10, TMP6, TMP6);
  1542. vis_st64(TMP6, dest[0]);
  1543. vis_padd16(TMP12, TMP8, TMP8);
  1544. vis_st64_2(TMP8, dest, 8);
  1545. dest += stride;
  1546. vis_faligndata(TMP0, TMP2, REF_0);
  1547. vis_faligndata(TMP2, TMP4, REF_2);
  1548. vis_xor(DST_0, REF_0, TMP20);
  1549. vis_and(TMP20, MASK_fe, TMP20);
  1550. vis_xor(DST_2, REF_2, TMP22);
  1551. vis_mul8x16(CONST_128, TMP20, TMP20);
  1552. vis_and(TMP22, MASK_fe, TMP22);
  1553. vis_and(DST_0, REF_0, TMP24);
  1554. vis_mul8x16(CONST_128, TMP22, TMP22);
  1555. vis_and(DST_2, REF_2, TMP26);
  1556. vis_and(TMP20, MASK_7f, TMP20);
  1557. vis_and(TMP22, MASK_7f, TMP22);
  1558. vis_padd16(TMP24, TMP20, TMP20);
  1559. vis_st64(TMP20, dest[0]);
  1560. vis_padd16(TMP26, TMP22, TMP22);
  1561. vis_st64_2(TMP22, dest, 8);
  1562. }
  1563. static void MC_put_no_round_x_16_vis (uint8_t * dest, const uint8_t * ref,
  1564. const ptrdiff_t stride, int height)
  1565. {
  1566. unsigned long off = (unsigned long) ref & 0x7;
  1567. unsigned long off_plus_1 = off + 1;
  1568. ref = vis_alignaddr(ref);
  1569. vis_ld64(ref[0], TMP0);
  1570. vis_ld64_2(ref, 8, TMP2);
  1571. vis_ld64_2(ref, 16, TMP4);
  1572. vis_ld64(constants_fe[0], MASK_fe);
  1573. vis_ld64(constants_7f[0], MASK_7f);
  1574. vis_faligndata(TMP0, TMP2, REF_0);
  1575. vis_ld64(constants128[0], CONST_128);
  1576. vis_faligndata(TMP2, TMP4, REF_4);
  1577. if (off != 0x7) {
  1578. vis_alignaddr_g0((void *)off_plus_1);
  1579. vis_faligndata(TMP0, TMP2, REF_2);
  1580. vis_faligndata(TMP2, TMP4, REF_6);
  1581. } else {
  1582. vis_src1(TMP2, REF_2);
  1583. vis_src1(TMP4, REF_6);
  1584. }
  1585. ref += stride;
  1586. height = (height >> 1) - 1;
  1587. do { /* 34 cycles */
  1588. vis_ld64(ref[0], TMP0);
  1589. vis_xor(REF_0, REF_2, TMP6);
  1590. vis_ld64_2(ref, 8, TMP2);
  1591. vis_xor(REF_4, REF_6, TMP8);
  1592. vis_ld64_2(ref, 16, TMP4);
  1593. vis_and(TMP6, MASK_fe, TMP6);
  1594. ref += stride;
  1595. vis_ld64(ref[0], TMP14);
  1596. vis_mul8x16(CONST_128, TMP6, TMP6);
  1597. vis_and(TMP8, MASK_fe, TMP8);
  1598. vis_ld64_2(ref, 8, TMP16);
  1599. vis_mul8x16(CONST_128, TMP8, TMP8);
  1600. vis_and(REF_0, REF_2, TMP10);
  1601. vis_ld64_2(ref, 16, TMP18);
  1602. ref += stride;
  1603. vis_and(REF_4, REF_6, TMP12);
  1604. vis_alignaddr_g0((void *)off);
  1605. vis_faligndata(TMP0, TMP2, REF_0);
  1606. vis_faligndata(TMP2, TMP4, REF_4);
  1607. if (off != 0x7) {
  1608. vis_alignaddr_g0((void *)off_plus_1);
  1609. vis_faligndata(TMP0, TMP2, REF_2);
  1610. vis_faligndata(TMP2, TMP4, REF_6);
  1611. } else {
  1612. vis_src1(TMP2, REF_2);
  1613. vis_src1(TMP4, REF_6);
  1614. }
  1615. vis_and(TMP6, MASK_7f, TMP6);
  1616. vis_and(TMP8, MASK_7f, TMP8);
  1617. vis_padd16(TMP10, TMP6, TMP6);
  1618. vis_st64(TMP6, dest[0]);
  1619. vis_padd16(TMP12, TMP8, TMP8);
  1620. vis_st64_2(TMP8, dest, 8);
  1621. dest += stride;
  1622. vis_xor(REF_0, REF_2, TMP6);
  1623. vis_xor(REF_4, REF_6, TMP8);
  1624. vis_and(TMP6, MASK_fe, TMP6);
  1625. vis_mul8x16(CONST_128, TMP6, TMP6);
  1626. vis_and(TMP8, MASK_fe, TMP8);
  1627. vis_mul8x16(CONST_128, TMP8, TMP8);
  1628. vis_and(REF_0, REF_2, TMP10);
  1629. vis_and(REF_4, REF_6, TMP12);
  1630. vis_alignaddr_g0((void *)off);
  1631. vis_faligndata(TMP14, TMP16, REF_0);
  1632. vis_faligndata(TMP16, TMP18, REF_4);
  1633. if (off != 0x7) {
  1634. vis_alignaddr_g0((void *)off_plus_1);
  1635. vis_faligndata(TMP14, TMP16, REF_2);
  1636. vis_faligndata(TMP16, TMP18, REF_6);
  1637. } else {
  1638. vis_src1(TMP16, REF_2);
  1639. vis_src1(TMP18, REF_6);
  1640. }
  1641. vis_and(TMP6, MASK_7f, TMP6);
  1642. vis_and(TMP8, MASK_7f, TMP8);
  1643. vis_padd16(TMP10, TMP6, TMP6);
  1644. vis_st64(TMP6, dest[0]);
  1645. vis_padd16(TMP12, TMP8, TMP8);
  1646. vis_st64_2(TMP8, dest, 8);
  1647. dest += stride;
  1648. } while (--height);
  1649. vis_ld64(ref[0], TMP0);
  1650. vis_xor(REF_0, REF_2, TMP6);
  1651. vis_ld64_2(ref, 8, TMP2);
  1652. vis_xor(REF_4, REF_6, TMP8);
  1653. vis_ld64_2(ref, 16, TMP4);
  1654. vis_and(TMP6, MASK_fe, TMP6);
  1655. vis_mul8x16(CONST_128, TMP6, TMP6);
  1656. vis_and(TMP8, MASK_fe, TMP8);
  1657. vis_mul8x16(CONST_128, TMP8, TMP8);
  1658. vis_and(REF_0, REF_2, TMP10);
  1659. vis_and(REF_4, REF_6, TMP12);
  1660. vis_alignaddr_g0((void *)off);
  1661. vis_faligndata(TMP0, TMP2, REF_0);
  1662. vis_faligndata(TMP2, TMP4, REF_4);
  1663. if (off != 0x7) {
  1664. vis_alignaddr_g0((void *)off_plus_1);
  1665. vis_faligndata(TMP0, TMP2, REF_2);
  1666. vis_faligndata(TMP2, TMP4, REF_6);
  1667. } else {
  1668. vis_src1(TMP2, REF_2);
  1669. vis_src1(TMP4, REF_6);
  1670. }
  1671. vis_and(TMP6, MASK_7f, TMP6);
  1672. vis_and(TMP8, MASK_7f, TMP8);
  1673. vis_padd16(TMP10, TMP6, TMP6);
  1674. vis_st64(TMP6, dest[0]);
  1675. vis_padd16(TMP12, TMP8, TMP8);
  1676. vis_st64_2(TMP8, dest, 8);
  1677. dest += stride;
  1678. vis_xor(REF_0, REF_2, TMP6);
  1679. vis_xor(REF_4, REF_6, TMP8);
  1680. vis_and(TMP6, MASK_fe, TMP6);
  1681. vis_mul8x16(CONST_128, TMP6, TMP6);
  1682. vis_and(TMP8, MASK_fe, TMP8);
  1683. vis_mul8x16(CONST_128, TMP8, TMP8);
  1684. vis_and(REF_0, REF_2, TMP10);
  1685. vis_and(REF_4, REF_6, TMP12);
  1686. vis_and(TMP6, MASK_7f, TMP6);
  1687. vis_and(TMP8, MASK_7f, TMP8);
  1688. vis_padd16(TMP10, TMP6, TMP6);
  1689. vis_st64(TMP6, dest[0]);
  1690. vis_padd16(TMP12, TMP8, TMP8);
  1691. vis_st64_2(TMP8, dest, 8);
  1692. }
  1693. static void MC_put_no_round_x_8_vis (uint8_t * dest, const uint8_t * ref,
  1694. const ptrdiff_t stride, int height)
  1695. {
  1696. unsigned long off = (unsigned long) ref & 0x7;
  1697. unsigned long off_plus_1 = off + 1;
  1698. ref = vis_alignaddr(ref);
  1699. vis_ld64(ref[0], TMP0);
  1700. vis_ld64(ref[8], TMP2);
  1701. vis_ld64(constants_fe[0], MASK_fe);
  1702. vis_ld64(constants_7f[0], MASK_7f);
  1703. vis_ld64(constants128[0], CONST_128);
  1704. vis_faligndata(TMP0, TMP2, REF_0);
  1705. if (off != 0x7) {
  1706. vis_alignaddr_g0((void *)off_plus_1);
  1707. vis_faligndata(TMP0, TMP2, REF_2);
  1708. } else {
  1709. vis_src1(TMP2, REF_2);
  1710. }
  1711. ref += stride;
  1712. height = (height >> 1) - 1;
  1713. do { /* 20 cycles */
  1714. vis_ld64(ref[0], TMP0);
  1715. vis_xor(REF_0, REF_2, TMP4);
  1716. vis_ld64_2(ref, 8, TMP2);
  1717. vis_and(TMP4, MASK_fe, TMP4);
  1718. ref += stride;
  1719. vis_ld64(ref[0], TMP8);
  1720. vis_and(REF_0, REF_2, TMP6);
  1721. vis_mul8x16(CONST_128, TMP4, TMP4);
  1722. vis_alignaddr_g0((void *)off);
  1723. vis_ld64_2(ref, 8, TMP10);
  1724. ref += stride;
  1725. vis_faligndata(TMP0, TMP2, REF_0);
  1726. if (off != 0x7) {
  1727. vis_alignaddr_g0((void *)off_plus_1);
  1728. vis_faligndata(TMP0, TMP2, REF_2);
  1729. } else {
  1730. vis_src1(TMP2, REF_2);
  1731. }
  1732. vis_and(TMP4, MASK_7f, TMP4);
  1733. vis_padd16(TMP6, TMP4, DST_0);
  1734. vis_st64(DST_0, dest[0]);
  1735. dest += stride;
  1736. vis_xor(REF_0, REF_2, TMP12);
  1737. vis_and(TMP12, MASK_fe, TMP12);
  1738. vis_and(REF_0, REF_2, TMP14);
  1739. vis_mul8x16(CONST_128, TMP12, TMP12);
  1740. vis_alignaddr_g0((void *)off);
  1741. vis_faligndata(TMP8, TMP10, REF_0);
  1742. if (off != 0x7) {
  1743. vis_alignaddr_g0((void *)off_plus_1);
  1744. vis_faligndata(TMP8, TMP10, REF_2);
  1745. } else {
  1746. vis_src1(TMP10, REF_2);
  1747. }
  1748. vis_and(TMP12, MASK_7f, TMP12);
  1749. vis_padd16(TMP14, TMP12, DST_0);
  1750. vis_st64(DST_0, dest[0]);
  1751. dest += stride;
  1752. } while (--height);
  1753. vis_ld64(ref[0], TMP0);
  1754. vis_xor(REF_0, REF_2, TMP4);
  1755. vis_ld64_2(ref, 8, TMP2);
  1756. vis_and(TMP4, MASK_fe, TMP4);
  1757. vis_and(REF_0, REF_2, TMP6);
  1758. vis_mul8x16(CONST_128, TMP4, TMP4);
  1759. vis_alignaddr_g0((void *)off);
  1760. vis_faligndata(TMP0, TMP2, REF_0);
  1761. if (off != 0x7) {
  1762. vis_alignaddr_g0((void *)off_plus_1);
  1763. vis_faligndata(TMP0, TMP2, REF_2);
  1764. } else {
  1765. vis_src1(TMP2, REF_2);
  1766. }
  1767. vis_and(TMP4, MASK_7f, TMP4);
  1768. vis_padd16(TMP6, TMP4, DST_0);
  1769. vis_st64(DST_0, dest[0]);
  1770. dest += stride;
  1771. vis_xor(REF_0, REF_2, TMP12);
  1772. vis_and(TMP12, MASK_fe, TMP12);
  1773. vis_and(REF_0, REF_2, TMP14);
  1774. vis_mul8x16(CONST_128, TMP12, TMP12);
  1775. vis_and(TMP12, MASK_7f, TMP12);
  1776. vis_padd16(TMP14, TMP12, DST_0);
  1777. vis_st64(DST_0, dest[0]);
  1778. dest += stride;
  1779. }
  1780. static void MC_avg_no_round_x_16_vis (uint8_t * dest, const uint8_t * ref,
  1781. const ptrdiff_t stride, int height)
  1782. {
  1783. unsigned long off = (unsigned long) ref & 0x7;
  1784. unsigned long off_plus_1 = off + 1;
  1785. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  1786. vis_ld64(constants3[0], CONST_3);
  1787. vis_fzero(ZERO);
  1788. vis_ld64(constants256_512[0], CONST_256);
  1789. ref = vis_alignaddr(ref);
  1790. do { /* 26 cycles */
  1791. vis_ld64(ref[0], TMP0);
  1792. vis_ld64(ref[8], TMP2);
  1793. vis_alignaddr_g0((void *)off);
  1794. vis_ld64(ref[16], TMP4);
  1795. vis_ld64(dest[0], DST_0);
  1796. vis_faligndata(TMP0, TMP2, REF_0);
  1797. vis_ld64(dest[8], DST_2);
  1798. vis_faligndata(TMP2, TMP4, REF_4);
  1799. if (off != 0x7) {
  1800. vis_alignaddr_g0((void *)off_plus_1);
  1801. vis_faligndata(TMP0, TMP2, REF_2);
  1802. vis_faligndata(TMP2, TMP4, REF_6);
  1803. } else {
  1804. vis_src1(TMP2, REF_2);
  1805. vis_src1(TMP4, REF_6);
  1806. }
  1807. vis_mul8x16au(REF_0, CONST_256, TMP0);
  1808. vis_pmerge(ZERO, REF_2, TMP4);
  1809. vis_mul8x16au(REF_0_1, CONST_256, TMP2);
  1810. vis_pmerge(ZERO, REF_2_1, TMP6);
  1811. vis_padd16(TMP0, TMP4, TMP0);
  1812. vis_mul8x16al(DST_0, CONST_512, TMP4);
  1813. vis_padd16(TMP2, TMP6, TMP2);
  1814. vis_mul8x16al(DST_1, CONST_512, TMP6);
  1815. vis_mul8x16au(REF_6, CONST_256, TMP12);
  1816. vis_padd16(TMP0, TMP4, TMP0);
  1817. vis_mul8x16au(REF_6_1, CONST_256, TMP14);
  1818. vis_padd16(TMP2, TMP6, TMP2);
  1819. vis_mul8x16au(REF_4, CONST_256, TMP16);
  1820. vis_padd16(TMP0, CONST_3, TMP8);
  1821. vis_mul8x16au(REF_4_1, CONST_256, TMP18);
  1822. vis_padd16(TMP2, CONST_3, TMP10);
  1823. vis_pack16(TMP8, DST_0);
  1824. vis_pack16(TMP10, DST_1);
  1825. vis_padd16(TMP16, TMP12, TMP0);
  1826. vis_st64(DST_0, dest[0]);
  1827. vis_mul8x16al(DST_2, CONST_512, TMP4);
  1828. vis_padd16(TMP18, TMP14, TMP2);
  1829. vis_mul8x16al(DST_3, CONST_512, TMP6);
  1830. vis_padd16(TMP0, CONST_3, TMP0);
  1831. vis_padd16(TMP2, CONST_3, TMP2);
  1832. vis_padd16(TMP0, TMP4, TMP0);
  1833. vis_padd16(TMP2, TMP6, TMP2);
  1834. vis_pack16(TMP0, DST_2);
  1835. vis_pack16(TMP2, DST_3);
  1836. vis_st64(DST_2, dest[8]);
  1837. ref += stride;
  1838. dest += stride;
  1839. } while (--height);
  1840. }
  1841. static void MC_put_no_round_y_16_vis (uint8_t * dest, const uint8_t * ref,
  1842. const ptrdiff_t stride, int height)
  1843. {
  1844. ref = vis_alignaddr(ref);
  1845. vis_ld64(ref[0], TMP0);
  1846. vis_ld64_2(ref, 8, TMP2);
  1847. vis_ld64_2(ref, 16, TMP4);
  1848. ref += stride;
  1849. vis_ld64(ref[0], TMP6);
  1850. vis_faligndata(TMP0, TMP2, REF_0);
  1851. vis_ld64_2(ref, 8, TMP8);
  1852. vis_faligndata(TMP2, TMP4, REF_4);
  1853. vis_ld64_2(ref, 16, TMP10);
  1854. ref += stride;
  1855. vis_ld64(constants_fe[0], MASK_fe);
  1856. vis_faligndata(TMP6, TMP8, REF_2);
  1857. vis_ld64(constants_7f[0], MASK_7f);
  1858. vis_faligndata(TMP8, TMP10, REF_6);
  1859. vis_ld64(constants128[0], CONST_128);
  1860. height = (height >> 1) - 1;
  1861. do { /* 24 cycles */
  1862. vis_ld64(ref[0], TMP0);
  1863. vis_xor(REF_0, REF_2, TMP12);
  1864. vis_ld64_2(ref, 8, TMP2);
  1865. vis_xor(REF_4, REF_6, TMP16);
  1866. vis_ld64_2(ref, 16, TMP4);
  1867. ref += stride;
  1868. vis_and(REF_0, REF_2, TMP14);
  1869. vis_ld64(ref[0], TMP6);
  1870. vis_and(REF_4, REF_6, TMP18);
  1871. vis_ld64_2(ref, 8, TMP8);
  1872. vis_faligndata(TMP0, TMP2, REF_0);
  1873. vis_ld64_2(ref, 16, TMP10);
  1874. ref += stride;
  1875. vis_faligndata(TMP2, TMP4, REF_4);
  1876. vis_and(TMP12, MASK_fe, TMP12);
  1877. vis_and(TMP16, MASK_fe, TMP16);
  1878. vis_mul8x16(CONST_128, TMP12, TMP12);
  1879. vis_mul8x16(CONST_128, TMP16, TMP16);
  1880. vis_xor(REF_0, REF_2, TMP0);
  1881. vis_xor(REF_4, REF_6, TMP2);
  1882. vis_and(REF_0, REF_2, TMP20);
  1883. vis_and(TMP12, MASK_7f, TMP12);
  1884. vis_and(TMP16, MASK_7f, TMP16);
  1885. vis_padd16(TMP14, TMP12, TMP12);
  1886. vis_st64(TMP12, dest[0]);
  1887. vis_padd16(TMP18, TMP16, TMP16);
  1888. vis_st64_2(TMP16, dest, 8);
  1889. dest += stride;
  1890. vis_and(REF_4, REF_6, TMP18);
  1891. vis_and(TMP0, MASK_fe, TMP0);
  1892. vis_and(TMP2, MASK_fe, TMP2);
  1893. vis_mul8x16(CONST_128, TMP0, TMP0);
  1894. vis_faligndata(TMP6, TMP8, REF_2);
  1895. vis_mul8x16(CONST_128, TMP2, TMP2);
  1896. vis_faligndata(TMP8, TMP10, REF_6);
  1897. vis_and(TMP0, MASK_7f, TMP0);
  1898. vis_and(TMP2, MASK_7f, TMP2);
  1899. vis_padd16(TMP20, TMP0, TMP0);
  1900. vis_st64(TMP0, dest[0]);
  1901. vis_padd16(TMP18, TMP2, TMP2);
  1902. vis_st64_2(TMP2, dest, 8);
  1903. dest += stride;
  1904. } while (--height);
  1905. vis_ld64(ref[0], TMP0);
  1906. vis_xor(REF_0, REF_2, TMP12);
  1907. vis_ld64_2(ref, 8, TMP2);
  1908. vis_xor(REF_4, REF_6, TMP16);
  1909. vis_ld64_2(ref, 16, TMP4);
  1910. vis_and(REF_0, REF_2, TMP14);
  1911. vis_and(REF_4, REF_6, TMP18);
  1912. vis_faligndata(TMP0, TMP2, REF_0);
  1913. vis_faligndata(TMP2, TMP4, REF_4);
  1914. vis_and(TMP12, MASK_fe, TMP12);
  1915. vis_and(TMP16, MASK_fe, TMP16);
  1916. vis_mul8x16(CONST_128, TMP12, TMP12);
  1917. vis_mul8x16(CONST_128, TMP16, TMP16);
  1918. vis_xor(REF_0, REF_2, TMP0);
  1919. vis_xor(REF_4, REF_6, TMP2);
  1920. vis_and(REF_0, REF_2, TMP20);
  1921. vis_and(TMP12, MASK_7f, TMP12);
  1922. vis_and(TMP16, MASK_7f, TMP16);
  1923. vis_padd16(TMP14, TMP12, TMP12);
  1924. vis_st64(TMP12, dest[0]);
  1925. vis_padd16(TMP18, TMP16, TMP16);
  1926. vis_st64_2(TMP16, dest, 8);
  1927. dest += stride;
  1928. vis_and(REF_4, REF_6, TMP18);
  1929. vis_and(TMP0, MASK_fe, TMP0);
  1930. vis_and(TMP2, MASK_fe, TMP2);
  1931. vis_mul8x16(CONST_128, TMP0, TMP0);
  1932. vis_mul8x16(CONST_128, TMP2, TMP2);
  1933. vis_and(TMP0, MASK_7f, TMP0);
  1934. vis_and(TMP2, MASK_7f, TMP2);
  1935. vis_padd16(TMP20, TMP0, TMP0);
  1936. vis_st64(TMP0, dest[0]);
  1937. vis_padd16(TMP18, TMP2, TMP2);
  1938. vis_st64_2(TMP2, dest, 8);
  1939. }
  1940. static void MC_put_no_round_y_8_vis (uint8_t * dest, const uint8_t * ref,
  1941. const ptrdiff_t stride, int height)
  1942. {
  1943. ref = vis_alignaddr(ref);
  1944. vis_ld64(ref[0], TMP0);
  1945. vis_ld64_2(ref, 8, TMP2);
  1946. ref += stride;
  1947. vis_ld64(ref[0], TMP4);
  1948. vis_ld64_2(ref, 8, TMP6);
  1949. ref += stride;
  1950. vis_ld64(constants_fe[0], MASK_fe);
  1951. vis_faligndata(TMP0, TMP2, REF_0);
  1952. vis_ld64(constants_7f[0], MASK_7f);
  1953. vis_faligndata(TMP4, TMP6, REF_2);
  1954. vis_ld64(constants128[0], CONST_128);
  1955. height = (height >> 1) - 1;
  1956. do { /* 12 cycles */
  1957. vis_ld64(ref[0], TMP0);
  1958. vis_xor(REF_0, REF_2, TMP4);
  1959. vis_ld64_2(ref, 8, TMP2);
  1960. ref += stride;
  1961. vis_and(TMP4, MASK_fe, TMP4);
  1962. vis_and(REF_0, REF_2, TMP6);
  1963. vis_mul8x16(CONST_128, TMP4, TMP4);
  1964. vis_faligndata(TMP0, TMP2, REF_0);
  1965. vis_ld64(ref[0], TMP0);
  1966. vis_ld64_2(ref, 8, TMP2);
  1967. ref += stride;
  1968. vis_xor(REF_0, REF_2, TMP12);
  1969. vis_and(TMP4, MASK_7f, TMP4);
  1970. vis_and(TMP12, MASK_fe, TMP12);
  1971. vis_mul8x16(CONST_128, TMP12, TMP12);
  1972. vis_and(REF_0, REF_2, TMP14);
  1973. vis_padd16(TMP6, TMP4, DST_0);
  1974. vis_st64(DST_0, dest[0]);
  1975. dest += stride;
  1976. vis_faligndata(TMP0, TMP2, REF_2);
  1977. vis_and(TMP12, MASK_7f, TMP12);
  1978. vis_padd16(TMP14, TMP12, DST_0);
  1979. vis_st64(DST_0, dest[0]);
  1980. dest += stride;
  1981. } while (--height);
  1982. vis_ld64(ref[0], TMP0);
  1983. vis_xor(REF_0, REF_2, TMP4);
  1984. vis_ld64_2(ref, 8, TMP2);
  1985. vis_and(TMP4, MASK_fe, TMP4);
  1986. vis_and(REF_0, REF_2, TMP6);
  1987. vis_mul8x16(CONST_128, TMP4, TMP4);
  1988. vis_faligndata(TMP0, TMP2, REF_0);
  1989. vis_xor(REF_0, REF_2, TMP12);
  1990. vis_and(TMP4, MASK_7f, TMP4);
  1991. vis_and(TMP12, MASK_fe, TMP12);
  1992. vis_mul8x16(CONST_128, TMP12, TMP12);
  1993. vis_and(REF_0, REF_2, TMP14);
  1994. vis_padd16(TMP6, TMP4, DST_0);
  1995. vis_st64(DST_0, dest[0]);
  1996. dest += stride;
  1997. vis_and(TMP12, MASK_7f, TMP12);
  1998. vis_padd16(TMP14, TMP12, DST_0);
  1999. vis_st64(DST_0, dest[0]);
  2000. }
  2001. static void MC_avg_no_round_y_16_vis (uint8_t * dest, const uint8_t * ref,
  2002. const ptrdiff_t stride, int height)
  2003. {
  2004. int stride_8 = stride + 8;
  2005. int stride_16 = stride + 16;
  2006. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  2007. ref = vis_alignaddr(ref);
  2008. vis_ld64(ref[ 0], TMP0);
  2009. vis_fzero(ZERO);
  2010. vis_ld64(ref[ 8], TMP2);
  2011. vis_ld64(ref[16], TMP4);
  2012. vis_ld64(constants3[0], CONST_3);
  2013. vis_faligndata(TMP0, TMP2, REF_2);
  2014. vis_ld64(constants256_512[0], CONST_256);
  2015. vis_faligndata(TMP2, TMP4, REF_6);
  2016. height >>= 1;
  2017. do { /* 31 cycles */
  2018. vis_ld64_2(ref, stride, TMP0);
  2019. vis_pmerge(ZERO, REF_2, TMP12);
  2020. vis_mul8x16au(REF_2_1, CONST_256, TMP14);
  2021. vis_ld64_2(ref, stride_8, TMP2);
  2022. vis_pmerge(ZERO, REF_6, TMP16);
  2023. vis_mul8x16au(REF_6_1, CONST_256, TMP18);
  2024. vis_ld64_2(ref, stride_16, TMP4);
  2025. ref += stride;
  2026. vis_ld64(dest[0], DST_0);
  2027. vis_faligndata(TMP0, TMP2, REF_0);
  2028. vis_ld64_2(dest, 8, DST_2);
  2029. vis_faligndata(TMP2, TMP4, REF_4);
  2030. vis_ld64_2(ref, stride, TMP6);
  2031. vis_pmerge(ZERO, REF_0, TMP0);
  2032. vis_mul8x16au(REF_0_1, CONST_256, TMP2);
  2033. vis_ld64_2(ref, stride_8, TMP8);
  2034. vis_pmerge(ZERO, REF_4, TMP4);
  2035. vis_ld64_2(ref, stride_16, TMP10);
  2036. ref += stride;
  2037. vis_ld64_2(dest, stride, REF_S0/*DST_4*/);
  2038. vis_faligndata(TMP6, TMP8, REF_2);
  2039. vis_mul8x16au(REF_4_1, CONST_256, TMP6);
  2040. vis_ld64_2(dest, stride_8, REF_S2/*DST_6*/);
  2041. vis_faligndata(TMP8, TMP10, REF_6);
  2042. vis_mul8x16al(DST_0, CONST_512, TMP20);
  2043. vis_padd16(TMP0, CONST_3, TMP0);
  2044. vis_mul8x16al(DST_1, CONST_512, TMP22);
  2045. vis_padd16(TMP2, CONST_3, TMP2);
  2046. vis_mul8x16al(DST_2, CONST_512, TMP24);
  2047. vis_padd16(TMP4, CONST_3, TMP4);
  2048. vis_mul8x16al(DST_3, CONST_512, TMP26);
  2049. vis_padd16(TMP6, CONST_3, TMP6);
  2050. vis_padd16(TMP12, TMP20, TMP12);
  2051. vis_mul8x16al(REF_S0, CONST_512, TMP20);
  2052. vis_padd16(TMP14, TMP22, TMP14);
  2053. vis_mul8x16al(REF_S0_1, CONST_512, TMP22);
  2054. vis_padd16(TMP16, TMP24, TMP16);
  2055. vis_mul8x16al(REF_S2, CONST_512, TMP24);
  2056. vis_padd16(TMP18, TMP26, TMP18);
  2057. vis_mul8x16al(REF_S2_1, CONST_512, TMP26);
  2058. vis_padd16(TMP12, TMP0, TMP12);
  2059. vis_mul8x16au(REF_2, CONST_256, TMP28);
  2060. vis_padd16(TMP14, TMP2, TMP14);
  2061. vis_mul8x16au(REF_2_1, CONST_256, TMP30);
  2062. vis_padd16(TMP16, TMP4, TMP16);
  2063. vis_mul8x16au(REF_6, CONST_256, REF_S4);
  2064. vis_padd16(TMP18, TMP6, TMP18);
  2065. vis_mul8x16au(REF_6_1, CONST_256, REF_S6);
  2066. vis_pack16(TMP12, DST_0);
  2067. vis_padd16(TMP28, TMP0, TMP12);
  2068. vis_pack16(TMP14, DST_1);
  2069. vis_st64(DST_0, dest[0]);
  2070. vis_padd16(TMP30, TMP2, TMP14);
  2071. vis_pack16(TMP16, DST_2);
  2072. vis_padd16(REF_S4, TMP4, TMP16);
  2073. vis_pack16(TMP18, DST_3);
  2074. vis_st64_2(DST_2, dest, 8);
  2075. dest += stride;
  2076. vis_padd16(REF_S6, TMP6, TMP18);
  2077. vis_padd16(TMP12, TMP20, TMP12);
  2078. vis_padd16(TMP14, TMP22, TMP14);
  2079. vis_pack16(TMP12, DST_0);
  2080. vis_padd16(TMP16, TMP24, TMP16);
  2081. vis_pack16(TMP14, DST_1);
  2082. vis_st64(DST_0, dest[0]);
  2083. vis_padd16(TMP18, TMP26, TMP18);
  2084. vis_pack16(TMP16, DST_2);
  2085. vis_pack16(TMP18, DST_3);
  2086. vis_st64_2(DST_2, dest, 8);
  2087. dest += stride;
  2088. } while (--height);
  2089. }
  2090. static void MC_put_no_round_xy_16_vis (uint8_t * dest, const uint8_t * ref,
  2091. const ptrdiff_t stride, int height)
  2092. {
  2093. unsigned long off = (unsigned long) ref & 0x7;
  2094. unsigned long off_plus_1 = off + 1;
  2095. int stride_8 = stride + 8;
  2096. int stride_16 = stride + 16;
  2097. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  2098. ref = vis_alignaddr(ref);
  2099. vis_ld64(ref[ 0], TMP0);
  2100. vis_fzero(ZERO);
  2101. vis_ld64(ref[ 8], TMP2);
  2102. vis_ld64(ref[16], TMP4);
  2103. vis_ld64(constants1[0], CONST_1);
  2104. vis_faligndata(TMP0, TMP2, REF_S0);
  2105. vis_ld64(constants256_512[0], CONST_256);
  2106. vis_faligndata(TMP2, TMP4, REF_S4);
  2107. if (off != 0x7) {
  2108. vis_alignaddr_g0((void *)off_plus_1);
  2109. vis_faligndata(TMP0, TMP2, REF_S2);
  2110. vis_faligndata(TMP2, TMP4, REF_S6);
  2111. } else {
  2112. vis_src1(TMP2, REF_S2);
  2113. vis_src1(TMP4, REF_S6);
  2114. }
  2115. height >>= 1;
  2116. do {
  2117. vis_ld64_2(ref, stride, TMP0);
  2118. vis_mul8x16au(REF_S0, CONST_256, TMP12);
  2119. vis_pmerge(ZERO, REF_S0_1, TMP14);
  2120. vis_alignaddr_g0((void *)off);
  2121. vis_ld64_2(ref, stride_8, TMP2);
  2122. vis_mul8x16au(REF_S2, CONST_256, TMP16);
  2123. vis_pmerge(ZERO, REF_S2_1, TMP18);
  2124. vis_ld64_2(ref, stride_16, TMP4);
  2125. ref += stride;
  2126. vis_mul8x16au(REF_S4, CONST_256, TMP20);
  2127. vis_pmerge(ZERO, REF_S4_1, TMP22);
  2128. vis_ld64_2(ref, stride, TMP6);
  2129. vis_mul8x16au(REF_S6, CONST_256, TMP24);
  2130. vis_pmerge(ZERO, REF_S6_1, TMP26);
  2131. vis_ld64_2(ref, stride_8, TMP8);
  2132. vis_faligndata(TMP0, TMP2, REF_0);
  2133. vis_ld64_2(ref, stride_16, TMP10);
  2134. ref += stride;
  2135. vis_faligndata(TMP2, TMP4, REF_4);
  2136. vis_faligndata(TMP6, TMP8, REF_S0);
  2137. vis_faligndata(TMP8, TMP10, REF_S4);
  2138. if (off != 0x7) {
  2139. vis_alignaddr_g0((void *)off_plus_1);
  2140. vis_faligndata(TMP0, TMP2, REF_2);
  2141. vis_faligndata(TMP2, TMP4, REF_6);
  2142. vis_faligndata(TMP6, TMP8, REF_S2);
  2143. vis_faligndata(TMP8, TMP10, REF_S6);
  2144. } else {
  2145. vis_src1(TMP2, REF_2);
  2146. vis_src1(TMP4, REF_6);
  2147. vis_src1(TMP8, REF_S2);
  2148. vis_src1(TMP10, REF_S6);
  2149. }
  2150. vis_mul8x16au(REF_0, CONST_256, TMP0);
  2151. vis_pmerge(ZERO, REF_0_1, TMP2);
  2152. vis_mul8x16au(REF_2, CONST_256, TMP4);
  2153. vis_pmerge(ZERO, REF_2_1, TMP6);
  2154. vis_padd16(TMP0, CONST_2, TMP8);
  2155. vis_mul8x16au(REF_4, CONST_256, TMP0);
  2156. vis_padd16(TMP2, CONST_1, TMP10);
  2157. vis_mul8x16au(REF_4_1, CONST_256, TMP2);
  2158. vis_padd16(TMP8, TMP4, TMP8);
  2159. vis_mul8x16au(REF_6, CONST_256, TMP4);
  2160. vis_padd16(TMP10, TMP6, TMP10);
  2161. vis_mul8x16au(REF_6_1, CONST_256, TMP6);
  2162. vis_padd16(TMP12, TMP8, TMP12);
  2163. vis_padd16(TMP14, TMP10, TMP14);
  2164. vis_padd16(TMP12, TMP16, TMP12);
  2165. vis_padd16(TMP14, TMP18, TMP14);
  2166. vis_pack16(TMP12, DST_0);
  2167. vis_pack16(TMP14, DST_1);
  2168. vis_st64(DST_0, dest[0]);
  2169. vis_padd16(TMP0, CONST_1, TMP12);
  2170. vis_mul8x16au(REF_S0, CONST_256, TMP0);
  2171. vis_padd16(TMP2, CONST_1, TMP14);
  2172. vis_mul8x16au(REF_S0_1, CONST_256, TMP2);
  2173. vis_padd16(TMP12, TMP4, TMP12);
  2174. vis_mul8x16au(REF_S2, CONST_256, TMP4);
  2175. vis_padd16(TMP14, TMP6, TMP14);
  2176. vis_mul8x16au(REF_S2_1, CONST_256, TMP6);
  2177. vis_padd16(TMP20, TMP12, TMP20);
  2178. vis_padd16(TMP22, TMP14, TMP22);
  2179. vis_padd16(TMP20, TMP24, TMP20);
  2180. vis_padd16(TMP22, TMP26, TMP22);
  2181. vis_pack16(TMP20, DST_2);
  2182. vis_pack16(TMP22, DST_3);
  2183. vis_st64_2(DST_2, dest, 8);
  2184. dest += stride;
  2185. vis_padd16(TMP0, TMP4, TMP24);
  2186. vis_mul8x16au(REF_S4, CONST_256, TMP0);
  2187. vis_padd16(TMP2, TMP6, TMP26);
  2188. vis_mul8x16au(REF_S4_1, CONST_256, TMP2);
  2189. vis_padd16(TMP24, TMP8, TMP24);
  2190. vis_padd16(TMP26, TMP10, TMP26);
  2191. vis_pack16(TMP24, DST_0);
  2192. vis_pack16(TMP26, DST_1);
  2193. vis_st64(DST_0, dest[0]);
  2194. vis_pmerge(ZERO, REF_S6, TMP4);
  2195. vis_pmerge(ZERO, REF_S6_1, TMP6);
  2196. vis_padd16(TMP0, TMP4, TMP0);
  2197. vis_padd16(TMP2, TMP6, TMP2);
  2198. vis_padd16(TMP0, TMP12, TMP0);
  2199. vis_padd16(TMP2, TMP14, TMP2);
  2200. vis_pack16(TMP0, DST_2);
  2201. vis_pack16(TMP2, DST_3);
  2202. vis_st64_2(DST_2, dest, 8);
  2203. dest += stride;
  2204. } while (--height);
  2205. }
  2206. static void MC_put_no_round_xy_8_vis (uint8_t * dest, const uint8_t * ref,
  2207. const ptrdiff_t stride, int height)
  2208. {
  2209. unsigned long off = (unsigned long) ref & 0x7;
  2210. unsigned long off_plus_1 = off + 1;
  2211. int stride_8 = stride + 8;
  2212. vis_set_gsr(5 << VIS_GSR_SCALEFACT_SHIFT);
  2213. ref = vis_alignaddr(ref);
  2214. vis_ld64(ref[ 0], TMP0);
  2215. vis_fzero(ZERO);
  2216. vis_ld64(ref[ 8], TMP2);
  2217. vis_ld64(constants1[0], CONST_1);
  2218. vis_ld64(constants256_512[0], CONST_256);
  2219. vis_faligndata(TMP0, TMP2, REF_S0);
  2220. if (off != 0x7) {
  2221. vis_alignaddr_g0((void *)off_plus_1);
  2222. vis_faligndata(TMP0, TMP2, REF_S2);
  2223. } else {
  2224. vis_src1(TMP2, REF_S2);
  2225. }
  2226. height >>= 1;
  2227. do { /* 26 cycles */
  2228. vis_ld64_2(ref, stride, TMP0);
  2229. vis_mul8x16au(REF_S0, CONST_256, TMP8);
  2230. vis_pmerge(ZERO, REF_S2, TMP12);
  2231. vis_alignaddr_g0((void *)off);
  2232. vis_ld64_2(ref, stride_8, TMP2);
  2233. ref += stride;
  2234. vis_mul8x16au(REF_S0_1, CONST_256, TMP10);
  2235. vis_pmerge(ZERO, REF_S2_1, TMP14);
  2236. vis_ld64_2(ref, stride, TMP4);
  2237. vis_ld64_2(ref, stride_8, TMP6);
  2238. ref += stride;
  2239. vis_faligndata(TMP0, TMP2, REF_S4);
  2240. vis_pmerge(ZERO, REF_S4, TMP18);
  2241. vis_pmerge(ZERO, REF_S4_1, TMP20);
  2242. vis_faligndata(TMP4, TMP6, REF_S0);
  2243. if (off != 0x7) {
  2244. vis_alignaddr_g0((void *)off_plus_1);
  2245. vis_faligndata(TMP0, TMP2, REF_S6);
  2246. vis_faligndata(TMP4, TMP6, REF_S2);
  2247. } else {
  2248. vis_src1(TMP2, REF_S6);
  2249. vis_src1(TMP6, REF_S2);
  2250. }
  2251. vis_padd16(TMP18, CONST_1, TMP18);
  2252. vis_mul8x16au(REF_S6, CONST_256, TMP22);
  2253. vis_padd16(TMP20, CONST_1, TMP20);
  2254. vis_mul8x16au(REF_S6_1, CONST_256, TMP24);
  2255. vis_mul8x16au(REF_S0, CONST_256, TMP26);
  2256. vis_pmerge(ZERO, REF_S0_1, TMP28);
  2257. vis_mul8x16au(REF_S2, CONST_256, TMP30);
  2258. vis_padd16(TMP18, TMP22, TMP18);
  2259. vis_mul8x16au(REF_S2_1, CONST_256, TMP32);
  2260. vis_padd16(TMP20, TMP24, TMP20);
  2261. vis_padd16(TMP8, TMP18, TMP8);
  2262. vis_padd16(TMP10, TMP20, TMP10);
  2263. vis_padd16(TMP8, TMP12, TMP8);
  2264. vis_padd16(TMP10, TMP14, TMP10);
  2265. vis_pack16(TMP8, DST_0);
  2266. vis_pack16(TMP10, DST_1);
  2267. vis_st64(DST_0, dest[0]);
  2268. dest += stride;
  2269. vis_padd16(TMP18, TMP26, TMP18);
  2270. vis_padd16(TMP20, TMP28, TMP20);
  2271. vis_padd16(TMP18, TMP30, TMP18);
  2272. vis_padd16(TMP20, TMP32, TMP20);
  2273. vis_pack16(TMP18, DST_2);
  2274. vis_pack16(TMP20, DST_3);
  2275. vis_st64(DST_2, dest[0]);
  2276. dest += stride;
  2277. } while (--height);
  2278. }
  2279. static void MC_avg_no_round_xy_16_vis (uint8_t * dest, const uint8_t * ref,
  2280. const ptrdiff_t stride, int height)
  2281. {
  2282. unsigned long off = (unsigned long) ref & 0x7;
  2283. unsigned long off_plus_1 = off + 1;
  2284. int stride_8 = stride + 8;
  2285. int stride_16 = stride + 16;
  2286. vis_set_gsr(4 << VIS_GSR_SCALEFACT_SHIFT);
  2287. ref = vis_alignaddr(ref);
  2288. vis_ld64(ref[ 0], TMP0);
  2289. vis_fzero(ZERO);
  2290. vis_ld64(ref[ 8], TMP2);
  2291. vis_ld64(ref[16], TMP4);
  2292. vis_ld64(constants6[0], CONST_6);
  2293. vis_faligndata(TMP0, TMP2, REF_S0);
  2294. vis_ld64(constants256_1024[0], CONST_256);
  2295. vis_faligndata(TMP2, TMP4, REF_S4);
  2296. if (off != 0x7) {
  2297. vis_alignaddr_g0((void *)off_plus_1);
  2298. vis_faligndata(TMP0, TMP2, REF_S2);
  2299. vis_faligndata(TMP2, TMP4, REF_S6);
  2300. } else {
  2301. vis_src1(TMP2, REF_S2);
  2302. vis_src1(TMP4, REF_S6);
  2303. }
  2304. height >>= 1;
  2305. do { /* 55 cycles */
  2306. vis_ld64_2(ref, stride, TMP0);
  2307. vis_mul8x16au(REF_S0, CONST_256, TMP12);
  2308. vis_pmerge(ZERO, REF_S0_1, TMP14);
  2309. vis_alignaddr_g0((void *)off);
  2310. vis_ld64_2(ref, stride_8, TMP2);
  2311. vis_mul8x16au(REF_S2, CONST_256, TMP16);
  2312. vis_pmerge(ZERO, REF_S2_1, TMP18);
  2313. vis_ld64_2(ref, stride_16, TMP4);
  2314. ref += stride;
  2315. vis_mul8x16au(REF_S4, CONST_256, TMP20);
  2316. vis_pmerge(ZERO, REF_S4_1, TMP22);
  2317. vis_ld64_2(ref, stride, TMP6);
  2318. vis_mul8x16au(REF_S6, CONST_256, TMP24);
  2319. vis_pmerge(ZERO, REF_S6_1, TMP26);
  2320. vis_ld64_2(ref, stride_8, TMP8);
  2321. vis_faligndata(TMP0, TMP2, REF_0);
  2322. vis_ld64_2(ref, stride_16, TMP10);
  2323. ref += stride;
  2324. vis_faligndata(TMP2, TMP4, REF_4);
  2325. vis_ld64(dest[0], DST_0);
  2326. vis_faligndata(TMP6, TMP8, REF_S0);
  2327. vis_ld64_2(dest, 8, DST_2);
  2328. vis_faligndata(TMP8, TMP10, REF_S4);
  2329. if (off != 0x7) {
  2330. vis_alignaddr_g0((void *)off_plus_1);
  2331. vis_faligndata(TMP0, TMP2, REF_2);
  2332. vis_faligndata(TMP2, TMP4, REF_6);
  2333. vis_faligndata(TMP6, TMP8, REF_S2);
  2334. vis_faligndata(TMP8, TMP10, REF_S6);
  2335. } else {
  2336. vis_src1(TMP2, REF_2);
  2337. vis_src1(TMP4, REF_6);
  2338. vis_src1(TMP8, REF_S2);
  2339. vis_src1(TMP10, REF_S6);
  2340. }
  2341. vis_mul8x16al(DST_0, CONST_1024, TMP30);
  2342. vis_pmerge(ZERO, REF_0, TMP0);
  2343. vis_mul8x16al(DST_1, CONST_1024, TMP32);
  2344. vis_pmerge(ZERO, REF_0_1, TMP2);
  2345. vis_mul8x16au(REF_2, CONST_256, TMP4);
  2346. vis_pmerge(ZERO, REF_2_1, TMP6);
  2347. vis_mul8x16al(DST_2, CONST_1024, REF_0);
  2348. vis_padd16(TMP0, CONST_6, TMP0);
  2349. vis_mul8x16al(DST_3, CONST_1024, REF_2);
  2350. vis_padd16(TMP2, CONST_6, TMP2);
  2351. vis_padd16(TMP0, TMP4, TMP0);
  2352. vis_mul8x16au(REF_4, CONST_256, TMP4);
  2353. vis_padd16(TMP2, TMP6, TMP2);
  2354. vis_mul8x16au(REF_4_1, CONST_256, TMP6);
  2355. vis_padd16(TMP12, TMP0, TMP12);
  2356. vis_mul8x16au(REF_6, CONST_256, TMP8);
  2357. vis_padd16(TMP14, TMP2, TMP14);
  2358. vis_mul8x16au(REF_6_1, CONST_256, TMP10);
  2359. vis_padd16(TMP12, TMP16, TMP12);
  2360. vis_mul8x16au(REF_S0, CONST_256, REF_4);
  2361. vis_padd16(TMP14, TMP18, TMP14);
  2362. vis_mul8x16au(REF_S0_1, CONST_256, REF_6);
  2363. vis_padd16(TMP12, TMP30, TMP12);
  2364. vis_padd16(TMP14, TMP32, TMP14);
  2365. vis_pack16(TMP12, DST_0);
  2366. vis_pack16(TMP14, DST_1);
  2367. vis_st64(DST_0, dest[0]);
  2368. vis_padd16(TMP4, CONST_6, TMP4);
  2369. vis_ld64_2(dest, stride, DST_0);
  2370. vis_padd16(TMP6, CONST_6, TMP6);
  2371. vis_mul8x16au(REF_S2, CONST_256, TMP12);
  2372. vis_padd16(TMP4, TMP8, TMP4);
  2373. vis_mul8x16au(REF_S2_1, CONST_256, TMP14);
  2374. vis_padd16(TMP6, TMP10, TMP6);
  2375. vis_padd16(TMP20, TMP4, TMP20);
  2376. vis_padd16(TMP22, TMP6, TMP22);
  2377. vis_padd16(TMP20, TMP24, TMP20);
  2378. vis_padd16(TMP22, TMP26, TMP22);
  2379. vis_padd16(TMP20, REF_0, TMP20);
  2380. vis_mul8x16au(REF_S4, CONST_256, REF_0);
  2381. vis_padd16(TMP22, REF_2, TMP22);
  2382. vis_pack16(TMP20, DST_2);
  2383. vis_pack16(TMP22, DST_3);
  2384. vis_st64_2(DST_2, dest, 8);
  2385. dest += stride;
  2386. vis_ld64_2(dest, 8, DST_2);
  2387. vis_mul8x16al(DST_0, CONST_1024, TMP30);
  2388. vis_pmerge(ZERO, REF_S4_1, REF_2);
  2389. vis_mul8x16al(DST_1, CONST_1024, TMP32);
  2390. vis_padd16(REF_4, TMP0, TMP8);
  2391. vis_mul8x16au(REF_S6, CONST_256, REF_4);
  2392. vis_padd16(REF_6, TMP2, TMP10);
  2393. vis_mul8x16au(REF_S6_1, CONST_256, REF_6);
  2394. vis_padd16(TMP8, TMP12, TMP8);
  2395. vis_padd16(TMP10, TMP14, TMP10);
  2396. vis_padd16(TMP8, TMP30, TMP8);
  2397. vis_padd16(TMP10, TMP32, TMP10);
  2398. vis_pack16(TMP8, DST_0);
  2399. vis_pack16(TMP10, DST_1);
  2400. vis_st64(DST_0, dest[0]);
  2401. vis_padd16(REF_0, TMP4, REF_0);
  2402. vis_mul8x16al(DST_2, CONST_1024, TMP30);
  2403. vis_padd16(REF_2, TMP6, REF_2);
  2404. vis_mul8x16al(DST_3, CONST_1024, TMP32);
  2405. vis_padd16(REF_0, REF_4, REF_0);
  2406. vis_padd16(REF_2, REF_6, REF_2);
  2407. vis_padd16(REF_0, TMP30, REF_0);
  2408. /* stall */
  2409. vis_padd16(REF_2, TMP32, REF_2);
  2410. vis_pack16(REF_0, DST_2);
  2411. vis_pack16(REF_2, DST_3);
  2412. vis_st64_2(DST_2, dest, 8);
  2413. dest += stride;
  2414. } while (--height);
  2415. }
  2416. /* End of no rounding code */
  2417. #define ACCEL_SPARC_VIS 1
  2418. #define ACCEL_SPARC_VIS2 2
  2419. static int vis_level(void)
  2420. {
  2421. int accel = 0;
  2422. accel |= ACCEL_SPARC_VIS;
  2423. accel |= ACCEL_SPARC_VIS2;
  2424. return accel;
  2425. }
  2426. /* libavcodec initialization code */
  2427. void ff_hpeldsp_init_vis(HpelDSPContext* c, int flags)
  2428. {
  2429. /* VIS-specific optimizations */
  2430. int accel = vis_level ();
  2431. if (accel & ACCEL_SPARC_VIS) {
  2432. c->put_pixels_tab[0][0] = MC_put_o_16_vis;
  2433. c->put_pixels_tab[0][1] = MC_put_x_16_vis;
  2434. c->put_pixels_tab[0][2] = MC_put_y_16_vis;
  2435. c->put_pixels_tab[0][3] = MC_put_xy_16_vis;
  2436. c->put_pixels_tab[1][0] = MC_put_o_8_vis;
  2437. c->put_pixels_tab[1][1] = MC_put_x_8_vis;
  2438. c->put_pixels_tab[1][2] = MC_put_y_8_vis;
  2439. c->put_pixels_tab[1][3] = MC_put_xy_8_vis;
  2440. c->avg_pixels_tab[0][0] = MC_avg_o_16_vis;
  2441. c->avg_pixels_tab[0][1] = MC_avg_x_16_vis;
  2442. c->avg_pixels_tab[0][2] = MC_avg_y_16_vis;
  2443. c->avg_pixels_tab[0][3] = MC_avg_xy_16_vis;
  2444. c->avg_pixels_tab[1][0] = MC_avg_o_8_vis;
  2445. c->avg_pixels_tab[1][1] = MC_avg_x_8_vis;
  2446. c->avg_pixels_tab[1][2] = MC_avg_y_8_vis;
  2447. c->avg_pixels_tab[1][3] = MC_avg_xy_8_vis;
  2448. c->put_no_rnd_pixels_tab[0][0] = MC_put_no_round_o_16_vis;
  2449. c->put_no_rnd_pixels_tab[0][1] = MC_put_no_round_x_16_vis;
  2450. c->put_no_rnd_pixels_tab[0][2] = MC_put_no_round_y_16_vis;
  2451. c->put_no_rnd_pixels_tab[0][3] = MC_put_no_round_xy_16_vis;
  2452. c->put_no_rnd_pixels_tab[1][0] = MC_put_no_round_o_8_vis;
  2453. c->put_no_rnd_pixels_tab[1][1] = MC_put_no_round_x_8_vis;
  2454. c->put_no_rnd_pixels_tab[1][2] = MC_put_no_round_y_8_vis;
  2455. c->put_no_rnd_pixels_tab[1][3] = MC_put_no_round_xy_8_vis;
  2456. c->avg_no_rnd_pixels_tab[0] = MC_avg_no_round_o_16_vis;
  2457. c->avg_no_rnd_pixels_tab[1] = MC_avg_no_round_x_16_vis;
  2458. c->avg_no_rnd_pixels_tab[2] = MC_avg_no_round_y_16_vis;
  2459. c->avg_no_rnd_pixels_tab[3] = MC_avg_no_round_xy_16_vis;
  2460. }
  2461. }