在 BLAS 中转置还是先自己做?
我正在用 Fortran 77 编写一些科学代码,并且正在讨论什么会更快。
基本上,我有一个 MxN 矩阵,我们称之为 A。M 大于 N。稍后在代码中,我需要将转置(A)乘以一堆向量。
我的问题是,采用 A,自行转置并存储它,或者当我调用 BLAS 时,只给它转置标志,会更快吗?
谢谢! -帕特里克
I'm putting together some scientific code in Fortran 77, and I am having a debate on what would be faster.
Basically, I have an MxN matrix, let's call it A. M is larger than N. Later on in the code, I need to multiply transpose(A) by a bunch of vectors.
My question is, would it be faster to take A, transpose it on my own and store that, or when I call BLAS, just give it the transpose flag?
Thanks!
-Patrick
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论
评论(1)
我的直觉告诉我使用转置标志。
在这种情况下,你会以 1 的步长进行大量的点积。
事实上,如果不实际运行代码,很难判断。
现代 blas 采用缓存阻塞技术,这使得简单的分析变得困难。
my gut feeling tells me to use transpose flag.
in that case you doing lots of dot products with stride of one.
In reality, it's very hard to tell without actually running codes.
modern blas employs cache blocking techniques which make simple analysis difficult at best.