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

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

DirextX编程

开发平台：

Visual C++

idctref.c：源码内容
  1. /* idctref_miha.c, Inverse Discrete Fourier Transform, double precision */ /*************************************************************/
  2. /*                                                           */
  3. /* x87 hand-optimized assembly by Miha Peternel              */
  4. /*                                     27.11. - 20.1.2001    */
  5. /*                                                           */
  6. /* You are free to use this code in your project if:         */
  7. /* - no changes are made to this message                     */
  8. /* - any changes to this code are publicly available         */
  9. /* - your project documentation contains the following text: */
  10. /*   "This software contains fast high-quality IDCT decoder  */
  11. /*    by Miha Peternel."                                     */
  12. /*                                                           */
  13. /*************************************************************/
  14. /*  Perform IEEE 1180 reference (64-bit floating point, separable 8x1  *  direct matrix multiply) Inverse Discrete Cosine Transform */
  15. #define ModelX 123 // enable C-level optimizations by Miha Peternel /* Here we use math.h to generate constants.  Compiler results may    vary a little */ #include <math.h> #define M_PI 3.1415926535897932384626433832795 const static double HALF = 0.5;
  17. /* private data */
  18. static short iclip[1024+1024]; /* clipping table */
  19. static short *iclp;
  20. /* cosine transform matrix for 8x1 IDCT */ static double c[8][8]; /* initialize DCT coefficient matrix */ void Initialize_REF_IDCT() {   int freq, time, i;   double scale;   for (freq=0; freq < 8; freq++)   {     scale = (freq == 0) ? sqrt(0.125) : 0.5;     for (time=0; time<8; time++)       c[freq][time] = scale*cos((M_PI/8.0)*freq*(time + 0.5));   }
  22. #ifdef ModelX   iclp = iclip+1024;
  23.   for (i= -1024; i<1024; i++)
  24.     iclp[i] = (i<-256) ? -256 : ((i>255) ? 255 : i);
  25. #endif
  26. }
  27. void REF_IDCT(short *block) { double tmp[64];
  28. double rnd[64];
  29. int int0, int1, int2, int3, int4, int5, int6, int7;
  30. unsigned short fpold;
  31. unsigned short fpnew;
  33. int *b = (int *) block;
  34.   //if( b[0]==0 && (b[31]==0x10000 || b[31]==0) )
  35.   //if( b[0]==0 && ((b[31]&~0x10000)==0) )
  36.   if( !(b[0]|(b[31]&~0x10000)) )
  37. {
  38.   if( b[ 1]|b[ 2]|b[ 3]|b[ 4]|b[ 5]|b[ 6] )
  39.   goto normal;
  40.   if( b[ 7]|b[ 8]|b[ 9]|b[10]|b[11]|b[12] )
  41.   goto normal;
  42.   if( b[13]|b[14]|b[15]|b[16]|b[17]|b[18] )
  43.   goto normal;
  44.   if( b[19]|b[20]|b[21]|b[22]|b[23]|b[24] )
  45.   goto normal;
  46.   if( b[25]|b[26]|b[27]|b[28]|b[29]|b[30] )
  47.   goto normal;
  48. b[31]=0;
  49. return;
  50. }
  51. normal:
  53. __asm
  54. {
  55. // do the IDCT
  56. mov esi,[block]
  57. lea eax,[c]
  58. lea edi,[tmp]
  59. //mov ebx,8
  60. mov ebx,8 // 0x77000000 // 8
  61. align 16
  62. __col1:
  63.     movzx edx,[esi+1*2]
  64. mov   ecx,[esi+2*2]
  65. or    edx,[esi+4*2]
  66. or    ecx,[esi+6*2]
  67. or edx,ecx
  68. //mov ecx,8
  69. mov ecx,8/2 // 0x77000000 // 8
  70. //fild  word ptr [esi+0*2]
  71. //fmul qword ptr [eax+0*8*8]
  72. jnz __row1
  73. fild  word ptr [esi+0*2]
  74. fmul qword ptr [eax+0*8*8]
  75. fst  qword ptr [edi+0*8]
  76. fst  qword ptr [edi+1*8]
  77. fst  qword ptr [edi+2*8]
  78. fst  qword ptr [edi+3*8]
  79. fst  qword ptr [edi+4*8]
  80. fst  qword ptr [edi+5*8]
  81. fst  qword ptr [edi+6*8]
  82. fstp qword ptr [edi+7*8]
  83. add edi,8*8
  84. jmp __next1
  85. align 16
  86. __row1:
  87. fild  word ptr [esi+0*2]
  88. fmul qword ptr [eax+0*8*8]
  89. fild  word ptr [esi+1*2]
  90. fmul qword ptr [eax+1*8*8]
  91. fadd
  92. fild  word ptr [esi+2*2]
  93. fmul qword ptr [eax+2*8*8]
  94. fadd
  95. fild  word ptr [esi+3*2]
  96. fmul qword ptr [eax+3*8*8]
  97. fadd
  98. fild  word ptr [esi+4*2]
  99. fmul qword ptr [eax+4*8*8]
  100. fadd
  101. fild  word ptr [esi+5*2]
  102. fmul qword ptr [eax+5*8*8]
  103. fadd
  104. fild  word ptr [esi+6*2]
  105. fmul qword ptr [eax+6*8*8]
  106. fadd
  107. fild  word ptr [esi+7*2]
  108. fmul qword ptr [eax+7*8*8]
  109. fadd
  110. //add eax,8
  111. //fstp qword ptr [edi]
  112. //add edi,8
  113. //
  114. fild  word ptr [esi+0*2]
  115. fmul qword ptr [eax+0*8*8+8]
  116. fild  word ptr [esi+1*2]
  117. fmul qword ptr [eax+1*8*8+8]
  118. fadd
  119. fild  word ptr [esi+2*2]
  120. fmul qword ptr [eax+2*8*8+8]
  121. fadd
  122. fild  word ptr [esi+3*2]
  123. fmul qword ptr [eax+3*8*8+8]
  124. fadd
  125. fild  word ptr [esi+4*2]
  126. fmul qword ptr [eax+4*8*8+8]
  127. fadd
  128. fild  word ptr [esi+5*2]
  129. fmul qword ptr [eax+5*8*8+8]
  130. fadd
  131. fild  word ptr [esi+6*2]
  132. fmul qword ptr [eax+6*8*8+8]
  133. fadd
  134. fild  word ptr [esi+7*2]
  135. fmul qword ptr [eax+7*8*8+8]
  136. fadd
  137. add eax,8*2
  138. fxch st(1)
  139. fstp qword ptr [edi]//
  140. fstp qword ptr [edi+8]
  141. add edi,8*2
  142. dec ecx
  143. // sub ecx,0x80000001 // add ecx,ecx 
  144. // js  __row1
  145.   //align 16
  146. //   test ecx,ecx // align jump &| redo flags
  147. jnz __row1
  148. add eax,-8*8
  149.   //align 16
  150. __next1:
  151. add esi,+8*2
  152. //dec ebx
  153. sub ebx,0x80000001 // add ebx,ebx 
  154. js  __col1
  155. //align 16
  156. test ebx,ebx // align jump &| redo flags
  157. jnz __col1
  159. lea esi,[tmp]
  160. lea eax,[c]
  161. lea edi,[rnd]
  162. //mov edi,[block]
  163.     fld qword ptr [HALF]
  164. mov ebx,8
  165. __row2:
  166. mov ecx,8/2
  167. align 16
  168. __col2:
  169. fld  qword ptr [esi+0*8*8]
  170. fmul qword ptr [eax+0*8*8]
  171. fld  qword ptr [esi+1*8*8]
  172. fmul qword ptr [eax+1*8*8]
  173. fadd
  174. fld  qword ptr [esi+2*8*8]
  175. fmul qword ptr [eax+2*8*8]
  176. fadd
  177. fld  qword ptr [esi+3*8*8]
  178. fmul qword ptr [eax+3*8*8]
  179. fadd
  180. fld  qword ptr [esi+4*8*8]
  181. fmul qword ptr [eax+4*8*8]
  182. fadd
  183. fld  qword ptr [esi+5*8*8]
  184. fmul qword ptr [eax+5*8*8]
  185. fadd
  186. fld  qword ptr [esi+6*8*8]
  187. fmul qword ptr [eax+6*8*8]
  188. fadd
  189. fld  qword ptr [esi+7*8*8]
  190. fmul qword ptr [eax+7*8*8]
  191. fadd
  192. fadd st(0),st(1)
  193. //add eax,8
  194. //fstp qword ptr [edi]
  195. //add edi,8*8
  196. fxch st(1)
  197. //
  198. fld  qword ptr [esi+0*8*8]
  199. fmul qword ptr [eax+0*8*8+8]
  200. fld  qword ptr [esi+1*8*8]
  201. fmul qword ptr [eax+1*8*8+8]
  202. fadd
  203. fld  qword ptr [esi+2*8*8]
  204. fmul qword ptr [eax+2*8*8+8]
  205. fadd
  206. fld  qword ptr [esi+3*8*8]
  207. fmul qword ptr [eax+3*8*8+8]
  208. fadd
  209. fld  qword ptr [esi+4*8*8]
  210. fmul qword ptr [eax+4*8*8+8]
  211. fadd
  212. fld  qword ptr [esi+5*8*8]
  213. fmul qword ptr [eax+5*8*8+8]
  214. fadd
  215. fld  qword ptr [esi+6*8*8]
  216. fmul qword ptr [eax+6*8*8+8]
  217. fadd
  218. fld  qword ptr [esi+7*8*8]
  219. fmul qword ptr [eax+7*8*8+8]
  220. fadd
  221. fadd st(0),st(1)
  222. add eax,8*2
  224. fxch st(2)
  225. fstp qword ptr [edi]
  226. fxch st(1)
  227. fstp qword ptr [edi+8*8]
  228. add edi,8*8*2
  229. /*
  230. fistp dword ptr [int0]
  231. mov edx,dword ptr [int0]
  232. mov dx,word ptr [iclip+2*1024+2*edx]
  233. mov word ptr [edi],dx
  234. add edi,8*2
  235. */
  236. dec ecx
  237. // sub ecx,0x80000001
  238. // js  __col2
  239.   //align 16
  240. // test ecx,ecx // align jump &| redo flags
  241. jnz __col2
  242. add eax,-8*8
  243. add esi,+8
  244. add edi,8-8*8*8
  245. //add edi,2-8*8*2
  246. //dec ebx
  247. sub ebx,0x80000001
  248. js  __row2
  249.   //align 16
  250. test ebx,ebx // align jump &| redo flags
  251. jnz __row2
  252. ffree st(0) // bye bye 0.5
  254.   // set x87 to floor mode
  255. fstcw [fpold]
  256. movzx eax, [fpold]
  257. //or ax, 0x0C00 // 3072
  258. //or ax, 0x0800 // round up
  259. or eax, 0x0400 // round down - floor
  260. mov [fpnew], ax
  261. fldcw [fpnew]
  263. // now floor the damn array
  264. lea esi, [rnd]
  265. mov edi, [block]
  266. mov ebx, -256 // clip min
  267. mov edx, +255 // clip max
  268. mov ecx, 8
  269. align 16
  270. __floor:
  271.   fld   qword ptr [esi+0*8]
  272. fistp dword ptr [int0]
  273.   mov eax,[int0]
  274. cmp   eax,ebx
  275. cmovl eax,ebx
  276. cmp   eax,edx
  277. cmovg eax,edx
  278.   fld   qword ptr [esi+1*8]
  279. fistp dword ptr [int1]
  280. mov word ptr [edi+0*2],ax
  281.   mov eax,[int1]
  282. cmp   eax,ebx
  283. cmovl eax,ebx
  284. cmp   eax,edx
  285. cmovg eax,edx
  286.   fld   qword ptr [esi+2*8]
  287. fistp dword ptr [int2]
  288. mov word ptr [edi+1*2],ax
  289.   mov eax,[int2]
  290. cmp   eax,ebx
  291. cmovl eax,ebx
  292. cmp   eax,edx
  293. cmovg eax,edx
  294.   fld   qword ptr [esi+3*8]
  295. fistp dword ptr [int3]
  296. mov word ptr [edi+2*2],ax
  297.   mov eax,[int3]
  298. cmp   eax,ebx
  299. cmovl eax,ebx
  300. cmp   eax,edx
  301. cmovg eax,edx
  302.   fld   qword ptr [esi+4*8]
  303. fistp dword ptr [int4]
  304. mov word ptr [edi+3*2],ax
  305.   mov eax,[int4]
  306. cmp   eax,ebx
  307. cmovl eax,ebx
  308. cmp   eax,edx
  309. cmovg eax,edx
  310.   fld   qword ptr [esi+5*8]
  311. fistp dword ptr [int5]
  312. mov word ptr [edi+4*2],ax
  313.   mov eax,[int5]
  314. cmp   eax,ebx
  315. cmovl eax,ebx
  316. cmp   eax,edx
  317. cmovg eax,edx
  318.   fld   qword ptr [esi+6*8]
  319. fistp dword ptr [int6]
  320. mov word ptr [edi+5*2],ax
  321.   mov eax,[int6]
  322. cmp   eax,ebx
  323. cmovl eax,ebx
  324. cmp   eax,edx
  325. cmovg eax,edx
  326.   fld   qword ptr [esi+7*8]
  327. fistp dword ptr [int7]
  328. mov word ptr [edi+6*2],ax
  329.   mov eax,[int7]
  330. cmp   eax,ebx
  331. cmovl eax,ebx
  332. cmp   eax,edx
  333. cmovg eax,edx
  334. mov word ptr [edi+7*2],ax
  336. add esi, 8*8
  337. add edi, 8*2
  338. //dec ecx
  339. sub ecx,0x80000001
  340. js  __floor
  341.   //align 16
  342. test ecx,ecx // align jump &| redo flags
  343. jnz __floor
  345. // set x87 to default mode
  346. fldcw [fpold]
  347. };
  348. }