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

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 多媒体编程 > 查看源码
idctref.cpp
资源名称：fullMPCcode.rar　[点击查看]
上传用户：xjjlds
上传日期：2015-12-05
资源大小：22823k
文件大小：8k
源码类别：

多媒体编程

开发平台：

Visual C++

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