fortran - 在 BLAS 中转置还是先自己做？

Question

我在 Fortran 77 中整理了一些科学代码，我正在争论什么会更快。

基本上，我有一个 MxN 矩阵，我们称之为 A。M 大于 N。稍后在代码中，我需要将 transpose(A) 乘以一堆向量。

我的问题是，取 A、自己转置并存储它会更快，还是当我调用 BLAS 时，只给它转置标志？

谢谢!
-帕特里克

Answer 1

我的直觉告诉我使用转置标志。
在这种情况下，您会以 1 的步幅进行大量的点积。

实际上，如果不实际运行代码，很难判断。
现代 blas 采用缓存阻塞技术，这使得简单的分析充其量是困难的。