36.cuBLAS开发指南中文版--cuBLAS中的Level-2函数hpmv()

2.6.19. cublashpmv()

cublasStatus_t cublasChpmv(cublasHandle_t handle, cublasFillMode_t uplo,
                           int n, const cuComplex       *alpha,
                           const cuComplex       *AP,
                           const cuComplex       *x, int incx,
                           const cuComplex       *beta,
                           cuComplex       *y, int incy)
cublasStatus_t cublasZhpmv(cublasHandle_t handle, cublasFillMode_t uplo,
                           int n, const cuDoubleComplex *alpha,
                           const cuDoubleComplex *AP,
                           const cuDoubleComplex *x, int incx,
                           const cuDoubleComplex *beta,
                           cuDoubleComplex *y, int incy)
1
2
3
4
5
6
7
8
9
10
11
12

$y=\alpha Ax+\beta y$

其中 A 是以压缩格式存储的 nxn Hermitian 矩阵，x 和 y 是向量，α 和 β 是标量。

如果 uplo == CUBLAS_FILL_MODE_LOWER 那么Hermitian矩阵 A 的下三角部分的元素被逐列无间隙地打包在一起，使得元素 A(i,j) 存储在内存位置 AP[i+((2* n-j+1)*j)/2] 对于 j=1,…,n 和 i>=j。 因此，打包格式只需要 n(n+1)/2 个元素进行存储。

如果 uplo == CUBLAS_FILL_MODE_UPPER 则将厄米矩阵 A 的上三角部分的元素逐列无间隙地打包在一起，使得元素 A(i,j) 存储在内存位置 AP[i+(j*( j+1))/2] 对于 j=1,…,n 和 i<=j 。 因此，打包格式只需要 n(n+1)/2 个元素进行存储。

该函数可能返回的错误值及其含义如下所列。

详细信息, 请参考:

chpmv, zhpmv