mkl库DCT函数实现问题

我想通过使用mkl.lib DCT函数替换ipp.lib DCT函数。我用下列基于伊普API做DCT：mkl库DCT函数实现问题

ippsDCTFwdInitAlloc_32f 
ippsDCTFwd_32f 
ippsDCTFwdFree_32f

PFB我MKL基于DCT的示例代码：

#include <time.h> 
#include <stdlib.h> 
#include "mkl.h" 
int main(int argc, char* argv[]){ 

    float *dpar; 
    float *out; 
    MKL_INT *ipar; 
    MKL_INT tt_type,stat,n_1,nn; 
    FILE *fp,*fw,*fonce; 
    fp = fopen("D:\\dump\\fileinput.txt","r"); 
    if(fp == NULL){ 
     cout<<"file not created properly"<<endl; 
    } 
    DFTI_DESCRIPTOR_HANDLE handle = 0; 
    int n = 150; /*Hardcoded to run for my code TODO: going to change after integrating into my main codebase.As I can't attach my full fileinput.txt .So just taken n value = 150 */ 
    nn = (MKL_INT)n; 
    tt_type = MKL_STAGGERED_COSINE_TRANSFORM; 

    n_1 = nn + 1 ; 
    out = (float*)malloc((n+1)*sizeof(float)); 
    dpar= (float*)malloc((5*n_1/2+2)*sizeof(float)); 
    ipar= (MKL_INT*)malloc((128)*sizeof(int)); 
    s_init_trig_transform(&n_1,&tt_type,ipar,dpar,&stat); 
    for (int srcSize =0 ;srcSize< n ; srcSize++) 
    { 
     fscanf(fp,"%f\n",&out[srcSize]); 
    } 
    fclose(fp); 
    if (stat != 0) 
    { 
     printf("\n============================================================================\n"); 
     printf("FFTW2MKL FATAL ERROR: MKL TT initialization has failed with status=%d\n",(MKL_INT)stat); 
     printf("Please refer to the Trigonometric Transform Routines Section of MKL Manual\n"); 
     printf("to find what went wrong...\n"); 
     printf("============================================================================\n"); 
     return NULL; 
    } 
    ipar[10] = 1; //nx, that is, the number of intervals along the x-axis, in the Cartesian case. 
    ipar[14] = n_1; //specifies the internal partitioning of the dpar array. 
    ipar[15] = 1; //value of ipar[14]+1,Specifies the internal partitioning of the dpar array. 
    s_commit_trig_transform(out,&handle,ipar,dpar,&stat); 
    if (stat != 0) 
    { 
     printf("\n============================================================================\n"); 
     printf("FFTW2MKL FATAL ERROR: MKL TT commit step has failed with status=%d\n",(MKL_INT)stat); 
     printf("Please refer to the Trigonometric Transform Routines Section of MKL Manual\n"); 
     printf("to find what went wrong...\n"); 
     printf("============================================================================\n"); 
     return NULL; 
    } 
    s_forward_trig_transform(out,&handle,ipar,dpar,&stat); 
    if (stat != 0) 
    { 
     printf("\n============================================================================\n"); 
     printf("FFTW2MKL FATAL ERROR: MKL TT commit step has failed with status=%d\n",(MKL_INT)stat); 
     printf("Please refer to the Trigonometric Transform Routines Section of MKL Manual\n"); 
     printf("to find what went wrong...\n"); 
     printf("============================================================================\n"); 
     return NULL; 
    } 
    free_trig_transform(&handle,ipar,&stat); 
    printf("\n===== DCT GOT OVER ======== \n"); 

    return 0; 

}

但我没有得到正确的输出。我感到困惑，我的执行是否正确。

因为我不知道如何附加我的fileinput.txt，所以我把文件内容放在这里。对不起，该:(

129.000000 128.000000 145.000000 140.000000 102.000000 93.000000 100.000000 97.000000 94.000000 88.000000 58.000000 67.000000 
79.000000 62.000000 62.000000 64.000000 63.000000 63.000000 70.000000 80.000000 61.000000 62.000000 63.000000 64.000000 
62.000000 78.000000 66.000000 65.000000 68.000000 74.000000 78.000000 66.000000 67.000000 76.000000 67.000000 65.000000 
76.000000 83.000000 66.000000 66.000000 65.000000 66.000000 72.000000 80.000000 65.000000 66.000000 69.000000 80.000000 
75.000000 66.000000 71.000000 80.000000 69.000000 80.000000 85.000000 70.000000 77.000000 68.000000 67.000000 80.000000 
83.000000 69.000000 71.000000 73.000000 69.000000 70.000000 71.000000 71.000000 70.000000 71.000000 88.000000 86.000000 
78.000000 78.000000 70.000000 73.000000 81.000000 83.000000 76.000000 75.000000 72.000000 72.000000 73.000000 75.000000 
79.000000 76.000000 80.000000 87.000000 84.000000 79.000000 73.000000 72.000000 72.000000 72.000000 73.000000 79.000000 
76.000000 74.000000 76.000000 78.000000 76.000000 74.000000 73.000000 76.000000 78.000000 78.000000 73.000000 74.000000 
73.000000 72.000000 72.000000 73.000000 75.000000 77.000000 76.000000 75.000000 74.000000 74.000000 75.000000 78.000000 
81.000000 78.000000 75.000000 74.000000 74.000000 74.000000 75.000000 77.000000 79.000000 80.000000 81.000000 78.000000 
76.000000 75.000000 74.000000 74.000000 74.000000 75.000000 76.000000 76.000000 75.000000 77.000000 78.000000 79.000000 
80.000000 79.000000 79.000000 79.000000 78.000000 76.000000 76.000000 76.000000 76.000000 77.000000 77.000000 77.000000 
78.000000 80.000000 80.000000 80.000000 80.000000 80.000000 78.000000 78.000000 78.000000 78.000000 78.000000 80.000000

来源

2014-12-23 Nilay Srivastava

我会检查这个。但至少运行你的代码，因为它是，您应该使用“COUT”或“printf”上。或者添加＃“的iostream” using namespace std; – Ashwin

我觉得你说的这些函数的参数DPAR。实际上，它们不被用于计算的结果。他们只是用来存储中间正弦/余弦数据在TT计算用于

对于运算功能： d_forward_trig_transform（双F []，DFTI_DESCRIPTOR_HANDLE *手柄，MKL_INT IPAR []，双DPAR []，MKL_INT * STAT）;

两者输入。和输出是F []，它是相同的大小以获取更多信息，您可以参考以下链接： link here

来源

2014-12-24 06:19:30 Ashwin

它帮助我理解了一些观点，但它仍然不能解决问题！ –

mkl库DCT函数实现问题

回答

相关问题