0

我在文件中有 C++ 代码,它在文件test-Q.cpp中调用 Fortran 子例程getqpf.F。在 filetest-Q.cpp中,我已将 Fortran 代码声明为外部代码,并且我正在使用getqpf_()名称修改约定调用该函数。和编译gccgfortran正在 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 代码函数。gfortran在使用和编译器的 GNU/Linux 上,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 代码之间正确传递内存。

更新

正如 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.

4

2 回答 2

3

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

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

于 2012-11-16T21:02:15.323 回答
2

我推荐使用 Fortran ISO C 绑定。Stackoverflow 和 gfortran 手册中有示例。它是 Fortran 2003 语言标准的一部分,在此之前是 Fortran 95 的技术报告。这使其编译器和平台可移植。您不必担心编译器特定的调用约定或名称修改。

于 2012-11-16T17:50:01.233 回答