开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > DirextX编程 > 查看源码
idctfpu.c
资源名称：DirectShowcode.rar　[点击查看]
上传用户：hhs829
上传日期：2022-06-17
资源大小：586k
文件大小：12k
源码类别：

DirextX编程

开发平台：

Visual C++

idctfpu.c：源码内容
  1. /* idct.c, inverse fast discrete cosine transform                           */
  4. /*************************************************************/
  5. /* inverse two dimensional DCT, Chen-Wang algorithm          */
  6. /* (cf. IEEE ASSP-32, pp. 803-816, Aug. 1984)                */
  7. /*                                                           */
  8. /* floating point conversion by Miha Peternel                */
  9. /* x87 hand-optimized assembly by Miha Peternel              */
  10. /*                                    27.11. - 11.12.2000    */
  11. /*                                                           */
  12. /* You are free to use this code in your project if:         */
  13. /* - no changes are made to this message                     */
  14. /* - any changes to this code are publicly available         */
  15. /* - your project documentation contains the following text: */
  16. /*   "This software contains fast high-quality IDCT decoder  */
  17. /*    by Miha Peternel."                                     */
  18. /*                                                           */
  19. /*************************************************************/
  22. /////////////////////////////////////////////////////
  23. //
  24. // TODO:
  25. // - loops can be easily vectorized for SIMD
  26. //
  27. /////////////////////////////////////////////////////
  29. #include <math.h>
  30. #  define PI 3.1415926535897932384626433832795
  32. #define FLOAT double
  34. const static double RC = 1.0*1024*1024*1024*1024*256*16 + 1024; // magic + clip center
  36. static FLOAT W1; // /* sqrt(2)*cos(1*pi/16) */
  37. static FLOAT W2; // /* sqrt(2)*cos(2*pi/16) */
  38. static FLOAT W5; // /* sqrt(2)*cos(5*pi/16) */
  40. static FLOAT W1_8;
  41. static FLOAT W2_8;
  42. static FLOAT W5_8;
  44. static FLOAT W7; // /* sqrt(2)*cos(7*pi/16) */
  45. static FLOAT W1mW7; // W1-W7
  46. static FLOAT W1pW7; // W1+W7
  48. static FLOAT W3; // /* sqrt(2)*cos(3*pi/16) */
  49. static FLOAT W3mW5; // W3-W5
  50. static FLOAT W3pW5; // W3+W5
  52. static FLOAT W6; // /* sqrt(2)*cos(6*pi/16) */
  53. static FLOAT W2mW6; // W2-W6
  54. static FLOAT W2pW6; // W2+W6
  56. static FLOAT S2; // 1/sqrt(2)
  57. static FLOAT D8 = 1.0/8;
  59. static FLOAT W7_8;
  60. static FLOAT W1mW7_8;
  61. static FLOAT W1pW7_8;
  63. static FLOAT W3_8;
  64. static FLOAT W3mW5_8;
  65. static FLOAT W3pW5_8;
  67. static FLOAT W6_8;
  68. static FLOAT W2mW6_8;
  69. static FLOAT W2pW6_8;
  71. /* private data */
  72. static short iclip[1024+1024]; /* clipping table */
  73. static short *iclp;
  75. void Initialize_FPU_IDCT()
  76. {
  77.   int i;
  79.   S2 = sqrt(0.5); // 1.0/sqrt(2);
  81.   W1 = sqrt(2)*cos(PI*(1.0/16)); 
  82. W1_8 = W1/8;
  83.   W2 = sqrt(2)*cos(PI*(2.0/16)); 
  84. W2_8 = W2/8;
  85.   W3 = sqrt(2)*cos(PI*(3.0/16)); 
  86. W3_8 = W3/8;
  87.   W5 = sqrt(2)*cos(PI*(5.0/16)); 
  88. W5_8 = W5/8;
  89.   W6 = sqrt(2)*cos(PI*(6.0/16)); 
  90. W6_8 = W6/8;
  91.   W7 = sqrt(2)*cos(PI*(7.0/16));
  92. W7_8 = W7/8;
  94.   W1mW7 = W1-W7;  W1mW7_8 = W1mW7/8;
  95.   W1pW7 = W1+W7;  W1pW7_8 = W1pW7/8;
  96.   W3mW5 = W3-W5;  W3mW5_8 = W3mW5/8;
  97.   W3pW5 = W3+W5;  W3pW5_8 = W3pW5/8;
  98.   W2mW6 = W2-W6;  W2mW6_8 = W2mW6/8;
  99.   W2pW6 = W2+W6;  W2pW6_8 = W2pW6/8;
  101.   iclp = iclip+1024;
  102.   for (i= -1024; i<1024; i++)
  103.     iclp[i] = (i<-256) ? -256 : ((i>255) ? 255 : i);
  104. }
  106. /*
  107. int all = 0;
  108. int turn = 0;
  109. int diff = 0;
  110. int empty = 0;
  111. int allused = 0;
  112. int used[65] = {0};
  113. int trows[9] = {0};
  114. int tcols[9] = {0};
  115. */
  117. void FPU_IDCT(short *block)
  118. {
  119. /*
  120.   int row[8] = {0};
  121.   int col[8] = {0};
  122. int dirty = 0;
  123. int x,y,nrows,ncols;
  124. */
  125. int *b = (int *) block;
  126.   if( b[0]==0 && (b[31]==0x10000 || b[31]==0) )
  127. {
  128.   if( b[ 1]|b[ 2]|b[ 3]|b[ 4]|b[ 5] )
  129.   goto normal;
  130.   if( b[ 6]|b[ 7]|b[ 8]|b[ 9]|b[10] )
  131.   goto normal;
  132.   if( b[11]|b[12]|b[13]|b[14]|b[15] )
  133.   goto normal;
  134.   if( b[16]|b[17]|b[18]|b[19]|b[20] )
  135.   goto normal;
  136.   if( b[21]|b[22]|b[23]|b[24]|b[25] )
  137.   goto normal;
  138.   if( b[26]|b[27]|b[28]|b[29]|b[30] )
  139.   goto normal;
  140. b[31]=0;
  141. ////empty++;
  142. return;
  143. }
  144. normal:
  145. /*
  146. for( y = 0; y < 8; y++ )
  147. for( x = 0; x < 8; x++ )
  148.   if( block[y*8+x] )
  149. {
  150.   row[y] = 1;
  151. col[x] = 1;
  152. dirty++;
  153. }
  154. allused += dirty;
  155. used[dirty]++;
  156. nrows = row[0]+row[1]+row[2]+row[3]+row[4]+row[5]+row[6]+row[7];
  157. trows[nrows]++;
  158. ncols = col[0]+col[1]+col[2]+col[3]+col[4]+col[5]+col[6]+col[7];
  159. tcols[ncols]++;
  160. if( ncols < nrows )
  161.   turn += nrows-ncols;
  163.   if( block[63] != 1 )
  164.   {
  165.   diff++;
  166. }
  167. all++;
  168. if( (all % 10000) == 0 )
  169. {
  170.   all *= 1;
  171. }
  172. if( (all % 100000) == 0 )
  173. {
  174.   all *= 1;
  175. }
  176. //*/
  177. //  FLOAT tmp[8*8];
  178. // FLOAT tmp1, tmp2, tmp3;
  179. // int int0, int1, int2, int3, int4, int5, int6, int7;
  180. #define tmp  ebx
  181. #define tmp1 ebx-1*8
  182. #define tmp2 ebx-2*8
  183. #define tmp3 ebx-3*8
  184. #define int0 ebx-3*8-1*4
  185. #define int1 ebx-3*8-2*4
  186. #define int2 ebx-3*8-3*4
  187. #define int3 ebx-3*8-4*4
  188. #define int4 ebx-3*8-5*4
  189. #define int5 ebx-3*8-6*4
  190. #define int6 ebx-3*8-7*4
  191. #define int7 ebx-3*8-8*4
  192. #define SIZE 8*8*8+3*8+8*4+16 // locals + 16-byte alignment area
  193. __asm
  194. {
  195.   lea ebx,[esp-8*8*8]
  196. sub esp,SIZE
  197. and ebx,-16 // force 16-byte alignment of locals
  199. // rows
  200. mov esi,[block]
  201. lea edi,[tmp]
  202. mov ecx,8
  204. align 16
  205. Lrows:
  206.     movsx eax,word ptr [esi+2]
  207. or    eax,         [esi+4]
  208. or    eax,         [esi+8]
  209. or    eax,         [esi+12]
  210. jnz L1
  212. fild word ptr [esi+0*2]
  213. fst  qword ptr [edi+7*8]
  214. fst  qword ptr [edi+6*8]
  215. fst  qword ptr [edi+5*8]
  216. fst  qword ptr [edi+4*8]
  217. fst  qword ptr [edi+3*8]
  218. fst  qword ptr [edi+2*8]
  219. fst  qword ptr [edi+1*8]
  220. fstp qword ptr [edi+0*8]
  221. jmp L2
  223. align 16
  224. L1:
  226. fild word ptr [esi+7*2]
  227. fld st(0)
  228. fild word ptr [esi+1*2]
  229. fadd st(1),st(0)
  230. fld qword ptr [W7]
  231. fxch st(1)
  232. fmul qword ptr [W1mW7]
  233. fxch st(1)
  234. fmulp st(2),st(0)
  235. fadd st(0),st(1)
  236. fstp qword ptr [tmp1]
  237. fild word ptr [esi+3*2]
  238. fld st(0)
  239. fxch st(3)
  240. fmul qword ptr [W1pW7]
  241. fild word ptr [esi+5*2]
  242. fadd st(4),st(0)
  243. fmul qword ptr [W3mW5]
  244. fxch st(1)
  245. fsubp st(3),st(0)//fsubrp
  246. fld qword ptr [W3]
  247. fmulp st(4),st(0)
  248. fsubr st(0),st(3)
  249. fstp qword ptr [tmp2]
  250. fmul qword ptr [W3pW5]
  251. fsubp st(2),st(0)//fsubrp
  252. fxch st(1)
  253. fstp qword ptr [tmp3]
  254. fild word ptr [esi+0*2]
  255. fild word ptr [esi+4*2]
  256. fild word ptr [esi+2*2]
  257. fld st(0)
  258. fmul qword ptr [W2mW6]
  259. fld st(3)
  260. fild word ptr [esi+6*2]
  261. fxch st(5)
  262. fsub st(0),st(4)
  263. fxch st(3)
  264. fadd st(0),st(5)
  265. fxch st(1)
  266. faddp st(4),st(0)
  267. fld qword ptr [W6]
  268. fmulp st(1),st(0)
  269. fxch st(4)
  270. fmul qword ptr [W2pW6]
  271. fld qword ptr [tmp1]
  272. fsub qword ptr [tmp2]
  273. fld st(5)
  274. fxch st(3)
  275. faddp st(6),st(0)
  276. fld qword ptr [tmp1]
  277. fxch st(1)
  278. fstp qword ptr [tmp1]
  279. fld st(6)
  280. fadd qword ptr [tmp3]
  281. fxch st(1)
  282. fadd qword ptr [tmp2]
  283. fxch st(7)
  284. fsub qword ptr [tmp3]
  285. fxch st(1)
  286. fstp qword ptr [tmp2]
  287. fld st(4)
  288. fxch st(3)
  289. fsubrp st(2),st(0)//fsubp
  290. fxch st(4)
  291. fsub st(0),st(5)
  292. fxch st(2)
  293. faddp st(5),st(0)
  294. fld st(2)
  295. fsub st(0),st(1)
  296. fxch st(5)
  297. fstp qword ptr [tmp3]
  298. fld qword ptr [tmp1]
  299. fld qword ptr [S2]
  300. fxch st(4)
  301. faddp st(2),st(0)
  302. fld st(3)
  303. fxch st(1)
  304. fadd st(0),st(5)
  305. fmulp st(1),st(0)
  306. //
  307. /*
  308. fld st(1)
  309. fxch st(5)
  310. fsubr qword ptr [tmp1]
  311. fxch st(5)
  312. fadd st(0),st(1)
  313. fxch st(4)
  314. fmulp st(5),st(0)
  315. fxch st(3)
  316. fstp qword ptr [edi+1*8]
  317. fld st(4)
  318. fxch st(1)
  319. fsubrp st(3),st(0)//fsubp
  320. fxch st(4)
  321. fsub st(0),st(3)
  322. fxch st(4)
  323. faddp st(3),st(0)
  324. fxch st(1)
  325. fstp qword ptr [edi+6*8]
  326. fxch st(2)
  327. fstp qword ptr [edi+5*8]
  328. fstp qword ptr [edi+2*8]
  329. fld qword ptr [tmp3]
  330. fadd st(0),st(2)
  331. fstp qword ptr [edi+0*8]
  332. fxch st(1)
  333. fsubr qword ptr [tmp3]
  334. fld st(1)
  335. fxch st(1)
  336. fstp qword ptr [edi+7*8]
  337. fadd qword ptr [tmp2]
  338. fstp qword ptr [edi+3*8]
  339. fsub qword ptr [tmp2]
  340. fstp qword ptr [edi+4*8]
  341. */
  342. //
  343. fld qword ptr [tmp3]
  344. fadd st(0),st(7)
  345. fxch st(5)
  346. fsubr qword ptr [tmp1]
  347. fxch st(5)
  348. fstp qword ptr [edi+0*8]
  349. fxch st(6)
  350. fsubr qword ptr [tmp3]
  351. fld st(2)
  352. fxch st(1)
  353. fstp qword ptr [edi+7*8]
  354. fadd qword ptr [tmp2]
  355. fxch st(3)
  356. fmulp st(4),st(0)
  357. fxch st(2)
  358. fstp qword ptr [edi+3*8]
  359. fld st(1)
  360. fadd st(0),st(5)
  361. fxch st(1)
  362. fsub qword ptr [tmp2]
  363. fxch st(2)
  364. fsubrp st(5),st(0)//fsubp
  365. fstp qword ptr [edi+1*8]
  366. fld st(2)
  367. fxch st(1)
  368. fstp qword ptr [edi+4*8]
  369. fxch st(2)
  370. fsub st(0),st(1)
  371. fxch st(2)
  372. faddp st(1),st(0)
  373. fxch st(2)
  374. fstp qword ptr [edi+6*8]
  375. fstp qword ptr [edi+5*8]
  376. fstp qword ptr [edi+2*8]
  377. L2:
  378.   add esi,8*2
  379. add edi,8*8
  380. dec ecx
  381. jnz Lrows
  383. // columns
  384.     lea esi,[tmp]
  385. mov edi,[block]
  386. lea edx,[iclip+1024*2]
  387. mov ecx,8
  389.     align 16
  390. Lcols:
  391. /*
  392.     mov eax,[esi+1*8*8]
  393. or  eax,[esi+2*8*8]
  394. or  eax,[esi+3*8*8]
  395. or  eax,[esi+4*8*8]
  396. or  eax,[esi+5*8*8]
  397. or  eax,[esi+6*8*8]
  398. or  eax,[esi+7*8*8]
  399. jnz L3
  401. fld qword ptr [esi+0*8*8]
  402. fmul qword ptr [D8]
  403. fistp dword ptr [int0]
  404. mov eax,[int0]
  405. mov ax,word ptr [edx+2*eax]
  406. mov [edi+0*8*2],ax
  407. mov [edi+1*8*2],ax
  408. mov [edi+2*8*2],ax
  409. mov [edi+3*8*2],ax
  410. mov [edi+4*8*2],ax
  411. mov [edi+5*8*2],ax
  412. mov [edi+6*8*2],ax
  413. mov [edi+7*8*2],ax
  414. jmp L4
  416. align 16
  417. L3:
  418. //*/
  419. fld qword ptr [esi+7*8*8]
  420. fld st(0)
  421. fld qword ptr [esi+1*8*8]
  422. fadd st(1),st(0)
  423. fld qword ptr [W7_8]
  424. fxch st(1)
  425. fmul qword ptr [W1mW7_8]
  426. fxch st(1)
  427. fmulp st(2),st(0)
  428. fadd st(0),st(1)
  429. fstp qword ptr [tmp2]
  430. fld qword ptr [esi+3*8*8]
  431. fld st(0)
  432. fxch st(3)
  433. fmul qword ptr [W1pW7_8]
  434. fld qword ptr [esi+5*8*8]
  435. fadd st(4),st(0)
  436. fmul qword ptr [W3mW5_8]
  437. fxch st(1)
  438. fsubp st(3),st(0)//fsubrp
  439. fld qword ptr [W3_8]
  440. fmulp st(4),st(0)
  441. fsubr st(0),st(3)
  442. fstp qword ptr [tmp3]
  443. fld qword ptr [D8]
  444. fld qword ptr [esi+0*8*8]
  445. fmul st(0),st(1)
  446. fxch st(2)
  447. fmul qword ptr [W3pW5_8]
  448. fld qword ptr [esi+4*8*8]
  449. fmulp st(2),st(0)
  450. fld qword ptr [esi+6*8*8]
  451. fld st(3)
  452. fxch st(6)
  453. fsubrp st(2),st(0)//fsubp
  454. fld qword ptr [esi+2*8*8]
  455. fld st(0)
  456. fxch st(5)
  457. fsub st(0),st(4)
  458. fxch st(7)
  459. faddp st(4),st(0)
  460. fxch st(4)
  461. fadd st(0),st(1)
  462. fld qword ptr [W6_8]
  463. fxch st(2)
  464. fmul qword ptr [W2pW6_8]
  465. fxch st(2)
  466. fmulp st(1),st(0)
  467. fxch st(4)
  468. fmul qword ptr [W2mW6_8]
  469. fld qword ptr [tmp2]
  470. fsub qword ptr [tmp3]
  471. fxch st(2)
  472. fsubr st(0),st(5)
  473. fxch st(1)
  474. faddp st(5),st(0)
  475. fld qword ptr [tmp2]
  476. fxch st(2)
  477. fstp qword ptr [tmp2]
  478. fld st(5)
  479. fxch st(2)
  480. fadd qword ptr [tmp3]
  481. fxch st(6)
  482. fsub st(0),st(3)
  483. fxch st(2)
  484. faddp st(3),st(0)
  485. fld st(3)
  486. fsub st(0),st(5)
  487. fxch st(3)
  488. fstp qword ptr [tmp3]
  489. fxch st(3)
  490. faddp st(4),st(0)
  491. fld st(5)
  492. fld qword ptr [tmp2]
  493. fxch st(7)
  494. fsub st(0),st(4)
  495. fxch st(7)
  496. fadd st(0),st(2)
  497. fxch st(1)
  498. faddp st(4),st(0)
  499. fld qword ptr [S2]
  500. fmul st(1),st(0)
  501. fxch st(1)
  502. fstp qword ptr [tmp1]
  503. fld st(4)
  504. fadd st(0),st(6)
  505. fxch st(2)
  506. fsubr qword ptr [tmp2]
  507. fxch st(5)
  508. fsubrp st(6),st(0)//fsubp
  509. fxch st(1)
  510. fistp dword ptr [int0]
  511. fxch st(4)
  512. mov eax,[int0]
  513. movsx eax,word ptr [edx+2*eax]
  514. mov [edi+0*8*2],ax
  515. fistp dword ptr [int7]
  516. mov eax,[int7]
  517. fld st(0)
  518. movsx eax,word ptr [edx+2*eax]
  519. mov [edi+7*8*2],ax
  520. fadd qword ptr [tmp3]
  521. fistp dword ptr [int3]
  522. mov eax,[int3]
  523. movsx eax,word ptr [edx+2*eax]
  524. mov [edi+3*8*2],ax
  525. fsub qword ptr [tmp3]
  526. fld st(1)
  527. fxch st(1)
  528. fistp dword ptr [int4]
  529. mov eax,[int4]
  530. movsx eax,word ptr [edx+2*eax]
  531. mov [edi+4*8*2],ax
  532. fadd qword ptr [tmp1]
  533. fxch st(3)
  534. fmulp st(2),st(0)
  535. fxch st(2)
  536. fistp dword ptr [int1]
  537. fxch st(1)
  538. mov eax,[int1]
  539. movsx eax,word ptr [edx+2*eax]
  540. mov [edi+1*8*2],ax
  541. fsub qword ptr [tmp1]
  542. fld st(2)
  543. fsub st(0),st(2)
  544. fxch st(1)
  545. fistp dword ptr [int6]
  546. fxch st(2)
  547. mov eax,[int6]
  548. faddp st(1),st(0)
  549. movsx eax,word ptr [edx+2*eax]
  550. mov [edi+6*8*2],ax
  551. fistp dword ptr [int2]
  552. mov eax,[int2]
  553. movsx eax,word ptr [edx+2*eax]
  554. mov [edi+2*8*2],ax
  555. fistp dword ptr [int5]
  556. mov eax,[int5]
  557. movsx eax,word ptr [edx+2*eax]
  558. mov [edi+5*8*2],ax
  560.   add esi,8
  561. add edi,2
  562. dec ecx
  563. jnz Lcols
  565. add esp,SIZE
  566.   }
  567. }