我在将 char 数组从 C++ 传递到 Fortran (f90) 时遇到问题。

这是我的 C++ 文件，“cmain.cxx”:

#include <iostream>

using namespace std;

extern "C" int ftest_( char (*string)[4] );

int main() {
    char string[2][4];

    strcpy(string[0],"abc");
    strcpy(string[1],"xyz");

    cout << "c++: string[0] = '" << string[0] << "'" << endl;
    cout << "c++: string[1] = '" << string[1] << "'" << endl;

    ftest_(string);

    return 0;
}

这是我的 Fortran 文件“ftest.f90”:

SUBROUTINE FTEST(string)

CHARACTER*3 string(2)
CHARACTER*3 expected(2)
data expected(1)/'abc'/
data expected(2)/'xyz'/

DO i=1,2
    WRITE(6,10) i,string(i)
10  FORMAT("fortran: string(",i1,") = '", a, "'" )

    IF(string(i).eq.expected(i)) THEN
        WRITE(6,20) string(i),expected(i)
20      FORMAT("'",a,"' equals '",a,"'")
    ELSE
        WRITE(6,30) string(i),expected(i)
30      FORMAT("'",a,"' does not equal '",a,"'")
    END IF
ENDDO

RETURN
END

构建过程是:

gfortran -c -m64   ftest.f90 
g++ -c  cmain.cxx
gfortran -m64 -lstdc++ -gnofor_main -o test ftest.o cmain.o

编辑:请注意，可执行文件也可以通过以下方式构建:

g++ -lgfortran -o test ftest.o cmain.o

此外，-m64 标志是必需的，因为我运行的是 OSX 10.6。

执行'test'的输出是:

c++: string[0] = 'abc'
c++: string[1] = 'xyz'
fortran: string(1) = 'abc'
'abc' equals 'abc'
fortran: string(2) = 'xy'
'xy' does not equal 'xyz'

在 ftest.f90 中声明大小为 4 的 'string' 和 'expected' 字符数组，即:

CHARACTER*4 string(2)
CHARACTER*4 expected(2)

重新编译后输出如下:

c++: string[0] = 'abc'
c++: string[1] = 'xyz'
fortran: string(1) = 'abc'
'abc' does not equal 'abc '
fortran: string(2) = 'xyz'
'xyz' does not equal 'xyz '

在“cmain.cxx”中声明大小为 3 的字符数组，即:

extern "C" int ftest_( char (*string)[3] );

int main() {
    char string[2][3];

并恢复到 fortran 文件 (3) 中的原始大小，即:

CHARACTER*3 string(2)
CHARACTER*3 expected(2)

重新编译后输出如下:

c++: string[0] = 'abcxyz'
c++: string[1] = 'xyz'
fortran: string(1) = 'abc'
'abc' equals 'abc'
fortran: string(2) = 'xyz'
'xyz' equals 'xyz'

所以最后一种情况是唯一可行的，但在这里我将 3 个字符分配给大小为 3 的字符数组，这意味着缺少终止符 '\0'，并导致输出 'abcxyz' - 这是不适合我的预期应用。

非常感谢任何帮助，这让我抓狂!

最佳答案

C 字符串以零结尾，而 Fortran 字符串按照惯例以空格填充但大小固定。您不应该期望能够在不进行某些转换的情况下将 C 字符串传递给 Fortran。

例如:

#include <algorithm>

void ConvertToFortran(char* fstring, std::size_t fstring_len,
                      const char* cstring)
{
    std::size_t inlen = std::strlen(cstring);
    std::size_t cpylen = std::min(inlen, fstring_len);

    if (inlen > fstring_len)
    {
        // TODO: truncation error or warning
    }

    std::copy(cstring, cstring + cpylen, fstring);
    std::fill(fstring + cpylen, fstring + fstring_len, ' ');
}

然后您可以将其与 ftest 的 3 或 4 长度版本一起使用:

#include <iostream>
#include <ostream>
extern "C" int ftest_( char string[][4] );

void ConvertToFortran(char* fstring, std::size_t fstring_len,
                      const char* cstring);

int main()
{
    char cstring[2][4] = { "abc", "xyz" };
    char string[2][4];

    ConvertToFortran(string[0], sizeof string[0], cstring[0]);
    ConvertToFortran(string[1], sizeof string[1], cstring[1]);

    std::cout << "c++: string[0] = '" << cstring[0] << "'" << std::endl;
    std::cout << "c++: string[1] = '" << cstring[1] << "'" << std::endl;

    ftest_(string);

    return 0;
}

关于c++ - 将 char 数组从 C++ 传递到 Fortran，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/10163485/