自定义3D(并行）FFTW变换

我有一个3D真实数据阵列，我想沿一个轴沿一个轴和正常的DFT在其他两个轴上转换。结果应该是一个3D复合阵列，该阵列保持此转换系数。

您知道FFTW3软件包是否提供了这样的例行程序 - 可能是并行的 - 开箱即用？FFTW3在所有三个方向上为简单的3D DFT提供了这样的例程。

，如果没有，您可能会暗示如何实现C/C 中的最佳方法？

我的天真想法：DST/DCT的组装，然后是沿着某些包装程序内的第一个轴的2D真实复合变换。然后，人们可以想到实现并行性的一维分解。2D会更好，但工作量更多。

ps：此转换用于求解Navier-Stokes方程的光谱方法。