为什么CuBlas不对三角形矩阵使用一维数组？

Question

这可能是对旧的BLAS库设计的一种回想，但我刚才感到惊讶的是，发现CuBlas通过对三角形矩阵使用常规2d数组来浪费内存。 我想这使与API其余部分的接口更加方便。

Answer 1

刚才我很惊讶地发现CuBlas通过对三角形矩阵使用常规2d数组浪费了内存

严格说来并非如此。

如果您查看2级BLAS例程，您会发现它们在以打包格式存储的三角或Hermitian矩阵上运行。

3级BLAS例程没有，但是有两个很好的理由将它们以完全密集的格式存储。

• BLAS就是这样

• 这些例程大部分已添加到BLAS中，以支持LAPACK求解器。 而且这些求解器通常在提供的全密集输入中原位存储分解的结果，因此在BLAS中使用该格式是合乎逻辑的

我想如果您不喜欢这种设计选择，可以随时尝试写信给Jack Dongarra投诉。