FFTW和智慧的签名问题

问题描述 投票：0回答：1

我在Centos 7 Linux目标上使用FFTW 3.3.8。

[第一次测试，我使用FFTW执行一些计算，结果很好，并且我使用fftw_export_wisdom_to_filename保存计划。

在第二个测试中，与第一个测试相同，我刚刚在计划的计算中添加了对fftw_import_wisdom_from_filename和属性FFTW_WISDOM的调用，它也可以正常工作。

[当我在智慧文件上使用fftw-wisdom-to-conf，编译生成的输出，并将其添加到测试二进制文件中时，就会出现问题（在这种情况下，二进制文件的大小会减小，这对我来说很有趣）。

在这种情况下，导入智慧功能由于签名错误而返回“ 0”。智慧文件中的前四个md5值不正确（与预期值不同）。

如果我在智慧文件中修补这些值，则一切工作正常，可以正常导入，并且还产生FFTW。没有激活gcc优化选项，

我不使用多线程，我的配置看起来很简单，您有什么想法可以解释这种行为吗？

用于生成FFTW库的命令如下：./configure CC=gcc CFLAGS="-g" --with-slow-timer --disable-threads --disable-fortran --host=x86_64-pc-linux-gnu --prefix=/lvcugen/workspace/workspace_ext6/GR712/TMP_FFTW/INSTALL_LINUX

用于执行此测试并从test_fftw_wisdom.c安装目录编译的fftw文件的内容如下：

/*
 * P.SERAUD : FFTW Wisdom test (define MODE_USE to generate wisdom load binary, default binary generate wisdom
 *
 * Gen wisdom build command : gcc -u main -o0 -g test_fftw_wisdom.c -I ./include -L ./lib -lfftw3 -lm -lc -o test_fftw_gen_wisdom
 * Use wisdom build command : gcc -u main -o0 -g -DMODE_USE test_fftw_wisdom.c conf_wisdom.o -I ./include -L ./lib -lfftw3 -lm -lc -o test_fftw_use_wisdom
 *
 * conf_wisdom file generated by : ./bin/fftw-wisdom-to-conf < wisdom.txt
 *
 */
#include <stdio.h>
#include <string.h>

#include "fftw3.h"


/* in : input image 
 * reOut : real part (output)
 * imOut : imaginary part (output)
 * largeur : largeur des images d'entrée et de sortie
 * hauteur : idem
 */
void fourierForward(const double* in,
                    double* reOut,
                    double* imOut,
                    unsigned int largeur,
                    unsigned int hauteur)
{

/* The data is an array of type 'fftw_complex', which is by default a
   'double[2]' composed of the real ('in[i][0]') and imaginary ('in[i][1]') 
     parts of a complex number.
     */

   fftw_complex* spatial_repr;
   fftw_complex* frequency_repr;
   unsigned int i;
   fftw_plan plan;
   fftw_plan plan2;
     int resWisdom=0;

     spatial_repr= fftw_malloc(sizeof(fftw_complex)*largeur*hauteur);
   frequency_repr= fftw_malloc(sizeof(fftw_complex)*largeur*hauteur);

     fftw_set_timelimit(FFTW_NO_TIMELIMIT);

   /* Use Wisdom mode called when MODE_USE is set  */
#ifdef MODE_USE
     printf ("Importing Wisdom %d\n",0);
     resWisdom = fftw_import_wisdom_from_filename("wisdom.txt");

     printf ("Import result %d\n",resWisdom);

     if(resWisdom == 0)
         return;
#endif

   printf("Launching fftw_plan_dft_2d FORWARD with %d / %d\n",largeur, hauteur);
     /*on calcule le plan d'exécution*/
#ifdef MODE_USE
   plan=fftw_plan_dft_2d(hauteur, largeur, spatial_repr, frequency_repr, FFTW_FORWARD, FFTW_WISDOM_ONLY | FFTW_ESTIMATE);
#else
   plan=fftw_plan_dft_2d(hauteur, largeur, spatial_repr, frequency_repr, FFTW_FORWARD, FFTW_ESTIMATE);
#endif

     if (plan == NULL)
        { 
    printf("NULL plan generated by fftw_plan_dft_2d with %d / %d !!!\n",largeur, hauteur);
        return;
        }

   /*On remplit la structure qui sera utilisée par fftw*/
   for(i=0;i<largeur*hauteur;i++)
   {
          spatial_repr[i][0] = in[i]; /* real part */
      spatial_repr[i][1] =  0.0;  /* imaginary part */
      frequency_repr[i][0] =  0.0;  /* real part */
      frequency_repr[i][1] =  0.0;  /* imaginary part */
   }

     printf("Executing fftw_execute with %d / %d\n\n",largeur, hauteur);

   /*on effectue la transformée de Fourier*/
   fftw_execute(plan);


   /*on retranscrit le résultat en 2 images, l'une représentant la partie réelle, l'autre
     la partie imaginaire*/
   for(i=0;i<largeur*hauteur;i++)
      {
          reOut[i]=frequency_repr[i][0];
          imOut[i]=frequency_repr[i][1];
                    printf("%d : in (%lf) / reOut (%lf) / imOut (%lf) \n",i, in[i],reOut[i],imOut[i]);
            }

   printf("Launching fftw_plan_dft_2d BACKWARD with %d / %d\n",largeur, hauteur);
#ifdef MODE_USE
   plan2=fftw_plan_dft_2d(hauteur, largeur, spatial_repr, frequency_repr, FFTW_BACKWARD, FFTW_WISDOM_ONLY | FFTW_ESTIMATE);
#else
   plan2=fftw_plan_dft_2d(hauteur, largeur, spatial_repr, frequency_repr, FFTW_BACKWARD, FFTW_ESTIMATE);
#endif

     if (plan2 == NULL)
        { 
    printf("NULL plan generated by fftw_plan_dft_2d with %d / %d !!!\n",largeur, hauteur);
        return;
        }

   /*On remplit la structure qui sera utilisée par fftw*/
   for(i=0;i<largeur*hauteur;i++)
   {
     spatial_repr[i][0] = reOut[i]; /* real part */
     spatial_repr[i][1] = imOut[i]; /* imaginary part */
     frequency_repr[i][0] =  0.0;  /* real part */
     frequency_repr[i][1] =  0.0;  /* imaginary part */
   }

     printf("Executing fftw_execute with %d / %d\n\n",largeur, hauteur);
   /*on effectue la transformée de Fourier*/
   fftw_execute(plan2);

   /*on retrouve la matrice d'origine dans la partie réelle en la divisant par "largeur * hauteur" */ 
   for(i=0;i<largeur*hauteur;i++)
      {
          reOut[i]=frequency_repr[i][0] /(largeur*hauteur);
          imOut[i]=frequency_repr[i][1];
                    printf("%d : in (%lf) / reOut (%lf) / imOut (%lf) \n",i, in[i],reOut[i],imOut[i]);
            }

#ifndef MODE_USE
     resWisdom = fftw_export_wisdom_to_filename("wisdom.txt"); 
     printf ("Wisdom export res = %d\n", resWisdom);
#endif

   /*on détruit les objets*/
   fftw_destroy_plan(plan);
   fftw_destroy_plan(plan2);
   fftw_free(spatial_repr);
   fftw_free(frequency_repr);

}

int main(void) 
{
        printf("Starting !!!!\n");

        double in [8][8] =    {{0.0, 0.0,   0.0,   0.0,   0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0,  70.0,  80.0,  90.0, 0.0, 0.0, 0.0},
            {0.0, 0.0,  90.0, 100.0, 110.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 110.0, 120.0, 130.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 130.0, 140.0, 150.0, 0.0, 0.0, 0.0},
            {0.0, 0.0,   0.0,   0.0,   0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0,   0.0,   0.0,   0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0,   0.0,   0.0,   0.0, 0.0, 0.0, 0.0}};

        double reOut [8][8] = {{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}};

        double imOut [8][8] = {{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
            {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0}};

        fourierForward(&in[0][0], &reOut[0][0], &imOut[0][0], 8, 8);

        return 0;
}

conf_wisdom.o可执行文件中包含的test_fftw_use_wisdom目标文件是通过conf_wisdom.c编译gcc -c conf_wisdom.c源文件的结果。该源文件的内容如下：

/* Automatically generated by fftw-wisdom-to-conf from fftw 3.3.8.
   DO NOT EDIT!  (Unless you really, really want to.  Then it's okay.) */
void fftw_configure_planner(void *plnr)
{
    struct solvtab_s { void (*reg)(void *); const char *reg_nam; };
    extern void fftw_solvtab_exec(const struct solvtab_s s[], void *);

#define DECLARE(name) extern void name(void *);
#define STRINGIZEx(x) #x
#define STRINGIZE(x) STRINGIZEx(x)
#define SOLVTAB(s) { s, STRINGIZE(s) },
#define DO(X) \
    X(fftw_codelet_n1_8)\
    X(fftw_dft_rank_geq2_register)\
    /* end DO(X) */

    DO(DECLARE)

    const struct solvtab_s s[] = {
        DO(SOLVTAB)
        { 0, 0 }
    };

    fftw_solvtab_exec(s, plnr);
}

通过调用函数fftw_export_wisdom_to_filename生成的witness.txt文件的内容如下：

(fftw-3.3.8 fftw_wisdom #x4be12fff #x7b2df9b2 #xa5975329 #x385b0041
  (fftw_codelet_n1_8 0 #x31bff #x31bff #x0 #x2874cdd9 #x63ab3b26 #xcb8c1374 #xedc69a82)
  (fftw_codelet_n1_8 0 #x31bff #x31bff #x0 #x4e583e34 #x14aedc27 #x54a2fad7 #xd124b261)
  (fftw_dft_rank_geq2_register 0 #x31bff #x31bff #x0 #xcd87e2bd #x8b6dd8ee #x2d70ed77 #xd20ea4af)
  (fftw_codelet_n1_8 0 #x31bff #x31bff #x0 #xece23a64 #x64207145 #x3f794c2c #x32ef2fbf)
  (fftw_dft_rank_geq2_register 0 #x31bff #x31bff #x0 #xd416e586 #xdc5a4e14 #x2b04907f #x69b22bd8)
  (fftw_codelet_n1_8 0 #x31bff #x31bff #x0 #x155aa639 #x5f99241d #x6351b1a0 #xc7ccb78d)
)

未成功检查的错误md5值是＃x4be12fff＃x7b2df9b2＃xa5975329＃x385b0041。如果将这些值替换为期望的值，则一切正常。

此外，当库的全部内容包含在test_fftw_use_wisdom二进制文件中（不使用conf_wisdom.o文件）时，它也可以正常工作。

实际上，问题出在导入和导出功能都使用的功能signature_of_configuration中。当包含conf_wisdom.o时，就完整列表而言，planner内容（与nslvdesc和slvdescsiz关联的解析器列表）的内容将减少，并且关联的md5值必须有所不同。

c fftw
1个回答
0
投票

测试表明，当从智慧数据生成的wisdom.c既包含在导入程序中，又包含在导出程序中时，导入就起作用了。当两个程序中均不包含wisdom.c时，它也适用。如果将它包含在导入程序中但不包含在导出程序中，则它将失败。

因此，解决方法是在不使用wisdom.c的情况下构建导出程序，导出智慧，生成wisdom.c，使用wisdom.c重建导出程序，重新导出智慧，然后使用wisdom.c生成导入程序。

这似乎是FFTW的缺陷，可能在校验和检查的设计中。

最新问题
© www.soinside.com 2019 - 2024. All rights reserved.