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.

666 lines
13KB

  1. ;*****************************************************************************
  2. ;* x86util.asm
  3. ;*****************************************************************************
  4. ;* Copyright (C) 2008-2010 x264 project
  5. ;*
  6. ;* Authors: Loren Merritt <lorenm@u.washington.edu>
  7. ;* Holger Lubitz <holger@lubitz.org>
  8. ;*
  9. ;* This file is part of Libav.
  10. ;*
  11. ;* Libav is free software; you can redistribute it and/or
  12. ;* modify it under the terms of the GNU Lesser General Public
  13. ;* License as published by the Free Software Foundation; either
  14. ;* version 2.1 of the License, or (at your option) any later version.
  15. ;*
  16. ;* Libav is distributed in the hope that it will be useful,
  17. ;* but WITHOUT ANY WARRANTY; without even the implied warranty of
  18. ;* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  19. ;* Lesser General Public License for more details.
  20. ;*
  21. ;* You should have received a copy of the GNU Lesser General Public
  22. ;* License along with Libav; if not, write to the Free Software
  23. ;* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  24. ;******************************************************************************
  25. %define cpuflags_mmxext cpuflags_mmx2
  26. %include "libavutil/x86/x86inc.asm"
  27. %macro SBUTTERFLY 4
  28. %if avx_enabled == 0
  29. mova m%4, m%2
  30. punpckl%1 m%2, m%3
  31. punpckh%1 m%4, m%3
  32. %else
  33. punpckh%1 m%4, m%2, m%3
  34. punpckl%1 m%2, m%3
  35. %endif
  36. SWAP %3, %4
  37. %endmacro
  38. %macro SBUTTERFLY2 4
  39. punpckl%1 m%4, m%2, m%3
  40. punpckh%1 m%2, m%2, m%3
  41. SWAP %2, %4, %3
  42. %endmacro
  43. %macro SBUTTERFLYPS 3
  44. unpcklps m%3, m%1, m%2
  45. unpckhps m%1, m%1, m%2
  46. SWAP %1, %3, %2
  47. %endmacro
  48. %macro TRANSPOSE4x4B 5
  49. SBUTTERFLY bw, %1, %2, %5
  50. SBUTTERFLY bw, %3, %4, %5
  51. SBUTTERFLY wd, %1, %3, %5
  52. SBUTTERFLY wd, %2, %4, %5
  53. SWAP %2, %3
  54. %endmacro
  55. %macro TRANSPOSE4x4W 5
  56. SBUTTERFLY wd, %1, %2, %5
  57. SBUTTERFLY wd, %3, %4, %5
  58. SBUTTERFLY dq, %1, %3, %5
  59. SBUTTERFLY dq, %2, %4, %5
  60. SWAP %2, %3
  61. %endmacro
  62. %macro TRANSPOSE2x4x4W 5
  63. SBUTTERFLY wd, %1, %2, %5
  64. SBUTTERFLY wd, %3, %4, %5
  65. SBUTTERFLY dq, %1, %3, %5
  66. SBUTTERFLY dq, %2, %4, %5
  67. SBUTTERFLY qdq, %1, %2, %5
  68. SBUTTERFLY qdq, %3, %4, %5
  69. %endmacro
  70. %macro TRANSPOSE4x4D 5
  71. SBUTTERFLY dq, %1, %2, %5
  72. SBUTTERFLY dq, %3, %4, %5
  73. SBUTTERFLY qdq, %1, %3, %5
  74. SBUTTERFLY qdq, %2, %4, %5
  75. SWAP %2, %3
  76. %endmacro
  77. ; identical behavior to TRANSPOSE4x4D, but using SSE1 float ops
  78. %macro TRANSPOSE4x4PS 5
  79. SBUTTERFLYPS %1, %2, %5
  80. SBUTTERFLYPS %3, %4, %5
  81. movlhps m%5, m%1, m%3
  82. movhlps m%3, m%1
  83. SWAP %5, %1
  84. movlhps m%5, m%2, m%4
  85. movhlps m%4, m%2
  86. SWAP %5, %2, %3
  87. %endmacro
  88. %macro TRANSPOSE8x8W 9-11
  89. %if ARCH_X86_64
  90. SBUTTERFLY wd, %1, %2, %9
  91. SBUTTERFLY wd, %3, %4, %9
  92. SBUTTERFLY wd, %5, %6, %9
  93. SBUTTERFLY wd, %7, %8, %9
  94. SBUTTERFLY dq, %1, %3, %9
  95. SBUTTERFLY dq, %2, %4, %9
  96. SBUTTERFLY dq, %5, %7, %9
  97. SBUTTERFLY dq, %6, %8, %9
  98. SBUTTERFLY qdq, %1, %5, %9
  99. SBUTTERFLY qdq, %2, %6, %9
  100. SBUTTERFLY qdq, %3, %7, %9
  101. SBUTTERFLY qdq, %4, %8, %9
  102. SWAP %2, %5
  103. SWAP %4, %7
  104. %else
  105. ; in: m0..m7, unless %11 in which case m6 is in %9
  106. ; out: m0..m7, unless %11 in which case m4 is in %10
  107. ; spills into %9 and %10
  108. %if %0<11
  109. movdqa %9, m%7
  110. %endif
  111. SBUTTERFLY wd, %1, %2, %7
  112. movdqa %10, m%2
  113. movdqa m%7, %9
  114. SBUTTERFLY wd, %3, %4, %2
  115. SBUTTERFLY wd, %5, %6, %2
  116. SBUTTERFLY wd, %7, %8, %2
  117. SBUTTERFLY dq, %1, %3, %2
  118. movdqa %9, m%3
  119. movdqa m%2, %10
  120. SBUTTERFLY dq, %2, %4, %3
  121. SBUTTERFLY dq, %5, %7, %3
  122. SBUTTERFLY dq, %6, %8, %3
  123. SBUTTERFLY qdq, %1, %5, %3
  124. SBUTTERFLY qdq, %2, %6, %3
  125. movdqa %10, m%2
  126. movdqa m%3, %9
  127. SBUTTERFLY qdq, %3, %7, %2
  128. SBUTTERFLY qdq, %4, %8, %2
  129. SWAP %2, %5
  130. SWAP %4, %7
  131. %if %0<11
  132. movdqa m%5, %10
  133. %endif
  134. %endif
  135. %endmacro
  136. ; PABSW macro assumes %1 != %2, while ABS1/2 macros work in-place
  137. %macro PABSW 2
  138. %if cpuflag(ssse3)
  139. pabsw %1, %2
  140. %elif cpuflag(mmxext)
  141. pxor %1, %1
  142. psubw %1, %2
  143. pmaxsw %1, %2
  144. %else
  145. pxor %1, %1
  146. pcmpgtw %1, %2
  147. pxor %2, %1
  148. psubw %2, %1
  149. SWAP %1, %2
  150. %endif
  151. %endmacro
  152. %macro PSIGNW_MMX 2
  153. pxor %1, %2
  154. psubw %1, %2
  155. %endmacro
  156. %macro PSIGNW_SSSE3 2
  157. psignw %1, %2
  158. %endmacro
  159. %macro ABS1_MMX 2 ; a, tmp
  160. pxor %2, %2
  161. pcmpgtw %2, %1
  162. pxor %1, %2
  163. psubw %1, %2
  164. %endmacro
  165. %macro ABS2_MMX 4 ; a, b, tmp0, tmp1
  166. pxor %3, %3
  167. pxor %4, %4
  168. pcmpgtw %3, %1
  169. pcmpgtw %4, %2
  170. pxor %1, %3
  171. pxor %2, %4
  172. psubw %1, %3
  173. psubw %2, %4
  174. %endmacro
  175. %macro ABS1_MMXEXT 2 ; a, tmp
  176. pxor %2, %2
  177. psubw %2, %1
  178. pmaxsw %1, %2
  179. %endmacro
  180. %macro ABS2_MMXEXT 4 ; a, b, tmp0, tmp1
  181. pxor %3, %3
  182. pxor %4, %4
  183. psubw %3, %1
  184. psubw %4, %2
  185. pmaxsw %1, %3
  186. pmaxsw %2, %4
  187. %endmacro
  188. %macro ABS1_SSSE3 2
  189. pabsw %1, %1
  190. %endmacro
  191. %macro ABS2_SSSE3 4
  192. pabsw %1, %1
  193. pabsw %2, %2
  194. %endmacro
  195. %macro ABSB_MMX 2
  196. pxor %2, %2
  197. psubb %2, %1
  198. pminub %1, %2
  199. %endmacro
  200. %macro ABSB2_MMX 4
  201. pxor %3, %3
  202. pxor %4, %4
  203. psubb %3, %1
  204. psubb %4, %2
  205. pminub %1, %3
  206. pminub %2, %4
  207. %endmacro
  208. %macro ABSD2_MMX 4
  209. pxor %3, %3
  210. pxor %4, %4
  211. pcmpgtd %3, %1
  212. pcmpgtd %4, %2
  213. pxor %1, %3
  214. pxor %2, %4
  215. psubd %1, %3
  216. psubd %2, %4
  217. %endmacro
  218. %macro ABSB_SSSE3 2
  219. pabsb %1, %1
  220. %endmacro
  221. %macro ABSB2_SSSE3 4
  222. pabsb %1, %1
  223. pabsb %2, %2
  224. %endmacro
  225. %macro ABS4 6
  226. ABS2 %1, %2, %5, %6
  227. ABS2 %3, %4, %5, %6
  228. %endmacro
  229. %define ABS1 ABS1_MMX
  230. %define ABS2 ABS2_MMX
  231. %define ABSB ABSB_MMX
  232. %define ABSB2 ABSB2_MMX
  233. %macro SPLATB_LOAD 3
  234. %if cpuflag(ssse3)
  235. movd %1, [%2-3]
  236. pshufb %1, %3
  237. %else
  238. movd %1, [%2-3] ;to avoid crossing a cacheline
  239. punpcklbw %1, %1
  240. SPLATW %1, %1, 3
  241. %endif
  242. %endmacro
  243. %macro SPLATB_REG 3
  244. %if cpuflag(ssse3)
  245. movd %1, %2d
  246. pshufb %1, %3
  247. %else
  248. movd %1, %2d
  249. punpcklbw %1, %1
  250. SPLATW %1, %1, 0
  251. %endif
  252. %endmacro
  253. %macro PALIGNR 4-5
  254. %if cpuflag(ssse3)
  255. %if %0==5
  256. palignr %1, %2, %3, %4
  257. %else
  258. palignr %1, %2, %3
  259. %endif
  260. %elif cpuflag(mmx) ; [dst,] src1, src2, imm, tmp
  261. %define %%dst %1
  262. %if %0==5
  263. %ifnidn %1, %2
  264. mova %%dst, %2
  265. %endif
  266. %rotate 1
  267. %endif
  268. %ifnidn %4, %2
  269. mova %4, %2
  270. %endif
  271. %if mmsize==8
  272. psllq %%dst, (8-%3)*8
  273. psrlq %4, %3*8
  274. %else
  275. pslldq %%dst, 16-%3
  276. psrldq %4, %3
  277. %endif
  278. por %%dst, %4
  279. %endif
  280. %endmacro
  281. %macro PSHUFLW 1+
  282. %if mmsize == 8
  283. pshufw %1
  284. %else
  285. pshuflw %1
  286. %endif
  287. %endmacro
  288. %macro PSWAPD 2
  289. %if cpuflag(mmxext)
  290. pshufw %1, %2, q1032
  291. %elif cpuflag(3dnowext)
  292. pswapd %1, %2
  293. %elif cpuflag(3dnow)
  294. movq %1, %2
  295. psrlq %1, 32
  296. punpckldq %1, %2
  297. %endif
  298. %endmacro
  299. %macro DEINTB 5 ; mask, reg1, mask, reg2, optional src to fill masks from
  300. %ifnum %5
  301. pand m%3, m%5, m%4 ; src .. y6 .. y4
  302. pand m%1, m%5, m%2 ; dst .. y6 .. y4
  303. %else
  304. mova m%1, %5
  305. pand m%3, m%1, m%4 ; src .. y6 .. y4
  306. pand m%1, m%1, m%2 ; dst .. y6 .. y4
  307. %endif
  308. psrlw m%2, 8 ; dst .. y7 .. y5
  309. psrlw m%4, 8 ; src .. y7 .. y5
  310. %endmacro
  311. %macro SUMSUB_BA 3-4
  312. %if %0==3
  313. padd%1 m%2, m%3
  314. padd%1 m%3, m%3
  315. psub%1 m%3, m%2
  316. %else
  317. %if avx_enabled == 0
  318. mova m%4, m%2
  319. padd%1 m%2, m%3
  320. psub%1 m%3, m%4
  321. %else
  322. padd%1 m%4, m%2, m%3
  323. psub%1 m%3, m%2
  324. SWAP %2, %4
  325. %endif
  326. %endif
  327. %endmacro
  328. %macro SUMSUB_BADC 5-6
  329. %if %0==6
  330. SUMSUB_BA %1, %2, %3, %6
  331. SUMSUB_BA %1, %4, %5, %6
  332. %else
  333. padd%1 m%2, m%3
  334. padd%1 m%4, m%5
  335. padd%1 m%3, m%3
  336. padd%1 m%5, m%5
  337. psub%1 m%3, m%2
  338. psub%1 m%5, m%4
  339. %endif
  340. %endmacro
  341. %macro SUMSUB2_AB 4
  342. %ifnum %3
  343. psub%1 m%4, m%2, m%3
  344. psub%1 m%4, m%3
  345. padd%1 m%2, m%2
  346. padd%1 m%2, m%3
  347. %else
  348. mova m%4, m%2
  349. padd%1 m%2, m%2
  350. padd%1 m%2, %3
  351. psub%1 m%4, %3
  352. psub%1 m%4, %3
  353. %endif
  354. %endmacro
  355. %macro SUMSUB2_BA 4
  356. %if avx_enabled == 0
  357. mova m%4, m%2
  358. padd%1 m%2, m%3
  359. padd%1 m%2, m%3
  360. psub%1 m%3, m%4
  361. psub%1 m%3, m%4
  362. %else
  363. padd%1 m%4, m%2, m%3
  364. padd%1 m%4, m%3
  365. psub%1 m%3, m%2
  366. psub%1 m%3, m%2
  367. SWAP %2, %4
  368. %endif
  369. %endmacro
  370. %macro SUMSUBD2_AB 5
  371. %ifnum %4
  372. psra%1 m%5, m%2, 1 ; %3: %3>>1
  373. psra%1 m%4, m%3, 1 ; %2: %2>>1
  374. padd%1 m%4, m%2 ; %3: %3>>1+%2
  375. psub%1 m%5, m%3 ; %2: %2>>1-%3
  376. SWAP %2, %5
  377. SWAP %3, %4
  378. %else
  379. mova %5, m%2
  380. mova %4, m%3
  381. psra%1 m%3, 1 ; %3: %3>>1
  382. psra%1 m%2, 1 ; %2: %2>>1
  383. padd%1 m%3, %5 ; %3: %3>>1+%2
  384. psub%1 m%2, %4 ; %2: %2>>1-%3
  385. %endif
  386. %endmacro
  387. %macro DCT4_1D 5
  388. %ifnum %5
  389. SUMSUB_BADC w, %4, %1, %3, %2, %5
  390. SUMSUB_BA w, %3, %4, %5
  391. SUMSUB2_AB w, %1, %2, %5
  392. SWAP %1, %3, %4, %5, %2
  393. %else
  394. SUMSUB_BADC w, %4, %1, %3, %2
  395. SUMSUB_BA w, %3, %4
  396. mova [%5], m%2
  397. SUMSUB2_AB w, %1, [%5], %2
  398. SWAP %1, %3, %4, %2
  399. %endif
  400. %endmacro
  401. %macro IDCT4_1D 6-7
  402. %ifnum %6
  403. SUMSUBD2_AB %1, %3, %5, %7, %6
  404. ; %3: %3>>1-%5 %5: %3+%5>>1
  405. SUMSUB_BA %1, %4, %2, %7
  406. ; %4: %2+%4 %2: %2-%4
  407. SUMSUB_BADC %1, %5, %4, %3, %2, %7
  408. ; %5: %2+%4 + (%3+%5>>1)
  409. ; %4: %2+%4 - (%3+%5>>1)
  410. ; %3: %2-%4 + (%3>>1-%5)
  411. ; %2: %2-%4 - (%3>>1-%5)
  412. %else
  413. %ifidn %1, w
  414. SUMSUBD2_AB %1, %3, %5, [%6], [%6+16]
  415. %else
  416. SUMSUBD2_AB %1, %3, %5, [%6], [%6+32]
  417. %endif
  418. SUMSUB_BA %1, %4, %2
  419. SUMSUB_BADC %1, %5, %4, %3, %2
  420. %endif
  421. SWAP %2, %5, %4
  422. ; %2: %2+%4 + (%3+%5>>1) row0
  423. ; %3: %2-%4 + (%3>>1-%5) row1
  424. ; %4: %2-%4 - (%3>>1-%5) row2
  425. ; %5: %2+%4 - (%3+%5>>1) row3
  426. %endmacro
  427. %macro LOAD_DIFF 5
  428. %ifidn %3, none
  429. movh %1, %4
  430. movh %2, %5
  431. punpcklbw %1, %2
  432. punpcklbw %2, %2
  433. psubw %1, %2
  434. %else
  435. movh %1, %4
  436. punpcklbw %1, %3
  437. movh %2, %5
  438. punpcklbw %2, %3
  439. psubw %1, %2
  440. %endif
  441. %endmacro
  442. %macro STORE_DCT 6
  443. movq [%5+%6+ 0], m%1
  444. movq [%5+%6+ 8], m%2
  445. movq [%5+%6+16], m%3
  446. movq [%5+%6+24], m%4
  447. movhps [%5+%6+32], m%1
  448. movhps [%5+%6+40], m%2
  449. movhps [%5+%6+48], m%3
  450. movhps [%5+%6+56], m%4
  451. %endmacro
  452. %macro LOAD_DIFF_8x4P 7-10 r0,r2,0 ; 4x dest, 2x temp, 2x pointer, increment?
  453. LOAD_DIFF m%1, m%5, m%7, [%8], [%9]
  454. LOAD_DIFF m%2, m%6, m%7, [%8+r1], [%9+r3]
  455. LOAD_DIFF m%3, m%5, m%7, [%8+2*r1], [%9+2*r3]
  456. LOAD_DIFF m%4, m%6, m%7, [%8+r4], [%9+r5]
  457. %if %10
  458. lea %8, [%8+4*r1]
  459. lea %9, [%9+4*r3]
  460. %endif
  461. %endmacro
  462. %macro DIFFx2 6-7
  463. movh %3, %5
  464. punpcklbw %3, %4
  465. psraw %1, 6
  466. paddsw %1, %3
  467. movh %3, %6
  468. punpcklbw %3, %4
  469. psraw %2, 6
  470. paddsw %2, %3
  471. packuswb %2, %1
  472. %endmacro
  473. %macro STORE_DIFF 4
  474. movh %2, %4
  475. punpcklbw %2, %3
  476. psraw %1, 6
  477. paddsw %1, %2
  478. packuswb %1, %1
  479. movh %4, %1
  480. %endmacro
  481. %macro STORE_DIFFx2 8 ; add1, add2, reg1, reg2, zero, shift, source, stride
  482. movh %3, [%7]
  483. movh %4, [%7+%8]
  484. psraw %1, %6
  485. psraw %2, %6
  486. punpcklbw %3, %5
  487. punpcklbw %4, %5
  488. paddw %3, %1
  489. paddw %4, %2
  490. packuswb %3, %5
  491. packuswb %4, %5
  492. movh [%7], %3
  493. movh [%7+%8], %4
  494. %endmacro
  495. %macro PMINUB 3 ; dst, src, ignored
  496. %if cpuflag(mmxext)
  497. pminub %1, %2
  498. %else ; dst, src, tmp
  499. mova %3, %1
  500. psubusb %3, %2
  501. psubb %1, %3
  502. %endif
  503. %endmacro
  504. %macro SPLATW 2-3 0
  505. %if mmsize == 16
  506. pshuflw %1, %2, (%3)*0x55
  507. punpcklqdq %1, %1
  508. %elif cpuflag(mmx2)
  509. pshufw %1, %2, (%3)*0x55
  510. %else
  511. %ifnidn %1, %2
  512. mova %1, %2
  513. %endif
  514. %if %3 & 2
  515. punpckhwd %1, %1
  516. %else
  517. punpcklwd %1, %1
  518. %endif
  519. %if %3 & 1
  520. punpckhwd %1, %1
  521. %else
  522. punpcklwd %1, %1
  523. %endif
  524. %endif
  525. %endmacro
  526. %macro SPLATD 2-3 0
  527. %if mmsize == 16
  528. pshufd %1, %2, (%3)*0x55
  529. %else
  530. pshufw %1, %2, (%3)*0x11 + ((%3)+1)*0x44
  531. %endif
  532. %endmacro
  533. %macro SPLATD_MMX 1
  534. punpckldq %1, %1
  535. %endmacro
  536. %macro SPLATD_SSE 1
  537. shufps %1, %1, 0
  538. %endmacro
  539. %macro SPLATD_SSE2 1
  540. pshufd %1, %1, 0
  541. %endmacro
  542. %macro CLIPW 3 ;(dst, min, max)
  543. pmaxsw %1, %2
  544. pminsw %1, %3
  545. %endmacro
  546. %macro PMINSD_MMX 3 ; dst, src, tmp
  547. mova %3, %2
  548. pcmpgtd %3, %1
  549. pxor %1, %2
  550. pand %1, %3
  551. pxor %1, %2
  552. %endmacro
  553. %macro PMAXSD_MMX 3 ; dst, src, tmp
  554. mova %3, %1
  555. pcmpgtd %3, %2
  556. pand %1, %3
  557. pandn %3, %2
  558. por %1, %3
  559. %endmacro
  560. %macro CLIPD_MMX 3-4 ; src/dst, min, max, tmp
  561. PMINSD_MMX %1, %3, %4
  562. PMAXSD_MMX %1, %2, %4
  563. %endmacro
  564. %macro CLIPD_SSE2 3-4 ; src/dst, min (float), max (float), unused
  565. cvtdq2ps %1, %1
  566. minps %1, %3
  567. maxps %1, %2
  568. cvtps2dq %1, %1
  569. %endmacro
  570. %macro CLIPD_SSE41 3-4 ; src/dst, min, max, unused
  571. pminsd %1, %3
  572. pmaxsd %1, %2
  573. %endmacro
  574. %macro VBROADCASTSS 2 ; dst xmm/ymm, src m32
  575. %if cpuflag(avx)
  576. vbroadcastss %1, %2
  577. %else ; sse
  578. movss %1, %2
  579. shufps %1, %1, 0
  580. %endif
  581. %endmacro
  582. %macro SHUFFLE_MASK_W 8
  583. %rep 8
  584. %if %1>=0x80
  585. db %1, %1
  586. %else
  587. db %1*2
  588. db %1*2+1
  589. %endif
  590. %rotate 1
  591. %endrep
  592. %endmacro
  593. %macro PMOVSXWD 2; dst, src
  594. %if cpuflag(sse4)
  595. pmovsxwd %1, %2
  596. %else
  597. %ifnidn %1, %2
  598. mova %1, %2
  599. %endif
  600. punpcklwd %1, %1
  601. psrad %1, 16
  602. %endif
  603. %endmacro