cuBLAS中的Level-1函数asum()和axpy()

2.5.3. cublasasum()

cublasStatus_t  cublasSasum(cublasHandle_t handle, int n,
                            const float           *x, int incx, float  *result)
cublasStatus_t  cublasDasum(cublasHandle_t handle, int n,
                            const double          *x, int incx, double *result)
cublasStatus_t cublasScasum(cublasHandle_t handle, int n,
                            const cuComplex       *x, int incx, float  *result)
cublasStatus_t cublasDzasum(cublasHandle_t handle, int n,
                            const cuDoubleComplex *x, int incx, double *result)
1
2
3
4
5
6
7
8

此函数计算向量 x 的元素的绝对值之和。 因此，结果是 ∑ i = 1 n | I m ( x [ j ] ) | + | R e ( x [ j ] ) | 其中 j = 1 + ( i - 1 ) * incx 。 请注意，最后一个等式反映了用于与 Fortran 兼容的基于 1 的索引。

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

请参考:

sasum, dasum, scasum, dzasum

2.5.4. cublasaxpy()

cublasStatus_t cublasSaxpy(cublasHandle_t handle, int n,
                           const float           *alpha,
                           const float           *x, int incx,
                           float                 *y, int incy)
cublasStatus_t cublasDaxpy(cublasHandle_t handle, int n,
                           const double          *alpha,
                           const double          *x, int incx,
                           double                *y, int incy)
cublasStatus_t cublasCaxpy(cublasHandle_t handle, int n,
                           const cuComplex       *alpha,
                           const cuComplex       *x, int incx,
                           cuComplex             *y, int incy)
cublasStatus_t cublasZaxpy(cublasHandle_t handle, int n,
                           const cuDoubleComplex *alpha,
                           const cuDoubleComplex *x, int incx,
                           cuDoubleComplex       *y, int incy)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

此函数将向量 x 乘以标量并将其添加到向量 y 中，并用结果覆盖最新的向量。 因此，对 i = 1 , … , n 、 k = 1 + ( i - 1 ) * incx 和 j = 1 + ( i - 1 ) * incy 执行的操作是 y [ j ] = α × x [ k ] + y [ j ]。 请注意，最后两个等式反映了用于与 Fortran 兼容的基于 1 的索引。

请参考:

saxpy, daxpy, caxpy, zaxpy