为什么 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 实现和Fortran API定义的。BLAS可以用C或汇编语言实现，但API设置为Fortran。

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

但是，您可能可以sdsdot直接从 C 调用该函数（以实际实现的任何语言）cblas_sdsdot。CBLAS的作者选择实现一个Fortran中间子例程sdsdotsub。我现在不知道为什么有必要这样做。区别非常小，实际上只是将函数更改为子例程。

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