为什么 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 编写的，但它不是最快的。

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

由于@jxh 正确地设置，调用函数与调用子例程（类似于 void 函数）存在更大的 ABI 不兼容风险。