为什么在 C 编程中需要 Fortran 包装器？

Question

我最近在阅读 cBLAS 中的一些源代码，有些地方让我不清楚。在很多函数中，.c 文件调用 Fortran Wrapper 而不是直接在 C 文件中编写代码，如下文件：

/*
 * cblas_sdsdot.c
 *
 * The program is a C interface to sdsdot.
 * It calls the fortran wrapper before calling sdsdot.
 *
 * Written by Keita Teranishi.  2/11/1998
 *
 */
#include "cblas.h"
#include "cblas_f77.h"
float cblas_sdsdot( const int N, const float alpha, const float *X,
                      const int incX, const float *Y, const int incY)
{
   float dot;
#ifdef F77_INT
   F77_INT F77_N=N, F77_incX=incX, F77_incY=incY;
#else 
   #define F77_N N
   #define F77_incX incX
   #define F77_incY incY
#endif
   F77_sdsdot_sub( &F77_N, &alpha, X, &F77_incX, Y, &F77_incY, &dot);
   return dot;
}

我完全糊涂了，为什么要这样做？是因为Fortran的计算效率更高吗？

Answer 1

"其实我想问的是为什么需要中间包装器，为什么不用C写呢？"

整个 CBLAS 是 BLAS 的包装器。 BLAS 是使用参考 Fortran 实现和 a Fortran API 定义的。 BLAS 可以用 C 或汇编实现，但 API 设置为 Fortran。

因此，CBLAS 实际上并不包含全部功能。该功能位于您安装的任何 BLAS 实现中。非常常见的参考实现是用 Fortran 编写的，但它不是最快的。

但是，您可能可以直接从 C cblas_sdsdot 调用 sdsdot 函数（无论它实际实现的语言是什么）。 CBLAS 的作者选择实现一个Fortran 中间子程序sdsdotsub。我不知道为什么现在有必要这样做。差别很小，真的只是把一个函数改成一个子程序而已。

正如@jxh 正确评论的那样，调用函数与调用子例程（类似于 void 函数）相比，ABI 不兼容的风险更大。

为什么在 C 编程中需要 Fortran 包装器？

Why is fortran wrapper needed in C programming?

c

fortran

wrapper

blas