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

2.6.22. cublashpr()

cublasStatus_t cublasChpr(cublasHandle_t handle, cublasFillMode_t uplo,
                          int n, const float *alpha,
                          const cuComplex       *x, int incx,
                          cuComplex       *AP)
cublasStatus_t cublasZhpr(cublasHandle_t handle, cublasFillMode_t uplo,
                          int n, const double *alpha,
                          const cuDoubleComplex *x, int incx,
                          cuDoubleComplex *AP)

1
2
3
4
5
6
7
8
9

此函数执行打包的 Hermitian rank-1 更新

$A=\alpha x{x}^H+A$

其中 A 是以列优先格式存储的 NxN Hermitian 矩阵，x 是向量，$\alpha$ 是标量。

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

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

详细信息, 请参考:

chpr, zhpr