---
title: ?spr2
description: "对称展开矩阵秩2更新。"
url: https://www.hikunpeng.com/document/detail/zh/kunpengboostkithistory/240RC2/accel/kunpengaccel_kml_16_0047.html
sourcePath: /source/zh/kunpengboostkithistory/240RC2/accel/kunpengaccel_kml_16_0047.html
indexId: 5b6171c424599be3b9265576d0e6c3aa9272494a0a7f0f34dc434fdd599a4fa776
---
# ?spr2

对称展开矩阵秩2更新。

即：。alpha为乘法系数，x和y为包含n个元素的向量，A为n*n的三角展开对称矩阵。

#### 接口定义

C interface：

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

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

Fortran interface：

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

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


#### 参数

| 参数名 | 类型 | 描述 | 输入/输出 |
| --- | --- | --- | --- |
| order | 枚举类型CBLAS\_ORDER | 表示矩阵是行主序或列主序。 | 输入 |
| Uplo | 枚举类型CBLAS\_UPLO | 对称矩阵展开存储方式（上三角或下三角展开）。 如果Uplo=CblasUpper，则使用A的上三角展开。 如果Uplo=CblasLower，则使用A的下三角展开。 | 输入 |
| N | 整型数 | 向量X的元素个数。 | 输入 |
| alpha | 在sspr2中是单精度复数类型。 在dspr2中是双精度复数类型。 | 乘法系数。 | 输入 |
| X | 在sspr2中是单精度复数类型。 在dspr2中是双精度复数类型。 | 矩阵X，长度至少 1+(n\-1)\*abs(incX)。 | 输入 |
| incX | 整型数 | 向量X的增长步长，不能为零。 | 输入 |
| Y | 在sspr2中是单精度复数类型。 在dspr2中是双精度复数类型。 | 矩阵Y，长度至少 1+(n\-1)\*abs(incY)。 | 输入 |
| incY | 整型数 | 向量Y的增长步长，不能为零。 | 输入 |
| A | 在sspr2中是单精度复数类型。 在dspr2中是双精度复数类型。 | 矩阵A三角存储。 | 输入/输出 |


#### 依赖

#include "kblas.h"


#### 示例

C interface：

```
float alpha = 1.0;
int n = 3, incx = 1, incy = 1;
float x[3] = {2.0, 2.0, 1.0};
float y[3] = {1.0, 2.0, 1.0};
/**
*             | 3.0  1.0  2.0 |
*         A = | 1.0  6.0  3.0 |
*             | 2.0  3.0  3.0 |
*/
float a[6] = {3.0, 1.0, 2.0, 6.0, 3.0, 3.0};
cblas_sspr2(CblasColMajor,CblasLower, n, alpha, x, incx, y, incy, a);
/**
*  Output a = |7.0, 7.0, 5.0, 14.0, 7.0, 5.0|
*/
```

Fortran interface：

```
REAL(4) :: ALPHA=1.0
INTEGER :: N=3
INTEGER :: INCX=1
INTEGER :: INCY=1
REAL(4) :: X(3), Y(3), A(6)
DATA X/2.0, 2.0, 1.0/
DATA Y/1.0, 2.0, 1.0/
DATA A/3.0, 1.0, 2.0, 6.0, 3.0, 3.0/
EXTERNAL SSPR2
CALL SSPR2('L', N, ALPHA, X, INCX, Y, INCX, A)
*     Output A = |7.0, 7.0, 5.0, 14.0, 7.0, 5.0|
```