vpx_convolve_copy_msa.c 8.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249
  1. /*
  2. * Copyright (c) 2015 The WebM project authors. All Rights Reserved.
  3. *
  4. * Use of this source code is governed by a BSD-style license
  5. * that can be found in the LICENSE file in the root of the source
  6. * tree. An additional intellectual property rights grant can be found
  7. * in the file PATENTS. All contributing project authors may
  8. * be found in the AUTHORS file in the root of the source tree.
  9. */
  10. #include <string.h>
  11. #include "./vpx_dsp_rtcd.h"
  12. #include "vpx_dsp/mips/macros_msa.h"
  13. static void copy_width8_msa(const uint8_t *src, int32_t src_stride,
  14. uint8_t *dst, int32_t dst_stride, int32_t height) {
  15. int32_t cnt;
  16. uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
  17. v16u8 src0, src1, src2, src3, src4, src5, src6, src7;
  18. if (0 == height % 12) {
  19. for (cnt = (height / 12); cnt--;) {
  20. LD_UB8(src, src_stride, src0, src1, src2, src3, src4, src5, src6, src7);
  21. src += (8 * src_stride);
  22. out0 = __msa_copy_u_d((v2i64)src0, 0);
  23. out1 = __msa_copy_u_d((v2i64)src1, 0);
  24. out2 = __msa_copy_u_d((v2i64)src2, 0);
  25. out3 = __msa_copy_u_d((v2i64)src3, 0);
  26. out4 = __msa_copy_u_d((v2i64)src4, 0);
  27. out5 = __msa_copy_u_d((v2i64)src5, 0);
  28. out6 = __msa_copy_u_d((v2i64)src6, 0);
  29. out7 = __msa_copy_u_d((v2i64)src7, 0);
  30. SD4(out0, out1, out2, out3, dst, dst_stride);
  31. dst += (4 * dst_stride);
  32. SD4(out4, out5, out6, out7, dst, dst_stride);
  33. dst += (4 * dst_stride);
  34. LD_UB4(src, src_stride, src0, src1, src2, src3);
  35. src += (4 * src_stride);
  36. out0 = __msa_copy_u_d((v2i64)src0, 0);
  37. out1 = __msa_copy_u_d((v2i64)src1, 0);
  38. out2 = __msa_copy_u_d((v2i64)src2, 0);
  39. out3 = __msa_copy_u_d((v2i64)src3, 0);
  40. SD4(out0, out1, out2, out3, dst, dst_stride);
  41. dst += (4 * dst_stride);
  42. }
  43. } else if (0 == height % 8) {
  44. for (cnt = height >> 3; cnt--;) {
  45. LD_UB8(src, src_stride, src0, src1, src2, src3, src4, src5, src6, src7);
  46. src += (8 * src_stride);
  47. out0 = __msa_copy_u_d((v2i64)src0, 0);
  48. out1 = __msa_copy_u_d((v2i64)src1, 0);
  49. out2 = __msa_copy_u_d((v2i64)src2, 0);
  50. out3 = __msa_copy_u_d((v2i64)src3, 0);
  51. out4 = __msa_copy_u_d((v2i64)src4, 0);
  52. out5 = __msa_copy_u_d((v2i64)src5, 0);
  53. out6 = __msa_copy_u_d((v2i64)src6, 0);
  54. out7 = __msa_copy_u_d((v2i64)src7, 0);
  55. SD4(out0, out1, out2, out3, dst, dst_stride);
  56. dst += (4 * dst_stride);
  57. SD4(out4, out5, out6, out7, dst, dst_stride);
  58. dst += (4 * dst_stride);
  59. }
  60. } else if (0 == height % 4) {
  61. for (cnt = (height / 4); cnt--;) {
  62. LD_UB4(src, src_stride, src0, src1, src2, src3);
  63. src += (4 * src_stride);
  64. out0 = __msa_copy_u_d((v2i64)src0, 0);
  65. out1 = __msa_copy_u_d((v2i64)src1, 0);
  66. out2 = __msa_copy_u_d((v2i64)src2, 0);
  67. out3 = __msa_copy_u_d((v2i64)src3, 0);
  68. SD4(out0, out1, out2, out3, dst, dst_stride);
  69. dst += (4 * dst_stride);
  70. }
  71. } else if (0 == height % 2) {
  72. for (cnt = (height / 2); cnt--;) {
  73. LD_UB2(src, src_stride, src0, src1);
  74. src += (2 * src_stride);
  75. out0 = __msa_copy_u_d((v2i64)src0, 0);
  76. out1 = __msa_copy_u_d((v2i64)src1, 0);
  77. SD(out0, dst);
  78. dst += dst_stride;
  79. SD(out1, dst);
  80. dst += dst_stride;
  81. }
  82. }
  83. }
  84. static void copy_16multx8mult_msa(const uint8_t *src, int32_t src_stride,
  85. uint8_t *dst, int32_t dst_stride,
  86. int32_t height, int32_t width) {
  87. int32_t cnt, loop_cnt;
  88. const uint8_t *src_tmp;
  89. uint8_t *dst_tmp;
  90. v16u8 src0, src1, src2, src3, src4, src5, src6, src7;
  91. for (cnt = (width >> 4); cnt--;) {
  92. src_tmp = src;
  93. dst_tmp = dst;
  94. for (loop_cnt = (height >> 3); loop_cnt--;) {
  95. LD_UB8(src_tmp, src_stride, src0, src1, src2, src3, src4, src5, src6,
  96. src7);
  97. src_tmp += (8 * src_stride);
  98. ST_UB8(src0, src1, src2, src3, src4, src5, src6, src7, dst_tmp,
  99. dst_stride);
  100. dst_tmp += (8 * dst_stride);
  101. }
  102. src += 16;
  103. dst += 16;
  104. }
  105. }
  106. static void copy_width16_msa(const uint8_t *src, int32_t src_stride,
  107. uint8_t *dst, int32_t dst_stride, int32_t height) {
  108. int32_t cnt;
  109. v16u8 src0, src1, src2, src3, src4, src5, src6, src7;
  110. if (0 == height % 12) {
  111. for (cnt = (height / 12); cnt--;) {
  112. LD_UB8(src, src_stride, src0, src1, src2, src3, src4, src5, src6, src7);
  113. src += (8 * src_stride);
  114. ST_UB8(src0, src1, src2, src3, src4, src5, src6, src7, dst, dst_stride);
  115. dst += (8 * dst_stride);
  116. LD_UB4(src, src_stride, src0, src1, src2, src3);
  117. src += (4 * src_stride);
  118. ST_UB4(src0, src1, src2, src3, dst, dst_stride);
  119. dst += (4 * dst_stride);
  120. }
  121. } else if (0 == height % 8) {
  122. copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 16);
  123. } else if (0 == height % 4) {
  124. for (cnt = (height >> 2); cnt--;) {
  125. LD_UB4(src, src_stride, src0, src1, src2, src3);
  126. src += (4 * src_stride);
  127. ST_UB4(src0, src1, src2, src3, dst, dst_stride);
  128. dst += (4 * dst_stride);
  129. }
  130. }
  131. }
  132. static void copy_width32_msa(const uint8_t *src, int32_t src_stride,
  133. uint8_t *dst, int32_t dst_stride, int32_t height) {
  134. int32_t cnt;
  135. v16u8 src0, src1, src2, src3, src4, src5, src6, src7;
  136. if (0 == height % 12) {
  137. for (cnt = (height / 12); cnt--;) {
  138. LD_UB4(src, src_stride, src0, src1, src2, src3);
  139. LD_UB4(src + 16, src_stride, src4, src5, src6, src7);
  140. src += (4 * src_stride);
  141. ST_UB4(src0, src1, src2, src3, dst, dst_stride);
  142. ST_UB4(src4, src5, src6, src7, dst + 16, dst_stride);
  143. dst += (4 * dst_stride);
  144. LD_UB4(src, src_stride, src0, src1, src2, src3);
  145. LD_UB4(src + 16, src_stride, src4, src5, src6, src7);
  146. src += (4 * src_stride);
  147. ST_UB4(src0, src1, src2, src3, dst, dst_stride);
  148. ST_UB4(src4, src5, src6, src7, dst + 16, dst_stride);
  149. dst += (4 * dst_stride);
  150. LD_UB4(src, src_stride, src0, src1, src2, src3);
  151. LD_UB4(src + 16, src_stride, src4, src5, src6, src7);
  152. src += (4 * src_stride);
  153. ST_UB4(src0, src1, src2, src3, dst, dst_stride);
  154. ST_UB4(src4, src5, src6, src7, dst + 16, dst_stride);
  155. dst += (4 * dst_stride);
  156. }
  157. } else if (0 == height % 8) {
  158. copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 32);
  159. } else if (0 == height % 4) {
  160. for (cnt = (height >> 2); cnt--;) {
  161. LD_UB4(src, src_stride, src0, src1, src2, src3);
  162. LD_UB4(src + 16, src_stride, src4, src5, src6, src7);
  163. src += (4 * src_stride);
  164. ST_UB4(src0, src1, src2, src3, dst, dst_stride);
  165. ST_UB4(src4, src5, src6, src7, dst + 16, dst_stride);
  166. dst += (4 * dst_stride);
  167. }
  168. }
  169. }
  170. static void copy_width64_msa(const uint8_t *src, int32_t src_stride,
  171. uint8_t *dst, int32_t dst_stride, int32_t height) {
  172. copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 64);
  173. }
  174. void vpx_convolve_copy_msa(const uint8_t *src, ptrdiff_t src_stride,
  175. uint8_t *dst, ptrdiff_t dst_stride,
  176. const InterpKernel *filter, int x0_q4,
  177. int32_t x_step_q4, int y0_q4, int32_t y_step_q4,
  178. int32_t w, int32_t h) {
  179. (void)filter;
  180. (void)x0_q4;
  181. (void)x_step_q4;
  182. (void)y0_q4;
  183. (void)y_step_q4;
  184. switch (w) {
  185. case 4: {
  186. uint32_t cnt, tmp;
  187. /* 1 word storage */
  188. for (cnt = h; cnt--;) {
  189. tmp = LW(src);
  190. SW(tmp, dst);
  191. src += src_stride;
  192. dst += dst_stride;
  193. }
  194. break;
  195. }
  196. case 8: {
  197. copy_width8_msa(src, src_stride, dst, dst_stride, h);
  198. break;
  199. }
  200. case 16: {
  201. copy_width16_msa(src, src_stride, dst, dst_stride, h);
  202. break;
  203. }
  204. case 32: {
  205. copy_width32_msa(src, src_stride, dst, dst_stride, h);
  206. break;
  207. }
  208. case 64: {
  209. copy_width64_msa(src, src_stride, dst, dst_stride, h);
  210. break;
  211. }
  212. default: {
  213. uint32_t cnt;
  214. for (cnt = h; cnt--;) {
  215. memcpy(dst, src, w);
  216. src += src_stride;
  217. dst += dst_stride;
  218. }
  219. break;
  220. }
  221. }
  222. }