C 代码调用 Fortran 子程序时，子程序顶部出现分段错误

我的文件中有 C++ 代码test-Q.cpp调用文件中的 Fortran 子例程getqpf.F。在文件中test-Q.cpp，我已将 Fortran 代码声明为外部代码，并且我正在使用getqpf_()名称修改约定。这gcc and gfortran编译器正在 GNU/Linux 上使用。

以下是 C++ 文件顶部的片段：

extern "C" {
            void  getqpf_  (double *tri, 
                    int nsamp, 
                    int lwin,
                    int nfreqfit, 
                    double dt, 
                    float null, 
                    int L2,
                    double df,
                    double *qq, 
                    double *pf, 
                    double *ampls, 
                    double *work1, 
                    double *work2, 
                    double *work3, 
                    double *work4,
                    int mem, 
                    int morder, 
                    int nfs, 
                    double *xReal, 
                    double *xImag, 
                    double *xAbs,
                    double *x1,
                    int cen,
                    int top,
                    int bot, 
                    float cut,
                    int nfst,
                    int raw);  

        } // end

以下是 Fortran 文件中的相应片段：

   subroutine getqpf (tri, nsamp, lwin, nfreqfit, dt, null, L2, df,
     1                   qq, pf, ampls, work1, work2, work3, work4,
     2                   mem, morder, nfs, xReal, xImag, xAbs, x1,
     3                   cen,top,bot, cut,nfst,raw)



      integer  morder, lwin, nsamp, nfreqfit, delay, nfs

      real     tri(*)
      real     qq(*), pf(*), ampls(*)

      real * 8 work1(*), work2(*), work3(*), work4(*)
      real * 8 xReal(*), xImag(*), xabs(*), x1(*)

      real * 8 dt8, cut8, df8
      real     null, cut
      integer  nfst
      logical  mem, L2, cen, top, bot, raw


      integer nf

C program logic code starts here
          nf = nfreqfit
          delay = 0
          dt8  = dt
          cut8 = cut

Fortran 代码调用其他 C 代码函数。在 GNU/Linux 上使用gfortran and gcc我已按以下方式编译并链接了所有文件：

 g++ -c test-Q.cpp -I./boost/boost_1_52_0/ -g
 gcc -c paul2.c -g
 gcc -c paul2_L1.c -g
 gcc -c paul6.c -g
 gcc -c paul6_L1.c -g 
 gcc -c fit_slope.c -g
 gfortran -c getqpf.F -g
 g++ -o test-Q test-Q.o paul2.o paul2_L1.o paul6.o paul6_L1.o fit_slope.o getqpf.o -g

尽管我能够成功构建二进制文件，但该行出现了段错误nf = nfreqfit。它位于 Fortran 文件的最顶部。跑步gdb二进制文件会产生以下输出：

Program received signal SIGSEGV, Segmentation fault.
0x0000000000406fd3 in getqpf (tri=..., nsamp=Cannot access memory at address 0x3e9
) at getqpf.F:44
44        nf = nfreqfit

这里发生了什么，为什么会出现段错误？内存似乎没有在 C++ 代码和 Fortran 代码之间正确传递。

UPDATE

正如 IanH 在下面的答案中提到的，问题是由于没有通过引用传递参数造成的。使用 C++，该函数必须声明为：

 extern"C" {
            void  getqpf_  (float *tri, 
                    int &nsamp, 
                    int &lwin,
                    int &nfreqfit, 
                    float &dt, 
                    float &null, 
                    int &L2,
                    float &df,
                    float *qq, 
                    float *pf, 
                    float *ampls, 
                    double *work1, 
                    double *work2, 
                    double *work3, 
                    double *work4,
                    int &mem, 
                    int &morder, 
                    int &nfs, 
                    double *xReal, 
                    double *xImag, 
                    double *xAbs,
                    double *x1,
                    int &cen,
                    int &top,
                    int &bot, 
                    float &cut,
                    int &nfst,
                    int &raw);  

        } // end 

请注意 & 符号的存在。然后，可以在代码中调用该函数，如下所示：

getqpf_ (tri,       
    nsamp, 
    lwin,
    nfreqfit, 
    dt, 
    null, 
    L2,
    df,
    qq, 
    pf, 
    ampls, 
    work1, 
    work2, 
    work3, 
    work4,
    mem, 
    morder, 
    nfs, 
    xReal, 
    xImag, 
    xAbs,
    x1,
    cen,
    top,
    bot, 
    cut,
    nfst,
    raw); 

请注意，变量如nsamp被声明为int nsamp = 1001.

在赞同 M.S.B. 关于使用 F2003 的 C 互操作性的建议时，请注意，您的具体问题是按引用传递/按值传递不匹配（即使在使用 C 互操作性时，这仍然是您必须考虑的问题）。典型的 Fortran 实现通过引用传递所有参数，而在 C(++) 中默认是通过值传递。在 C++ 方面，请注意所有 int 和 float 参数以及一些 double 参数缺少指针说明符 (*）。这些参数是按值传递的 - 但 Fortran 方面没有任何内容表明这一点。在 F2003 之前，这通常是使用 Fortran 代码中的编译器特定指令来完成的。

使用 F2003 的 C 互操作，具有 BIND(C) 属性的过程参数的默认传递约定是通过引用。按值传递的参数需要在其声明中具有 VALUE 属性。