/* fdctref.c, 前向離散余弦變換，雙精度*/#include <math.h>#include "config.h"#ifndef PI# ifdef M_PI#  define PI M_PI# else#  define PI 3.14159265358979323846# endif#endif/* 全局聲明 */void init_fdct _ANSI_ARGS_((void));void fdct _ANSI_ARGS_((short *block));static double c[8][8]; /* 變換系數 */void init_fdct(){  int i, j;  double s;  for (i=0; i<8; i++)  {    s = (i==0) ? sqrt(0.125) : 0.5;    for (j=0; j<8; j++)      c[i][j] = s * cos((PI/8.0)*i*(j+0.5));  }}void fdct(block)short *block;{  int i, j, k;  double s;  double tmp[64];  for (i=0; i<8; i++)    for (j=0; j<8; j++)    {      s = 0.0;      for (k=0; k<8; k++)        s += c[j][k] * block[8*i+k];      tmp[8*i+j] = s;    }  for (j=0; j<8; j++)    for (i=0; i<8; i++)    {      s = 0.0;      for (k=0; k<8; k++)        s += c[i][k] * tmp[8*k+j];      block[8*i+j] = (int)floor(s+0.499999);      /*       這里需要說明一下，為什么采用0.499999而不采用0.5？
	   這是因為s的取值經常是以.5結尾的（至少在i，j等于
	   0或4的時候是這樣），如果將閾值設為0.5將導致非常
	   高的算術誤差，而如果采用0.499999則產生誤差的概率
	   也就是s以.5和.500001之間的某個數值結尾的概率，顯
	   然，這樣的幾率要小很多。       */    }}