---
title: ?spmv
description: "向量与压缩的对称矩阵乘积"
url: https://www.hikunpeng.com/document/detail/zh/kunpengboostkithistory/240RC1/accel/kunpengaccel_kml_16_0045.html
sourcePath: /source/zh/kunpengboostkithistory/240RC1/accel/kunpengaccel_kml_16_0045.html
indexId: 2db74f01c5153defbf0b773a0f3f62a9ba5ccbdbbf742ce70efdf3ef65b515c576
---
# ?spmv

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

即。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|
```