Я пытаюсь заменить функцию DCT ipp.lib с помощью функции DCT mkl.lib.
Я использовал следующий API на основе Ipp для выполнения DCT:
ippsDCTFwdInitAlloc_32f
ippsDCTFwd_32f
ippsDCTFwdFree_32f
Пример кода DCT на основе MKL для моего PFB:
#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
Я думаю, что вы говорите параметр dpar для этих функций.
На самом деле, они не используются для вычисленного результата. Они просто используются для хранения промежуточных данных sin / cos, которые используются в вычислениях TT.
Для функции вычисления:
d_forward_trig_transform (двойной f [], DFTI_DESCRIPTOR_HANDLE * handle, MKL_INT ipar [], двойной dpar [], MKL_INT * stat);
И вход и выход — F [], который имеет одинаковый размер.
Для получения дополнительной информации вы можете обратиться по ссылке ниже:
ссылка здесь