Question

我正在努力为FFTW操作实现更好的性能。出于这个原因，我决定使用智慧文件来创建计划，但问题是它无法从智慧文件中加载计划（导出智慧很有效）。我试图将智慧导出到一个文件，并在下一个程序运行时将其加载回智慧文件，但该函数总是为任何文件名返回0（即使对于非现有文件名）。我也尝试从字符串加载但它也不起作用。

这里似乎有类似的问题，但没有一个问题得到解答或者问题出在其他地方。那么这是图书馆的一个错误，还是我做错了什么？

修改仅用于演示并查看导出功能是否真的有效：下面的代码显示了使用文件进行内容导入从字符串加载智慧的操作（仅显示内容在文件中可用）：

FILE * pfile;
vector<char> buffer;
pfile=fopen("WisdomFile.txt","r");

if(pfile==0)
    cout<<"Could not open file"<<endl;
else
    cout<<"Could open file successfully"<<endl;

long length;

if(pfile)
{
    fseek(pfile,0,SEEK_END);
    length=ftell(pfile);
    fseek(pfile,0,SEEK_SET);
    buffer.assign(length+1,'\0'); //allocate space with the same length as the file
    int n=fread(&buffer[0],1,length,pfile); //read whole file to the buffer
    assert(n==length);
    fclose(pfile);
}
string show_wisdom(buffer);
cout<<show_wisdom<<endl;    //content could be read

int ret=fftwf_import_wisdom_from_string(reinterpret_cast<const char*>(&buffer[0])); // returns 0 for every filename
buffer.clear();
//... further FFTW code trying to use `FFTW_PATIENT | FFTW_USE_WISDOM_ONLY` ->but never uses the wisdom file

输出如下：

    (fftw-3.3.4 fftwf_wisdom #xbedb7e38 #x1ac524dc #x7a69378e #x21629161
  (fftwf_codelet_t2_8 3 #x11048 #x11048 #x0 #xa75017ef #xb6eb4747 #x4bef8a59 #xb03d9427)
  (fftwf_dft_thr_vrank_geq1_register 1 #x11048 #x11048 #x0 #x27a0c32d #x4e3441f9 #xb3fb3f2d #x90ae8374)
  (fftwf_dft_vrank_geq1_register 0 #x11048 #x11048 #x0 #x9be02645 #x53c7643d #xf6cf9608 #xed5460b7)
  (fftwf_dft_r2hc_register 0 #x11048 #x11048 #x0 #x52a71bc4 #x3c83e70d #x942dd977 #xf047f7e9)
  (fftwf_codelet_n1_64 0 #x11448 #x11448 #x0 #x11559ac4 #xea86db86 #xad6ae8e4 #x97f477c6)
  (fftwf_codelet_t1_16 0 #x11048 #x11048 #x0 #x8811820f #xea00b698 #x861ae7ed #x109ec45a)
  (fftwf_rdft_rank0_register 2 #x11048 #x11048 #x0 #x0095ff64 #x86e47338 #x76e9cf55 #x6cde6434)
  (fftwf_codelet_t1_16 1 #x11048 #x11048 #x0 #x29eda2bf #x97038fb2 #x0eddb089 #xafc2b57e)
  (fftwf_dft_indirect_register 0 #x11048 #x11048 #x0 #x1bea55f5 #x48417896 #x04bc4c58 #x571ce0b9)
  (fftwf_dft_thr_vrank_geq1_register 0 #x11048 #x11048 #x0 #x7b53c8cd #xda17faa2 #x220c1322 #x7c207bbd)
)

因此，如上所示，导出函数似乎有效，但是如FFTW教程中所示导入它似乎没有效果（程序仍尝试使用FFTW_PATIENT选项创建一个新计划，大约需要5分钟。

Answer 1

根据FFTW's documentation：

智慧被自动用于任何适用的大小，只要规划者标志不比那些创造智慧的人更“耐心”。

换句话说，如果您想为使用FFTW_PATIENT标志创建的计划导入智慧，那么必须为使用相同大小和FFTW_PATIENT或FFTW_EXHAUSTIVE旗帜。

然后，序列将首先提前创建一个计划，然后通过以下方式输出智慧：

// Might as well use FFTW_EXHAUSTIVE if we can afford it ahead of time
fftw_plan plan = fftw_plan_dft_r2c_1d(N, in, out, FFTW_EXHAUSTIVE);
fftw_export_wisdom_to_filename("WisdomFile.txt");

最后，导出的智慧可以导入并用于：

if (!fftw_import_wisdom_from_filename("WisdomFile.txt"))
{
  std::cout << "Warning: could not import wisdom file" << std::endl;
}
// plan can now be created with any flags less or equal to 
// FFTW_EXHAUSTIVE used for the exported wisdom.
fftw_plan plan = fftw_plan_dft_r2c_1d(N, in, out, FFTW_PATIENT);

Answer 2

我遇到了这个问题，直到我注意到我使用函数void fftw_set_timelimit(double seconds)来限制总计划创建时间。无法加载使用此时间限制生成的智慧，而可以加载在禁用时间限制时生成的智慧。即使没有超过时限，情况似乎也是如此。

由于您的代码已经过编辑，我不确定这是否适用于您的具体问题，但其他人可能会认为这有用。

FFTW导出智慧文件但无法加载它

2 个答案: