Blas中是否有任何功能可以执行偏斜的矩阵矢量产物？

Question

我正在考虑使用Intel-MKL进行一些计算，特别是用于Fortran中程序的矩阵矢量稀疏BLAS函数。

矩阵表达我的计算，这些矩阵恰好是稀疏和偏斜的

我可以用我所看到的 ，稀疏的blas具有稀疏的一般和对称矩阵的功能，所以我想知道我是否有一种方法可以使用稀疏的偏斜工作 - 合成矩阵，因为我认为这会减少内存足迹。

Answer 1

tldr; Mkl稀疏Blas可以用 mkl_scsrmv subroutine subroutine subroutine and subsing 'a'a'a'在矩阵描述符阵列。

好的，当我开始测试一般MKL稀疏Blas矩阵 - 矢量乘法CSR格式时，我设法找到了问题的答案（ mkl_？csrmv ），

我了解到有一个字符数组用于描述输入矩阵（ matdescra ）。该数组中的第一个字符可以设置为'a'，该字符导致子例程将输入矩阵解释为偏斜 - 对称。例如（不一定是一个好的），

给定矩阵，a和vector，x，：

A = [ 0  1  2      x = [ 1
     -1  0  3            2
     -2 -3  0 ]           3 ]

a 的上三角形可以

val = [1, 2, 3]
col = [2, 3, 3]
rowstart = [1, 3, 3]
rowend = [3 3 4]

用字符阵列 matdescra =表示为和 代表。 ['a'，'u'，'n'，'f'] ，
矩阵向量产品是通过

call mkl_scsrmv('n', 3, 3, 1., matdescra, val, rowstart, rowend, x, 1., y

将输出（向量）添加到向量阵列的地方获得的， y 。