C-Fortran混合编程

Question

我正在进行一个项目，其中我必须使用C中的Fortran库。在Fortran库中，有一个包含复杂*16，4x4数组的公共块。在C中，复变量只是包含两个元素的结构，由于它是复数*16，元素应该是长双，这是Fortran中对应的C数据类型。所以我有两个长双打的结构。

我将访问这个数组的元素。好的是，我已经可以和库中的所有其他公共变量一起使用它了。问题是，我从数组中导入的元素是，

1)不应按顺序排列，“即使考虑了C和数组结构的差异”。

2)虽然大多数要素都是正确的，但两个要素与它们应有的性质有很大不同。

( 3)只有在使用double而不是long double的情况下，我才能得到正确的元素(除了这两个元素)。当我使用长双(以及正确的字符转换)时，我得到了完全不同的东西，这清楚地指出了转换的问题。

我已经用尽了所有的解释，但都没有用。我在C中使用priting数组的代码如下所示：

for (j=0;j<=3;j++){
    printf("%s", "\n");
    for(k=0;k<=3;k++){            
        printf("%s %d %s %d %s %s %LE %s %LE %s",
          "(", k+1, "," ,j+1, ")", "{",
          (long double)mssmmixing_.neunmx[k][j].dr,
          " ",
          (long double)mssmmixing_.neunmx[k][j].di,
          "}\n");           
    }
}

附加信息:由于我必须混合Fortran对象文件，所以我使用gfortran来编译C文件。如果我使用GNU编译器，它会引发关于不识别gfortran例程的错误。这也可能是一个问题的根源，可能是gfortran不承认长双打在C。

任何帮助都是有用的。

Answer 1

对于混合Fortran和C，我建议使用ISO_C_Binding。它甚至有一个Fortran类型的C_LONG_DOUBLE_COMPLEX，它与C类型的long double _Complex相匹配--参见http://gcc.gnu.org/onlinedocs/gfortran/ISO_005fC_005fBINDING.html。作为Fortran语言标准的一部分，这些类型保证匹配(当您使用兼容的编译器时)。很可能C_LONG_DOUBLE_COMPLEX实际上与复杂*16相同，但在Fortran中，您可以尝试在这两种类型之间进行复制，以防内存布局不同。您应该用gfortran编译Fortran源文件，用gcc编译C。这是最容易链接到gfortran。(或者使用其他“品牌”编译器。)