?spmv

向量与压缩的对称矩阵乘积

即 $\text{[math]}$ 。alpha和beta是缩放系数，x和y是含有n个元素的向量，A是n阶压缩的对称矩阵。

接口定义

C interface：

void cblas_sspmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const float alpha, const float *Ap, const float *X, const BLASINT incX, const float beta, float *Y, const BLASINT incY);

void cblas_dspmv(const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo, const BLASINT N, const double alpha, const double *Ap, const double *X, const BLASINT incX, const double beta, double *Y, const BLASINT incY);

Fortran interface：

CALL SSPMV(UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY)

CALL DSPMV(UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY)

参数

参数名	类型	描述	输入/输出
order	枚举类型CBLAS_ORDER	表示矩阵是行主序或列主序。	输入
Uplo	枚举类型CBLAS_UPLO	表示使用矩阵A存储展开方式（上三角或下三角）。如果Uplo=CblasUpper，则使用A的上三角展开。如果Uplo=CblasLower，则使用A的下三角展开。	输入
N	整型数	表示矩阵A的阶数，要求N大于或等于零。	输入
alpha	在sspmv中是单精度复数类型。在dspmv中是双精度复数类型。	表示乘法系数。	输入
Ap	在sspmv中是单精度复数类型。在dspmv中是双精度复数类型。	压缩的对称矩阵，矩阵规模至少是(n*(n+1)/2)。	输入
X	在sspmv中是单精度复数类型。在dspmv中是双精度复数类型。	向量X，向量规模至少是(1+(N-1)*abs(incX))。	输入
incX	整型数	表示X中向量增长步长，不能为零。	输入
beta	在sspmv中是单精度复数类型。在dspmv中是双精度复数类型。	乘法系数。	输入
Y	在sspmv中是单精度复数类型。在dspmv中是双精度复数类型。	向量Y，向量规模至少是(1+(N-1)*abs(incY))。	输入/输出
incY	整型数	表示Y中向量增长步长，不能为零。	输入

依赖

#include "kblas.h"

示例

C interface：

    int n = 3; 
    float alpha = 1.0, beta = 1.0; 
    int incx = 1, incy = 1; 
    /* 
     * A =  | 2.0  1.0  3.0 | 
     *      | 1.0  6.0  9.0 | 
     *      | 3.0  9.0  2.0 | 
     */ 
    float a[6] = {2.0, 1.0, 3.0, 6.0, 9.0, 2.0}; 
    float x[3] = {1.0, 1.0, 1.0}; 
    float y[3] = {3.0, 2.0, 2.0}; 
 
    cblas_sspmv(CblasColMajor,CblasLower, n, alpha, a, x, incx, beta, y, incy); 
    /* 
     * Output y = |9.0, 18.0, 16.0| 
     */

Fortran interface：

      INTEGER :: N=3 
      REAL(4) :: ALPHA=2.0 
      REAL(4) :: BETA=1.0 
      REAL(4) :: A(6), X(3), Y(3) 
      DATA A/2.0, 1.0, 3.0, 6.0, 9.0, 2.0/ 
      DATA X/1.0, 1.0, 1.0/ 
      DATA Y/3.0, 2.0, 2.0/ 
      EXTERNAL SSPMV 
      CALL SSPMV('L', N, ALPHA, A, X, INCX, BETA, Y, INCY) 
 
*     Output Y = |9.0, 18.0, 16.0|

父主题： KML_BLAS Level 2函数