如何在Fortran代码中引入cufft.h头文件并解决nvfortran编译报错问题

首先，你遇到的问题本质是直接在Fortran代码中包含C头文件导致的语法冲突——nvfortran虽然支持预处理，但C头文件的语法（//注释、续行方式、文件引用逻辑）和Fortran自由格式的规则不匹配，手动修改头文件显然不是长久之计，这里给你几个更优雅的解决方案：

1. 不要直接包含cufft.h，改用Fortran接口封装 #

cuFFT本身没有官方提供Fortran模块，但你可以通过Fortran的iso_c_binding手动封装需要的cuFFT函数，彻底绕开C头文件的语法问题。

示例：封装基础cuFFT接口 #

创建一个cufft_interface.f90模块文件，只封装你代码中用到的cuFFT函数（比如计划创建、执行FFT、销毁计划等）：

module cufft_interface
  use iso_c_binding
  implicit none

  ! 对应C中的cufftHandle类型
  type, bind(c) :: cufftHandle
    integer(c_int) :: handle
  end type cufftHandle

  ! 定义cuFFT的常量（从cufft.h中提取对应值）
  integer(c_int), parameter :: CUFFT_C2C = 1
  integer(c_int), parameter :: CUFFT_FORWARD = -1
  integer(c_int), parameter :: CUFFT_INVERSE = 1

  ! 封装cufftPlan1d函数
  interface
    integer(c_int) function cufftPlan1d(plan, nx, type, batch) bind(c, name='cufftPlan1d')
      import :: cufftHandle, c_int
      type(cufftHandle), intent(out) :: plan
      integer(c_int), value :: nx, type, batch
    end function cufftPlan1d

    ! 封装cufftExecC2C
    integer(c_int) function cufftExecC2C(plan, idata, odata, direction) bind(c, name='cufftExecC2C')
      import :: cufftHandle, c_int
      type(cufftHandle), value :: plan
      complex(c_float_complex), intent(inout) :: idata(*), odata(*)
      integer(c_int), value :: direction
    end function cufftExecC2C

    ! 封装cufftDestroy
    integer(c_int) function cufftDestroy(plan) bind(c, name='cufftDestroy')
      import :: cufftHandle, c_int
      type(cufftHandle), value :: plan
    end function cufftDestroy
  end interface

end module cufft_interface

之后在你的主代码中，只需要use cufft_interface就可以调用cuFFT函数，完全不用接触C头文件：

program fft_gpu_demo
  use cufft_interface
  use openacc
  implicit none

  type(cufftHandle) :: plan
  complex(c_float_complex), allocatable :: data(:)
  integer(c_int) :: nx, status

  nx = 1024
  allocate(data(nx))

  ! 创建FFT计划
  status = cufftPlan1d(plan, nx, CUFFT_C2C, 1)
  if (status /= 0) error stop "Failed to create cuFFT plan"

  ! 用OpenACC将数据拷贝到GPU
  !$acc enter data copyin(data)

  ! 执行正向FFT
  status = cufftExecC2C(plan, data, data, CUFFT_FORWARD)
  if (status /= 0) error stop "FFT execution failed"

  ! 处理结果...

  ! 清理资源
  !$acc exit data delete(data)
  status = cufftDestroy(plan)

end program fft_gpu_demo

2. 利用nvfortran自带的cuFFT兼容模块 #

部分版本的nvfortran提供了预封装的cuFFT Fortran模块（比如cufft.f90），你可以检查CUDA安装目录下的include文件夹（通常是/usr/local/cuda/include），如果有这个文件，直接在代码中use cufft即可，编译时不需要手动包含头文件，只需链接-lcufft。

3. 解决头文件路径问题（如果必须包含C头文件） #

如果你因为某些原因一定要包含cufft.h，首先要确保编译时通过-I指定CUDA的include目录，比如：

nvfortran -cpp -Mfree -I/usr/local/cuda/include -lcufft your_code.f90 -o your_code

这样nvfortran就能正确找到cufft.h及其依赖的头文件，不需要手动修改头文件中的#include路径。

对于//注释和续行错误，可以添加-CC选项让nvfortran保留C风格注释，同时用-xc强制预处理阶段按C语法解析头文件，但这种方法还是容易出现格式冲突，不如接口封装可靠。

为什么会出现"续行标签字段不为空"错误？ #

即使你用了-Mfree（自由格式），nvfortran在处理C头文件时，会把C的续行符（\）或者行首的非空白字符误判为Fortran固定格式的行标签（固定格式下第1-5列是标签区），这就是为什么你添加7个空格能解决第1行，但后面的行还是出错——C头文件的行结构完全不符合Fortran自由格式的规则，所以这个问题无法通过调整格式彻底解决，只能避免直接包含C头文件。

最后，建议你移除-lcufftw链接选项，这个是cuFFTW兼容库，用于让FFT代码无缝切换到cuFFT，但你已经在重构代码用原生cuFFT，这个库没必要再链接了。

内容的提问来源于stack exchange，提问作者riky_cv