고속 DCT 문제 해결

[한창호]

고속 DCT를 실행하면 다음과 같이 오히려 fast가 더 늦게 나온다.

문제점은 
1. cos 계산은 많이 줄었는데, 현재 컴퓨터가 산술 연산을 월등히 빨리 처리해서 효과가 적다.
2. at 함수 사용이 좀더 늘었는데 그로 인해 오히려 시간이 더걸렸다.

따라서 at -> ptr로 변경해서 시도했더니 겨우 빨라졌다.

at 함수를 ptr 함수로 변경한 소스코드

Mat DCT_block(Mat& g, Mat& cos_mat, int dir)
{
  Mat dst(g.size(), g.type(), Scalar(0));
  int N = g.rows, M = g.cols;
  Vec2f CK(sqrt(1.0f / N), sqrt(2.0f / N));
  Vec2f CL(sqrt(1.0f / M), sqrt(2.0f / M));

  for (int k = 0; k < N; k++)
  {
    for (int l = 0; l < M; l++)
    {
      float ck = (k) ? CK[1] : CK[0];
      float cl = (l) ? CL[1] : CL[0];
      Point idx(l * M, k * N);

      for (int n = 0; n < N; n++)
      {
        for (int m = 0; m < M; m++)
        {
          //float cos_value = cos_mat.at<float>(idx + Point(m, n));
          //float coef = (dir == 0) ? g.at<float>(n, m) : g.at<float>(k, l);
          //if (dir == 0) dst.at<float>(k, l) += ck * cl * coef * cos_value;//정변환
          //else dst.at<float>(n, m) += ck * cl * coef * cos_value;//역변환

          float cos_value = *(cos_mat.ptr<float>(idx.y + n) + idx.x + m);
          float coef = (dir == 0) ? *(g.ptr<float>(m) + n) : *(g.ptr<float>(l) + k) ;
          if (dir == 0) *(dst.ptr<float>(l) + k) += ck * cl * coef * cos_value;//정변환
          else *(dst.ptr<float>(m) +n) += ck * cl * coef * cos_value;//역변환          
        }
      }
    }
  }
  return dst;
}