试图在 C++ 中调用 fortran 子例程

Question

我对 C++ 的了解非常有限，对 fortran 的了解更少，我目前正在尝试从 c++ 主程序调用 fortran 子例程。按照一些示例，我能够想出以下代码来调用

subroutine fireballess(ear,ne,parames,ifl,photar,photer)

这是我的 C++ 代码：

#include <stdlib.h>
#include <stdio.h>
#include <iostream>

using namespace std;

extern "C" void fireballess_( double *fear, int fne,double* fparames, int fifl, double *fphotar, double *fphoter);

int main(int argc, char ** argv)
{
    int ne,ifl;
    double *ear;
    double *parames;

    double *photar;
    double *photer;


    parames = new double[9];

// parames=[4.3,0.23,0.5,0.5,1.5,1.,1000.,2.15,3.]
        parames[0]=4.3;
        parames[1]=0.23;
        parames[2]=0.5;
        parames[3]=0.5;
        parames[4]=1.5;
        parames[5]=1.;
        parames[6]=1000.;
        parames[7]=2.15;
        parames[8]=3.;

ne = 2;

    ear = new double[ne];
 ear[0] = 0.;
 ear[1] = 20.;
 ear[2] = 40.;
 ifl=2;

    photar = new double[ne];
    photer = new double[ne];

    // Call a Fortran subroutine
    //subroutine_sum_(&size,vec,&sum);
    fireballess_(&ear,ne,&parames,ifl,&photar,&photer);

    cout << "Calling a Fortran subroutine" << endl;
    cout << "===============================" << endl;

for (int i=0;i<=ne;i++){
    cout << "ear = " <<ear[i-1]<< " - "<<ear[i] << endl;
    cout << "photar = " << photar[i] << endl;
    cout << "photer = " << photer[i] << endl << endl;
}

    delete[] ear;
    delete[] parames;
    delete[] photar;
    delete[] photer;
}

但是，当我尝试编译时，出现以下错误：

call_fortran.cpp: In function ‘int main(int, char**)’:
call_fortran.cpp:45:53: error: cannot convert ‘double**’ to ‘double*’ for argument ‘1’ to ‘void fireballess_(double*, int, double*, int, double*, double)’
     fireballess_(&ear,ne,&parames,ifl,&photar,photer);

我不明白，因为 photer 变量遵循与 photar 变量完全相同的路径，因此不会出现错误。 我希望比我更了解指针的人可以在这里帮助我。 谢谢

Answer 1

ear 是 double* 类型，因此 &ear 是 double** 类型，这对于 fortran 函数原型是不合适的。摆脱&s 可能解决问题（&parames, &photar 和 &photer 也是双**）：

改变

fireballess_(&ear,ne,&parames,ifl,&photar,&photer);

到

fireballess_(ear,ne,parames,ifl,photar,photer);

fortran 函数原型似乎是void fireballess_(double*, int, double*, int, double*, double)。如果确实如此，那么photer 仍然必须更改为photer[i]，其中（i 是一个数组索引，如 0、1、2...）。

如果你真的想了解发生了什么，我建议阅读一本好的 C++ 指针教程。