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

2.6.7. cublasspr2()

cublasStatus_t cublasSspr2(cublasHandle_t handle, cublasFillMode_t uplo,      int n, const float  *alpha,      const float  *x, int incx,      const float  *y, int incy, float  *AP)cublasStatus_t cublasDspr2(cublasHandle_t handle, cublasFillMode_t uplo,      int n, const double *alpha,      const double *x, int incx,      const double *y, int incy, double *AP)

此函数执行打包的对称 rank-2 更新

y = α ( x y T + y x T ) + A y = \alpha (xy^T + yx^T)+A y=α(xyT+yxT)+A

其中 A 是以压缩格式存储的 n*n 对称矩阵，x 是向量，而 $\alpha$ 是标量。

如果 uplo == CUBLAS_FILL_MODE_LOWER 则对称矩阵 a 的下三角部分中的元素被逐列填充在一起，没有间隙. 因此，一般来说，元素 A(i,j) 存储在内存位置 AP[i+((2*n-j+1)*j)/2] 中 j=1,…,n 和 $i>=j$ 。 因此，打包格式只需要存储n(n+1) 2 \frac{n(n+1)}{2} 2n(n+1)​个元素。

如果 uplo == CUBLAS_FILL_MODE_UPPER 则对称矩阵 a 的下三角部分中的元素被逐列填充在一起，没有间隙. 因此，一般来说，元素 A(i,j) 存储在内存位置 AP[i+(j*(j+1))/2] 中 j=1,…,n 和 $i<=j$ 。 因此，打包格式只需要存储n(n+1) 2 \frac{n(n+1)}{2} 2n(n+1)​个元素。

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

请参考:
sspr2,