计算科学 - 从“fftw_plan”中提取旋转因子的实部 - 吾爱随笔录

使用 FFTW 计算 II 型 DCT 时，您创建一个fftw_plan, 通过

fftw_plan plan = fftw_plan_r2r_1d(size, in_ptr, out_ptr, FFTW_REDFT10, FFTW_ESTIMATE);

据我了解，该计划预先计算了旋转因子，以便该计划可以有效地重复用于多次调用fftw_execute_r2r(plan, in_ptr, out_ptr).

因为我试图计算

\begin{aligned} c_{k} = \frac{1}{N} \sum_{j = 0}^{N - 1} f (\cos (\frac{(j + 1 / 2) π}{N})) \cos (\frac{(j + 1 / 2) k π}{N}) \end{aligned}

$\begin{align*} c_{k} = \frac{1}{N} \sum_{j=0}^{N-1} f\left(\cos\left(\frac{(j+1/2)\pi}{N}\right)\right) \cos\left(\frac{(j+1/2)k\pi}{N}\right) \end{align*}$ 我认为fftw_plan已经计算了参数

f

$f$ ，所以我不想重新计算它。

你如何从中提取旋转因子fftw_plan？